My augmented Lagrangian local minimization code gives slightly different results for the same inputs for a problem involving packing 50 integer circles.
Past work has shown that circle packing is chaotic in that slightly different inputs give different results but I don't understand why the same input gives different results. For the first run, I get an objective function value of 232.568 and 232.132 for the second. The starting point is generated randomly but is used for both runs.
The attached notebook has the code. It was run with Mathematica Version 10.1 on a PC running Windows 8.
Attachments: