Skip to content

Commit

Permalink
sqrt() -> std::sqrt() and exp() -> std::exp() (#1330)
Browse files Browse the repository at this point in the history
  • Loading branch information
zingale authored Sep 10, 2023
1 parent 3c75538 commit 2c408a9
Show file tree
Hide file tree
Showing 2 changed files with 4 additions and 4 deletions.
2 changes: 1 addition & 1 deletion EOS/breakout/actual_eos.H
Original file line number Diff line number Diff line change
Expand Up @@ -123,7 +123,7 @@ void actual_eos (I input, T& state)
state.gam1 = gamma_const;

// sound speed
state.cs = sqrt(gamma_const * poverrho);
state.cs = std::sqrt(gamma_const * poverrho);

state.dpdr_e = poverrho;
state.dpde = (gamma_const - 1.0_rt) * state.rho;
Expand Down
6 changes: 3 additions & 3 deletions neutrinos/sneut5.H
Original file line number Diff line number Diff line change
Expand Up @@ -347,7 +347,7 @@ void sneut5(const Real temp, const Real den,
t9 = temp * 1.0e-9_rt;
xl = t9 * con1;
xldt = 1.0e-9_rt * con1;
xlp5 = sqrt(xl);
xlp5 = std::sqrt(xl);
xl2 = xl*xl;
xl3 = xl2*xl;
xl4 = xl3*xl;
Expand Down Expand Up @@ -391,10 +391,10 @@ void sneut5(const Real temp, const Real den,
a2 = 2.084e20_rt + 2.0e0_rt*1.872e21_rt*zeta;

if (t9 < 10.0_rt) {
b1 = exp(-5.5924e0_rt*zeta);
b1 = std::exp(-5.5924e0_rt*zeta);
b2 = -b1*5.5924e0_rt;
} else {
b1 = exp(-4.9924e0_rt*zeta);
b1 = std::exp(-4.9924e0_rt*zeta);
b2 = -b1*4.9924e0_rt;
}

Expand Down

0 comments on commit 2c408a9

Please sign in to comment.