Continuous Optimization for Satisfiability Modulo Theories on Linear Real Arithmetic

arXiv cs.AI / 3/25/2026

💬 OpinionIdeas & Deep AnalysisTools & Practical UsageModels & Research

Key Points

  • The paper proposes FourierSMT, a continuous, highly parallelizable optimization-based approach to satisfiability modulo theories (SMT) that is easier to scale than conflict-driven clause learning in parallel settings.
  • It extends the Walsh-Fourier expansion (WFE) to a mixed Boolean-real domain via an extended WFE (xWFE), enabling gradient methods and local updates for high-arity SMT constraints.
  • To make xWFE efficient, the authors introduce an extended binary decision diagram (xBDD) and reformulate xWFE constraint evaluation using randomized rounding and the circuit-output probability (COP), which matches the expected value of xWFE.
  • The reduced optimization problem is proven to converge while preserving satisfiability, supporting soundness of the computed solutions.
  • Experiments on large scheduling and placement instances (up to 10,000 variables and 700,000 constraints) report up to 8× speedups over state-of-the-art SMT solvers and suggest GPU-based optimization for continuous systems.

Abstract

Efficient solutions for satisfiability modulo theories (SMT) are integral in industrial applications such as hardware verification and design automation. Existing approaches are predominantly based on conflict-driven clause learning, which is structurally difficult to parallelize and therefore scales poorly. In this work, we introduce FourierSMT as a scalable and highly parallelizable continuous-variable optimization framework for SMT. We generalize the Walsh-Fourier expansion (WFE), called extended WFE (xWFE), from the Boolean domain to a mixed Boolean-real domain, which allows the use of gradient methods for SMT. This addresses the challenge of finding satisfying variable assignments to high-arity constraints by local updates of discrete variables. To reduce the evaluation complexity of xWFE, we present the extended binary decision diagram (xBDD) and map the constraints from xWFE to xBDDs. We then show that sampling the circuit-output probability (COP) of xBDDs under randomized rounding is equivalent to the expectation value of the xWFEs. This allows for efficient computation of the constraints. We show that the reduced problem is guaranteed to converge and preserves satisfiability, ensuring the soundness of the solutions. The framework is benchmarked for large-scale scheduling and placement problems with up to 10,000 variables and 700,000 constraints, achieving 8-fold speedups compared to state-of-the-art SMT solvers. These results pave the way for GPU-based optimization of SMTs with continuous systems.