a little more factoring

AMReX-Astro · Oct 29, 2023 · 0f7b90f · 0f7b90f
1 parent 0dff547
commit 0f7b90f
Showing 1 changed file with 2 additions and 2 deletions.
diff --git 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);