  We numerically compute features in the power-spectrum that originate from the
decay of fields during inflation. Using a simple, phenomenological, multi-field
setup, we increase the number of fields from a few to thousands. Whenever a
field decays, its associated potential energy is transferred into radiation,
causing a jump in the equation of state parameter and mode mixing at the
perturbed level. We observe discrete steps in the power-spectrum if the number
of fields is low, in agreement with analytic arguments in the literature. These
features become increasingly smeared out once many fields decay within a given
Hubble time. In this regime we confirm the validity of the analytic approach to
staggered inflation, which is based on a coarse-graining procedure. Our
numerical approach bridges the aforementioned analytic treatments, and can be
used in more complicated scenarios.
