Skip to content

Commit

Permalink
fix some comments in sneut5
Browse files Browse the repository at this point in the history
  • Loading branch information
zingale committed Nov 9, 2023
1 parent 39a0f40 commit 63e1fda
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions neutrinos/sneut5.H
Original file line number Diff line number Diff line change
Expand Up @@ -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);

Expand Down

0 comments on commit 63e1fda

Please sign in to comment.