Coreset-Induced Conditional Velocity Flow Matching

Department of Statistics, Purdue University
NeurIPS 2026
Three panels: (a) two Gaussian clusters summarised by two weighted coreset atoms; (b) samples from noise connected by straight velocity segments to endpoints near the clusters; (c) Stage II endpoints corrected onto the clusters.
CCVFM in one picture. (a) Stage I summarises the data by weighted coreset atoms. (b) Stage II draws velocities from the induced conditional law in closed form, with no network. (c) Stage III learns only the small correction that remains.
Abstract

Hierarchical flow matching models the full conditional law of the velocity, but its inner flow starts from pure Gaussian noise and must learn an entire noise-to-data transport. CCVFM starts that inner flow from a closed-form coreset of the data instead, so the network only has to learn a small correction.

We propose Coreset-Induced Conditional Velocity Flow Matching (CCVFM), a generative model that augments hierarchical rectified flow with a data-informed source distribution. CCVFM first compresses the target into weighted atoms using an entropic Sinkhorn coreset and lifts them to a Gaussian mixture. The induced conditional velocity law is then a closed-form Gaussian mixture that can be sampled without a learned neural sampler. A lightweight correction flow, trained from this surrogate source, refines the remaining surrogate-to-target residual rather than learning an entire noise-to-data map. We prove that the surrogate transport cost equals the target–surrogate Wasserstein gap under an explicit compression assumption, whereas the noise-source analogue has a dimension-scale lower bound. On MNIST, CIFAR-10, ImageNet-32 and CelebA-HQ 256, CCVFM reaches competitive few-step generation under matched architectures.

Method

Three stages, and only the last one is learned

STAGE I · CORESET

Compress the data, then lift to a mixture

An entropic-Sinkhorn coreset summarises the data by K weighted atoms. Closed-form probabilistic PCA lifts each atom to a Gaussian N(μk, LkLk⊤+σ²I), giving the surrogate ρ̃1.

No network · 95 s for CIFAR-10 at K = 10,000
STAGE II · CLOSED-FORM VELOCITY LAW

Read off the velocity law, no training

Plugging ρ̃1 into the conditional velocity identity gives an explicit K-component mixture π̃(v | xt, t). At t = 0 it is ρ̃1(x0+v): one categorical draw and one Gaussian draw.

No network · one closed-form draw
STAGE III · LEARNED CORRECTION

Learn only the residual

A correction network is trained by flow matching in velocity space, with the source drawn from π̃ and anchored to the mixture component of each training image. It refines a surrogate-to-target residual, not a noise-to-data map.

The only learned part · L Euler steps at sampling
Results

Strong few-step generation from a better source

0.75MNIST · FID50k
6.35CIFAR-10 · FID50k
8.76ImageNet-32 · FID50k
4.17CelebA-HQ 256 · FID

All four numbers use 51 function evaluations (50 correction steps plus the closed-form draw) and no classifier-free guidance.

Uncurated MNIST samples
MNIST 28×28
Uncurated CIFAR-10 samples
CIFAR-10 32×32
Uncurated ImageNet-32 samples
ImageNet-32 32×32
Uncurated CelebA-HQ 256 samples
CelebA-HQ 256×256
FID versus number of function evaluations. Left: MNIST, CCVFM reaches 0.75 at 51 NFE while HRF2 reaches 2.57 at 500 NFE. Right: ImageNet-32, CCVFM is 38 percent and 24 percent lower than HRF2 at 11 and 21 NFE.
Fewer steps for the same quality. On MNIST, CCVFM reaches FID50k 0.75 at 51 NFE; the best HRF2 result is 2.57 at 500 NFE. On ImageNet-32 the U-Net, data and protocol match HRF2 and only the source changes: FID drops by 38% at 11 NFE and 24% at 21 NFE.
Ablation

The source is what matters

Same correction U-Net, same Stage III budget (80k steps) on MNIST. The only change is where the inner flow starts.

FID versus NFE on MNIST for the coreset surrogate source and the Gaussian source
Inner-flow source6 NFE11 NFE21 NFE
Gaussian N(0, I) (HRF2), 2 seeds17.1710.749.94
Coreset surrogate (CCVFM), 5 seeds3.25 ± 0.312.60 ± 0.222.86 ± 0.36

FID50k, lower is better. Theorem 3 explains the gap: the surrogate source leaves only a surrogate-to-target residual to learn (see the checks below).

Three 2D targets (two moons, spiral, checkerboard): coreset atoms, Stage II samples, CCVFM samples at 2 steps, Gaussian-source samples at 2 steps, and sliced Wasserstein distance versus steps.
The same story in 2D. With 1–2 network evaluations, CCVFM is 7–28× closer to the target (sliced W2) than the same network trained from a Gaussian source. Reproduce with python examples/toy2d.py.
From the review

Checks added during peer review

Experiments from our public responses to the NeurIPS 2026 reviewers. Lightweight budgets, so absolute values are not comparable to the main tables. Code: experiments/rebuttal/.

Both identities of Theorem 3 hold numerically

Dataset (d)SurrogateThm 3(i)GaussianThm 3(ii)Ratio
ring-6 (2)8.1718.1588.0888.0741.0
moons (2)9.4009.4269.6189.6181.0
pinwheel (2)4.9934.9886.5086.4861.3
checkerboard (2)5.3075.3186.6836.6721.3
helix (3)4.6684.6778.3448.3441.8
MNIST (784)105.5105.41655.21655.615.7
CIFAR-10 (3072)1528.41523.26929.26929.04.5

Monte-Carlo estimates of E‖V1−V0‖² at t = 0 agree with both closed forms to 0.03–0.3%. On images the surrogate source shrinks the regression target 16× (MNIST) and 4.5× (CIFAR-10).

The low-rank covariance saturates near r = 50

r6 NFE11 NFE21 NFE
103.183.474.35
303.252.602.86
503.142.382.60
2003.222.362.64
7842.902.402.82

MNIST FID50k, mean over 3 training seeds (5 for r = 30).

Compression follows an effective dimension

TargetSlope in KReading
thin circle−0.503K−1/2
two moons−0.467K−1/2
helix (3D)−0.44steeper than −1/3
MNIST−0.162deff ≈ 12
CIFAR-10−0.136deff ≈ 15

Sliced W2 on low-dimensional targets; transport-cost proxy on images (nominal d = 784 / 3072).

Theory

Why a coreset source helps

  1. Transport-task reduction. At the generation boundary, the task left to the correction network is exactly the surrogate gap W2(ρ1, ρ̃1), which is O(K−1/d) under the Stage I compression assumption. With a Gaussian source the same task is bounded below by √d(√(σ1²+1)−1), a constant however large K or n is.
  2. Fewer Euler steps. The number of Euler steps needed for a given discretisation error scales with the square of that task size, so a smaller surrogate gap directly buys lower NFE.
  3. A smaller regression target. The conditional second moment of V1−V0 under the surrogate source depends on within-mode quantities; the Gaussian source pays the cross-mode diameter plus d.
Citation

BibTeX

@inproceedings{she2026ccvfm,
  title     = {Coreset-Induced Conditional Velocity Flow Matching},
  author    = {She, Zihua and Su, Jianxi and Wang, Xiao},
  booktitle = {Advances in Neural Information Processing Systems (NeurIPS)},
  year      = {2026},
  note      = {arXiv:2605.12951}
}