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
