From 3fb8c0ad75366639b85f3e06a3cf1260d130a8b2 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Sun, 29 Oct 2023 14:57:52 -0400 Subject: [PATCH] fix zeta --- neutrinos/sneut5.H | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/neutrinos/sneut5.H b/neutrinos/sneut5.H index 2734855984..862b259025 100644 --- a/neutrinos/sneut5.H +++ b/neutrinos/sneut5.H @@ -779,16 +779,16 @@ void nu_recomb(const sneutf_t& sf, } // equation 6.5 - z = std::exp(sf.zeta + nu); + z = std::exp(zeta + nu); dum = 1.0e0_rt + z; a1 = 1.0e0_rt/dum; a2 = 1.0e0_rt/bigj; sreco = nu_constants::tfac6 * 2.649e-18_rt * sf.ye * std::pow(sf.zbar, 13.0_rt) * sf.den * bigj * a1; if constexpr (do_derivatives) { - srecodt = sreco * (bigjdt * a2 - z * (sf.zetadt + nudt) * a1); - srecoda = sreco * (-1.0e0_rt * sf.abari + bigjda * a2 - z * (sf.zetada + nuda) * a1); - srecodz = sreco * (14.0e0_rt * sf.zbari + bigjdz * a2 - z * (sf.zetadz + nudz) * a1); + srecodt = sreco * (bigjdt * a2 - z * (zetadt + nudt) * a1); + srecoda = sreco * (-1.0e0_rt * sf.abari + bigjda * a2 - z * (zetada + nuda) * a1); + srecodz = sreco * (14.0e0_rt * sf.zbari + bigjdz * a2 - z * (zetadz + nudz) * a1); } }