Skip to content

Commit

Permalink
update He-C-Fe-group to include the positron parts of table rates
Browse files Browse the repository at this point in the history 
now we include positron-capture and positron-decay
we also include n -> p and p -> n from Langanke.
  • Loading branch information
@zingale
zingale committed Oct 18, 2023
1 parent 8b2c485 commit 427bb0f
Show file tree
Hide file tree
Showing 20 changed files with 2,167 additions and 1,553 deletions.
244 changes: 122 additions & 122 deletions networks/He-C-Fe-group/55co-55fe_electroncapture.dat
View file

Large diffs are not rendered by default.

286 changes: 143 additions & 143 deletions networks/He-C-Fe-group/55fe-55co_betadecay.dat
View file

Large diffs are not rendered by default.

148 changes: 148 additions & 0 deletions networks/He-C-Fe-group/55fe-55mn_electroncapture.dat
View file

Large diffs are not rendered by default.

148 changes: 148 additions & 0 deletions networks/He-C-Fe-group/55mn-55fe_betadecay.dat
View file

Large diffs are not rendered by default.

234 changes: 117 additions & 117 deletions networks/He-C-Fe-group/56co-56fe_electroncapture.dat
View file

Large diffs are not rendered by default.

264 changes: 132 additions & 132 deletions networks/He-C-Fe-group/56co-56ni_betadecay.dat
View file

Large diffs are not rendered by default.

286 changes: 143 additions & 143 deletions networks/He-C-Fe-group/56fe-56co_betadecay.dat
View file

Large diffs are not rendered by default.

122 changes: 61 additions & 61 deletions networks/He-C-Fe-group/56ni-56co_electroncapture.dat
View file

Large diffs are not rendered by default.

286 changes: 143 additions & 143 deletions networks/He-C-Fe-group/57co-57ni_betadecay.dat
View file

Large diffs are not rendered by default.

230 changes: 115 additions & 115 deletions networks/He-C-Fe-group/57ni-57co_electroncapture.dat
View file

Large diffs are not rendered by default.

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.
542 changes: 272 additions & 270 deletions networks/He-C-Fe-group/actual_network.H
View file

Large diffs are not rendered by default.

141 changes: 71 additions & 70 deletions networks/He-C-Fe-group/actual_network_data.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -11,76 +11,75 @@ namespace network
namespace NSE_INDEX
{
AMREX_GPU_MANAGED amrex::Array2D<int, 1, Rates::NumRates, 1, 7, Order::C> rate_indices {
-1, 1, 3, -1, -1, 4, 101,
-1, 2, 3, -1, -1, 6, 102,
-1, 2, 5, -1, -1, 7, 103,
-1, 2, 6, -1, -1, 8, 104,
-1, 2, 7, -1, -1, 10, 106,
-1, 2, 8, -1, -1, 12, 108,
-1, 1, 9, -1, -1, 10, 105,
-1, 1, 11, -1, -1, 12, 107,
-1, 2, 12, -1, -1, 14, 110,
-1, 1, 13, -1, -1, 14, 109,
-1, 2, 14, -1, -1, 16, 112,
-1, 1, 15, -1, -1, 16, 111,
-1, 2, 20, -1, -1, 23, 114,
-1, 1, 21, -1, -1, 23, 113,
-1, 2, 21, -1, -1, 28, 115,
-1, 2, 23, -1, -1, 31, 117,
-1, 1, 28, -1, -1, 31, 116,
-1, 2, 28, -1, -1, 34, 118,
-1, 2, 31, -1, -1, 35, 120,
-1, 1, 34, -1, -1, 35, 119,
-1, 1, 3, -1, -1, 4, 102,
-1, 2, 3, -1, -1, 6, 103,
-1, 2, 5, -1, -1, 7, 104,
-1, 2, 6, -1, -1, 8, 105,
-1, 2, 7, -1, -1, 10, 107,
-1, 2, 8, -1, -1, 12, 109,
-1, 1, 9, -1, -1, 10, 106,
-1, 1, 11, -1, -1, 12, 108,
-1, 2, 12, -1, -1, 14, 111,
-1, 1, 13, -1, -1, 14, 110,
-1, 2, 14, -1, -1, 16, 113,
-1, 1, 15, -1, -1, 16, 112,
-1, 2, 20, -1, -1, 23, 115,
-1, 1, 21, -1, -1, 23, 114,
-1, 2, 21, -1, -1, 28, 116,
-1, 2, 23, -1, -1, 31, 118,
-1, 1, 28, -1, -1, 31, 117,
-1, 2, 28, -1, -1, 34, 119,
-1, 2, 31, -1, -1, 35, 121,
-1, 1, 34, -1, -1, 35, 120,
-1, 3, 3, -1, 1, 11, -1,
-1, 3, 3, -1, 2, 8, -1,
-1, 2, 4, -1, 1, 6, 122,
-1, 2, 4, -1, 1, 6, 123,
-1, 3, 6, -1, 1, 13, -1,
-1, 3, 6, -1, 2, 12, -1,
-1, 6, 6, -1, 1, 15, -1,
-1, 6, 6, -1, 2, 14, -1,
-1, 2, 7, -1, 1, 9, 124,
-1, 1, 11, -1, 2, 8, 123,
-1, 1, 13, -1, 2, 12, 125,
-1, 1, 15, -1, 2, 14, 126,
-1, 2, 20, -1, 1, 21, 127,
-1, 2, 23, -1, 1, 28, 128,
-1, 1, 34, -1, 2, 31, 129,
2, 2, 2, -1, -1, 3, 121,
-1, 2, 7, -1, 1, 9, 125,
-1, 1, 11, -1, 2, 8, 124,
-1, 1, 13, -1, 2, 12, 126,
-1, 1, 15, -1, 2, 14, 127,
-1, 2, 20, -1, 1, 21, 128,
-1, 2, 23, -1, 1, 28, 129,
-1, 1, 34, -1, 2, 31, 130,
2, 2, 2, -1, -1, 3, 122,
-1, 3, 3, -1, -1, 12, -1,
-1, 6, 6, -1, -1, 16, -1,
-1, 3, 6, -1, -1, 14, -1,
-1, -1, 0, -1, -1, 1, -1,
-1, 1, 22, -1, -1, 27, 134,
-1, 0, 23, -1, -1, 24, 130,
-1, 0, 24, -1, -1, 25, 131,
-1, 2, 24, -1, -1, 32, 142,
-1, 0, 25, -1, -1, 26, 132,
-1, 1, 25, -1, -1, 28, 135,
-1, 2, 25, -1, -1, 33, 145,
-1, 0, 26, -1, -1, 27, 133,
-1, 1, 26, -1, -1, 29, 137,
-1, 1, 27, -1, -1, 30, 139,
-1, 0, 28, -1, -1, 29, 136,
-1, 0, 29, -1, -1, 30, 138,
-1, 1, 29, -1, -1, 32, 141,
-1, 1, 30, -1, -1, 33, 144,
-1, 0, 31, -1, -1, 32, 140,
-1, 0, 32, -1, -1, 33, 143,
-1, 1, 33, -1, -1, 34, 146,
-1, 2, 21, -1, 1, 25, 149,
-1, 2, 24, -1, 1, 29, 156,
-1, 0, 26, -1, 1, 22, 147,
-1, 0, 28, -1, 1, 26, 152,
-1, 2, 28, -1, 1, 33, 160,
-1, 0, 29, -1, 1, 27, 154,
-1, 1, 30, -1, 2, 25, 151,
-1, 0, 31, -1, 1, 29, 155,
-1, 0, 31, -1, 2, 24, 148,
-1, 0, 32, -1, 1, 30, 158,
-1, 0, 32, -1, 2, 25, 150,
-1, 0, 33, -1, 2, 26, 153,
-1, 0, 34, -1, 2, 29, 157,
-1, 0, 35, -1, 2, 32, 159,
-1, 1, 22, -1, -1, 27, 135,
-1, 0, 23, -1, -1, 24, 131,
-1, 0, 24, -1, -1, 25, 132,
-1, 2, 24, -1, -1, 32, 143,
-1, 0, 25, -1, -1, 26, 133,
-1, 1, 25, -1, -1, 28, 136,
-1, 2, 25, -1, -1, 33, 146,
-1, 0, 26, -1, -1, 27, 134,
-1, 1, 26, -1, -1, 29, 138,
-1, 1, 27, -1, -1, 30, 140,
-1, 0, 28, -1, -1, 29, 137,
-1, 0, 29, -1, -1, 30, 139,
-1, 1, 29, -1, -1, 32, 142,
-1, 1, 30, -1, -1, 33, 145,
-1, 0, 31, -1, -1, 32, 141,
-1, 0, 32, -1, -1, 33, 144,
-1, 1, 33, -1, -1, 34, 147,
-1, 2, 21, -1, 1, 25, 150,
-1, 2, 24, -1, 1, 29, 157,
-1, 0, 26, -1, 1, 22, 148,
-1, 0, 28, -1, 1, 26, 153,
-1, 2, 28, -1, 1, 33, 161,
-1, 0, 29, -1, 1, 27, 155,
-1, 1, 30, -1, 2, 25, 152,
-1, 0, 31, -1, 1, 29, 156,
-1, 0, 31, -1, 2, 24, 149,
-1, 0, 32, -1, 1, 30, 159,
-1, 0, 32, -1, 2, 25, 151,
-1, 0, 33, -1, 2, 26, 154,
-1, 0, 34, -1, 2, 29, 158,
-1, 0, 35, -1, 2, 32, 160,
-1, -1, -1, -1, -1, -1, -1,
-1, -1, -1, -1, -1, -1, -1,
-1, -1, -1, -1, -1, -1, -1,
Expand All @@ -93,23 +92,25 @@ namespace NSE_INDEX
-1, -1, -1, -1, -1, -1, -1,
-1, -1, -1, -1, -1, -1, -1,
-1, -1, -1, -1, -1, -1, -1,
-1, -1, 28, -1, -1, 26, 87,
-1, -1, 29, -1, -1, 27, 89,
-1, -1, 28, -1, -1, 26, 86,
-1, -1, 29, -1, -1, 27, 88,
-1, -1, 29, -1, -1, 31, -1,
-1, -1, 30, -1, -1, 32, -1,
-1, -1, 26, -1, -1, 28, -1,
-1, -1, 26, -1, -1, 22, 90,
-1, -1, 26, -1, -1, 22, 89,
-1, -1, 27, -1, -1, 29, -1,
-1, -1, 22, -1, -1, 26, -1,
-1, -1, 31, -1, -1, 29, 85,
-1, -1, 32, -1, -1, 30, 86,
-1, 2, 16, -1, -1, 17, 94,
-1, -1, 31, -1, -1, 29, 84,
-1, -1, 32, -1, -1, 30, 85,
-1, -1, 0, -1, -1, 1, -1,
-1, -1, 1, -1, -1, 0, -1,
-1, 2, 16, -1, -1, 17, 95,
-1, -1, 17, -1, 2, 16, -1,
-1, 2, 17, -1, -1, 18, 96,
-1, 2, 17, -1, -1, 18, 97,
-1, -1, 18, -1, 2, 17, -1,
-1, 2, 18, -1, -1, 19, 98,
-1, 2, 18, -1, -1, 19, 99,
-1, -1, 19, -1, 2, 18, -1,
-1, 2, 19, -1, -1, 20, 100,
-1, 2, 19, -1, -1, 20, 101,
-1, -1, 20, -1, 2, 19, -1,
-1, -1, 4, -1, 1, 3, -1,
-1, -1, 6, -1, 2, 3, -1,
Expand Down
36 changes: 28 additions & 8 deletions networks/He-C-Fe-group/actual_rhs.H
View file
Original file line number Diff line number Diff line change
Expand Up @@ -1349,6 +1349,22 @@ void evaluate_rates(const burn_t& state, T& rate_eval) {
}
rate_eval.add_energy_rate(k_ni57_to_co57) = edot_nu + edot_gamma;

tabular_evaluate(j_n_p_meta, j_n_p_rhoy, j_n_p_temp, j_n_p_data,
rhoy, state.T, rate, drate_dt, edot_nu, edot_gamma);
rate_eval.screened_rates(k_n_to_p) = rate;
if constexpr (std::is_same<T, rate_derivs_t>::value) {
rate_eval.dscreened_rates_dT(k_n_to_p) = drate_dt;
}
rate_eval.add_energy_rate(k_n_to_p) = edot_nu + edot_gamma;

tabular_evaluate(j_p_n_meta, j_p_n_rhoy, j_p_n_temp, j_p_n_data,
rhoy, state.T, rate, drate_dt, edot_nu, edot_gamma);
rate_eval.screened_rates(k_p_to_n) = rate;
if constexpr (std::is_same<T, rate_derivs_t>::value) {
rate_eval.dscreened_rates_dT(k_p_to_n) = drate_dt;
}
rate_eval.add_energy_rate(k_p_to_n) = edot_nu + edot_gamma;


}

Expand All @@ -1361,7 +1377,6 @@ void rhs_nuc(const burn_t& state,
using namespace Rates;

ydot_nuc(N) =
-screened_rates(k_n_to_p_weak_wc12)*Y(N) +
(-screened_rates(k_n_fe52_to_fe53)*Y(Fe52)*Y(N)*state.rho + screened_rates(k_fe53_to_n_fe52_derived)*Y(Fe53)) +
(-screened_rates(k_n_fe53_to_fe54)*Y(Fe53)*Y(N)*state.rho + screened_rates(k_fe54_to_n_fe53_derived)*Y(Fe54)) +
(-screened_rates(k_n_fe54_to_fe55)*Y(Fe54)*Y(N)*state.rho + screened_rates(k_fe55_to_n_fe54_derived)*Y(Fe55)) +
Expand All @@ -1379,7 +1394,9 @@ void rhs_nuc(const burn_t& state,
(-screened_rates(k_n_ni57_to_he4_fe54)*Y(N)*Y(Ni57)*state.rho + screened_rates(k_he4_fe54_to_n_ni57_derived)*Y(Fe54)*Y(He4)*state.rho) +
(-screened_rates(k_n_ni58_to_he4_fe55)*Y(N)*Y(Ni58)*state.rho + screened_rates(k_he4_fe55_to_n_ni58_derived)*Y(Fe55)*Y(He4)*state.rho) +
(-screened_rates(k_n_cu59_to_he4_co56)*Y(Cu59)*Y(N)*state.rho + screened_rates(k_he4_co56_to_n_cu59_derived)*Y(Co56)*Y(He4)*state.rho) +
(-screened_rates(k_n_zn60_to_he4_ni57)*Y(N)*Y(Zn60)*state.rho + screened_rates(k_he4_ni57_to_n_zn60_derived)*Y(He4)*Y(Ni57)*state.rho);
(-screened_rates(k_n_zn60_to_he4_ni57)*Y(N)*Y(Zn60)*state.rho + screened_rates(k_he4_ni57_to_n_zn60_derived)*Y(He4)*Y(Ni57)*state.rho) +
-screened_rates(k_n_to_p)*Y(N) +
screened_rates(k_p_to_n)*Y(H1);

ydot_nuc(H1) =
(-screened_rates(k_p_c12_to_n13)*Y(C12)*Y(H1)*state.rho + screened_rates(k_n13_to_p_c12_derived)*Y(N13)) +
Expand All @@ -1401,7 +1418,6 @@ void rhs_nuc(const burn_t& state,
(screened_rates(k_he4_cr48_to_p_mn51)*Y(Cr48)*Y(He4)*state.rho + -screened_rates(k_p_mn51_to_he4_cr48_derived)*Y(Mn51)*Y(H1)*state.rho) +
(screened_rates(k_he4_fe52_to_p_co55)*Y(Fe52)*Y(He4)*state.rho + -screened_rates(k_p_co55_to_he4_fe52_derived)*Y(Co55)*Y(H1)*state.rho) +
(-screened_rates(k_p_cu59_to_he4_ni56)*Y(Cu59)*Y(H1)*state.rho + screened_rates(k_he4_ni56_to_p_cu59_derived)*Y(He4)*Y(Ni56)*state.rho) +
screened_rates(k_n_to_p_weak_wc12)*Y(N) +
(-screened_rates(k_p_mn55_to_fe56)*Y(Mn55)*Y(H1)*state.rho + screened_rates(k_fe56_to_p_mn55_derived)*Y(Fe56)) +
(-screened_rates(k_p_fe54_to_co55)*Y(Fe54)*Y(H1)*state.rho + screened_rates(k_co55_to_p_fe54_derived)*Y(Co55)) +
(-screened_rates(k_p_fe55_to_co56)*Y(Fe55)*Y(H1)*state.rho + screened_rates(k_co56_to_p_fe55_derived)*Y(Co56)) +
Expand All @@ -1417,7 +1433,9 @@ void rhs_nuc(const burn_t& state,
(screened_rates(k_n_co56_to_p_fe56)*Y(Co56)*Y(N)*state.rho + -screened_rates(k_p_fe56_to_n_co56_derived)*Y(Fe56)*Y(H1)*state.rho) +
(-screened_rates(k_p_co57_to_he4_fe54)*Y(Co57)*Y(H1)*state.rho + screened_rates(k_he4_fe54_to_p_co57_derived)*Y(Fe54)*Y(He4)*state.rho) +
(screened_rates(k_n_ni56_to_p_co56)*Y(N)*Y(Ni56)*state.rho + -screened_rates(k_p_co56_to_n_ni56_derived)*Y(Co56)*Y(H1)*state.rho) +
(screened_rates(k_n_ni57_to_p_co57)*Y(N)*Y(Ni57)*state.rho + -screened_rates(k_p_co57_to_n_ni57_derived)*Y(Co57)*Y(H1)*state.rho);
(screened_rates(k_n_ni57_to_p_co57)*Y(N)*Y(Ni57)*state.rho + -screened_rates(k_p_co57_to_n_ni57_derived)*Y(Co57)*Y(H1)*state.rho) +
screened_rates(k_n_to_p)*Y(N) +
-screened_rates(k_p_to_n)*Y(H1);

ydot_nuc(He4) =
(-screened_rates(k_he4_c12_to_o16)*Y(C12)*Y(He4)*state.rho + screened_rates(k_o16_to_he4_c12_derived)*Y(O16)) +
Expand Down Expand Up @@ -1726,6 +1744,8 @@ void actual_rhs (burn_t& state, Array1D<Real, 1, neqs>& ydot)
enuc += C::Legacy::n_A * Y(Mn55) * rate_eval.add_energy_rate(k_mn55_to_fe55);
enuc += C::Legacy::n_A * Y(Ni56) * rate_eval.add_energy_rate(k_ni56_to_co56);
enuc += C::Legacy::n_A * Y(Ni57) * rate_eval.add_energy_rate(k_ni57_to_co57);
enuc += C::Legacy::n_A * Y(N) * rate_eval.add_energy_rate(k_n_to_p);
enuc += C::Legacy::n_A * Y(H1) * rate_eval.add_energy_rate(k_p_to_n);

// Get the thermal neutrino losses

Expand All @@ -1750,10 +1770,10 @@ void jac_nuc(const burn_t& state,

Real scratch;

scratch = -screened_rates(k_n_co55_to_co56)*Y(Co55)*state.rho - screened_rates(k_n_co55_to_p_fe55)*Y(Co55)*state.rho - screened_rates(k_n_co56_to_co57)*Y(Co56)*state.rho - screened_rates(k_n_co56_to_p_fe56)*Y(Co56)*state.rho - screened_rates(k_n_cu59_to_he4_co56)*Y(Cu59)*state.rho - screened_rates(k_n_fe52_to_fe53)*Y(Fe52)*state.rho - screened_rates(k_n_fe53_to_fe54)*Y(Fe53)*state.rho - screened_rates(k_n_fe54_to_fe55)*Y(Fe54)*state.rho - screened_rates(k_n_fe55_to_fe56)*Y(Fe55)*state.rho - screened_rates(k_n_fe55_to_p_mn55)*Y(Fe55)*state.rho - screened_rates(k_n_ni56_to_he4_fe53)*Y(Ni56)*state.rho - screened_rates(k_n_ni56_to_ni57)*Y(Ni56)*state.rho - screened_rates(k_n_ni56_to_p_co56)*Y(Ni56)*state.rho - screened_rates(k_n_ni57_to_he4_fe54)*Y(Ni57)*state.rho - screened_rates(k_n_ni57_to_ni58)*Y(Ni57)*state.rho - screened_rates(k_n_ni57_to_p_co57)*Y(Ni57)*state.rho - screened_rates(k_n_ni58_to_he4_fe55)*Y(Ni58)*state.rho - screened_rates(k_n_to_p_weak_wc12) - screened_rates(k_n_zn60_to_he4_ni57)*Y(Zn60)*state.rho;
scratch = -screened_rates(k_n_co55_to_co56)*Y(Co55)*state.rho - screened_rates(k_n_co55_to_p_fe55)*Y(Co55)*state.rho - screened_rates(k_n_co56_to_co57)*Y(Co56)*state.rho - screened_rates(k_n_co56_to_p_fe56)*Y(Co56)*state.rho - screened_rates(k_n_cu59_to_he4_co56)*Y(Cu59)*state.rho - screened_rates(k_n_fe52_to_fe53)*Y(Fe52)*state.rho - screened_rates(k_n_fe53_to_fe54)*Y(Fe53)*state.rho - screened_rates(k_n_fe54_to_fe55)*Y(Fe54)*state.rho - screened_rates(k_n_fe55_to_fe56)*Y(Fe55)*state.rho - screened_rates(k_n_fe55_to_p_mn55)*Y(Fe55)*state.rho - screened_rates(k_n_ni56_to_he4_fe53)*Y(Ni56)*state.rho - screened_rates(k_n_ni56_to_ni57)*Y(Ni56)*state.rho - screened_rates(k_n_ni56_to_p_co56)*Y(Ni56)*state.rho - screened_rates(k_n_ni57_to_he4_fe54)*Y(Ni57)*state.rho - screened_rates(k_n_ni57_to_ni58)*Y(Ni57)*state.rho - screened_rates(k_n_ni57_to_p_co57)*Y(Ni57)*state.rho - screened_rates(k_n_ni58_to_he4_fe55)*Y(Ni58)*state.rho - screened_rates(k_n_to_p) - screened_rates(k_n_zn60_to_he4_ni57)*Y(Zn60)*state.rho;
jac.set(N, N, scratch);

scratch = screened_rates(k_p_co56_to_n_ni56_derived)*Y(Co56)*state.rho + screened_rates(k_p_co57_to_n_ni57_derived)*Y(Co57)*state.rho + screened_rates(k_p_fe55_to_n_co55_derived)*Y(Fe55)*state.rho + screened_rates(k_p_fe56_to_n_co56_derived)*Y(Fe56)*state.rho + screened_rates(k_p_mn55_to_n_fe55_derived)*Y(Mn55)*state.rho;
scratch = screened_rates(k_p_co56_to_n_ni56_derived)*Y(Co56)*state.rho + screened_rates(k_p_co57_to_n_ni57_derived)*Y(Co57)*state.rho + screened_rates(k_p_fe55_to_n_co55_derived)*Y(Fe55)*state.rho + screened_rates(k_p_fe56_to_n_co56_derived)*Y(Fe56)*state.rho + screened_rates(k_p_mn55_to_n_fe55_derived)*Y(Mn55)*state.rho + screened_rates(k_p_to_n);
jac.set(N, H1, scratch);

scratch = screened_rates(k_he4_co56_to_n_cu59_derived)*Y(Co56)*state.rho + screened_rates(k_he4_fe53_to_n_ni56_derived)*Y(Fe53)*state.rho + screened_rates(k_he4_fe54_to_n_ni57_derived)*Y(Fe54)*state.rho + screened_rates(k_he4_fe55_to_n_ni58_derived)*Y(Fe55)*state.rho + screened_rates(k_he4_ni57_to_n_zn60_derived)*Y(Ni57)*state.rho;
Expand Down Expand Up @@ -1801,10 +1821,10 @@ void jac_nuc(const burn_t& state,
scratch = -screened_rates(k_n_zn60_to_he4_ni57)*Y(N)*state.rho;
jac.set(N, Zn60, scratch);

scratch = screened_rates(k_n_co55_to_p_fe55)*Y(Co55)*state.rho + screened_rates(k_n_co56_to_p_fe56)*Y(Co56)*state.rho + screened_rates(k_n_fe55_to_p_mn55)*Y(Fe55)*state.rho + screened_rates(k_n_ni56_to_p_co56)*Y(Ni56)*state.rho + screened_rates(k_n_ni57_to_p_co57)*Y(Ni57)*state.rho + screened_rates(k_n_to_p_weak_wc12);
scratch = screened_rates(k_n_co55_to_p_fe55)*Y(Co55)*state.rho + screened_rates(k_n_co56_to_p_fe56)*Y(Co56)*state.rho + screened_rates(k_n_fe55_to_p_mn55)*Y(Fe55)*state.rho + screened_rates(k_n_ni56_to_p_co56)*Y(Ni56)*state.rho + screened_rates(k_n_ni57_to_p_co57)*Y(Ni57)*state.rho + screened_rates(k_n_to_p);
jac.set(H1, N, scratch);

scratch = -screened_rates(k_p_al27_to_he4_mg24)*Y(Al27)*state.rho - screened_rates(k_p_al27_to_si28)*Y(Al27)*state.rho - screened_rates(k_p_c12_to_n13)*Y(C12)*state.rho - screened_rates(k_p_co55_to_he4_fe52_derived)*Y(Co55)*state.rho - screened_rates(k_p_co55_to_ni56)*Y(Co55)*state.rho - screened_rates(k_p_co56_to_he4_fe53_derived)*Y(Co56)*state.rho - screened_rates(k_p_co56_to_n_ni56_derived)*Y(Co56)*state.rho - screened_rates(k_p_co56_to_ni57)*Y(Co56)*state.rho - screened_rates(k_p_co57_to_he4_fe54)*Y(Co57)*state.rho - screened_rates(k_p_co57_to_n_ni57_derived)*Y(Co57)*state.rho - screened_rates(k_p_co57_to_ni58)*Y(Co57)*state.rho - screened_rates(k_p_cu59_to_he4_ni56)*Y(Cu59)*state.rho - screened_rates(k_p_cu59_to_zn60)*Y(Cu59)*state.rho - screened_rates(k_p_fe54_to_co55)*Y(Fe54)*state.rho - screened_rates(k_p_fe54_to_he4_mn51_derived)*Y(Fe54)*state.rho - screened_rates(k_p_fe55_to_co56)*Y(Fe55)*state.rho - screened_rates(k_p_fe55_to_n_co55_derived)*Y(Fe55)*state.rho - screened_rates(k_p_fe56_to_co57)*Y(Fe56)*state.rho - screened_rates(k_p_fe56_to_n_co56_derived)*Y(Fe56)*state.rho - screened_rates(k_p_mn51_to_fe52)*Y(Mn51)*state.rho - screened_rates(k_p_mn51_to_he4_cr48_derived)*Y(Mn51)*state.rho - screened_rates(k_p_mn55_to_fe56)*Y(Mn55)*state.rho - screened_rates(k_p_mn55_to_n_fe55_derived)*Y(Mn55)*state.rho - screened_rates(k_p_na23_to_he4_ne20)*Y(Na23)*state.rho - screened_rates(k_p_na23_to_mg24)*Y(Na23)*state.rho - screened_rates(k_p_ne21_to_he4_f18_derived)*Y(Ne21)*state.rho - screened_rates(k_p_ne21_to_na22)*Y(Ne21)*state.rho - screened_rates(k_p_ni58_to_cu59)*Y(Ni58)*state.rho - screened_rates(k_p_ni58_to_he4_co55_derived)*Y(Ni58)*state.rho - screened_rates(k_p_o16_to_he4_n13_derived)*Y(O16)*state.rho - screened_rates(k_p_p31_to_he4_si28)*Y(P31)*state.rho - screened_rates(k_p_p31_to_s32)*Y(P31)*state.rho;
scratch = -screened_rates(k_p_al27_to_he4_mg24)*Y(Al27)*state.rho - screened_rates(k_p_al27_to_si28)*Y(Al27)*state.rho - screened_rates(k_p_c12_to_n13)*Y(C12)*state.rho - screened_rates(k_p_co55_to_he4_fe52_derived)*Y(Co55)*state.rho - screened_rates(k_p_co55_to_ni56)*Y(Co55)*state.rho - screened_rates(k_p_co56_to_he4_fe53_derived)*Y(Co56)*state.rho - screened_rates(k_p_co56_to_n_ni56_derived)*Y(Co56)*state.rho - screened_rates(k_p_co56_to_ni57)*Y(Co56)*state.rho - screened_rates(k_p_co57_to_he4_fe54)*Y(Co57)*state.rho - screened_rates(k_p_co57_to_n_ni57_derived)*Y(Co57)*state.rho - screened_rates(k_p_co57_to_ni58)*Y(Co57)*state.rho - screened_rates(k_p_cu59_to_he4_ni56)*Y(Cu59)*state.rho - screened_rates(k_p_cu59_to_zn60)*Y(Cu59)*state.rho - screened_rates(k_p_fe54_to_co55)*Y(Fe54)*state.rho - screened_rates(k_p_fe54_to_he4_mn51_derived)*Y(Fe54)*state.rho - screened_rates(k_p_fe55_to_co56)*Y(Fe55)*state.rho - screened_rates(k_p_fe55_to_n_co55_derived)*Y(Fe55)*state.rho - screened_rates(k_p_fe56_to_co57)*Y(Fe56)*state.rho - screened_rates(k_p_fe56_to_n_co56_derived)*Y(Fe56)*state.rho - screened_rates(k_p_mn51_to_fe52)*Y(Mn51)*state.rho - screened_rates(k_p_mn51_to_he4_cr48_derived)*Y(Mn51)*state.rho - screened_rates(k_p_mn55_to_fe56)*Y(Mn55)*state.rho - screened_rates(k_p_mn55_to_n_fe55_derived)*Y(Mn55)*state.rho - screened_rates(k_p_na23_to_he4_ne20)*Y(Na23)*state.rho - screened_rates(k_p_na23_to_mg24)*Y(Na23)*state.rho - screened_rates(k_p_ne21_to_he4_f18_derived)*Y(Ne21)*state.rho - screened_rates(k_p_ne21_to_na22)*Y(Ne21)*state.rho - screened_rates(k_p_ni58_to_cu59)*Y(Ni58)*state.rho - screened_rates(k_p_ni58_to_he4_co55_derived)*Y(Ni58)*state.rho - screened_rates(k_p_o16_to_he4_n13_derived)*Y(O16)*state.rho - screened_rates(k_p_p31_to_he4_si28)*Y(P31)*state.rho - screened_rates(k_p_p31_to_s32)*Y(P31)*state.rho - screened_rates(k_p_to_n);
jac.set(H1, H1, scratch);

scratch = screened_rates(k_he4_co55_to_p_ni58)*Y(Co55)*state.rho + screened_rates(k_he4_cr48_to_p_mn51)*Y(Cr48)*state.rho + screened_rates(k_he4_f18_to_p_ne21)*Y(F18)*state.rho + screened_rates(k_he4_fe52_to_p_co55)*Y(Fe52)*state.rho + screened_rates(k_he4_fe53_to_p_co56)*Y(Fe53)*state.rho + screened_rates(k_he4_fe54_to_p_co57_derived)*Y(Fe54)*state.rho + screened_rates(k_he4_mg24_to_p_al27_derived)*Y(Mg24)*state.rho + screened_rates(k_he4_mn51_to_p_fe54)*Y(Mn51)*state.rho + screened_rates(k_he4_n13_to_p_o16)*Y(N13)*state.rho + screened_rates(k_he4_ne20_to_p_na23_derived)*Y(Ne20)*state.rho + screened_rates(k_he4_ni56_to_p_cu59_derived)*Y(Ni56)*state.rho + screened_rates(k_he4_si28_to_p_p31_derived)*Y(Si28)*state.rho;
Expand Down
Loading

0 comments on commit 427bb0f

Please sign in to comment.