diff --git a/neutrinos/sneut5.H b/neutrinos/sneut5.H index 91206a73bc..de06839d3d 100644 --- a/neutrinos/sneut5.H +++ b/neutrinos/sneut5.H @@ -1178,19 +1178,19 @@ void nu_brem(const sneutf_t& sf, const auto [sin1, cos1] = amrex::Math::sincos(u); // double, triple, etc. angle formulas - // sin/cos (2 fac1 tau) + // sin/cos (2 u) const Real sin2 = 2.0_rt * sin1 * cos1; const Real cos2 = 2.0_rt * cos1 * cos1 - 1.0_rt; - // sin/cos (3 fac1 tau) + // sin/cos (3 u) const Real sin3 = sin1 * (3.0_rt - 4.0_rt * sin1 * sin1); const Real cos3 = cos1 * (4.0_rt * cos1 * cos1 - 3.0_rt); - // sin/cos (4 fac1 tau) -- use double angle on sin2/cos2 + // sin/cos (4 u) -- use double angle on sin2/cos2 const Real sin4 = 2.0_rt * sin2 * cos2; const Real cos4 = 2.0_rt * cos2 * cos2 - 1.0_rt; - // sin/cos (5 fac1 tau) + // sin/cos (5 u) //const Real sin5 = sin1 * (5.0_rt - sin1 * sin1 * (20.0_rt - 16.0_rt * sin1 * sin1)); const Real cos5 = cos1 * (cos1 * cos1 * (16.0_rt * cos1 * cos1 - 20.0_rt) + 5.0_rt);