Packing of Circles on Square Flat Torus as Global Optimization of Mixed Integer Nonlinear problem

The article demonstrates rather general approach to problems of discrete geometry: treat them as global optimization problems to be solved by one of general purpose solver implementing branch-and-bound algorithm. This approach may be used for various types of problems, i.e. Tammes problems, Thomson problems, search of minimal potential energy of micro-clusters, etc. Here we consider a problem of densest packing of equal circles in special geometrical object, so called square flat torus $\mathbb{R}^2/\mathbb{Z}^2$ with the induced metric. It is formulated as Mixed-Integer Nonlinear Problem with linear and non-convex quadratic constraints. The open-source solver SCIP, this http URL, and its parallel implementation ParaSCIP, this http URL, had been used in computing experiments. The main result is "computer aided" proof of the conjecture on optimal packing for N=9 that was published in 2012. To do that, ParaSCIP took about 2059 CPU*hours (16 hours x 128 CPUs) of cluster HPC4/HPC5, National Research Centre "Kurchatov Institute".

[1]  Thorsten Koch,et al.  Solving Open MIP Instances with ParaSCIP on Supercomputers Using up to 80,000 Cores , 2016, 2016 IEEE International Parallel and Distributed Processing Symposium (IPDPS).

[2]  Brian W. Kernighan,et al.  AMPL: A Modeling Language for Mathematical Programming , 1993 .

[3]  János D. Pintér,et al.  Solving circle packing problems by global optimization: Numerical results and industrial applications , 2008, Eur. J. Oper. Res..

[4]  Oleg R. Musin,et al.  Optimal Packings of Congruent Circles on a Square Flat Torus , 2012, Discret. Comput. Geom..

[5]  Thorsten Koch,et al.  ParaSCIP: A Parallel Extension of SCIP , 2010, CHPC.

[6]  S. V. Smirnov,et al.  On Domain Decomposition Strategies to Parallelize Branch-and-Bound Method for Global Optimization in Everest Distributed Environment , 2018 .

[7]  David L. Woodruff,et al.  Pyomo — Optimization Modeling in Python , 2012, Springer Optimization and Its Applications.

[8]  Ambros M. Gleixner,et al.  SCIP: global optimization of mixed-integer nonlinear programs in a branch-and-cut framework , 2018, Optim. Methods Softw..

[9]  Benjamin Müller,et al.  The SCIP Optimization Suite 3.2 , 2016 .

[10]  S. V. Smirnov,et al.  Implementation of Concurrent Parallelization of Branch-and-bound algorithm in Everest Distributed Environment , 2017 .

[11]  Oleg Sukhoroslov,et al.  Implementation and Use of Coarse-grained Parallel Branch-and-bound in Everest Distributed Environment , 2017, ICCS.

[12]  G. Dantzig ON THE SIGNIFICANCE OF SOLVING LINEAR PROGRAMMING PROBLEMS WITH SOME INTEGER VARIABLES , 1960 .

[13]  Benjamin Müller,et al.  The SCIP Optimization Suite 5.0 , 2017, 2112.08872.