def smooth_decoupling_test(params): neutrino = Particle(params=params, **SMP.leptons.neutrino_e) energy_density = REGIMES.RADIATION.energy_density(neutrino) density = REGIMES.RADIATION.density(neutrino) pressure = REGIMES.RADIATION.pressure(neutrino) numerator = REGIMES.RADIATION.numerator(neutrino) denominator = REGIMES.RADIATION.denominator(neutrino) assert neutrino.energy_density() - energy_density < eps assert neutrino.density() - density < eps assert neutrino.pressure() - pressure < eps assert neutrino.numerator() - numerator < eps assert neutrino.denominator() - denominator < eps