Skip to content

Commit

Permalink
update pynucastro nets to have bion/mion constexpr
Browse files Browse the repository at this point in the history
  • Loading branch information
@zingale
zingale committed Jan 3, 2024
1 parent 6dc51ad commit 8ee600e
Show file tree
Hide file tree
Showing 32 changed files with 731 additions and 731 deletions.
14 changes: 7 additions & 7 deletions networks/CNO_extras/actual_rhs.H
View file
Original file line number Diff line number Diff line change
Expand Up @@ -635,7 +635,7 @@ void rhs_nuc(const burn_t& state,
-screened_rates(k_He4_O16_to_Ne20)*Y(He4)*Y(O16)*state.rho +
-screened_rates(k_He4_Ne18_to_Mg22)*Y(He4)*Y(Ne18)*state.rho +
-screened_rates(k_He4_Ne20_to_Mg24)*Y(He4)*Y(Ne20)*state.rho +
0.5*screened_rates(k_C12_C12_to_He4_Ne20)*amrex::Math::powi<2>(Y(C12))*state.rho +
0.5*screened_rates(k_C12_C12_to_He4_Ne20)*std::pow(Y(C12), 2)*state.rho +
-screened_rates(k_He4_N13_to_p_O16)*Y(He4)*Y(N13)*state.rho +
screened_rates(k_p_N15_to_He4_C12)*Y(N15)*Y(H1)*state.rho +
-screened_rates(k_He4_O14_to_p_F17)*Y(He4)*Y(O14)*state.rho +
Expand All @@ -644,16 +644,16 @@ void rhs_nuc(const burn_t& state,
screened_rates(k_p_O18_to_He4_N15)*Y(O18)*Y(H1)*state.rho +
screened_rates(k_p_F18_to_He4_O15)*Y(F18)*Y(H1)*state.rho +
screened_rates(k_p_F19_to_He4_O16)*Y(F19)*Y(H1)*state.rho +
-0.5*screened_rates(k_He4_He4_He4_to_C12)*amrex::Math::powi<3>(Y(He4))*amrex::Math::powi<2>(state.rho) +
-0.5*screened_rates(k_He4_He4_He4_to_C12)*std::pow(Y(He4), 3)*std::pow(state.rho, 2) +
screened_rates(k_p_Ne20_to_He4_F17)*Y(Ne20)*Y(H1)*state.rho;

ydot_nuc(C12) =
-screened_rates(k_p_C12_to_N13)*Y(C12)*Y(H1)*state.rho +
-screened_rates(k_He4_C12_to_O16)*Y(C12)*Y(He4)*state.rho +
-screened_rates(k_C12_C12_to_He4_Ne20)*amrex::Math::powi<2>(Y(C12))*state.rho +
-screened_rates(k_C12_C12_to_He4_Ne20)*std::pow(Y(C12), 2)*state.rho +
screened_rates(k_p_N15_to_He4_C12)*Y(N15)*Y(H1)*state.rho +
-screened_rates(k_C12_O16_to_He4_Mg24)*Y(C12)*Y(O16)*state.rho +
0.16666666666666667*screened_rates(k_He4_He4_He4_to_C12)*amrex::Math::powi<3>(Y(He4))*amrex::Math::powi<2>(state.rho);
0.16666666666666667*screened_rates(k_He4_He4_He4_to_C12)*std::pow(Y(He4), 3)*std::pow(state.rho, 2);

ydot_nuc(C13) =
screened_rates(k_N13_to_C13_weak_wc12)*Y(N13) +
Expand Down Expand Up @@ -747,7 +747,7 @@ void rhs_nuc(const burn_t& state,
screened_rates(k_He4_O16_to_Ne20)*Y(He4)*Y(O16)*state.rho +
screened_rates(k_p_F19_to_Ne20)*Y(F19)*Y(H1)*state.rho +
-screened_rates(k_He4_Ne20_to_Mg24)*Y(He4)*Y(Ne20)*state.rho +
0.5*screened_rates(k_C12_C12_to_He4_Ne20)*amrex::Math::powi<2>(Y(C12))*state.rho +
0.5*screened_rates(k_C12_C12_to_He4_Ne20)*std::pow(Y(C12), 2)*state.rho +
-screened_rates(k_p_Ne20_to_He4_F17)*Y(Ne20)*Y(H1)*state.rho;

ydot_nuc(Mg22) =
Expand Down Expand Up @@ -865,7 +865,7 @@ void jac_nuc(const burn_t& state,
scratch = screened_rates(k_p_F18_to_He4_O15)*Y(F18)*state.rho + screened_rates(k_p_F19_to_He4_O16)*Y(F19)*state.rho + screened_rates(k_p_N15_to_He4_C12)*Y(N15)*state.rho + screened_rates(k_p_Ne20_to_He4_F17)*Y(Ne20)*state.rho + screened_rates(k_p_O17_to_He4_N14)*Y(O17)*state.rho + screened_rates(k_p_O18_to_He4_N15)*Y(O18)*state.rho;
jac.set(He4, H1, scratch);

scratch = -screened_rates(k_He4_C12_to_O16)*Y(C12)*state.rho - 1.5*screened_rates(k_He4_He4_He4_to_C12)*amrex::Math::powi<2>(Y(He4))*amrex::Math::powi<2>(state.rho) - screened_rates(k_He4_N13_to_p_O16)*Y(N13)*state.rho - screened_rates(k_He4_N14_to_F18)*Y(N14)*state.rho - screened_rates(k_He4_N15_to_F19)*Y(N15)*state.rho - screened_rates(k_He4_Ne18_to_Mg22)*Y(Ne18)*state.rho - screened_rates(k_He4_Ne20_to_Mg24)*Y(Ne20)*state.rho - screened_rates(k_He4_O14_to_Ne18)*Y(O14)*state.rho - screened_rates(k_He4_O14_to_p_F17)*Y(O14)*state.rho - screened_rates(k_He4_O15_to_Ne19)*Y(O15)*state.rho - screened_rates(k_He4_O16_to_Ne20)*Y(O16)*state.rho;
scratch = -screened_rates(k_He4_C12_to_O16)*Y(C12)*state.rho - 1.5*screened_rates(k_He4_He4_He4_to_C12)*std::pow(Y(He4), 2)*std::pow(state.rho, 2) - screened_rates(k_He4_N13_to_p_O16)*Y(N13)*state.rho - screened_rates(k_He4_N14_to_F18)*Y(N14)*state.rho - screened_rates(k_He4_N15_to_F19)*Y(N15)*state.rho - screened_rates(k_He4_Ne18_to_Mg22)*Y(Ne18)*state.rho - screened_rates(k_He4_Ne20_to_Mg24)*Y(Ne20)*state.rho - screened_rates(k_He4_O14_to_Ne18)*Y(O14)*state.rho - screened_rates(k_He4_O14_to_p_F17)*Y(O14)*state.rho - screened_rates(k_He4_O15_to_Ne19)*Y(O15)*state.rho - screened_rates(k_He4_O16_to_Ne20)*Y(O16)*state.rho;
jac.set(He4, He4, scratch);

scratch = 1.0*screened_rates(k_C12_C12_to_He4_Ne20)*Y(C12)*state.rho + screened_rates(k_C12_O16_to_He4_Mg24)*Y(O16)*state.rho - screened_rates(k_He4_C12_to_O16)*Y(He4)*state.rho;
Expand Down Expand Up @@ -910,7 +910,7 @@ void jac_nuc(const burn_t& state,
scratch = -screened_rates(k_p_C12_to_N13)*Y(C12)*state.rho + screened_rates(k_p_N15_to_He4_C12)*Y(N15)*state.rho;
jac.set(C12, H1, scratch);

scratch = -screened_rates(k_He4_C12_to_O16)*Y(C12)*state.rho + 0.5*screened_rates(k_He4_He4_He4_to_C12)*amrex::Math::powi<2>(Y(He4))*amrex::Math::powi<2>(state.rho);
scratch = -screened_rates(k_He4_C12_to_O16)*Y(C12)*state.rho + 0.5*screened_rates(k_He4_He4_He4_to_C12)*std::pow(Y(He4), 2)*std::pow(state.rho, 2);
jac.set(C12, He4, scratch);

scratch = -2.0*screened_rates(k_C12_C12_to_He4_Ne20)*Y(C12)*state.rho - screened_rates(k_C12_O16_to_He4_Mg24)*Y(O16)*state.rho - screened_rates(k_He4_C12_to_O16)*Y(He4)*state.rho - screened_rates(k_p_C12_to_N13)*Y(H1)*state.rho;
Expand Down
Binary file modified networks/CNO_extras/cno_extras.png
View file
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified networks/CNO_extras/cno_extras_hide_alpha.png
View file
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
12 changes: 6 additions & 6 deletions networks/ECSN/actual_rhs.H
View file
Original file line number Diff line number Diff line change
Expand Up @@ -316,7 +316,7 @@ void rhs_nuc(const burn_t& state,
ydot_nuc(H1) =
-screened_rates(k_p_Al27_to_Si28)*Y(Al27)*Y(H1)*state.rho +
-screened_rates(k_p_P31_to_S32)*Y(P31)*Y(H1)*state.rho +
0.5*screened_rates(k_O16_O16_to_p_P31)*amrex::Math::powi<2>(Y(O16))*state.rho +
0.5*screened_rates(k_O16_O16_to_p_P31)*std::pow(Y(O16), 2)*state.rho +
(-screened_rates(k_p_Al27_to_He4_Mg24)*Y(Al27)*Y(H1)*state.rho + screened_rates(k_He4_Mg24_to_p_Al27)*Y(He4)*Y(Mg24)*state.rho) +
(-screened_rates(k_p_P31_to_He4_Si28)*Y(P31)*Y(H1)*state.rho + screened_rates(k_He4_Si28_to_p_P31)*Y(He4)*Y(Si28)*state.rho);

Expand All @@ -326,14 +326,14 @@ void rhs_nuc(const burn_t& state,
-screened_rates(k_He4_Mg24_to_Si28)*Y(He4)*Y(Mg24)*state.rho +
-screened_rates(k_He4_Al27_to_P31)*Y(Al27)*Y(He4)*state.rho +
-screened_rates(k_He4_Si28_to_S32)*Y(He4)*Y(Si28)*state.rho +
0.5*screened_rates(k_O16_O16_to_He4_Si28)*amrex::Math::powi<2>(Y(O16))*state.rho +
0.5*screened_rates(k_O16_O16_to_He4_Si28)*std::pow(Y(O16), 2)*state.rho +
(screened_rates(k_p_Al27_to_He4_Mg24)*Y(Al27)*Y(H1)*state.rho + -screened_rates(k_He4_Mg24_to_p_Al27)*Y(He4)*Y(Mg24)*state.rho) +
(screened_rates(k_p_P31_to_He4_Si28)*Y(P31)*Y(H1)*state.rho + -screened_rates(k_He4_Si28_to_p_P31)*Y(He4)*Y(Si28)*state.rho);

ydot_nuc(O16) =
(-screened_rates(k_He4_O16_to_Ne20)*Y(He4)*Y(O16)*state.rho + screened_rates(k_Ne20_to_He4_O16)*Y(Ne20)) +
-screened_rates(k_O16_O16_to_p_P31)*amrex::Math::powi<2>(Y(O16))*state.rho +
-screened_rates(k_O16_O16_to_He4_Si28)*amrex::Math::powi<2>(Y(O16))*state.rho;
-screened_rates(k_O16_O16_to_p_P31)*std::pow(Y(O16), 2)*state.rho +
-screened_rates(k_O16_O16_to_He4_Si28)*std::pow(Y(O16), 2)*state.rho;

ydot_nuc(O20) =
(-screened_rates(k_O20_to_F20)*Y(O20) + screened_rates(k_F20_to_O20)*Y(F20));
Expand Down Expand Up @@ -361,13 +361,13 @@ void rhs_nuc(const burn_t& state,
screened_rates(k_He4_Mg24_to_Si28)*Y(He4)*Y(Mg24)*state.rho +
screened_rates(k_p_Al27_to_Si28)*Y(Al27)*Y(H1)*state.rho +
-screened_rates(k_He4_Si28_to_S32)*Y(He4)*Y(Si28)*state.rho +
0.5*screened_rates(k_O16_O16_to_He4_Si28)*amrex::Math::powi<2>(Y(O16))*state.rho +
0.5*screened_rates(k_O16_O16_to_He4_Si28)*std::pow(Y(O16), 2)*state.rho +
(screened_rates(k_p_P31_to_He4_Si28)*Y(P31)*Y(H1)*state.rho + -screened_rates(k_He4_Si28_to_p_P31)*Y(He4)*Y(Si28)*state.rho);

ydot_nuc(P31) =
screened_rates(k_He4_Al27_to_P31)*Y(Al27)*Y(He4)*state.rho +
-screened_rates(k_p_P31_to_S32)*Y(P31)*Y(H1)*state.rho +
0.5*screened_rates(k_O16_O16_to_p_P31)*amrex::Math::powi<2>(Y(O16))*state.rho +
0.5*screened_rates(k_O16_O16_to_p_P31)*std::pow(Y(O16), 2)*state.rho +
(-screened_rates(k_p_P31_to_He4_Si28)*Y(P31)*Y(H1)*state.rho + screened_rates(k_He4_Si28_to_p_P31)*Y(He4)*Y(Si28)*state.rho);

ydot_nuc(S32) =
Expand Down
Binary file modified networks/He-C-Fe-group/He-C-Fe-group.png
View file
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Loading

0 comments on commit 8ee600e

Please sign in to comment.