From 0f7b90fbbad4f45452f1908ef1dc7b755d178ba9 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Sun, 29 Oct 2023 19:43:18 -0400 Subject: [PATCH] a little more factoring --- neutrinos/sneut5.H | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/neutrinos/sneut5.H b/neutrinos/sneut5.H index fdf11e160c..ec8281bf97 100644 --- a/neutrinos/sneut5.H +++ b/neutrinos/sneut5.H @@ -843,8 +843,8 @@ void nu_photo(const sneutf_t& sf, const Real cos4 = 2.0_rt * cos2 * cos2 - 1.0_rt; // sin/cos (5 fac1 tau) - const Real sin5 = 5.0_rt * sin1 - sin1 * 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 * cos1; + 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); const auto [xast, last] = amrex::Math::sincos(nu_constants::fac2 * tau);