Consistency Models

Yang Song, Prafulla Dhariwal, Mark Chen, Ilya Sutskever

2023 · ICML

Consistency Models

Problem

Framing

Diffusion models deliver strong sample quality but need 10–2000 denoising steps at inference. This paper replaces iterative denoising with a self-consistent map from any PF-ODE state to the trajectory origin, enabling one-step generation while retaining multistep refinement. One-step FID reaches 3.55 on CIFAR-10 and 6.20 on ImageNet 64×6464 \times 64.

Currently Used Methods

Direct antecedents

Proposed Method

Architecture

The model learns fθ(x,t)f_{\theta}(\mathbf{x}, t), a map from any PF-ODE point to the same origin xϵ\mathbf{x}_{\epsilon}.It enforces the boundary condition with a skip-style parameterization

fθ(x,t)=cskip(t)x+cout(t)Fθ(x,t),f_{\theta}(\mathbf{x}, t) = c_{\mathrm{skip}}(t)\mathbf{x} + c_{\mathrm{out}}(t)F_{\theta}(\mathbf{x}, t),

where cskip(ϵ)=1c_{\mathrm{skip}}(\epsilon)=1 and cout(ϵ)=0c_{\mathrm{out}}(\epsilon)=0. Experiments reuse EDM backbones with modified skip connections.

Architecture diagram: PF-ODE trajectories from data to noise, with the consistency model mapping multiple points on one trajectory to the same origin.

Loss / Objective

Distillation trains neighboring PF-ODE states to produce identical outputs.

LCDN(θ,θ;ϕ)=E[λ(tn)d ⁣(fθ(x^tnϕ,tn),fθ(xtn+1,tn+1))]\mathcal{L}^{N}_{\mathrm{CD}}(\theta, \theta^{-}; \phi) = \mathbb{E}\left[ \lambda(t_n) \, d\!\left(f_{\theta}(\hat{\mathbf{x}}^{\phi}_{t_n}, t_n), \, f_{\theta^{-}}(\mathbf{x}_{t_{n+1}}, t_{n+1})\right) \right]

Standalone training replaces the solver pair with two Gaussian-corrupted views of the same data point.

LCTN(θ,θ)=E[λ(tn)d ⁣(fθ(x+tn+1z,tn+1),fθ(x+tnz,tn))]\mathcal{L}^{N}_{\mathrm{CT}}(\theta, \theta^{-}) = \mathbb{E}\left[ \lambda(t_n) \, d\!\left(f_{\theta}(\mathbf{x}+t_{n+1}\mathbf{z}, t_{n+1}), \, f_{\theta^{-}}(\mathbf{x}+t_n\mathbf{z}, t_n)\right) \right]

Sampling Rule

One-step generation evaluates the model once at terminal noise. Multistep sampling alternates noise injection and consistency projection.

x^ϵ=fθ(x^T,T),x^τn=x+τn2ϵ2z,xfθ(x^τn,τn)\hat{\mathbf{x}}_{\epsilon} = f_{\theta}(\hat{\mathbf{x}}_{T}, T), \qquad \hat{\mathbf{x}}_{\tau_n} = \mathbf{x} + \sqrt{\tau_n^2 - \epsilon^2}\,\mathbf{z}, \qquad \mathbf{x} \leftarrow f_{\theta}(\hat{\mathbf{x}}_{\tau_n}, \tau_n)

For distillation, adjacent PF-ODE states are linked by a one-step solver update.

x^tnϕ:=xtn+1+(tntn+1)Φ(xtn+1,tn+1;ϕ)\hat{\mathbf{x}}^{\phi}_{t_n} := \mathbf{x}_{t_{n+1}} + (t_n - t_{n+1})\,\Phi(\mathbf{x}_{t_{n+1}}, t_{n+1}; \phi)

Training Procedure

Evaluation

Datasets

Metrics

Headline results

Ablations

Results figure

Figure 3 studies training choices on CIFAR-10: metric, solver, discretization steps, and adaptive schedules.

Results plots: CIFAR-10 FID versus training iterations for CD and CT under different metrics, ODE solvers, discretization steps, and adaptive schedules.

Table 1: Hyperparameters used for training CD and CT models

HyperparameterCIFAR-10 CDCIFAR-10 CTImageNet 64 × 64 CDImageNet 64 × 64 CTLSUN 256 × 256 CDLSUN 256 × 256 CT
Learning rate4e-44e-48e-68e-61e-51e-5
Batch size5125122048204820482048
μ\mu00.950.95
μ0\mu_00.90.950.95
s0s_0222
s1s_1150200150
NN184040
ODE solverHeunHeunHeun
EMA decay rate0.99990.99990.9999430.9999430.9999430.999943
Training iterations800k800k600k800k600k1000k
Mixed-Precision (FP16)NoNoYesYesYesYes
Dropout probability0.00.00.00.00.00.0
Number of GPUs8864646464

Method Strengths and Weaknesses

Strengths

Weaknesses

Suggestions from the authors

Links

Prior Papers

Further Papers

No vault papers identified as further work yet.