Recognition: unknown
Out-of-the-Box Global Optimization for Packing Problems: New Models and Improved Solutions
Pith reviewed 2026-05-08 15:56 UTC · model grok-4.3
The pith
Off-the-shelf global solvers find new packings for circles, polygons, and solids using direct nonlinear models.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The paper establishes that direct nonlinear formulations of packing problems, using the S-lemma for ellipse containment and Farkas-lemma-based non-overlap constraints for polygons and solids, enable off-the-shelf global optimization solvers to obtain numerous new incumbent solutions as well as first solutions for previously unstudied variants, without problem-specific modifications beyond writing down the models. Computational results also examine the impact of formulation choices, and the work supplies further evidence that global nonlinear optimization has matured into a practical technology for highly nonconvex problems.
What carries the argument
Nonlinear programming models of the packing constraints, with an S-lemma-based containment formulation for ellipses and compact Farkas-lemma-based non-overlap formulations for polygons and Platonic solids.
If this is right
- New best-known solutions for several circle-packing instances in squares and rectangles.
- First feasible solutions obtained for some previously unstudied packing variants of polygons and solids.
- Direct evidence that formulation variants affect solver success on these nonconvex problems.
- Support for treating global nonlinear optimization as a ready-to-use tool for geometric packing without custom algorithmic work.
Where Pith is reading between the lines
- The same modeling pattern could be tested on related problems such as packing with rotation or irregular containers.
- Success on these instances suggests global solvers may now scale to modestly larger packing tasks where only nonlinear models were previously considered intractable.
- If the Farkas-lemma non-overlap constraints generalize cleanly, they could simplify modeling for other rigid-body packing problems in three dimensions.
Load-bearing premise
The nonlinear models must accurately capture the geometric rules for non-overlap and containment, and the global solvers must be able to locate improved or optimal solutions in the resulting nonconvex problems.
What would settle it
For a small, well-studied packing instance whose optimal value is known independently, if the solver applied to one of these models returns either a feasible solution that violates non-overlap or a worse objective value than the known optimum, the central claim would be falsified.
read the original abstract
Recent LLM-driven discoveries have renewed interest in geometric packing problems. In this paper, we study several classes of such packing problems through the lens of modern global nonlinear optimization. Starting from comparatively direct nonlinear formulations, we consider packing circles in squares and fixed-perimeter rectangles, packing circles into minimum-area ellipses, packing regular polygons into regular polygons, and packing Platonic solids into Platonic solids. For ellipse packing, we derive a novel containment formulation based on the S-lemma. For polygon and Platonic solid packing, we develop compact non-overlap formulations based on the Farkas lemma and study several natural modeling variants computationally. Using off-the-shelf global optimization solvers, namely FICO Xpress and SCIP, we obtain numerous new incumbent solutions as well as first solutions for previously unstudied variants without any problem-specific solver modifications beyond writing down the models. Our computational results analyze the impact of various formulation choices. Beyond the individual packing results, the paper illustrates a broader point. It provides further evidence that global nonlinear optimization has matured into an increasingly practical model-and-solve technology for highly nonconvex problems.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops direct nonlinear programming formulations for several geometric packing problems: circles in squares and fixed-perimeter rectangles, circles in minimum-area ellipses (via a novel S-lemma containment model), and regular polygons/Platonic solids into like solids (via compact Farkas-lemma non-overlap constraints). Using unmodified off-the-shelf global solvers (FICO Xpress and SCIP), the authors report numerous new incumbent solutions and first solutions for previously unstudied variants, together with computational comparisons of formulation variants.
Significance. If the reformulations are shown to be exactly equivalent to the geometric constraints and the reported solutions are independently verified, the work supplies concrete evidence that modern global NLP solvers can be applied out-of-the-box to nonconvex packing problems to produce improved or first solutions. The explicit analysis of modeling variants and the absence of problem-specific solver tuning are strengths that support the broader claim about the maturing practicality of model-and-solve technology for such problems.
major comments (3)
- [ellipse-packing section] The S-lemma containment formulation for ellipse packing (described in the ellipse-packing section) is invoked to certify that one quadratic set lies inside another. The equivalence holds only when a strict-feasibility condition is satisfied (there exists a point at which the constraining quadratic is negative). The manuscript does not explicitly verify or prove that this condition holds for all solved instances; without it, the NLP may admit geometrically invalid points, undermining the claimed new solutions.
- [polygon and Platonic solid packing sections] The Farkas-lemma non-overlap formulations for polygons and Platonic solids (polygon-packing and Platonic-solid sections) introduce auxiliary variables for separation certificates. Equivalence to non-intersection requires that the original inequality systems satisfy the hypotheses of Farkas' lemma (e.g., no recession directions that render the alternative non-strict). The paper does not discuss or check these conditions for the considered shapes; violation could allow invalid configurations or exclude feasible ones.
- [computational results section] The computational results section reports new incumbents and first solutions but provides neither the complete model files, exact solver parameter settings, nor certificates (e.g., primal/dual bounds or geometric feasibility checks) that would allow independent verification that the NLP solutions satisfy the original non-overlap and containment constraints.
minor comments (2)
- [abstract and results summary] A summary table listing the number of new solutions per problem class, the previous best-known values, and the improvement margins would make the impact of the results clearer.
- [throughout] Notation for decision variables (centers, radii, orientations) is introduced separately for each problem class; a unified notation table would improve readability across sections.
Simulated Author's Rebuttal
We thank the referee for the careful review and constructive comments on the equivalence conditions of our formulations and on reproducibility. We address each major comment below.
read point-by-point responses
-
Referee: [ellipse-packing section] The S-lemma containment formulation for ellipse packing (described in the ellipse-packing section) is invoked to certify that one quadratic set lies inside another. The equivalence holds only when a strict-feasibility condition is satisfied (there exists a point at which the constraining quadratic is negative). The manuscript does not explicitly verify or prove that this condition holds for all solved instances; without it, the NLP may admit geometrically invalid points, undermining the claimed new solutions.
Authors: The referee correctly notes the strict-feasibility hypothesis required for equivalence in the S-lemma. In the ellipse-packing models, the inner quadratic inequality for each disk is strictly feasible at the disk center, where it evaluates to -r² < 0. All instances use non-degenerate ellipses of positive area containing circles of positive radius, so the condition holds for every solved instance. We will add an explicit remark in the ellipse-packing section confirming this fact for the computational study. This strengthens the justification but does not alter any reported solutions or models. revision: partial
-
Referee: [polygon and Platonic solid packing sections] The Farkas-lemma non-overlap formulations for polygons and Platonic solids (polygon-packing and Platonic-solid sections) introduce auxiliary variables for separation certificates. Equivalence to non-intersection requires that the original inequality systems satisfy the hypotheses of Farkas' lemma (e.g., no recession directions that render the alternative non-strict). The paper does not discuss or check these conditions for the considered shapes; violation could allow invalid configurations or exclude feasible ones.
Authors: We agree that Farkas' lemma requires the inequality systems to have no nonzero recession directions for full equivalence. The regular polygons and Platonic solids are compact convex bodies, so their inequality descriptions are bounded and possess no recession directions. This guarantees the required equivalence. We will insert a short clarifying discussion in the polygon-packing and Platonic-solid sections to state these facts and confirm they hold for the shapes studied. No changes to the computational results or models are needed. revision: partial
-
Referee: [computational results section] The computational results section reports new incumbents and first solutions but provides neither the complete model files, exact solver parameter settings, nor certificates (e.g., primal/dual bounds or geometric feasibility checks) that would allow independent verification that the NLP solutions satisfy the original non-overlap and containment constraints.
Authors: We accept that additional details improve verifiability. The revised manuscript will list the precise solver parameter settings used with FICO Xpress and SCIP. Complete model files for all instances will be supplied as supplementary material on the arXiv posting. The global solvers return primal-feasible points to the NLP; once the equivalence conditions (now to be stated explicitly) are satisfied, these points meet the original geometric constraints. We will also tabulate the final primal and dual bounds to support verification of optimality gaps. revision: yes
Circularity Check
No circularity: models derived from standard lemmas and solved with external solvers
full rationale
The derivation chain consists of applying known theorems (S-lemma for ellipse containment, Farkas lemma for polygon/Platonic non-overlap) to produce explicit nonlinear constraints, then feeding the resulting NLPs to independent commercial solvers (Xpress, SCIP). No parameter is fitted to the target packing outcomes, no result is renamed as a prediction, and no load-bearing premise reduces to a self-citation or self-definition. The reported incumbents are direct solver outputs on the externally formulated models.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math S-lemma
- standard math Farkas lemma
Reference graph
Works this paper leans on
-
[1]
Technical report, Optimization Online (2025)
Hojny, C., Besan¸ con, M., Bestuzheva, K., Borst, S., Chmiela, A., Dion´ ısio, J., Eifler, L., Ghannam, M., Gleixner, A., G¨ oß, A., Hoen, A., Hulst, R., Kamp, D., Koch, T., Kofler, K., Lentz, J., Maher, S.J., Mexi, G., M¨ uhmer, E., Pfetsch, M.E., Pokutta, S., Serrano, F., Shinano, Y., Turner, M., Vigerske, S., Wal- ter, M., Weninger, D., Xu, L.: The SCI...
2025
-
[2]
Optimization, 1–19 (2025) https: //doi.org/10.1080/02331934.2025.2595437
Belotti, P., Berthold, T., Gally, T., Gottwald, L., P´ olik, I.: Solving MINLPs to global optimality with FICO Xpress Global. Optimization, 1–19 (2025) https: //doi.org/10.1080/02331934.2025.2595437
-
[3]
PhD thesis, Technis- che Universit¨ at Berlin (2014)
Berthold, T.: Heuristic algorithms in global MINLP solvers. PhD thesis, Technis- che Universit¨ at Berlin (2014)
2014
-
[4]
Journal of Global Optimization70(1), 189–206 (2018) https://doi.org/10
Berthold, T.: A computational study of primal heuristics inside an MI(NL)P solver. Journal of Global Optimization70(1), 189–206 (2018) https://doi.org/10. 1007/s10898-017-0600-3
2018
-
[5]
Cambridge University Press, Cambridge, UK (2025)
Berthold, T., Lodi, A., Salvagnin, D.: Primal Heuristics in Integer Programming. Cambridge University Press, Cambridge, UK (2025)
2025
-
[6]
Georgiev, B., G´ omez-Serrano, J., Tao, T., Wagner, A.Z.: Mathematical explo- ration and discovery at scale. arXiv preprint arXiv:2511.02864 (2025)
-
[7]
AlphaEvolve: A coding agent for scientific and algorithmic discovery
Novikov, A., V˜ u, N., Eisenberger, M., Dupont, E., Huang, P.-S., Wagner, A.Z., Shirobokov, S., Kozlovskii, B., Ruiz, F.J., Mehrabian, A., et al.: Alphae- volve: A coding agent for scientific and algorithmic discovery. arXiv preprint arXiv:2506.13131 (2025)
work page internal anchor Pith review arXiv 2025
-
[8]
https: //github.com/algorithmicsuperintelligence/openevolve 28
Sharma, A.: OpenEvolve: an Open-source Evolutionary Coding Agent. https: //github.com/algorithmicsuperintelligence/openevolve 28
-
[9]
GitHub (2025)
Sharma, A.: Advanced circle packing for n=26 circles in a unit square. GitHub (2025). https://github.com/algorithmicsuperintelligence/openevolve/ blob/main/examples/circle packing/best program.py
2025
-
[10]
https://docs.scipy.org/doc/ scipy/reference/optimize.minimize-slsqp.html
SciPy API: Optimization and root finding (2025). https://docs.scipy.org/doc/ scipy/reference/optimize.minimize-slsqp.html
2025
-
[11]
Nature methods 17(3), 261–272 (2020)
Virtanen, P., Gommers, R., Oliphant, T.E., Haberland, M., Reddy, T., Cour- napeau, D., Burovski, E., Peterson, P., Weckesser, W., Bright, J.,et al.: SciPy 1.0: fundamental algorithms for scientific computing in Python. Nature methods 17(3), 261–272 (2020)
2020
-
[12]
https: //erich-friedman.github.io/packing/
Friedman, E.: Packing Center — Circle Packings and Related Problems. https: //erich-friedman.github.io/packing/. Online Resource
-
[13]
SIAM Review49(3), 371–418 (2007) https://doi.org/10.1137/S003614450444614X
P´ olik, I., Terlaky, T.: A survey of the S-lemma. SIAM Review49(3), 371–418 (2007) https://doi.org/10.1137/S003614450444614X
-
[14]
Journal f¨ ur die reine und angewandte Mathematik (Crelle’s Journal)124, 1–27 (1902)
Farkas, J.: Theorie der einfachen Ungleichungen. Journal f¨ ur die reine und angewandte Mathematik (Crelle’s Journal)124, 1–27 (1902)
1902
-
[15]
Berthold, T., Kamp, D., Mexi, G., Pokutta, S., P´ olik, I.: Global optimization for combinatorial geometry problems. arXiv preprint arXiv:2601.05943. Accepted for publication in ISCO2026 Proceedings. (2026) https://doi.org/10.48550/arXiv. 2601.05943
work page internal anchor Pith review doi:10.48550/arxiv 2026
-
[16]
arxiv preprint arxiv:2603.11107 (2026)
Sudermann-Merx, N.: From computational certification to exact coordinates: Heilbronn’s triangle problem on the unit square using mixed-integer optimization. arxiv preprint arxiv:2603.11107 (2026)
-
[17]
Annals of Mathematics162(3), 1065–1185 (2005) https://doi.org/10.4007/annals.2005.162.1065
Hales, T.C.: A proof of the Kepler conjecture. Annals of Mathematics162(3), 1065–1185 (2005) https://doi.org/10.4007/annals.2005.162.1065
-
[18]
Typis Collegii Reformatorum, Marosv´ as´ arhely (1832)
Bolyai, W.: Tentamen Juventutem Studiosam in Elementa Matheseos Purae, Ele- mentaris Ac Sublimioris. Typis Collegii Reformatorum, Marosv´ as´ arhely (1832). Two volumes
-
[19]
Pergamon Press, London, UK (1964)
Fejes T´ oth, L.: Regular Figures. Pergamon Press, London, UK (1964)
1964
-
[20]
European Journal of Operational Research183(3), 1109–1130 (2007)
W¨ ascher, G., Haußner, H., Schumann, H.: An improved typology of cutting and packing problems. European Journal of Operational Research183(3), 1109–1130 (2007)
2007
-
[21]
European Journal of Operational Research183(3), 1131–1135 (2007)
Fekete, S.P., Veen, J.C.: PackLib2: An integrated library of multi-dimensional packing problems. European Journal of Operational Research183(3), 1131–1135 (2007)
2007
-
[22]
Optimization Letters16(2), 471–480 (2022)
Iori, M., Lima, V.L., Martello, S., Monaci, M.: 2DPackLib: a two-dimensional 29 cutting and packing library. Optimization Letters16(2), 471–480 (2022)
2022
-
[23]
PhD thesis, Georgia Institute of Technology (2024)
Kuznetsov, A.: Convexification and global optimization of problems involving the Euclidean norm. PhD thesis, Georgia Institute of Technology (2024)
2024
-
[24]
http://www.gamsworld.org/minlp/minlplib.htm
MINLP Library. http://www.gamsworld.org/minlp/minlplib.htm
-
[25]
SIAM Review56(1), 3–69 (2014) https://doi.org/10.1137/ 120875909
Liberti, L., Lavor, C., Maculan, N., Mucherino, A.: Euclidean distance geome- try and applications. SIAM Review56(1), 3–69 (2014) https://doi.org/10.1137/ 120875909
2014
-
[26]
Discrete Applied Mathematics122(1-3), 139–166 (2002)
Locatelli, M., Raber, U.: Packing equal circles in a square: a deterministic global optimization approach. Discrete Applied Mathematics122(1-3), 139–166 (2002)
2002
-
[27]
European Journal of Operational Research191(3), 786–802 (2008)
Castillo, I., Kampas, F.J., Pint´ er, J.D.: Solving circle packing problems by global optimization: numerical results and industrial applications. European Journal of Operational Research191(3), 786–802 (2008)
2008
-
[28]
https://plus.maths.org/content/ prize-specimens
Wainwright, M.: Prize specimens (2001). https://plus.maths.org/content/ prize-specimens
2001
-
[29]
Operations Research Letters57(2024)
Khajavirad, A.: The circle packing problem: A theoretical comparison of various convexification techniques. Operations Research Letters57(2024)
2024
-
[30]
Blog post (2025)
Berthold, T.: FICO Xpress Optimization Surpasses AlphaEvolve’s Achievements. Blog post (2025). https://www.fico.com/blogs/best-global-optimization-solver
2025
-
[31]
SIAM Journal on Computing3(4), 262–279 (1974)
Sahni, S.: Computationally related problems. SIAM Journal on Computing3(4), 262–279 (1974)
1974
-
[32]
Soviet Mathematics
Matijaseviˇ c, J.V.: Enumerable sets are Diophantine. Soviet Mathematics. Dok- lady11(2), 354–358 (1970)
1970
-
[33]
Mathematical P rogramming 10, 147–175 (1976) https://doi.org/10.1007/BF01580665
McCormick, G.P.: Computability of global solutions to factorable nonconvex pro- grams: Part I — Convex underestimating problems. Mathematical Programming B10(1), 147–175 (1976) https://doi.org/10.1007/BF01580665
-
[34]
SIAM Journal on Discrete Mathematics3(3), 411–430 (1990)
Sherali, H.D., Adams, W.P.: A hierarchy of relaxations between the continu- ous and convex hull representations for zero-one programming problems. SIAM Journal on Discrete Mathematics3(3), 411–430 (1990)
1990
-
[35]
Journal of Global Optimization43(2-3), 471–484 (2009) https://doi.org/10.1007/ s10898-008-9372-0
Anstreicher, K.M.: Semidefinite programming versus the reformulation-lineariza- tion technique for nonconvex quadratically constrained quadratic programming. Journal of Global Optimization43(2-3), 471–484 (2009) https://doi.org/10.1007/ s10898-008-9372-0
2009
-
[36]
Optimization Methods & Software 30 24, 597–634 (2009)
Belotti, P., Lee, J., Liberti, L., Margot, F., W¨ achter, A.: Branching and bounds tightening techniques for non-convex MINLP. Optimization Methods & Software 30 24, 597–634 (2009)
2009
-
[37]
Journal of Global Optimizatio n 67, 731– 757 (2017) https://doi.org/10.1007/s10898-016-0450-4
Gleixner, A., Berthold, T., M¨ uller, B., Weltge, S.: Three enhancements for optimization-based bound tightening. Journal of Global Optimization67, 731– 757 (2017) https://doi.org/10.1007/s10898-016-0450-4
-
[38]
PhD thesis, Humboldt-Universit¨ at zu Berlin (2012)
Vigerske, S.: Decomposition in multistage stochastic programming and a con- straint integer programming approach to mixed-integer nonlinear programming. PhD thesis, Humboldt-Universit¨ at zu Berlin (2012)
2012
-
[39]
arXiv preprint arXiv:2602.09996 (2026)
Berthold, T., Geis, F.: Learning to choose branching rules for nonconvex MINLPs. arXiv preprint arXiv:2602.09996 (2026)
-
[40]
INFORMS Journal on Computing33, 421–435 (2021) https://doi.org/10.1287/ijoc.2020.1050
Berthold, T., Witzig, J.: Conflict analysis for MINLP. INFORMS Journal on Computing33, 421–435 (2021) https://doi.org/10.1287/ijoc.2020.1050
-
[41]
Berthold, T., Gleixner, A.M., Heinz, S., Vigerske, S.: Analyzing the computa- tional impact of MIQCP solver components. Numerical Algebra, Control and Optimization2(4), 739–748 (2012) https://doi.org/10.3934/naco.2012.2.739
-
[42]
Nonconvex Optimization and Its Applications, vol
Tawarmalani, M., Sahinidis, N.V.: Convexification and Global Optimization in Continuous and Mixed-integer Nonlinear Programming: Theory, Algorithms, Software, and Applications. Nonconvex Optimization and Its Applications, vol
-
[43]
Kluwer Academic Publishers, The Netherlands (2002)
2002
-
[44]
Mathematical Programming 99, 563–591 (2004)
Tawarmalani, M., Sahinidis, N.V.: Global optimization of mixed-integer nonlinear programs: A theoretical and computational study. Mathematical Programming 99, 563–591 (2004)
2004
-
[45]
Journal of Global Optimization59(2–3), 503–526 (2014) https://doi.org/10.1007/s10898-014-0166-2
Misener, R., Floudas, C.A.: ANTIGONE: Algorithms for coNTinuous / Integer Global Optimization of Nonlinear Equations. Journal of Global Optimization 59(2-3), 503–526 (2014) https://doi.org/10.1007/s10898-014-0166-2
-
[46]
Mathemati- cal Programming Computation1(1), 1–41 (2009) https://doi.org/10.1007/ s12532-008-0001-1
Achterberg, T.: SCIP: Solving Constraint Integer Programs. Mathemati- cal Programming Computation1(1), 1–41 (2009) https://doi.org/10.1007/ s12532-008-0001-1
2009
-
[47]
In: System Modelling and Optimization: Proceedings of the 15th IFIP Conference Zurich, Switzerland, September 2–6, 1991, pp
Peikert, R., W¨ urtz, D., Monagan, M., Groot, C.: Packing circles in a square: a review and new results. In: System Modelling and Optimization: Proceedings of the 15th IFIP Conference Zurich, Switzerland, September 2–6, 1991, pp. 45–54 (2007). Springer
1991
-
[48]
ACM SIGSAM Bulletin10(3), 19–29 (1976)
Buchberger, B.: A theoretical basis for the reduction of polynomials to canonical forms. ACM SIGSAM Bulletin10(3), 19–29 (1976)
1976
-
[49]
https://erich-friedman.github.io/packing/ cirRsqu/ Accessed 2025-12-15 31
Friedman, E.: Circles in Squares. https://erich-friedman.github.io/packing/ cirRsqu/ Accessed 2025-12-15 31
2025
-
[50]
Vestnik Leningrad University1, 62–77 (1971)
Yakubovich, V.A.: S-procedure in nonlinear control theory. Vestnik Leningrad University1, 62–77 (1971). In Russian
1971
-
[51]
Mathematical Programming Computation 11, 37–93 (2019)
Pfetsch, M.E., Rehn, T.: A computational comparison of symmetry handling methods for mixed integer programs. Mathematical Programming Computation 11, 37–93 (2019)
2019
-
[52]
Journal of Global Optimization43(2), 299–328 (2009)
Kallrath, J.: Cutting circles and polygons from area-minimizing rectangles. Journal of Global Optimization43(2), 299–328 (2009)
2009
-
[53]
Journal of Global Optimization65(2), 283–307 (2016)
Stoyan, Y., Pankratov, A., Romanova, T.: Quasi-phi-functions and optimal packing of ellipses. Journal of Global Optimization65(2), 283–307 (2016)
2016
-
[54]
European Journal of Operational Research268(1), 37–53 (2018)
Romanova, T., Bennell, J., Stoyan, Y., Pankratov, A.: Packing of concave poly- hedra with continuous rotations using nonlinear optimisation. European Journal of Operational Research268(1), 37–53 (2018)
2018
-
[55]
Bertsimas, D., Tsitsiklis, J.N.: Introduction to Linear Optimization vol. 6. Athena scientific, Belmont, MA (1997)
1997
-
[56]
https://erich-friedman.github.io/packing/ hexinhex/ Accessed 2025-12-15
Friedman, E.: Hexagons in Hexagons. https://erich-friedman.github.io/packing/ hexinhex/ Accessed 2025-12-15
2025
-
[57]
ORMS Today53(1) (2026) https://doi.org/10.1287/orms.2026.01
Berthold, T., P´ olik, I.: What GenAI taught us about the future of global optimization. ORMS Today53(1) (2026) https://doi.org/10.1287/orms.2026.01. 01
-
[58]
Diploma thesis, Technische Universit¨ at Berlin (2006)
Berthold, T.: Primal Heuristics for Mixed Integer Programs. Diploma thesis, Technische Universit¨ at Berlin (2006)
2006
-
[59]
Technical report (2025) 32
Monji, A., Modir, A., Kocuk, B.: Solving the Heilbronn triangle problem using global optimization methods. Technical report (2025) 32
2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.