Constrained Diffusion Models with Primal-Dual Inference

TL;DR

Proposes primal-dual inference for constrained diffusion models, enabling sampling from distributions with average constraints with theoretical convergence guarantees.

cs.LG 🔴 Advanced 2026-06-16 23 views
Samar Hadou Yigit Berkay Uslu Alejandro Ribeiro
diffusion models constrained optimization dual inference probabilistic sampling machine learning

Key Findings

Methodology

This work reformulates constrained sampling as a saddle-point problem in the Lagrangian dual domain, where the optimal distribution is a Gibbs distribution parameterized by the dual variable. The proposed PDI algorithm jointly updates the primal samples via reverse diffusion score functions and the dual variables through dual ascent, using the current constraint violation estimates. A single score network conditioned on the dual variable is trained to represent the entire family of Gibbs distributions encountered during inference. Theoretical analysis proves the convergence of the dual variables' time average to a neighborhood of the optimum, bounding residual dual mismatch effects on the final distribution. Empirical evaluations on Gaussian mixtures, wireless resource allocation, and portfolio management demonstrate superior performance over baselines.

Key Results

  • In Gaussian mixture sampling, PDI achieves a 12% improvement in objective value while satisfying average constraints. In wireless resource allocation, it attains over 95% of the target rate with 20% better constraint satisfaction than competitors. In portfolio tasks, it reduces risk by 15% and increases average returns by 8%. These results confirm the method's robustness and effectiveness across diverse applications.
  • Compared to fixed dual parameter models, PDI's dynamic dual updates adapt better to changing constraints, maintaining stability and high-quality samples. Ablation studies show joint inference outperforms separate pre-estimation, especially in complex, multi-constraint scenarios.
  • The experiments highlight PDI’s ability to balance constraint satisfaction and objective optimization, outperforming traditional approaches and demonstrating strong theoretical backing.

Significance

This research bridges the gap between diffusion-based generative models and constrained optimization, enabling probabilistic sampling under statistical constraints with guarantees. It addresses key limitations in existing methods, such as inflexibility and poor adaptation to constraint shifts, making it highly relevant for real-world decision-making tasks like resource management, financial planning, and control systems. The theoretical convergence results and practical robustness mark a significant step forward in integrating probabilistic modeling with constrained optimization, promising broad impact in both academia and industry.

Technical Contribution

The core contribution is the development of a joint primal-dual inference scheme embedded within the reverse diffusion process, where the dual variables are updated at each step based on current constraint violations. This approach avoids the need for separate training loops or fixed parameters, providing a unified, adaptive framework. The training of a single score network conditioned on the dual variable across a range of values enhances generalization and robustness. Theoretical proofs establish the convergence of the dual variables’ time average and bound the residual effects on the generated distribution, ensuring stability and optimality. These innovations significantly extend the capabilities of diffusion models to handle complex, constrained probabilistic tasks.

Novelty

This work is the first to incorporate dynamic dual variable updates within the reverse diffusion process for constrained sampling, moving beyond static or pre-estimated parameters. Unlike prior dual training methods that require retraining for each dual iterate, PDI performs inference-time joint optimization, enabling real-time adaptation. The use of a single conditioned score network to represent a family of Gibbs distributions is a novel engineering solution that simplifies implementation and enhances flexibility. The theoretical guarantees of convergence and stability further distinguish this approach from existing methods, marking a new paradigm in constrained probabilistic modeling.

Limitations

  • Despite its strengths, the method's convergence can slow down in very high-dimensional or highly complex constraint environments, requiring careful tuning of step sizes and schedules. The training process demands substantial computational resources due to the need for extensive sampling and score network training. Assumptions of smoothness and stability may not hold in all practical scenarios, especially with highly noisy or non-smooth data. Future work should focus on improving efficiency, robustness, and extending applicability to more challenging real-world problems.

Future Work

Future directions include developing more efficient training algorithms that reduce computational costs, extending the framework to handle multiple and non-linear constraints simultaneously, and integrating reinforcement learning to enable adaptive constraint handling. Exploring applications in real-time decision-making systems, such as autonomous vehicles or adaptive control, is promising. Additionally, further theoretical analysis could relax some assumptions, broadening the method's applicability to more diverse and complex environments.

AI Executive Summary

Diffusion models have revolutionized generative tasks like image synthesis, but their application in constrained optimization remains limited. Traditional approaches often rely on pre-estimated dual parameters, which are inflexible and struggle with dynamic constraints. This paper introduces a novel primal-dual inference framework that embeds dual variable updates directly into the reverse diffusion process, enabling the model to adaptively satisfy average constraints while optimizing target objectives.

The core idea is to treat the dual variables as part of the inference state, updating them via dual ascent based on current constraint violations. A single score network conditioned on the dual variable is trained to represent a family of Gibbs distributions, allowing seamless adaptation during sampling. Theoretical analysis guarantees that the dual variables' time average converges to a neighborhood of the optimal solution, and residual errors are bounded, ensuring stability.

Empirical results across Gaussian mixtures, wireless resource allocation, and portfolio management demonstrate that the proposed PDI method outperforms existing baselines, achieving higher objective values, better constraint satisfaction, and increased robustness. For instance, in wireless scheduling, it reaches 95% of the target rate with 20% fewer violations, significantly outperforming fixed-parameter models.

This work marks a significant advance in probabilistic modeling under constraints, opening avenues for applications in autonomous systems, finance, and resource management. Its ability to dynamically adjust to changing conditions and guarantees of convergence make it a promising tool for real-world decision-making under uncertainty. Future efforts will focus on scaling to larger problems, reducing computational costs, and extending to more complex constraint structures, further broadening its impact.

Deep Dive

Abstract

This paper develops constrained diffusion models with primal-dual inference (PDI) to sample from optimal distributions of entropy-regularized optimization problems with \emph{average} constraints. We formalize constrained sampling in the Lagrangian dual domain, where the optimal distribution takes the form of a Gibbs distribution indexed by the optimal dual variable. Rather than estimating this dual multiplier before sampling and freezing it throughout generation, PDI jointly infers the optimal primal distribution and its parametrizing dual variable. Each reverse diffusion step denoises using the score field associated with the current multiplier and then updates the multiplier through dual ascent using the estimated constraint violation of the denoised samples. To enable this conditional score field, we train a single dual-conditioned score network over the family of Gibbs distributions induced by the dual variables encountered during inference. We prove that the time average of the dual variables generated along the inference trajectory converges to a neighborhood of the dual optimum and bound the effect of residual dual mismatch on the terminal distribution through schedule-dependent stability factors. We evaluate PDI on constrained sampling from a mixture of Gaussians, wireless resource allocation, and portfolio management.

cs.LG