Which parts of modern generative modelling actually matter? Starting from IMLE and adding back only the essentials: we achieve FID 2.56 on ImageNet 256 in a single step.
Generative models have undergone many generations of evolution, from VAEs/GANs to diffusion/flow matching. Along the way, the underlying techniques have become more complicated and various beliefs about what drives strong empirical performance have taken hold. Due to the success of diffusion models and flow matching, one of the more common beliefs is the importance of transforming the noise distribution to the data distribution gradually through many small transformations.
We ask whether this is truly necessary, and take a minimalist approach to designing a competitive generative model. We start with the bare-bones essentials, namely just a training objective and a model. We purposefully make both simple. For the training objective, we choose Implicit Maximum Likelihood Estimation (IMLE), and eschew more complicated alternatives such as variational inference, adversarial training and numerical integration. For the model, we eschew transformers and instead choose a moderately sized convolutional network. Then we judiciously added elements that are truly essential, which surprisingly do not include iterative denoising. The result is a single-step parameter-efficient generative model that produces high quality samples at fast speed: it achieves an FID of 2.56 on ImageNet 256 and simultaneously attains good precision and recall.
Diffusion and flow matching models, both members of a broader family known as stochastic interpolant models
Their advantage is commonly attributed to how they decompose generation
During training, each prior sample is randomly coupled with a data sample. Many different couplings may intersect at a given point, so the network can only learn to predict their average. Thus during inference, a single step regresses to the data mean rather than a real sample. Multi-step iterative sampling is required to yield an accurate sample.
Implicit Maximum Likelihood Estimation
Training such a generator directly is generally difficult because the likelihood of the data under it is intractable. IMLE provides a workaround, implicitly maximizing that likelihood without requiring it to be known explicitly. IMLE optimizes a likelihood-based objective akin to regression, without the variational bounds and encoders of VAEs, the adversarial training of GANs, or the iterative numerical integration used by stochastic interpolant models.
The algorithm decomposes into two phases. In the matching phase, latent codes are sampled from the prior and each real image is matched to its nearest generated sample, establishing a coupling between the prior and the data. In the optimization phase, the parameters are updated to bring each generated sample closer to its matched real image. The two phases repeat until convergence.
This paper asks a simple question: is all of the machinery that modern generative models have accumulated, from elaborate architectures
Our two contributions improve the optimization and matching phases of IMLE respectively.
Stochastic interpolant methods are commonly understood to perform iterative denoising. At training time, a denoising network is trained to convert a noisy image into a clean one. At test time, that network is applied repeatedly to convert pure noise into a generated sample. This iterative application is their key distinguishing feature, and it is what they are commonly believed to owe their success to. It also creates a mismatch between training and test: the network is applied once per sample during training, but iteratively during generation.
We consider an alternative view of this mismatch. Rather than viewing the model as one instance of the denoising network, we unroll the iterations of the sampler and view the model as their composition, \(h_\theta = g^1_\theta \circ \cdots \circ g^T_\theta\), where \(g^t_\theta\) is the update applied at iteration \(t\). Under this view the model is applied just once at test time, and it is training that deviates from the norm, evaluating only part of the model at a time. Each update is now an intermediate layer of \(h_\theta\). We term the conventional view the training-centric view, and this one the testing-centric view.
Read this way, the training objective reveals something specific. These methods supervise the denoising network’s prediction at every time step, with equal weight applied to all steps, so every intermediate layer of \(h_\theta\) receives a direct training signal. We call this per-stage supervision, and we hypothesize that the sample quality of many-step models stems from it rather than from iterative sampling itself.
The unrolled view applies naturally to IMLE. Its generator can also be seen as a prior-to-data network composed of stages, only executed in a single forward pass rather than through iterative sampling. We build it as a stack of upsampling blocks, \(f_\theta = f_{\theta^{L}} \circ \cdots \circ f_{\theta^{1}}\), in which each stage operates at a progressively higher resolution. This recasts IMLE in a form similar to the diffusion construction, with each stage now analogous to \(g^t_\theta\), only indexed by resolution rather than noise level.
The architecture of existing IMLE methods
For the IMLE algorithm, the model distribution changes during training, so the matching between model samples and ground truth data points needs to be recomputed. Due to randomness in the sampling, there may be no samples drawn from a mode of the model distribution, and as a result, a ground truth data point that can otherwise be explained by that mode would not be well matched. Note that this is different from mode collapse: the poor matching results from randomness in the sampling, not a deficiency in the model distribution.
The standard squared $l_2$ loss lets these mismatches dominate the mean loss. In essence this is a problem of outliers: a small number of large-residual mismatched pairs contribute disproportionately to the mean loss compared to a much larger pool of well-matched pairs.
The classical remedy in robust statistics
We call the resulting method ROMS-IMLE, short for Robust Multi-Stage Implicit Maximum Likelihood Estimation. We ablate our contributions on Oxford Flowers
| Configuration | FID ↓ | Precision ↑ | Recall ↑ |
|---|---|---|---|
| Vanilla IMLE | 62.61 ± 2.28 | 0.45 ± 0.03 | 0.46 ± 0.01 |
| IMLE + Multi-Stage | 31.24 ± 1.08 | 0.79 ± 0.01 | 0.72 ± 0.02 |
| IMLE + Multi-Stage + Geman-McClure (ours) | 13.90 ± 0.08 | 0.93 ± 0.00 | 0.73 ± 0.02 |
Our approach works well in both pixel and latent space.
@article{vashist2026romsimle,
title={ROMS-IMLE: A Minimalist Approach to Competitive Single-Step Generative Modelling},
author={Vashist, Chirag and Li, Ke},
journal={arXiv preprint arXiv:2607.19332},
year={2026}
}