REVIEW 3 major objections 4 minor 25 references
Smooth Logic Constraints in Nonlinear Optimization and Optimal Control Problems
T0 review · 3 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read A single reformulation turns arbitrary logic constraints into smooth math, removing the need for binary variables and mixed-integer solvers.
desk verdict The central exactness claim is false: Lemma 1's epsilon is ill-defined and the practical epsilon=0 changes the feasible set; the numerical examples also implement different logic than stated. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is the max-min reading of logic: an OR of propositions $p_j \le 0$ is equivalent to $\min_j p_j \le 0$, and an AND of clauses is equivalent to $\max_i \left(\min_j p_{i,j}\right) \le 0$. The paper then applies a known smoothing result to replace each inner min over a finite set by a weighted sum $\sum_j \lambda_j p_j$ with $\lambda$ on the standard simplex, making the whole logic constraint differentiable while preserving the solution set under the stated regularity assumption. Equality propositions are first converted into paired inequalities, and negations are pushed to the proposition level and removed by flipping the sign of the function, so the logic program reduces to conjunctive normal form in inequality propositions only.
What would settle it
Solve a one-variable problem whose only logic constraint is NOT(p <= 0) with p(x) = x, using the paper's replacement p_bar = -p. At x = 0 the original constraint is false, because p=0 satisfies p<=0 and its negation is false, while the smoothed constraint -0 <= 0 is true, so the smoothed problem accepts x=0 as feasible. If a run of the method returns such a boundary point that the original logic rejects, the claimed exactness of the solution set fails in that instance.
Extended reading notes
Core claim
At the center is Theorem 1: for any logic program $L$ satisfying the paper's assumptions, the logical constraint $L(x,u)=true$ can be replaced by the family of smooth constraints $\sum_j \lambda_{i,j} p_{i,j}(x,u) \le 0$ for each conjunctive clause $i$, with each $\lambda$ vector on the simplex, without affecting the solution set. In plain terms, the paper claims an exact bridge from discrete logic to continuous differentiability: no binary variables, no relaxation penalty, and no loss of feasible solutions relative to the original logic-constrained problem. The price is the regularity assumption that the functions defining the propositions are twice differentiable and satisfy a non-degeneracy condition when gradients are tied together by a min operation.
Load-bearing premise
The load-bearing premise is that a strict inequality can be converted into a non-strict one by subtracting an arbitrarily small positive epsilon, but since no real epsilon is smaller than every positive delta, the negation step is exact only in the limit and can be wrong exactly on constraint boundaries.
Editorial extensions
If this is right
- Any optimal control problem whose logic can be written with equality and inequality propositions plus NOT, AND, and OR can be handed to a standard continuous solver without binary variables.
- For the two quadrotor benchmarks, the smoothed formulation finds the optimum in 81.3% of runs versus 66.0% for the best baseline (big-M), with lower average computation time and fewer infeasible runs.
- Because the reformulation preserves the solution set exactly under the stated assumptions, feasible trajectories found by the smoothed problem satisfy the original logical specification, not merely an approximation.
- The conversion to conjunctive normal form is not mandatory; direct min and max conversion may be numerically better in some problems, so the theorem provides a family of equivalent smooth formulations rather than a single fixed encoding.
Reading between the lines
- Boundary cases need care beyond the paper's proof: at points where a negated proposition function vanishes, flipping the sign can make both the proposition and its negation appear true, so a practical implementation should treat boundary satisfaction as approximate or add a safety margin.
- The smoothing variables are simplex weights, so their initialization may influence which local solution a nonconvex solver reaches; varying these initializations could be a practical tuning knob that the paper does not explore.
- The same reduction applies to any mathematical program with logic constraints, not only dynamic optimal control, and could be tested on general mixed-integer nonlinear benchmark problems to compare against binary encodings.
- The choice of representation, such as conjunctive normal form versus direct max-min conversion, changes the number of auxiliary variables and the curvature of the constraints; finding an optimal representation is left open and could be approached heuristically or empirically.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a method for reformulating logic constraints in optimal control problems as smooth continuous constraints, avoiding binary variables or mixed-integer solvers. The approach first rewrites any logic program built from equality and inequality propositions with NOT, AND, and OR into conjunctive normal form, then converts the resulting AND/OR structure into max-min constraints, and finally smooths the max-min constraints using auxiliary simplex variables, citing a result of Kirjner-Neto and Polak. The central claim, stated in Theorem 1 and Section III-D, is that the reformulation has an identical solution set to the original logic-constrained problem. The method is demonstrated on two planar quadrotor optimal control problems with conditional logic constraints, with numerical comparisons against complementarity and Big-M formulations.
Significance. If the central exactness claim were correct, the paper would offer a practically valuable bridge between discrete logic and continuous, derivative-based optimization, and the numerical comparisons on the quadrotor examples are informative. The paper is clearly written and the numerical setup is reproducible in principle, with randomized restarts and timing statistics. However, the load-bearing theoretical claim, namely that the smooth reformulation preserves the solution set exactly, rests on the negation-elimination step of Lemma 1, and that step is not mathematically sound. Because the exactness claim is the stated contribution, the flaw affects the main result of the paper.
major comments (3)
- [III-A, Lemma 1] The elimination of negated inequalities via Eq. (10b) requires an epsilon with 0 < epsilon < delta for every delta > 0; no real number has this property. The paper then states that in practice setting ̄p_i = -p_i is sufficient, but with epsilon = 0 the replacement turns ¬(p ≤ 0) into p ≥ 0. At p = 0, both p ≤ 0 and its replacement hold, so the claimed logical equivalence fails on constraint boundaries. This is not a harmless measure-zero issue: for the problem minimize x subject to ¬(x ≤ 0) and x ∈ [-1,1], the original feasible set is (0,1] and no minimum is attained, while the smoothed replacement p ≥ 0 admits x = 0, a point that violates the original logic. Since every negated inequality in the CNF of Section III-B passes through this step, the error propagates into Theorem 1 and the 'identical solution set' claim of Section III-D.
- [III-C, Theorem 1] The proof of Theorem 1 is deferred to [19] and does not address the negation-elimination step that precedes Lemma 2. Even accepting the max-min-to-simplex conversion from [19], that conversion can only preserve the solution set of the max-min constraint (12); it cannot repair the fact that the max-min constraint itself is derived using Lemma 1's inexact treatment of negated inequalities. The manuscript also asserts Assumption 3.9 of [19] without checking it for the two quadrotor examples, so the numerical demonstrations do not validate the hypotheses under which the external theorem is invoked.
- [III-A, paragraph after Eq. (10)] The statement that solvers cannot distinguish strict and non-strict inequalities explains why an epsilon-perturbation might be useful numerically, but it does not justify the exactness claim. The paper's contribution is explicitly stated as replacing the logic constraint 'without affecting the solution set' and producing an 'identical solution set' in Section III-D. As the counterexample above shows, boundary points where a proposition and its negation overlap can be selected by a continuous optimizer, so the practical recipe with epsilon = 0 changes the feasible set. A correct treatment would need to handle strict inequalities explicitly, which standard NLP solvers cannot do, or would require weakening the exactness claim substantially.
minor comments (4)
- [IV-C, first paragraph] The text says the optimal solution to Problem 1 is 'as seen in Figure 3 (c)', but the corresponding trajectory is shown in Figure 2 (c); this cross-reference should be corrected.
- [Eq. (21)] The displayed smoothed constraint for Problem 1 is ambiguous: it is not clear whether a separate constraint is generated for each i ∈ {5,...,9} or whether the expression is a single constraint, and the λ variables are not indexed by i even though they should belong to different simplex sets for different i.
- [II, Eq. (3)] The domain notation f : R^{n/(N+1)} × R^{m/N} → R^{n/(N+1)} is dimensionally opaque; the standard notation f : R^n × R^m → R^n with time indices would be clearer.
- [Abstract and Introduction] The paper repeatedly claims to handle 'arbitrary logic constraints', but Assumption 2 restricts to twice-differentiable constituent functions with a non-degeneracy condition. This restriction should be stated in the abstract and introduction so that the scope of the claim is accurate.
Circularity Check
No significant circularity: the reformulation chain rests on external standard results, not on self-citation or fitted inputs.
full rationale
The derivation chain is: Lemma 1 reduces equality propositions to inequality conjunctions using explicit algebraic constructions (4)-(5), then eliminates negations via De Morgan equivalences (6)-(8) and the replacement (10). Lemma 2 converts AND/OR to max/min using the standard logical identities quoted from [22] and CNF from [21]. Theorem 1 then invokes an external smoothing theorem, Kirjner-Neto and Polak [19], to replace each finite max-min constraint by smooth constraints (13); the proof of Theorem 1 is explicitly 'apply the smoothing approach of [19]', which is not authored by the present paper and is not itself derived from the paper's claims. Assumption 2 likewise defers to Assumption 3.9 in [19]. No parameter is fitted to data and no quantity that is predicted is an input by construction. The self-citations [20] and [23] are provenance and illustration only: [20] states the method 'was previously utilized in our work', and [23] is only the source of the quadrotor figure. Consequently there is no exhibited reduction of a claimed theorem to its own assumptions. The one genuinely questionable step is Lemma 1's epsilon condition (0<epsilon<delta for every delta>0 in R), which no real epsilon satisfies, and the practical choice epsilon=0 makes the replacement of negated inequalities approximate on the boundary p=0. That is a correctness/technical flaw in the claimed exactness, not a circular definition: it is not defined in terms of the target solution set, fitted to outputs, or justified exclusively by a self-citation. Under the hard rule that circularity must be exhibited as Eq. X = Eq. Y by construction or a self-citation chain, that flaw is outside the circularity count. The central result therefore has independent content and is benchmarked against external methods [24,25]; 0 is the appropriate score.
Assumptions & free parameters
free parameters (1)
- epsilon =
0 (practice); no real value satisfies 0<epsilon<delta for all delta>0
assumptions (6)
- standard math Functional completeness of NOT, AND, OR for finite propositional logic
- standard math Every propositional formula can be put in conjunctive normal form without changing truth
- standard math For real numbers, AND of inequalities is max and OR of inequalities is min
- domain assumption Assumption 1: logic consists of finitely many equality/inequality propositions linked by NOT, AND, OR
- domain assumption Assumption 2: functions twice differentiable and gradients satisfy Assumption 3.9 of [19]
- ad hoc to paper Solvers cannot distinguish strict and non-strict inequalities, so epsilon can be dropped
Cite this review
Pith. "Pith review of Smooth Logic Constraints in Nonlinear Optimization and Optimal Control Problems." pith.science (2026). https://pith.science/paper/QXTG6D56
@misc{pith2026250601742,
author = {Pith},
title = {Pith review of: Smooth Logic Constraints in Nonlinear Optimization and Optimal Control Problems},
year = {2026},
howpublished = {\url{https://pith.science/paper/QXTG6D56}},
note = {Machine review of arXiv:2506.01742}
}
read the original abstract
In some optimal control problems, complex relationships between states and inputs cannot be easily represented using continuous constraints, necessitating the use of discrete logic instead. This paper presents a method for incorporating such logic constraints directly within continuous optimization frameworks, eliminating the need for binary variables or specialized solvers. Our approach reformulates arbitrary logic constraints under minimal assumptions as max-min constraints, which are then smoothed by introducing auxiliary variables into the optimization problem. When these reformulated constraints are satisfied, they guarantee that the original logical conditions hold, ensuring correctness in the optimization process. We demonstrate the effectiveness of this method on two planar quadrotor control tasks with complex logic constraints. Compared to existing techniques for encoding logic in continuous optimization, our approach achieves faster computational performance and improved convergence to feasible solutions.
Figures
Reference graph
Works this paper leans on
-
[19]
C. Kirjner-Neto and E. Polak, “On the conversion of optimization problems with max-min constraints to standard optimization prob- lems,” SIAM Journal on Optimization , vol. 8, no. 4, pp. 887–915, 1998
work page 1998
-
[1]
R. B. Vinter, Optimal control. Springer, 2010, vol. 2, no. 1
work page 2010
-
[2]
J. Nocedal and S. J. Wright, Numerical optimization. Springer, 2006. (a) Going through the green circle allows passing through the red obstacle. (b) Sample sub-optimal trajectory. (c) Optimal solution trajectory. Fig. 2. Obstacle (red), logic trigger (green), target (black), and solution trajectory (blue) for Problem 1. (a) Going through the green circle ...
work page 2006
-
[3]
Optimization-based control for dynamic legged robots,
P. M. Wensing, M. Posa, Y . Hu, A. Escande, N. Mansard, and A. Del Prete, “Optimization-based control for dynamic legged robots,” IEEE Transactions on Robotics , vol. 40, pp. 43–63, 2023
2023
-
[4]
Robust optimal control for demand side management of multi- carrier microgrids,
R. Carli, G. Cavone, T. Pippia, B. De Schutter, and M. Dotoli, “Robust optimal control for demand side management of multi- carrier microgrids,” IEEE Transactions on Automation Science and Engineering, vol. 19, no. 3, pp. 1338–1351, 2022
work page 2022
-
[5]
Fast homotopy for spacecraft ren- dezvous trajectory optimization with discrete logic,
D. Malyuta and B. Ac ¸ıkmes ¸e, “Fast homotopy for spacecraft ren- dezvous trajectory optimization with discrete logic,” Journal of Guid- ance, Control, and Dynamics , vol. 46, no. 7, pp. 1262–1279, 2023
work page 2023
-
[6]
A review and comparison of solvers for convex MINLP,
J. Kronqvist, D. E. Bernal, A. Lundell, and I. E. Grossmann, “A review and comparison of solvers for convex MINLP,” Optimization and Engineering, vol. 20, pp. 397–455, 2019
work page 2019
-
[7]
Undecidability and hardness in mixed-integer nonlinear programming,
L. Liberti, “Undecidability and hardness in mixed-integer nonlinear programming,” RAIRO-Operations Research, vol. 53, no. 1, pp. 81– 109, 2019
work page 2019
Show all 25 references
-
[8]
Approximation properties and tight bounds for constrained mixed-integer optimal control,
C. Kirches, F. Lenders, and P. Manns, “Approximation properties and tight bounds for constrained mixed-integer optimal control,” SIAM Journal on Control and Optimization , vol. 58, no. 3, pp. 1371–1402, 2020
2020
-
[9]
Non-convex mixed-integer nonlinear programming: A survey,
S. Burer and A. N. Letchford, “Non-convex mixed-integer nonlinear programming: A survey,” Surveys in Operations Research and Man- agement Science, vol. 17, no. 2, pp. 97–106, 2012
2012
-
[10]
A unified approach to mixed-integer optimization problems with logical constraints,
D. Bertsimas, R. Cory-Wright, and J. Pauphilet, “A unified approach to mixed-integer optimization problems with logical constraints,” SIAM Journal on Optimization , vol. 31, no. 3, pp. 2340–2367, 2021
2021
-
[11]
M. S. Bazaraa, J. J. Jarvis, and H. D. Sherali, Linear programming and network flows . John Wiley & Sons, 2011
2011
-
[12]
On handling indicator constraints in mixed integer programming,
P. Belotti, P. Bonami, M. Fischetti, A. Lodi, M. Monaci, A. Nogales- G´omez, and D. Salvagnin, “On handling indicator constraints in mixed integer programming,” Computational Optimization and Applications , vol. 65, pp. 545–566, 2016
2016
-
[13]
S. P. Boyd and L. Vandenberghe, Convex optimization. Cambridge university press, 2004
2004
-
[14]
On the approximation of the step function by some sigmoid functions,
A. Iliev, N. Kyurkchiev, and S. Markov, “On the approximation of the step function by some sigmoid functions,” Mathematics and Computers in Simulation , vol. 133, pp. 223–234, 2017
2017
-
[15]
Perspective reformulations of mixed integer nonlinear programs with indicator variables,
O. G ¨unl¨uk and J. Linderoth, “Perspective reformulations of mixed integer nonlinear programs with indicator variables,” Mathematical Programming, vol. 124, pp. 183–205, 2010
2010
-
[16]
Optimization-based control of nonlinear systems with linear temporal logic specifications,
E. M. Wolff, U. Topcu, and R. M. Murray, “Optimization-based control of nonlinear systems with linear temporal logic specifications,” in Proc. of Int. Conf. on Robotics and Automation , 2014, pp. 5319–5325
2014
-
[17]
Formal methods for control synthesis: An optimization perspective,
C. Belta and S. Sadraddini, “Formal methods for control synthesis: An optimization perspective,” Annual Review of Control, Robotics, and Autonomous Systems , vol. 2, no. 1, pp. 115–140, 2019
2019
-
[18]
The continuous quadrant penalty formulation of logical constraints,
S. Cafieri, A. R. Conn, and M. Mongeau, “The continuous quadrant penalty formulation of logical constraints,” Open Journal of Mathe- matical Optimization, pp. 1–12, 2023
2023
-
[20]
Semi-infinite programs for robust con- trol and optimization: Efficient solutions and extensions to existence constraints,
J. Wehbeh and E. C. Kerrigan, “Semi-infinite programs for robust con- trol and optimization: Efficient solutions and extensions to existence constraints,” in 8th IFAC Conference on Nonlinear Model Predictive Control NMPC 2024 . IFAC, 2024, pp. 317–322
2024
-
[21]
H. B. Enderton, A mathematical introduction to logic. Elsevier, 2001
2001
-
[22]
L. A. Wolsey and G. L. Nemhauser, Integer and combinatorial optimization. John Wiley & Sons, 1999
1999
-
[23]
Robust output feedback of nonlin- ear systems through the efficient solution of min-max optimization problems,
J. Wehbeh and E. C. Kerrigan, “Robust output feedback of nonlin- ear systems through the efficient solution of min-max optimization problems,” in 2024 IEEE 63rd Conference on Decision and Control (CDC). IEEE, 2024, pp. 8870–8875
2024
-
[24]
JuMP 1.0: Recent improvements to a modeling language for mathematical optimization,
M. Lubin, O. Dowson, J. Dias Garcia, J. Huchette, B. Legat, and J. P. Vielma, “JuMP 1.0: Recent improvements to a modeling language for mathematical optimization,” Mathematical Programming Compu- tation, 2023
2023
-
[25]
On the implementation of an interior- point filter line-search algorithm for large-scale nonlinear program- ming,
A. W ¨achter and L. T. Biegler, “On the implementation of an interior- point filter line-search algorithm for large-scale nonlinear program- ming,” Mathematical Programming, vol. 106, pp. 25–57, 2006
2006
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.