Description
Warm inflation (WI) provides an alternative to the standard cold inflation paradigm, in which a subdominant radiation bath coexists with the inflaton field and continuously extracts energy through dissipative interactions. This framework allows for a smooth transition to a radiation-dominated universe and can accommodate steeper inflationary potentials due to the presence of additional friction. The inclusion of dissipative effects and radiation fluctuations makes the computation of primordial perturbations in WI significantly more involved, requiring a fully numerical treatment. In this work, we present the Stochastic Warm Inflation Module (SWIM), a numerical framework designed to solve the coupled perturbation equations and compute inflationary observables. Existing numerical tools, such as WarmSPy and WI2Easy, typically compute correction factors to approximate analytical expressions for the power spectrum. In contrast, SWIM computes the full primordial power spectrum directly as a function of the comoving wavenumber $k$, without relying on simplifying assumptions inherent in semi-analytical approaches while also efficiently computing correction factors. This enables a more accurate treatment of model-dependent effects and captures features that may be missed in approximate methods. Our results demonstrate that full numerical treatments can differ significantly from semi-analytical approaches in certain regimes, highlighting limitations of previous approximation-based approaches.