From a9f3f1617934b55272c1a69c6e50cfe9a9db6a28 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Tue, 10 Dec 2024 20:36:27 -0500 Subject: [PATCH] use nse protons --- .../he-burn/he-burn-31anp/actual_network.H | 101 ++-- .../he-burn-31anp/actual_network_data.cpp | 194 ++++---- networks/he-burn/he-burn-31anp/actual_rhs.H | 443 +++++++++--------- .../he-burn/he-burn-31anp/he-burn-31anp.py | 3 + .../he-burn-31anp/inputs.burn_cell.VODE | 1 + .../he-burn-31anp/partition_functions.H | 1 + networks/he-burn/he-burn-31anp/pynucastro.net | 1 + .../he-burn/he-burn-31anp/reaclib_rates.H | 384 +++++++-------- 8 files changed, 578 insertions(+), 550 deletions(-) diff --git a/networks/he-burn/he-burn-31anp/actual_network.H b/networks/he-burn/he-burn-31anp/actual_network.H index a19a20f736..2db369d324 100644 --- a/networks/he-burn/he-burn-31anp/actual_network.H +++ b/networks/he-burn/he-burn-31anp/actual_network.H @@ -32,6 +32,9 @@ namespace network else if constexpr (spec == H1) { return 0.0_rt; } + else if constexpr (spec == P_nse) { + return 0.0_rt; + } else if constexpr (spec == He4) { return 28.295662457999697_rt; } @@ -185,10 +188,10 @@ namespace Rates k_He4_Si28_to_S32 = 11, k_p_P31_to_S32 = 12, k_He4_Cr48_to_Fe52 = 13, - k_p_Mn51_to_Fe52 = 14, + k_p_nse_Mn51_to_Fe52 = 14, k_He4_Mn51_to_Co55 = 15, k_He4_Fe52_to_Ni56 = 16, - k_p_Co55_to_Ni56 = 17, + k_p_nse_Co55_to_Ni56 = 17, k_C12_C12_to_p_Na23 = 18, k_C12_C12_to_He4_Ne20 = 19, k_He4_N13_to_p_O16 = 20, @@ -200,23 +203,23 @@ namespace Rates k_p_Na23_to_He4_Ne20 = 26, k_p_Al27_to_He4_Mg24 = 27, k_p_P31_to_He4_Si28 = 28, - k_He4_Cr48_to_p_Mn51 = 29, - k_He4_Fe52_to_p_Co55 = 30, + k_He4_Cr48_to_p_nse_Mn51 = 29, + k_He4_Fe52_to_p_nse_Co55 = 30, k_He4_He4_He4_to_C12 = 31, k_C12_C12_to_Mg24_modified = 32, k_O16_O16_to_S32_modified = 33, k_C12_O16_to_Si28_modified = 34, - k_p_Fe54_to_Co55 = 35, + k_p_nse_Fe54_to_Co55 = 35, k_He4_Fe54_to_Ni58 = 36, - k_p_Fe56_to_Co57 = 37, + k_p_nse_Fe56_to_Co57 = 37, k_n_Co55_to_Co56 = 38, k_n_Co56_to_Co57 = 39, - k_p_Co57_to_Ni58 = 40, - k_He4_Mn51_to_p_Fe54 = 41, - k_He4_Co55_to_p_Ni58 = 42, - k_n_Co56_to_p_Fe56 = 43, - k_p_Co57_to_He4_Fe54 = 44, - k_n_Ni56_to_p_Co56 = 45, + k_p_nse_Co57_to_Ni58 = 40, + k_He4_Mn51_to_p_nse_Fe54 = 41, + k_He4_Co55_to_p_nse_Ni58 = 42, + k_n_Co56_to_p_nse_Fe56 = 43, + k_p_nse_Co57_to_He4_Fe54 = 44, + k_n_Ni56_to_p_nse_Co56 = 45, k_He4_S32_to_Ar36_removed = 46, k_p_Cl35_to_Ar36_removed = 47, k_p_Cl35_to_He4_S32_removed = 48, @@ -267,10 +270,10 @@ namespace Rates k_Si28_to_He4_Mg24_derived = 93, k_S32_to_p_P31_derived = 94, k_S32_to_He4_Si28_derived = 95, - k_Fe52_to_p_Mn51_derived = 96, + k_Fe52_to_p_nse_Mn51_derived = 96, k_Fe52_to_He4_Cr48_derived = 97, k_Co55_to_He4_Mn51_derived = 98, - k_Ni56_to_p_Co55_derived = 99, + k_Ni56_to_p_nse_Co55_derived = 99, k_Ni56_to_He4_Fe52_derived = 100, k_C12_to_He4_He4_He4_derived = 101, k_p_O16_to_He4_N13_derived = 102, @@ -278,19 +281,19 @@ namespace Rates k_p_Ne21_to_He4_F18_derived = 104, k_He4_Mg24_to_p_Al27_derived = 105, k_He4_Si28_to_p_P31_derived = 106, - k_p_Mn51_to_He4_Cr48_derived = 107, - k_p_Co55_to_He4_Fe52_derived = 108, - k_Co55_to_p_Fe54_derived = 109, + k_p_nse_Mn51_to_He4_Cr48_derived = 107, + k_p_nse_Co55_to_He4_Fe52_derived = 108, + k_Co55_to_p_nse_Fe54_derived = 109, k_Co56_to_n_Co55_derived = 110, k_Co57_to_n_Co56_derived = 111, - k_Co57_to_p_Fe56_derived = 112, - k_Ni58_to_p_Co57_derived = 113, + k_Co57_to_p_nse_Fe56_derived = 112, + k_Ni58_to_p_nse_Co57_derived = 113, k_Ni58_to_He4_Fe54_derived = 114, - k_p_Fe54_to_He4_Mn51_derived = 115, - k_He4_Fe54_to_p_Co57_derived = 116, - k_p_Fe56_to_n_Co56_derived = 117, - k_p_Co56_to_n_Ni56_derived = 118, - k_p_Ni58_to_He4_Co55_derived = 119, + k_p_nse_Fe54_to_He4_Mn51_derived = 115, + k_He4_Fe54_to_p_nse_Co57_derived = 116, + k_p_nse_Fe56_to_n_Co56_derived = 117, + k_p_nse_Co56_to_n_Ni56_derived = 118, + k_p_nse_Ni58_to_He4_Co55_derived = 119, k_He4_S32_to_p_Cl35_derived_removed = 120, k_Ar36_to_He4_S32_derived_removed = 121, k_Ar36_to_p_Cl35_derived_removed = 122, @@ -339,10 +342,10 @@ namespace Rates "He4_Si28_to_S32", // 11, "p_P31_to_S32", // 12, "He4_Cr48_to_Fe52", // 13, - "p_Mn51_to_Fe52", // 14, + "p_nse_Mn51_to_Fe52", // 14, "He4_Mn51_to_Co55", // 15, "He4_Fe52_to_Ni56", // 16, - "p_Co55_to_Ni56", // 17, + "p_nse_Co55_to_Ni56", // 17, "C12_C12_to_p_Na23", // 18, "C12_C12_to_He4_Ne20", // 19, "He4_N13_to_p_O16", // 20, @@ -354,23 +357,23 @@ namespace Rates "p_Na23_to_He4_Ne20", // 26, "p_Al27_to_He4_Mg24", // 27, "p_P31_to_He4_Si28", // 28, - "He4_Cr48_to_p_Mn51", // 29, - "He4_Fe52_to_p_Co55", // 30, + "He4_Cr48_to_p_nse_Mn51", // 29, + "He4_Fe52_to_p_nse_Co55", // 30, "He4_He4_He4_to_C12", // 31, "C12_C12_to_Mg24_modified", // 32, "O16_O16_to_S32_modified", // 33, "C12_O16_to_Si28_modified", // 34, - "p_Fe54_to_Co55", // 35, + "p_nse_Fe54_to_Co55", // 35, "He4_Fe54_to_Ni58", // 36, - "p_Fe56_to_Co57", // 37, + "p_nse_Fe56_to_Co57", // 37, "n_Co55_to_Co56", // 38, "n_Co56_to_Co57", // 39, - "p_Co57_to_Ni58", // 40, - "He4_Mn51_to_p_Fe54", // 41, - "He4_Co55_to_p_Ni58", // 42, - "n_Co56_to_p_Fe56", // 43, - "p_Co57_to_He4_Fe54", // 44, - "n_Ni56_to_p_Co56", // 45, + "p_nse_Co57_to_Ni58", // 40, + "He4_Mn51_to_p_nse_Fe54", // 41, + "He4_Co55_to_p_nse_Ni58", // 42, + "n_Co56_to_p_nse_Fe56", // 43, + "p_nse_Co57_to_He4_Fe54", // 44, + "n_Ni56_to_p_nse_Co56", // 45, "He4_S32_to_Ar36_removed", // 46, "p_Cl35_to_Ar36_removed", // 47, "p_Cl35_to_He4_S32_removed", // 48, @@ -421,10 +424,10 @@ namespace Rates "Si28_to_He4_Mg24_derived", // 93, "S32_to_p_P31_derived", // 94, "S32_to_He4_Si28_derived", // 95, - "Fe52_to_p_Mn51_derived", // 96, + "Fe52_to_p_nse_Mn51_derived", // 96, "Fe52_to_He4_Cr48_derived", // 97, "Co55_to_He4_Mn51_derived", // 98, - "Ni56_to_p_Co55_derived", // 99, + "Ni56_to_p_nse_Co55_derived", // 99, "Ni56_to_He4_Fe52_derived", // 100, "C12_to_He4_He4_He4_derived", // 101, "p_O16_to_He4_N13_derived", // 102, @@ -432,19 +435,19 @@ namespace Rates "p_Ne21_to_He4_F18_derived", // 104, "He4_Mg24_to_p_Al27_derived", // 105, "He4_Si28_to_p_P31_derived", // 106, - "p_Mn51_to_He4_Cr48_derived", // 107, - "p_Co55_to_He4_Fe52_derived", // 108, - "Co55_to_p_Fe54_derived", // 109, + "p_nse_Mn51_to_He4_Cr48_derived", // 107, + "p_nse_Co55_to_He4_Fe52_derived", // 108, + "Co55_to_p_nse_Fe54_derived", // 109, "Co56_to_n_Co55_derived", // 110, "Co57_to_n_Co56_derived", // 111, - "Co57_to_p_Fe56_derived", // 112, - "Ni58_to_p_Co57_derived", // 113, + "Co57_to_p_nse_Fe56_derived", // 112, + "Ni58_to_p_nse_Co57_derived", // 113, "Ni58_to_He4_Fe54_derived", // 114, - "p_Fe54_to_He4_Mn51_derived", // 115, - "He4_Fe54_to_p_Co57_derived", // 116, - "p_Fe56_to_n_Co56_derived", // 117, - "p_Co56_to_n_Ni56_derived", // 118, - "p_Ni58_to_He4_Co55_derived", // 119, + "p_nse_Fe54_to_He4_Mn51_derived", // 115, + "He4_Fe54_to_p_nse_Co57_derived", // 116, + "p_nse_Fe56_to_n_Co56_derived", // 117, + "p_nse_Co56_to_n_Ni56_derived", // 118, + "p_nse_Ni58_to_He4_Co55_derived", // 119, "He4_S32_to_p_Cl35_derived_removed", // 120, "Ar36_to_He4_S32_derived_removed", // 121, "Ar36_to_p_Cl35_derived_removed", // 122, @@ -472,7 +475,7 @@ namespace NSE_INDEX { constexpr int H1_index = 1; constexpr int N_index = 0; - constexpr int He4_index = 2; + constexpr int He4_index = 3; // Each row corresponds to the rate in NetworkRates enum // First 3 row indices for reactants, followed by 3 product indices diff --git a/networks/he-burn/he-burn-31anp/actual_network_data.cpp b/networks/he-burn/he-burn-31anp/actual_network_data.cpp index cd7f186007..602a8fcf18 100644 --- a/networks/he-burn/he-burn-31anp/actual_network_data.cpp +++ b/networks/he-burn/he-burn-31anp/actual_network_data.cpp @@ -5,51 +5,51 @@ namespace NSE_INDEX { AMREX_GPU_MANAGED amrex::Array2D rate_indices { - -1, 1, 3, -1, -1, 4, 84, - -1, 2, 3, -1, -1, 6, 85, - -1, 2, 5, -1, -1, 7, 86, - -1, 2, 6, -1, -1, 8, 87, - -1, 2, 7, -1, -1, 10, 89, - -1, 2, 8, -1, -1, 12, 91, - -1, 1, 9, -1, -1, 10, 88, - -1, 1, 11, -1, -1, 12, 90, - -1, 2, 12, -1, -1, 14, 93, - -1, 1, 13, -1, -1, 14, 92, - -1, 2, 14, -1, -1, 16, 95, - -1, 1, 15, -1, -1, 16, 94, - -1, 2, 20, -1, -1, 22, 97, - -1, 1, 21, -1, -1, 22, 96, - -1, 2, 21, -1, -1, 25, 98, - -1, 2, 22, -1, -1, 28, 100, - -1, 1, 25, -1, -1, 28, 99, - -1, 3, 3, -1, 1, 11, -1, - -1, 3, 3, -1, 2, 8, -1, - -1, 2, 4, -1, 1, 6, 102, - -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, 104, - -1, 1, 11, -1, 2, 8, 103, - -1, 1, 13, -1, 2, 12, 105, - -1, 1, 15, -1, 2, 14, 106, - -1, 2, 20, -1, 1, 21, 107, - -1, 2, 22, -1, 1, 25, 108, - 2, 2, 2, -1, -1, 3, 101, - -1, 3, 3, -1, -1, 12, -1, - -1, 6, 6, -1, -1, 16, -1, - -1, 3, 6, -1, -1, 14, -1, - -1, 1, 23, -1, -1, 25, 109, - -1, 2, 23, -1, -1, 29, 114, - -1, 1, 24, -1, -1, 27, 112, - -1, 0, 25, -1, -1, 26, 110, - -1, 0, 26, -1, -1, 27, 111, - -1, 1, 27, -1, -1, 29, 113, - -1, 2, 21, -1, 1, 23, 115, - -1, 2, 25, -1, 1, 29, 119, - -1, 0, 26, -1, 1, 24, 117, - -1, 1, 27, -1, 2, 23, 116, - -1, 0, 28, -1, 1, 26, 118, + -1, 1, 4, -1, -1, 5, 84, + -1, 3, 4, -1, -1, 7, 85, + -1, 3, 6, -1, -1, 8, 86, + -1, 3, 7, -1, -1, 9, 87, + -1, 3, 8, -1, -1, 11, 89, + -1, 3, 9, -1, -1, 13, 91, + -1, 1, 10, -1, -1, 11, 88, + -1, 1, 12, -1, -1, 13, 90, + -1, 3, 13, -1, -1, 15, 93, + -1, 1, 14, -1, -1, 15, 92, + -1, 3, 15, -1, -1, 17, 95, + -1, 1, 16, -1, -1, 17, 94, + -1, 3, 21, -1, -1, 23, 97, + -1, 2, 22, -1, -1, 23, 96, + -1, 3, 22, -1, -1, 26, 98, + -1, 3, 23, -1, -1, 29, 100, + -1, 2, 26, -1, -1, 29, 99, + -1, 4, 4, -1, 1, 12, -1, + -1, 4, 4, -1, 3, 9, -1, + -1, 3, 5, -1, 1, 7, 102, + -1, 4, 7, -1, 1, 14, -1, + -1, 4, 7, -1, 3, 13, -1, + -1, 7, 7, -1, 1, 16, -1, + -1, 7, 7, -1, 3, 15, -1, + -1, 3, 8, -1, 1, 10, 104, + -1, 1, 12, -1, 3, 9, 103, + -1, 1, 14, -1, 3, 13, 105, + -1, 1, 16, -1, 3, 15, 106, + -1, 3, 21, -1, 2, 22, 107, + -1, 3, 23, -1, 2, 26, 108, + 3, 3, 3, -1, -1, 4, 101, + -1, 4, 4, -1, -1, 13, -1, + -1, 7, 7, -1, -1, 17, -1, + -1, 4, 7, -1, -1, 15, -1, + -1, 2, 24, -1, -1, 26, 109, + -1, 3, 24, -1, -1, 30, 114, + -1, 2, 25, -1, -1, 28, 112, + -1, 0, 26, -1, -1, 27, 110, + -1, 0, 27, -1, -1, 28, 111, + -1, 2, 28, -1, -1, 30, 113, + -1, 3, 22, -1, 2, 24, 115, + -1, 3, 26, -1, 2, 30, 119, + -1, 0, 27, -1, 2, 25, 117, + -1, 2, 28, -1, 3, 24, 116, + -1, 0, 29, -1, 2, 27, 118, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, @@ -68,62 +68,62 @@ 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, 26, -1, -1, 24, 66, - -1, -1, 26, -1, -1, 28, -1, - -1, -1, 24, -1, -1, 26, -1, + -1, -1, 27, -1, -1, 25, 66, + -1, -1, 27, -1, -1, 29, -1, + -1, -1, 25, -1, -1, 27, -1, -1, -1, 0, -1, -1, 1, 69, - -1, -1, 28, -1, -1, 26, 65, + -1, -1, 29, -1, -1, 27, 65, -1, -1, 1, -1, -1, 0, -1, - -1, 2, 16, -1, -1, 17, 71, - -1, -1, 17, -1, 2, 16, -1, - -1, 2, 17, -1, -1, 18, 73, - -1, -1, 18, -1, 2, 17, -1, - -1, 2, 18, -1, -1, 19, 75, - -1, -1, 19, -1, 2, 18, -1, - -1, 2, 19, -1, -1, 20, 77, - -1, -1, 20, -1, 2, 19, -1, - 0, 0, 22, -1, -1, 23, 79, - -1, -1, 23, 0, 0, 22, -1, - 0, 0, 23, -1, -1, 24, 81, + -1, 3, 17, -1, -1, 18, 71, + -1, -1, 18, -1, 3, 17, -1, + -1, 3, 18, -1, -1, 19, 73, + -1, -1, 19, -1, 3, 18, -1, + -1, 3, 19, -1, -1, 20, 75, + -1, -1, 20, -1, 3, 19, -1, + -1, 3, 20, -1, -1, 21, 77, + -1, -1, 21, -1, 3, 20, -1, + 0, 0, 23, -1, -1, 24, 79, -1, -1, 24, 0, 0, 23, -1, - 0, 0, 28, -1, -1, 29, 83, - -1, -1, 29, 0, 0, 28, -1, - -1, -1, 4, -1, 1, 3, -1, - -1, -1, 6, -1, 2, 3, -1, - -1, -1, 7, -1, 2, 5, -1, - -1, -1, 8, -1, 2, 6, -1, - -1, -1, 10, -1, 1, 9, -1, - -1, -1, 10, -1, 2, 7, -1, - -1, -1, 12, -1, 1, 11, -1, - -1, -1, 12, -1, 2, 8, -1, - -1, -1, 14, -1, 1, 13, -1, - -1, -1, 14, -1, 2, 12, -1, - -1, -1, 16, -1, 1, 15, -1, - -1, -1, 16, -1, 2, 14, -1, - -1, -1, 22, -1, 1, 21, -1, - -1, -1, 22, -1, 2, 20, -1, - -1, -1, 25, -1, 2, 21, -1, - -1, -1, 28, -1, 1, 25, -1, - -1, -1, 28, -1, 2, 22, -1, - -1, -1, 3, 2, 2, 2, -1, - -1, 1, 6, -1, 2, 4, -1, - -1, 2, 8, -1, 1, 11, -1, - -1, 1, 9, -1, 2, 7, -1, - -1, 2, 12, -1, 1, 13, -1, - -1, 2, 14, -1, 1, 15, -1, - -1, 1, 21, -1, 2, 20, -1, - -1, 1, 25, -1, 2, 22, -1, - -1, -1, 25, -1, 1, 23, -1, - -1, -1, 26, -1, 0, 25, -1, + 0, 0, 24, -1, -1, 25, 81, + -1, -1, 25, 0, 0, 24, -1, + 0, 0, 29, -1, -1, 30, 83, + -1, -1, 30, 0, 0, 29, -1, + -1, -1, 5, -1, 1, 4, -1, + -1, -1, 7, -1, 3, 4, -1, + -1, -1, 8, -1, 3, 6, -1, + -1, -1, 9, -1, 3, 7, -1, + -1, -1, 11, -1, 1, 10, -1, + -1, -1, 11, -1, 3, 8, -1, + -1, -1, 13, -1, 1, 12, -1, + -1, -1, 13, -1, 3, 9, -1, + -1, -1, 15, -1, 1, 14, -1, + -1, -1, 15, -1, 3, 13, -1, + -1, -1, 17, -1, 1, 16, -1, + -1, -1, 17, -1, 3, 15, -1, + -1, -1, 23, -1, 2, 22, -1, + -1, -1, 23, -1, 3, 21, -1, + -1, -1, 26, -1, 3, 22, -1, + -1, -1, 29, -1, 2, 26, -1, + -1, -1, 29, -1, 3, 23, -1, + -1, -1, 4, 3, 3, 3, -1, + -1, 1, 7, -1, 3, 5, -1, + -1, 3, 9, -1, 1, 12, -1, + -1, 1, 10, -1, 3, 8, -1, + -1, 3, 13, -1, 1, 14, -1, + -1, 3, 15, -1, 1, 16, -1, + -1, 2, 22, -1, 3, 21, -1, + -1, 2, 26, -1, 3, 23, -1, + -1, -1, 26, -1, 2, 24, -1, -1, -1, 27, -1, 0, 26, -1, - -1, -1, 27, -1, 1, 24, -1, - -1, -1, 29, -1, 1, 27, -1, - -1, -1, 29, -1, 2, 23, -1, - -1, 1, 23, -1, 2, 21, -1, - -1, 2, 23, -1, 1, 27, -1, - -1, 1, 24, -1, 0, 26, -1, - -1, 1, 26, -1, 0, 28, -1, - -1, 1, 29, -1, 2, 25, -1, + -1, -1, 28, -1, 0, 27, -1, + -1, -1, 28, -1, 2, 25, -1, + -1, -1, 30, -1, 2, 28, -1, + -1, -1, 30, -1, 3, 24, -1, + -1, 2, 24, -1, 3, 22, -1, + -1, 3, 24, -1, 2, 28, -1, + -1, 2, 25, -1, 0, 27, -1, + -1, 2, 27, -1, 0, 29, -1, + -1, 2, 30, -1, 3, 26, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, diff --git a/networks/he-burn/he-burn-31anp/actual_rhs.H b/networks/he-burn/he-burn-31anp/actual_rhs.H index ef91969561..4917e7e02f 100644 --- a/networks/he-burn/he-burn-31anp/actual_rhs.H +++ b/networks/he-burn/he-burn-31anp/actual_rhs.H @@ -382,11 +382,11 @@ void evaluate_rates(const burn_t& state, T& rate_eval) { rate_eval.dscreened_rates_dT(k_He4_Cr48_to_Fe52) = ratraw * dscor_dt + dratraw_dT * scor; } - ratraw = rate_eval.screened_rates(k_He4_Cr48_to_p_Mn51); - rate_eval.screened_rates(k_He4_Cr48_to_p_Mn51) *= scor; + ratraw = rate_eval.screened_rates(k_He4_Cr48_to_p_nse_Mn51); + rate_eval.screened_rates(k_He4_Cr48_to_p_nse_Mn51) *= scor; if constexpr (std::is_same_v) { - dratraw_dT = rate_eval.dscreened_rates_dT(k_He4_Cr48_to_p_Mn51); - rate_eval.dscreened_rates_dT(k_He4_Cr48_to_p_Mn51) = ratraw * dscor_dt + dratraw_dT * scor; + dratraw_dT = rate_eval.dscreened_rates_dT(k_He4_Cr48_to_p_nse_Mn51); + rate_eval.dscreened_rates_dT(k_He4_Cr48_to_p_nse_Mn51) = ratraw * dscor_dt + dratraw_dT * scor; } @@ -401,18 +401,18 @@ void evaluate_rates(const burn_t& state, T& rate_eval) { } - ratraw = rate_eval.screened_rates(k_p_Mn51_to_Fe52); - rate_eval.screened_rates(k_p_Mn51_to_Fe52) *= scor; + ratraw = rate_eval.screened_rates(k_p_nse_Mn51_to_Fe52); + rate_eval.screened_rates(k_p_nse_Mn51_to_Fe52) *= scor; if constexpr (std::is_same_v) { - dratraw_dT = rate_eval.dscreened_rates_dT(k_p_Mn51_to_Fe52); - rate_eval.dscreened_rates_dT(k_p_Mn51_to_Fe52) = ratraw * dscor_dt + dratraw_dT * scor; + dratraw_dT = rate_eval.dscreened_rates_dT(k_p_nse_Mn51_to_Fe52); + rate_eval.dscreened_rates_dT(k_p_nse_Mn51_to_Fe52) = ratraw * dscor_dt + dratraw_dT * scor; } - ratraw = rate_eval.screened_rates(k_p_Mn51_to_He4_Cr48_derived); - rate_eval.screened_rates(k_p_Mn51_to_He4_Cr48_derived) *= scor; + ratraw = rate_eval.screened_rates(k_p_nse_Mn51_to_He4_Cr48_derived); + rate_eval.screened_rates(k_p_nse_Mn51_to_He4_Cr48_derived) *= scor; if constexpr (std::is_same_v) { - dratraw_dT = rate_eval.dscreened_rates_dT(k_p_Mn51_to_He4_Cr48_derived); - rate_eval.dscreened_rates_dT(k_p_Mn51_to_He4_Cr48_derived) = ratraw * dscor_dt + dratraw_dT * scor; + dratraw_dT = rate_eval.dscreened_rates_dT(k_p_nse_Mn51_to_He4_Cr48_derived); + rate_eval.dscreened_rates_dT(k_p_nse_Mn51_to_He4_Cr48_derived) = ratraw * dscor_dt + dratraw_dT * scor; } @@ -434,11 +434,11 @@ void evaluate_rates(const burn_t& state, T& rate_eval) { rate_eval.dscreened_rates_dT(k_He4_Mn51_to_Co55) = ratraw * dscor_dt + dratraw_dT * scor; } - ratraw = rate_eval.screened_rates(k_He4_Mn51_to_p_Fe54); - rate_eval.screened_rates(k_He4_Mn51_to_p_Fe54) *= scor; + ratraw = rate_eval.screened_rates(k_He4_Mn51_to_p_nse_Fe54); + rate_eval.screened_rates(k_He4_Mn51_to_p_nse_Fe54) *= scor; if constexpr (std::is_same_v) { - dratraw_dT = rate_eval.dscreened_rates_dT(k_He4_Mn51_to_p_Fe54); - rate_eval.dscreened_rates_dT(k_He4_Mn51_to_p_Fe54) = ratraw * dscor_dt + dratraw_dT * scor; + dratraw_dT = rate_eval.dscreened_rates_dT(k_He4_Mn51_to_p_nse_Fe54); + rate_eval.dscreened_rates_dT(k_He4_Mn51_to_p_nse_Fe54) = ratraw * dscor_dt + dratraw_dT * scor; } @@ -460,11 +460,11 @@ void evaluate_rates(const burn_t& state, T& rate_eval) { rate_eval.dscreened_rates_dT(k_He4_Fe52_to_Ni56) = ratraw * dscor_dt + dratraw_dT * scor; } - ratraw = rate_eval.screened_rates(k_He4_Fe52_to_p_Co55); - rate_eval.screened_rates(k_He4_Fe52_to_p_Co55) *= scor; + ratraw = rate_eval.screened_rates(k_He4_Fe52_to_p_nse_Co55); + rate_eval.screened_rates(k_He4_Fe52_to_p_nse_Co55) *= scor; if constexpr (std::is_same_v) { - dratraw_dT = rate_eval.dscreened_rates_dT(k_He4_Fe52_to_p_Co55); - rate_eval.dscreened_rates_dT(k_He4_Fe52_to_p_Co55) = ratraw * dscor_dt + dratraw_dT * scor; + dratraw_dT = rate_eval.dscreened_rates_dT(k_He4_Fe52_to_p_nse_Co55); + rate_eval.dscreened_rates_dT(k_He4_Fe52_to_p_nse_Co55) = ratraw * dscor_dt + dratraw_dT * scor; } @@ -479,18 +479,18 @@ void evaluate_rates(const burn_t& state, T& rate_eval) { } - ratraw = rate_eval.screened_rates(k_p_Co55_to_Ni56); - rate_eval.screened_rates(k_p_Co55_to_Ni56) *= scor; + ratraw = rate_eval.screened_rates(k_p_nse_Co55_to_Ni56); + rate_eval.screened_rates(k_p_nse_Co55_to_Ni56) *= scor; if constexpr (std::is_same_v) { - dratraw_dT = rate_eval.dscreened_rates_dT(k_p_Co55_to_Ni56); - rate_eval.dscreened_rates_dT(k_p_Co55_to_Ni56) = ratraw * dscor_dt + dratraw_dT * scor; + dratraw_dT = rate_eval.dscreened_rates_dT(k_p_nse_Co55_to_Ni56); + rate_eval.dscreened_rates_dT(k_p_nse_Co55_to_Ni56) = ratraw * dscor_dt + dratraw_dT * scor; } - ratraw = rate_eval.screened_rates(k_p_Co55_to_He4_Fe52_derived); - rate_eval.screened_rates(k_p_Co55_to_He4_Fe52_derived) *= scor; + ratraw = rate_eval.screened_rates(k_p_nse_Co55_to_He4_Fe52_derived); + rate_eval.screened_rates(k_p_nse_Co55_to_He4_Fe52_derived) *= scor; if constexpr (std::is_same_v) { - dratraw_dT = rate_eval.dscreened_rates_dT(k_p_Co55_to_He4_Fe52_derived); - rate_eval.dscreened_rates_dT(k_p_Co55_to_He4_Fe52_derived) = ratraw * dscor_dt + dratraw_dT * scor; + dratraw_dT = rate_eval.dscreened_rates_dT(k_p_nse_Co55_to_He4_Fe52_derived); + rate_eval.dscreened_rates_dT(k_p_nse_Co55_to_He4_Fe52_derived) = ratraw * dscor_dt + dratraw_dT * scor; } @@ -655,18 +655,18 @@ void evaluate_rates(const burn_t& state, T& rate_eval) { } - ratraw = rate_eval.screened_rates(k_p_Fe54_to_Co55); - rate_eval.screened_rates(k_p_Fe54_to_Co55) *= scor; + ratraw = rate_eval.screened_rates(k_p_nse_Fe54_to_Co55); + rate_eval.screened_rates(k_p_nse_Fe54_to_Co55) *= scor; if constexpr (std::is_same_v) { - dratraw_dT = rate_eval.dscreened_rates_dT(k_p_Fe54_to_Co55); - rate_eval.dscreened_rates_dT(k_p_Fe54_to_Co55) = ratraw * dscor_dt + dratraw_dT * scor; + dratraw_dT = rate_eval.dscreened_rates_dT(k_p_nse_Fe54_to_Co55); + rate_eval.dscreened_rates_dT(k_p_nse_Fe54_to_Co55) = ratraw * dscor_dt + dratraw_dT * scor; } - ratraw = rate_eval.screened_rates(k_p_Fe54_to_He4_Mn51_derived); - rate_eval.screened_rates(k_p_Fe54_to_He4_Mn51_derived) *= scor; + ratraw = rate_eval.screened_rates(k_p_nse_Fe54_to_He4_Mn51_derived); + rate_eval.screened_rates(k_p_nse_Fe54_to_He4_Mn51_derived) *= scor; if constexpr (std::is_same_v) { - dratraw_dT = rate_eval.dscreened_rates_dT(k_p_Fe54_to_He4_Mn51_derived); - rate_eval.dscreened_rates_dT(k_p_Fe54_to_He4_Mn51_derived) = ratraw * dscor_dt + dratraw_dT * scor; + dratraw_dT = rate_eval.dscreened_rates_dT(k_p_nse_Fe54_to_He4_Mn51_derived); + rate_eval.dscreened_rates_dT(k_p_nse_Fe54_to_He4_Mn51_derived) = ratraw * dscor_dt + dratraw_dT * scor; } @@ -688,11 +688,11 @@ void evaluate_rates(const burn_t& state, T& rate_eval) { rate_eval.dscreened_rates_dT(k_He4_Fe54_to_Ni58) = ratraw * dscor_dt + dratraw_dT * scor; } - ratraw = rate_eval.screened_rates(k_He4_Fe54_to_p_Co57_derived); - rate_eval.screened_rates(k_He4_Fe54_to_p_Co57_derived) *= scor; + ratraw = rate_eval.screened_rates(k_He4_Fe54_to_p_nse_Co57_derived); + rate_eval.screened_rates(k_He4_Fe54_to_p_nse_Co57_derived) *= scor; if constexpr (std::is_same_v) { - dratraw_dT = rate_eval.dscreened_rates_dT(k_He4_Fe54_to_p_Co57_derived); - rate_eval.dscreened_rates_dT(k_He4_Fe54_to_p_Co57_derived) = ratraw * dscor_dt + dratraw_dT * scor; + dratraw_dT = rate_eval.dscreened_rates_dT(k_He4_Fe54_to_p_nse_Co57_derived); + rate_eval.dscreened_rates_dT(k_He4_Fe54_to_p_nse_Co57_derived) = ratraw * dscor_dt + dratraw_dT * scor; } @@ -707,18 +707,18 @@ void evaluate_rates(const burn_t& state, T& rate_eval) { } - ratraw = rate_eval.screened_rates(k_p_Fe56_to_Co57); - rate_eval.screened_rates(k_p_Fe56_to_Co57) *= scor; + ratraw = rate_eval.screened_rates(k_p_nse_Fe56_to_Co57); + rate_eval.screened_rates(k_p_nse_Fe56_to_Co57) *= scor; if constexpr (std::is_same_v) { - dratraw_dT = rate_eval.dscreened_rates_dT(k_p_Fe56_to_Co57); - rate_eval.dscreened_rates_dT(k_p_Fe56_to_Co57) = ratraw * dscor_dt + dratraw_dT * scor; + dratraw_dT = rate_eval.dscreened_rates_dT(k_p_nse_Fe56_to_Co57); + rate_eval.dscreened_rates_dT(k_p_nse_Fe56_to_Co57) = ratraw * dscor_dt + dratraw_dT * scor; } - ratraw = rate_eval.screened_rates(k_p_Fe56_to_n_Co56_derived); - rate_eval.screened_rates(k_p_Fe56_to_n_Co56_derived) *= scor; + ratraw = rate_eval.screened_rates(k_p_nse_Fe56_to_n_Co56_derived); + rate_eval.screened_rates(k_p_nse_Fe56_to_n_Co56_derived) *= scor; if constexpr (std::is_same_v) { - dratraw_dT = rate_eval.dscreened_rates_dT(k_p_Fe56_to_n_Co56_derived); - rate_eval.dscreened_rates_dT(k_p_Fe56_to_n_Co56_derived) = ratraw * dscor_dt + dratraw_dT * scor; + dratraw_dT = rate_eval.dscreened_rates_dT(k_p_nse_Fe56_to_n_Co56_derived); + rate_eval.dscreened_rates_dT(k_p_nse_Fe56_to_n_Co56_derived) = ratraw * dscor_dt + dratraw_dT * scor; } @@ -733,18 +733,18 @@ void evaluate_rates(const burn_t& state, T& rate_eval) { } - ratraw = rate_eval.screened_rates(k_p_Co57_to_Ni58); - rate_eval.screened_rates(k_p_Co57_to_Ni58) *= scor; + ratraw = rate_eval.screened_rates(k_p_nse_Co57_to_Ni58); + rate_eval.screened_rates(k_p_nse_Co57_to_Ni58) *= scor; if constexpr (std::is_same_v) { - dratraw_dT = rate_eval.dscreened_rates_dT(k_p_Co57_to_Ni58); - rate_eval.dscreened_rates_dT(k_p_Co57_to_Ni58) = ratraw * dscor_dt + dratraw_dT * scor; + dratraw_dT = rate_eval.dscreened_rates_dT(k_p_nse_Co57_to_Ni58); + rate_eval.dscreened_rates_dT(k_p_nse_Co57_to_Ni58) = ratraw * dscor_dt + dratraw_dT * scor; } - ratraw = rate_eval.screened_rates(k_p_Co57_to_He4_Fe54); - rate_eval.screened_rates(k_p_Co57_to_He4_Fe54) *= scor; + ratraw = rate_eval.screened_rates(k_p_nse_Co57_to_He4_Fe54); + rate_eval.screened_rates(k_p_nse_Co57_to_He4_Fe54) *= scor; if constexpr (std::is_same_v) { - dratraw_dT = rate_eval.dscreened_rates_dT(k_p_Co57_to_He4_Fe54); - rate_eval.dscreened_rates_dT(k_p_Co57_to_He4_Fe54) = ratraw * dscor_dt + dratraw_dT * scor; + dratraw_dT = rate_eval.dscreened_rates_dT(k_p_nse_Co57_to_He4_Fe54); + rate_eval.dscreened_rates_dT(k_p_nse_Co57_to_He4_Fe54) = ratraw * dscor_dt + dratraw_dT * scor; } @@ -759,11 +759,11 @@ void evaluate_rates(const burn_t& state, T& rate_eval) { } - ratraw = rate_eval.screened_rates(k_He4_Co55_to_p_Ni58); - rate_eval.screened_rates(k_He4_Co55_to_p_Ni58) *= scor; + ratraw = rate_eval.screened_rates(k_He4_Co55_to_p_nse_Ni58); + rate_eval.screened_rates(k_He4_Co55_to_p_nse_Ni58) *= scor; if constexpr (std::is_same_v) { - dratraw_dT = rate_eval.dscreened_rates_dT(k_He4_Co55_to_p_Ni58); - rate_eval.dscreened_rates_dT(k_He4_Co55_to_p_Ni58) = ratraw * dscor_dt + dratraw_dT * scor; + dratraw_dT = rate_eval.dscreened_rates_dT(k_He4_Co55_to_p_nse_Ni58); + rate_eval.dscreened_rates_dT(k_He4_Co55_to_p_nse_Ni58) = ratraw * dscor_dt + dratraw_dT * scor; } @@ -797,11 +797,11 @@ void evaluate_rates(const burn_t& state, T& rate_eval) { } - ratraw = rate_eval.screened_rates(k_p_Co56_to_n_Ni56_derived); - rate_eval.screened_rates(k_p_Co56_to_n_Ni56_derived) *= scor; + ratraw = rate_eval.screened_rates(k_p_nse_Co56_to_n_Ni56_derived); + rate_eval.screened_rates(k_p_nse_Co56_to_n_Ni56_derived) *= scor; if constexpr (std::is_same_v) { - dratraw_dT = rate_eval.dscreened_rates_dT(k_p_Co56_to_n_Ni56_derived); - rate_eval.dscreened_rates_dT(k_p_Co56_to_n_Ni56_derived) = ratraw * dscor_dt + dratraw_dT * scor; + dratraw_dT = rate_eval.dscreened_rates_dT(k_p_nse_Co56_to_n_Ni56_derived); + rate_eval.dscreened_rates_dT(k_p_nse_Co56_to_n_Ni56_derived) = ratraw * dscor_dt + dratraw_dT * scor; } @@ -816,11 +816,11 @@ void evaluate_rates(const burn_t& state, T& rate_eval) { } - ratraw = rate_eval.screened_rates(k_p_Ni58_to_He4_Co55_derived); - rate_eval.screened_rates(k_p_Ni58_to_He4_Co55_derived) *= scor; + ratraw = rate_eval.screened_rates(k_p_nse_Ni58_to_He4_Co55_derived); + rate_eval.screened_rates(k_p_nse_Ni58_to_He4_Co55_derived) *= scor; if constexpr (std::is_same_v) { - dratraw_dT = rate_eval.dscreened_rates_dT(k_p_Ni58_to_He4_Co55_derived); - rate_eval.dscreened_rates_dT(k_p_Ni58_to_He4_Co55_derived) = ratraw * dscor_dt + dratraw_dT * scor; + dratraw_dT = rate_eval.dscreened_rates_dT(k_p_nse_Ni58_to_He4_Co55_derived); + rate_eval.dscreened_rates_dT(k_p_nse_Ni58_to_He4_Co55_derived) = ratraw * dscor_dt + dratraw_dT * scor; } @@ -1159,6 +1159,8 @@ void get_ydot_weak(const burn_t& state, ydot_nuc(H1) = (screened_rates(k_n_to_p)*Y(N) + -screened_rates(k_p_to_n)*Y(H1)); + ydot_nuc(P_nse) = 0.0_rt; + ydot_nuc(He4) = 0.0_rt; ydot_nuc(C12) = 0.0_rt; @@ -1235,8 +1237,8 @@ void rhs_nuc(const burn_t& state, ydot_nuc(N) = (-screened_rates(k_n_Co55_to_Co56)*Y(Co55)*Y(N)*state.rho + screened_rates(k_Co56_to_n_Co55_derived)*Y(Co56)) + (-screened_rates(k_n_Co56_to_Co57)*Y(Co56)*Y(N)*state.rho + screened_rates(k_Co57_to_n_Co56_derived)*Y(Co57)) + - (-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_n_Ni56_to_p_Co56)*Y(Ni56)*Y(N)*state.rho + screened_rates(k_p_Co56_to_n_Ni56_derived)*Y(Co56)*Y(H1)*state.rho) + + (-screened_rates(k_n_Co56_to_p_nse_Fe56)*Y(Co56)*Y(N)*state.rho + screened_rates(k_p_nse_Fe56_to_n_Co56_derived)*Y(Fe56)*Y(P_nse)*state.rho) + + (-screened_rates(k_n_Ni56_to_p_nse_Co56)*Y(Ni56)*Y(N)*state.rho + screened_rates(k_p_nse_Co56_to_n_Ni56_derived)*Y(Co56)*Y(P_nse)*state.rho) + (-2.0*screened_rates(k_Fe52_n_n_to_Fe54_approx)*Y(Fe52)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho) + 2.0*screened_rates(k_Fe54_to_Fe52_n_n_approx)*Y(Fe54)) + (-2.0*screened_rates(k_Fe54_n_n_to_Fe56_approx)*Y(Fe54)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho) + 2.0*screened_rates(k_Fe56_to_Fe54_n_n_approx)*Y(Fe56)) + (-2.0*screened_rates(k_Ni56_n_n_to_Ni58_approx)*Y(Ni56)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho) + 2.0*screened_rates(k_Ni58_to_Ni56_n_n_approx)*Y(Ni58)) + @@ -1248,8 +1250,6 @@ void rhs_nuc(const burn_t& state, (-screened_rates(k_p_Na23_to_Mg24)*Y(Na23)*Y(H1)*state.rho + screened_rates(k_Mg24_to_p_Na23_derived)*Y(Mg24)) + (-screened_rates(k_p_Al27_to_Si28)*Y(Al27)*Y(H1)*state.rho + screened_rates(k_Si28_to_p_Al27_derived)*Y(Si28)) + (-screened_rates(k_p_P31_to_S32)*Y(P31)*Y(H1)*state.rho + screened_rates(k_S32_to_p_P31_derived)*Y(S32)) + - (-screened_rates(k_p_Mn51_to_Fe52)*Y(Mn51)*Y(H1)*state.rho + screened_rates(k_Fe52_to_p_Mn51_derived)*Y(Fe52)) + - (-screened_rates(k_p_Co55_to_Ni56)*Y(Co55)*Y(H1)*state.rho + screened_rates(k_Ni56_to_p_Co55_derived)*Y(Ni56)) + 0.5*screened_rates(k_C12_C12_to_p_Na23)*amrex::Math::powi<2>(Y(C12))*state.rho + (screened_rates(k_He4_N13_to_p_O16)*Y(He4)*Y(N13)*state.rho + -screened_rates(k_p_O16_to_He4_N13_derived)*Y(O16)*Y(H1)*state.rho) + screened_rates(k_C12_O16_to_p_Al27)*Y(C12)*Y(O16)*state.rho + @@ -1258,18 +1258,22 @@ void rhs_nuc(const burn_t& state, (-screened_rates(k_p_Na23_to_He4_Ne20)*Y(Na23)*Y(H1)*state.rho + screened_rates(k_He4_Ne20_to_p_Na23_derived)*Y(He4)*Y(Ne20)*state.rho) + (-screened_rates(k_p_Al27_to_He4_Mg24)*Y(Al27)*Y(H1)*state.rho + screened_rates(k_He4_Mg24_to_p_Al27_derived)*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_derived)*Y(He4)*Y(Si28)*state.rho) + - (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_Fe54_to_Co55)*Y(Fe54)*Y(H1)*state.rho + screened_rates(k_Co55_to_p_Fe54_derived)*Y(Co55)) + - (-screened_rates(k_p_Fe56_to_Co57)*Y(Fe56)*Y(H1)*state.rho + screened_rates(k_Co57_to_p_Fe56_derived)*Y(Co57)) + - (-screened_rates(k_p_Co57_to_Ni58)*Y(Co57)*Y(H1)*state.rho + screened_rates(k_Ni58_to_p_Co57_derived)*Y(Ni58)) + - (screened_rates(k_He4_Mn51_to_p_Fe54)*Y(He4)*Y(Mn51)*state.rho + -screened_rates(k_p_Fe54_to_He4_Mn51_derived)*Y(Fe54)*Y(H1)*state.rho) + - (screened_rates(k_He4_Co55_to_p_Ni58)*Y(Co55)*Y(He4)*state.rho + -screened_rates(k_p_Ni58_to_He4_Co55_derived)*Y(Ni58)*Y(H1)*state.rho) + - (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(Ni56)*Y(N)*state.rho + -screened_rates(k_p_Co56_to_n_Ni56_derived)*Y(Co56)*Y(H1)*state.rho) + (screened_rates(k_n_to_p)*Y(N) + -screened_rates(k_p_to_n)*Y(H1)); + ydot_nuc(P_nse) = + (-screened_rates(k_p_nse_Mn51_to_Fe52)*Y(Mn51)*Y(P_nse)*state.rho + screened_rates(k_Fe52_to_p_nse_Mn51_derived)*Y(Fe52)) + + (-screened_rates(k_p_nse_Co55_to_Ni56)*Y(Co55)*Y(P_nse)*state.rho + screened_rates(k_Ni56_to_p_nse_Co55_derived)*Y(Ni56)) + + (screened_rates(k_He4_Cr48_to_p_nse_Mn51)*Y(Cr48)*Y(He4)*state.rho + -screened_rates(k_p_nse_Mn51_to_He4_Cr48_derived)*Y(Mn51)*Y(P_nse)*state.rho) + + (screened_rates(k_He4_Fe52_to_p_nse_Co55)*Y(Fe52)*Y(He4)*state.rho + -screened_rates(k_p_nse_Co55_to_He4_Fe52_derived)*Y(Co55)*Y(P_nse)*state.rho) + + (-screened_rates(k_p_nse_Fe54_to_Co55)*Y(Fe54)*Y(P_nse)*state.rho + screened_rates(k_Co55_to_p_nse_Fe54_derived)*Y(Co55)) + + (-screened_rates(k_p_nse_Fe56_to_Co57)*Y(Fe56)*Y(P_nse)*state.rho + screened_rates(k_Co57_to_p_nse_Fe56_derived)*Y(Co57)) + + (-screened_rates(k_p_nse_Co57_to_Ni58)*Y(Co57)*Y(P_nse)*state.rho + screened_rates(k_Ni58_to_p_nse_Co57_derived)*Y(Ni58)) + + (screened_rates(k_He4_Mn51_to_p_nse_Fe54)*Y(He4)*Y(Mn51)*state.rho + -screened_rates(k_p_nse_Fe54_to_He4_Mn51_derived)*Y(Fe54)*Y(P_nse)*state.rho) + + (screened_rates(k_He4_Co55_to_p_nse_Ni58)*Y(Co55)*Y(He4)*state.rho + -screened_rates(k_p_nse_Ni58_to_He4_Co55_derived)*Y(Ni58)*Y(P_nse)*state.rho) + + (screened_rates(k_n_Co56_to_p_nse_Fe56)*Y(Co56)*Y(N)*state.rho + -screened_rates(k_p_nse_Fe56_to_n_Co56_derived)*Y(Fe56)*Y(P_nse)*state.rho) + + (-screened_rates(k_p_nse_Co57_to_He4_Fe54)*Y(Co57)*Y(P_nse)*state.rho + screened_rates(k_He4_Fe54_to_p_nse_Co57_derived)*Y(Fe54)*Y(He4)*state.rho) + + (screened_rates(k_n_Ni56_to_p_nse_Co56)*Y(Ni56)*Y(N)*state.rho + -screened_rates(k_p_nse_Co56_to_n_Ni56_derived)*Y(Co56)*Y(P_nse)*state.rho); + 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)) + (-screened_rates(k_He4_N14_to_F18)*Y(He4)*Y(N14)*state.rho + screened_rates(k_F18_to_He4_N14_derived)*Y(F18)) + @@ -1289,13 +1293,13 @@ void rhs_nuc(const burn_t& state, (screened_rates(k_p_Na23_to_He4_Ne20)*Y(Na23)*Y(H1)*state.rho + -screened_rates(k_He4_Ne20_to_p_Na23_derived)*Y(He4)*Y(Ne20)*state.rho) + (screened_rates(k_p_Al27_to_He4_Mg24)*Y(Al27)*Y(H1)*state.rho + -screened_rates(k_He4_Mg24_to_p_Al27_derived)*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_derived)*Y(He4)*Y(Si28)*state.rho) + - (-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_He4_Cr48_to_p_nse_Mn51)*Y(Cr48)*Y(He4)*state.rho + screened_rates(k_p_nse_Mn51_to_He4_Cr48_derived)*Y(Mn51)*Y(P_nse)*state.rho) + + (-screened_rates(k_He4_Fe52_to_p_nse_Co55)*Y(Fe52)*Y(He4)*state.rho + screened_rates(k_p_nse_Co55_to_He4_Fe52_derived)*Y(Co55)*Y(P_nse)*state.rho) + (-0.5*screened_rates(k_He4_He4_He4_to_C12)*amrex::Math::powi<3>(Y(He4))*amrex::Math::powi<2>(state.rho) + 3.0*screened_rates(k_C12_to_He4_He4_He4_derived)*Y(C12)) + (-screened_rates(k_He4_Fe54_to_Ni58)*Y(Fe54)*Y(He4)*state.rho + screened_rates(k_Ni58_to_He4_Fe54_derived)*Y(Ni58)) + - (-screened_rates(k_He4_Mn51_to_p_Fe54)*Y(He4)*Y(Mn51)*state.rho + screened_rates(k_p_Fe54_to_He4_Mn51_derived)*Y(Fe54)*Y(H1)*state.rho) + - (-screened_rates(k_He4_Co55_to_p_Ni58)*Y(Co55)*Y(He4)*state.rho + screened_rates(k_p_Ni58_to_He4_Co55_derived)*Y(Ni58)*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_He4_Mn51_to_p_nse_Fe54)*Y(He4)*Y(Mn51)*state.rho + screened_rates(k_p_nse_Fe54_to_He4_Mn51_derived)*Y(Fe54)*Y(P_nse)*state.rho) + + (-screened_rates(k_He4_Co55_to_p_nse_Ni58)*Y(Co55)*Y(He4)*state.rho + screened_rates(k_p_nse_Ni58_to_He4_Co55_derived)*Y(Ni58)*Y(P_nse)*state.rho) + + (screened_rates(k_p_nse_Co57_to_He4_Fe54)*Y(Co57)*Y(P_nse)*state.rho + -screened_rates(k_He4_Fe54_to_p_nse_Co57_derived)*Y(Fe54)*Y(He4)*state.rho) + (-screened_rates(k_S32_He4_to_Ar36_approx)*Y(He4)*Y(S32)*state.rho + screened_rates(k_Ar36_to_S32_He4_approx)*Y(Ar36)) + (-screened_rates(k_Ar36_He4_to_Ca40_approx)*Y(Ar36)*Y(He4)*state.rho + screened_rates(k_Ca40_to_Ar36_He4_approx)*Y(Ca40)) + (-screened_rates(k_Ca40_He4_to_Ti44_approx)*Y(Ca40)*Y(He4)*state.rho + screened_rates(k_Ti44_to_Ca40_He4_approx)*Y(Ti44)) + @@ -1400,69 +1404,69 @@ void rhs_nuc(const burn_t& state, ydot_nuc(Cr48) = (-screened_rates(k_He4_Cr48_to_Fe52)*Y(Cr48)*Y(He4)*state.rho + screened_rates(k_Fe52_to_He4_Cr48_derived)*Y(Fe52)) + - (-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_Cr48_to_p_nse_Mn51)*Y(Cr48)*Y(He4)*state.rho + screened_rates(k_p_nse_Mn51_to_He4_Cr48_derived)*Y(Mn51)*Y(P_nse)*state.rho) + (screened_rates(k_Ti44_He4_to_Cr48_approx)*Y(He4)*Y(Ti44)*state.rho + -screened_rates(k_Cr48_to_Ti44_He4_approx)*Y(Cr48)); ydot_nuc(Mn51) = - (-screened_rates(k_p_Mn51_to_Fe52)*Y(Mn51)*Y(H1)*state.rho + screened_rates(k_Fe52_to_p_Mn51_derived)*Y(Fe52)) + + (-screened_rates(k_p_nse_Mn51_to_Fe52)*Y(Mn51)*Y(P_nse)*state.rho + screened_rates(k_Fe52_to_p_nse_Mn51_derived)*Y(Fe52)) + (-screened_rates(k_He4_Mn51_to_Co55)*Y(He4)*Y(Mn51)*state.rho + screened_rates(k_Co55_to_He4_Mn51_derived)*Y(Co55)) + - (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_Mn51_to_p_Fe54)*Y(He4)*Y(Mn51)*state.rho + screened_rates(k_p_Fe54_to_He4_Mn51_derived)*Y(Fe54)*Y(H1)*state.rho); + (screened_rates(k_He4_Cr48_to_p_nse_Mn51)*Y(Cr48)*Y(He4)*state.rho + -screened_rates(k_p_nse_Mn51_to_He4_Cr48_derived)*Y(Mn51)*Y(P_nse)*state.rho) + + (-screened_rates(k_He4_Mn51_to_p_nse_Fe54)*Y(He4)*Y(Mn51)*state.rho + screened_rates(k_p_nse_Fe54_to_He4_Mn51_derived)*Y(Fe54)*Y(P_nse)*state.rho); ydot_nuc(Fe52) = (screened_rates(k_He4_Cr48_to_Fe52)*Y(Cr48)*Y(He4)*state.rho + -screened_rates(k_Fe52_to_He4_Cr48_derived)*Y(Fe52)) + - (screened_rates(k_p_Mn51_to_Fe52)*Y(Mn51)*Y(H1)*state.rho + -screened_rates(k_Fe52_to_p_Mn51_derived)*Y(Fe52)) + + (screened_rates(k_p_nse_Mn51_to_Fe52)*Y(Mn51)*Y(P_nse)*state.rho + -screened_rates(k_Fe52_to_p_nse_Mn51_derived)*Y(Fe52)) + (-screened_rates(k_He4_Fe52_to_Ni56)*Y(Fe52)*Y(He4)*state.rho + screened_rates(k_Ni56_to_He4_Fe52_derived)*Y(Ni56)) + - (-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_He4_Fe52_to_p_nse_Co55)*Y(Fe52)*Y(He4)*state.rho + screened_rates(k_p_nse_Co55_to_He4_Fe52_derived)*Y(Co55)*Y(P_nse)*state.rho) + (-screened_rates(k_Fe52_n_n_to_Fe54_approx)*Y(Fe52)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho) + screened_rates(k_Fe54_to_Fe52_n_n_approx)*Y(Fe54)); ydot_nuc(Fe54) = - (-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_nse_Fe54_to_Co55)*Y(Fe54)*Y(P_nse)*state.rho + screened_rates(k_Co55_to_p_nse_Fe54_derived)*Y(Co55)) + (-screened_rates(k_He4_Fe54_to_Ni58)*Y(Fe54)*Y(He4)*state.rho + screened_rates(k_Ni58_to_He4_Fe54_derived)*Y(Ni58)) + - (screened_rates(k_He4_Mn51_to_p_Fe54)*Y(He4)*Y(Mn51)*state.rho + -screened_rates(k_p_Fe54_to_He4_Mn51_derived)*Y(Fe54)*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_He4_Mn51_to_p_nse_Fe54)*Y(He4)*Y(Mn51)*state.rho + -screened_rates(k_p_nse_Fe54_to_He4_Mn51_derived)*Y(Fe54)*Y(P_nse)*state.rho) + + (screened_rates(k_p_nse_Co57_to_He4_Fe54)*Y(Co57)*Y(P_nse)*state.rho + -screened_rates(k_He4_Fe54_to_p_nse_Co57_derived)*Y(Fe54)*Y(He4)*state.rho) + (screened_rates(k_Fe52_n_n_to_Fe54_approx)*Y(Fe52)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho) + -screened_rates(k_Fe54_to_Fe52_n_n_approx)*Y(Fe54)) + (-screened_rates(k_Fe54_n_n_to_Fe56_approx)*Y(Fe54)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho) + screened_rates(k_Fe56_to_Fe54_n_n_approx)*Y(Fe56)); ydot_nuc(Fe56) = - (-screened_rates(k_p_Fe56_to_Co57)*Y(Fe56)*Y(H1)*state.rho + screened_rates(k_Co57_to_p_Fe56_derived)*Y(Co57)) + - (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_nse_Fe56_to_Co57)*Y(Fe56)*Y(P_nse)*state.rho + screened_rates(k_Co57_to_p_nse_Fe56_derived)*Y(Co57)) + + (screened_rates(k_n_Co56_to_p_nse_Fe56)*Y(Co56)*Y(N)*state.rho + -screened_rates(k_p_nse_Fe56_to_n_Co56_derived)*Y(Fe56)*Y(P_nse)*state.rho) + (screened_rates(k_Fe54_n_n_to_Fe56_approx)*Y(Fe54)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho) + -screened_rates(k_Fe56_to_Fe54_n_n_approx)*Y(Fe56)) + (screened_rates(k_Co56_to_Fe56)*Y(Co56) + -screened_rates(k_Fe56_to_Co56)*Y(Fe56)); ydot_nuc(Co55) = (screened_rates(k_He4_Mn51_to_Co55)*Y(He4)*Y(Mn51)*state.rho + -screened_rates(k_Co55_to_He4_Mn51_derived)*Y(Co55)) + - (-screened_rates(k_p_Co55_to_Ni56)*Y(Co55)*Y(H1)*state.rho + screened_rates(k_Ni56_to_p_Co55_derived)*Y(Ni56)) + - (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_Fe54_to_Co55)*Y(Fe54)*Y(H1)*state.rho + -screened_rates(k_Co55_to_p_Fe54_derived)*Y(Co55)) + + (-screened_rates(k_p_nse_Co55_to_Ni56)*Y(Co55)*Y(P_nse)*state.rho + screened_rates(k_Ni56_to_p_nse_Co55_derived)*Y(Ni56)) + + (screened_rates(k_He4_Fe52_to_p_nse_Co55)*Y(Fe52)*Y(He4)*state.rho + -screened_rates(k_p_nse_Co55_to_He4_Fe52_derived)*Y(Co55)*Y(P_nse)*state.rho) + + (screened_rates(k_p_nse_Fe54_to_Co55)*Y(Fe54)*Y(P_nse)*state.rho + -screened_rates(k_Co55_to_p_nse_Fe54_derived)*Y(Co55)) + (-screened_rates(k_n_Co55_to_Co56)*Y(Co55)*Y(N)*state.rho + screened_rates(k_Co56_to_n_Co55_derived)*Y(Co56)) + - (-screened_rates(k_He4_Co55_to_p_Ni58)*Y(Co55)*Y(He4)*state.rho + screened_rates(k_p_Ni58_to_He4_Co55_derived)*Y(Ni58)*Y(H1)*state.rho); + (-screened_rates(k_He4_Co55_to_p_nse_Ni58)*Y(Co55)*Y(He4)*state.rho + screened_rates(k_p_nse_Ni58_to_He4_Co55_derived)*Y(Ni58)*Y(P_nse)*state.rho); ydot_nuc(Co56) = (screened_rates(k_n_Co55_to_Co56)*Y(Co55)*Y(N)*state.rho + -screened_rates(k_Co56_to_n_Co55_derived)*Y(Co56)) + (-screened_rates(k_n_Co56_to_Co57)*Y(Co56)*Y(N)*state.rho + screened_rates(k_Co57_to_n_Co56_derived)*Y(Co57)) + - (-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_n_Ni56_to_p_Co56)*Y(Ni56)*Y(N)*state.rho + -screened_rates(k_p_Co56_to_n_Ni56_derived)*Y(Co56)*Y(H1)*state.rho) + + (-screened_rates(k_n_Co56_to_p_nse_Fe56)*Y(Co56)*Y(N)*state.rho + screened_rates(k_p_nse_Fe56_to_n_Co56_derived)*Y(Fe56)*Y(P_nse)*state.rho) + + (screened_rates(k_n_Ni56_to_p_nse_Co56)*Y(Ni56)*Y(N)*state.rho + -screened_rates(k_p_nse_Co56_to_n_Ni56_derived)*Y(Co56)*Y(P_nse)*state.rho) + (-screened_rates(k_Co56_to_Fe56)*Y(Co56) + screened_rates(k_Fe56_to_Co56)*Y(Fe56)) + (screened_rates(k_Ni56_to_Co56)*Y(Ni56) + -screened_rates(k_Co56_to_Ni56)*Y(Co56)); ydot_nuc(Co57) = - (screened_rates(k_p_Fe56_to_Co57)*Y(Fe56)*Y(H1)*state.rho + -screened_rates(k_Co57_to_p_Fe56_derived)*Y(Co57)) + + (screened_rates(k_p_nse_Fe56_to_Co57)*Y(Fe56)*Y(P_nse)*state.rho + -screened_rates(k_Co57_to_p_nse_Fe56_derived)*Y(Co57)) + (screened_rates(k_n_Co56_to_Co57)*Y(Co56)*Y(N)*state.rho + -screened_rates(k_Co57_to_n_Co56_derived)*Y(Co57)) + - (-screened_rates(k_p_Co57_to_Ni58)*Y(Co57)*Y(H1)*state.rho + screened_rates(k_Ni58_to_p_Co57_derived)*Y(Ni58)) + - (-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_p_nse_Co57_to_Ni58)*Y(Co57)*Y(P_nse)*state.rho + screened_rates(k_Ni58_to_p_nse_Co57_derived)*Y(Ni58)) + + (-screened_rates(k_p_nse_Co57_to_He4_Fe54)*Y(Co57)*Y(P_nse)*state.rho + screened_rates(k_He4_Fe54_to_p_nse_Co57_derived)*Y(Fe54)*Y(He4)*state.rho); ydot_nuc(Ni56) = (screened_rates(k_He4_Fe52_to_Ni56)*Y(Fe52)*Y(He4)*state.rho + -screened_rates(k_Ni56_to_He4_Fe52_derived)*Y(Ni56)) + - (screened_rates(k_p_Co55_to_Ni56)*Y(Co55)*Y(H1)*state.rho + -screened_rates(k_Ni56_to_p_Co55_derived)*Y(Ni56)) + - (-screened_rates(k_n_Ni56_to_p_Co56)*Y(Ni56)*Y(N)*state.rho + screened_rates(k_p_Co56_to_n_Ni56_derived)*Y(Co56)*Y(H1)*state.rho) + + (screened_rates(k_p_nse_Co55_to_Ni56)*Y(Co55)*Y(P_nse)*state.rho + -screened_rates(k_Ni56_to_p_nse_Co55_derived)*Y(Ni56)) + + (-screened_rates(k_n_Ni56_to_p_nse_Co56)*Y(Ni56)*Y(N)*state.rho + screened_rates(k_p_nse_Co56_to_n_Ni56_derived)*Y(Co56)*Y(P_nse)*state.rho) + (-screened_rates(k_Ni56_n_n_to_Ni58_approx)*Y(Ni56)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho) + screened_rates(k_Ni58_to_Ni56_n_n_approx)*Y(Ni58)) + (-screened_rates(k_Ni56_to_Co56)*Y(Ni56) + screened_rates(k_Co56_to_Ni56)*Y(Co56)); ydot_nuc(Ni58) = (screened_rates(k_He4_Fe54_to_Ni58)*Y(Fe54)*Y(He4)*state.rho + -screened_rates(k_Ni58_to_He4_Fe54_derived)*Y(Ni58)) + - (screened_rates(k_p_Co57_to_Ni58)*Y(Co57)*Y(H1)*state.rho + -screened_rates(k_Ni58_to_p_Co57_derived)*Y(Ni58)) + - (screened_rates(k_He4_Co55_to_p_Ni58)*Y(Co55)*Y(He4)*state.rho + -screened_rates(k_p_Ni58_to_He4_Co55_derived)*Y(Ni58)*Y(H1)*state.rho) + + (screened_rates(k_p_nse_Co57_to_Ni58)*Y(Co57)*Y(P_nse)*state.rho + -screened_rates(k_Ni58_to_p_nse_Co57_derived)*Y(Ni58)) + + (screened_rates(k_He4_Co55_to_p_nse_Ni58)*Y(Co55)*Y(He4)*state.rho + -screened_rates(k_p_nse_Ni58_to_He4_Co55_derived)*Y(Ni58)*Y(P_nse)*state.rho) + (screened_rates(k_Ni56_n_n_to_Ni58_approx)*Y(Ni56)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho) + -screened_rates(k_Ni58_to_Ni56_n_n_approx)*Y(Ni58)); } @@ -1523,43 +1527,46 @@ void jac_nuc(const burn_t& state, amrex::Real scratch; - scratch = -4.0*screened_rates(k_Fe52_n_n_to_Fe54_approx)*Y(Fe52)*Y(N)*amrex::Math::powi<2>(state.rho) - 4.0*screened_rates(k_Fe54_n_n_to_Fe56_approx)*Y(Fe54)*Y(N)*amrex::Math::powi<2>(state.rho) - 4.0*screened_rates(k_Ni56_n_n_to_Ni58_approx)*Y(Ni56)*Y(N)*amrex::Math::powi<2>(state.rho) - screened_rates(k_n_Co55_to_Co56)*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_Ni56_to_p_Co56)*Y(Ni56)*state.rho - screened_rates(k_n_to_p); + scratch = -4.0*screened_rates(k_Fe52_n_n_to_Fe54_approx)*Y(Fe52)*Y(N)*amrex::Math::powi<2>(state.rho) - 4.0*screened_rates(k_Fe54_n_n_to_Fe56_approx)*Y(Fe54)*Y(N)*amrex::Math::powi<2>(state.rho) - 4.0*screened_rates(k_Ni56_n_n_to_Ni58_approx)*Y(Ni56)*Y(N)*amrex::Math::powi<2>(state.rho) - screened_rates(k_n_Co55_to_Co56)*Y(Co55)*state.rho - screened_rates(k_n_Co56_to_Co57)*Y(Co56)*state.rho - screened_rates(k_n_Co56_to_p_nse_Fe56)*Y(Co56)*state.rho - screened_rates(k_n_Ni56_to_p_nse_Co56)*Y(Ni56)*state.rho - screened_rates(k_n_to_p); jac.set(N, N, scratch); - scratch = screened_rates(k_p_Co56_to_n_Ni56_derived)*Y(Co56)*state.rho + screened_rates(k_p_Fe56_to_n_Co56_derived)*Y(Fe56)*state.rho + screened_rates(k_p_to_n); + scratch = screened_rates(k_p_to_n); jac.set(N, H1, scratch); + scratch = screened_rates(k_p_nse_Co56_to_n_Ni56_derived)*Y(Co56)*state.rho + screened_rates(k_p_nse_Fe56_to_n_Co56_derived)*Y(Fe56)*state.rho; + jac.set(N, P_nse, scratch); + scratch = -2.0*screened_rates(k_Fe52_n_n_to_Fe54_approx)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho); jac.set(N, Fe52, scratch); scratch = -2.0*screened_rates(k_Fe54_n_n_to_Fe56_approx)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho) + 2.0*screened_rates(k_Fe54_to_Fe52_n_n_approx); jac.set(N, Fe54, scratch); - scratch = 2.0*screened_rates(k_Fe56_to_Fe54_n_n_approx) + screened_rates(k_p_Fe56_to_n_Co56_derived)*Y(H1)*state.rho; + scratch = 2.0*screened_rates(k_Fe56_to_Fe54_n_n_approx) + screened_rates(k_p_nse_Fe56_to_n_Co56_derived)*Y(P_nse)*state.rho; jac.set(N, Fe56, scratch); scratch = -screened_rates(k_n_Co55_to_Co56)*Y(N)*state.rho; jac.set(N, Co55, scratch); - scratch = screened_rates(k_Co56_to_n_Co55_derived) - screened_rates(k_n_Co56_to_Co57)*Y(N)*state.rho - screened_rates(k_n_Co56_to_p_Fe56)*Y(N)*state.rho + screened_rates(k_p_Co56_to_n_Ni56_derived)*Y(H1)*state.rho; + scratch = screened_rates(k_Co56_to_n_Co55_derived) - screened_rates(k_n_Co56_to_Co57)*Y(N)*state.rho - screened_rates(k_n_Co56_to_p_nse_Fe56)*Y(N)*state.rho + screened_rates(k_p_nse_Co56_to_n_Ni56_derived)*Y(P_nse)*state.rho; jac.set(N, Co56, scratch); scratch = screened_rates(k_Co57_to_n_Co56_derived); jac.set(N, Co57, scratch); - scratch = -2.0*screened_rates(k_Ni56_n_n_to_Ni58_approx)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho) - screened_rates(k_n_Ni56_to_p_Co56)*Y(N)*state.rho; + scratch = -2.0*screened_rates(k_Ni56_n_n_to_Ni58_approx)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho) - screened_rates(k_n_Ni56_to_p_nse_Co56)*Y(N)*state.rho; jac.set(N, Ni56, scratch); scratch = 2.0*screened_rates(k_Ni58_to_Ni56_n_n_approx); jac.set(N, Ni58, scratch); - scratch = screened_rates(k_n_Co56_to_p_Fe56)*Y(Co56)*state.rho + screened_rates(k_n_Ni56_to_p_Co56)*Y(Ni56)*state.rho + screened_rates(k_n_to_p); + scratch = 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_n_Ni56_derived)*Y(Co56)*state.rho - screened_rates(k_p_Co57_to_He4_Fe54)*Y(Co57)*state.rho - screened_rates(k_p_Co57_to_Ni58)*Y(Co57)*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_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_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_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); + 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_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_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_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_Si28_to_p_P31_derived)*Y(Si28)*state.rho; + scratch = screened_rates(k_He4_F18_to_p_Ne21)*Y(F18)*state.rho + screened_rates(k_He4_Mg24_to_p_Al27_derived)*Y(Mg24)*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_Si28_to_p_P31_derived)*Y(Si28)*state.rho; jac.set(H1, He4, scratch); scratch = 1.0*screened_rates(k_C12_C12_to_p_Na23)*Y(C12)*state.rho + screened_rates(k_C12_O16_to_p_Al27)*Y(O16)*state.rho - screened_rates(k_p_C12_to_N13)*Y(H1)*state.rho; @@ -1601,40 +1608,52 @@ void jac_nuc(const burn_t& state, scratch = screened_rates(k_S32_to_p_P31_derived); jac.set(H1, S32, scratch); - scratch = screened_rates(k_He4_Cr48_to_p_Mn51)*Y(He4)*state.rho; - jac.set(H1, Cr48, scratch); + scratch = screened_rates(k_n_Co56_to_p_nse_Fe56)*Y(Co56)*state.rho + screened_rates(k_n_Ni56_to_p_nse_Co56)*Y(Ni56)*state.rho; + jac.set(P_nse, N, scratch); + + scratch = -screened_rates(k_p_nse_Co55_to_He4_Fe52_derived)*Y(Co55)*state.rho - screened_rates(k_p_nse_Co55_to_Ni56)*Y(Co55)*state.rho - screened_rates(k_p_nse_Co56_to_n_Ni56_derived)*Y(Co56)*state.rho - screened_rates(k_p_nse_Co57_to_He4_Fe54)*Y(Co57)*state.rho - screened_rates(k_p_nse_Co57_to_Ni58)*Y(Co57)*state.rho - screened_rates(k_p_nse_Fe54_to_Co55)*Y(Fe54)*state.rho - screened_rates(k_p_nse_Fe54_to_He4_Mn51_derived)*Y(Fe54)*state.rho - screened_rates(k_p_nse_Fe56_to_Co57)*Y(Fe56)*state.rho - screened_rates(k_p_nse_Fe56_to_n_Co56_derived)*Y(Fe56)*state.rho - screened_rates(k_p_nse_Mn51_to_Fe52)*Y(Mn51)*state.rho - screened_rates(k_p_nse_Mn51_to_He4_Cr48_derived)*Y(Mn51)*state.rho - screened_rates(k_p_nse_Ni58_to_He4_Co55_derived)*Y(Ni58)*state.rho; + jac.set(P_nse, P_nse, scratch); - scratch = screened_rates(k_He4_Mn51_to_p_Fe54)*Y(He4)*state.rho - screened_rates(k_p_Mn51_to_Fe52)*Y(H1)*state.rho - screened_rates(k_p_Mn51_to_He4_Cr48_derived)*Y(H1)*state.rho; - jac.set(H1, Mn51, scratch); + scratch = screened_rates(k_He4_Co55_to_p_nse_Ni58)*Y(Co55)*state.rho + screened_rates(k_He4_Cr48_to_p_nse_Mn51)*Y(Cr48)*state.rho + screened_rates(k_He4_Fe52_to_p_nse_Co55)*Y(Fe52)*state.rho + screened_rates(k_He4_Fe54_to_p_nse_Co57_derived)*Y(Fe54)*state.rho + screened_rates(k_He4_Mn51_to_p_nse_Fe54)*Y(Mn51)*state.rho; + jac.set(P_nse, He4, scratch); - scratch = screened_rates(k_Fe52_to_p_Mn51_derived) + screened_rates(k_He4_Fe52_to_p_Co55)*Y(He4)*state.rho; - jac.set(H1, Fe52, scratch); + scratch = screened_rates(k_He4_Cr48_to_p_nse_Mn51)*Y(He4)*state.rho; + jac.set(P_nse, Cr48, scratch); - scratch = screened_rates(k_He4_Fe54_to_p_Co57_derived)*Y(He4)*state.rho - screened_rates(k_p_Fe54_to_Co55)*Y(H1)*state.rho - screened_rates(k_p_Fe54_to_He4_Mn51_derived)*Y(H1)*state.rho; - jac.set(H1, Fe54, scratch); + scratch = screened_rates(k_He4_Mn51_to_p_nse_Fe54)*Y(He4)*state.rho - screened_rates(k_p_nse_Mn51_to_Fe52)*Y(P_nse)*state.rho - screened_rates(k_p_nse_Mn51_to_He4_Cr48_derived)*Y(P_nse)*state.rho; + jac.set(P_nse, Mn51, scratch); - scratch = -screened_rates(k_p_Fe56_to_Co57)*Y(H1)*state.rho - screened_rates(k_p_Fe56_to_n_Co56_derived)*Y(H1)*state.rho; - jac.set(H1, Fe56, scratch); + scratch = screened_rates(k_Fe52_to_p_nse_Mn51_derived) + screened_rates(k_He4_Fe52_to_p_nse_Co55)*Y(He4)*state.rho; + jac.set(P_nse, Fe52, scratch); - scratch = screened_rates(k_Co55_to_p_Fe54_derived) + screened_rates(k_He4_Co55_to_p_Ni58)*Y(He4)*state.rho - screened_rates(k_p_Co55_to_He4_Fe52_derived)*Y(H1)*state.rho - screened_rates(k_p_Co55_to_Ni56)*Y(H1)*state.rho; - jac.set(H1, Co55, scratch); + scratch = screened_rates(k_He4_Fe54_to_p_nse_Co57_derived)*Y(He4)*state.rho - screened_rates(k_p_nse_Fe54_to_Co55)*Y(P_nse)*state.rho - screened_rates(k_p_nse_Fe54_to_He4_Mn51_derived)*Y(P_nse)*state.rho; + jac.set(P_nse, Fe54, scratch); - scratch = screened_rates(k_n_Co56_to_p_Fe56)*Y(N)*state.rho - screened_rates(k_p_Co56_to_n_Ni56_derived)*Y(H1)*state.rho; - jac.set(H1, Co56, scratch); + scratch = -screened_rates(k_p_nse_Fe56_to_Co57)*Y(P_nse)*state.rho - screened_rates(k_p_nse_Fe56_to_n_Co56_derived)*Y(P_nse)*state.rho; + jac.set(P_nse, Fe56, scratch); - scratch = screened_rates(k_Co57_to_p_Fe56_derived) - screened_rates(k_p_Co57_to_He4_Fe54)*Y(H1)*state.rho - screened_rates(k_p_Co57_to_Ni58)*Y(H1)*state.rho; - jac.set(H1, Co57, scratch); + scratch = screened_rates(k_Co55_to_p_nse_Fe54_derived) + screened_rates(k_He4_Co55_to_p_nse_Ni58)*Y(He4)*state.rho - screened_rates(k_p_nse_Co55_to_He4_Fe52_derived)*Y(P_nse)*state.rho - screened_rates(k_p_nse_Co55_to_Ni56)*Y(P_nse)*state.rho; + jac.set(P_nse, Co55, scratch); - scratch = screened_rates(k_Ni56_to_p_Co55_derived) + screened_rates(k_n_Ni56_to_p_Co56)*Y(N)*state.rho; - jac.set(H1, Ni56, scratch); + scratch = screened_rates(k_n_Co56_to_p_nse_Fe56)*Y(N)*state.rho - screened_rates(k_p_nse_Co56_to_n_Ni56_derived)*Y(P_nse)*state.rho; + jac.set(P_nse, Co56, scratch); - scratch = screened_rates(k_Ni58_to_p_Co57_derived) - screened_rates(k_p_Ni58_to_He4_Co55_derived)*Y(H1)*state.rho; - jac.set(H1, Ni58, scratch); + scratch = screened_rates(k_Co57_to_p_nse_Fe56_derived) - screened_rates(k_p_nse_Co57_to_He4_Fe54)*Y(P_nse)*state.rho - screened_rates(k_p_nse_Co57_to_Ni58)*Y(P_nse)*state.rho; + jac.set(P_nse, Co57, scratch); - scratch = screened_rates(k_p_Al27_to_He4_Mg24)*Y(Al27)*state.rho + screened_rates(k_p_Co55_to_He4_Fe52_derived)*Y(Co55)*state.rho + screened_rates(k_p_Co57_to_He4_Fe54)*Y(Co57)*state.rho + screened_rates(k_p_Fe54_to_He4_Mn51_derived)*Y(Fe54)*state.rho + screened_rates(k_p_Mn51_to_He4_Cr48_derived)*Y(Mn51)*state.rho + screened_rates(k_p_Na23_to_He4_Ne20)*Y(Na23)*state.rho + screened_rates(k_p_Ne21_to_He4_F18_derived)*Y(Ne21)*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; + scratch = screened_rates(k_Ni56_to_p_nse_Co55_derived) + screened_rates(k_n_Ni56_to_p_nse_Co56)*Y(N)*state.rho; + jac.set(P_nse, Ni56, scratch); + + scratch = screened_rates(k_Ni58_to_p_nse_Co57_derived) - screened_rates(k_p_nse_Ni58_to_He4_Co55_derived)*Y(P_nse)*state.rho; + jac.set(P_nse, Ni58, scratch); + + scratch = screened_rates(k_p_Al27_to_He4_Mg24)*Y(Al27)*state.rho + screened_rates(k_p_Na23_to_He4_Ne20)*Y(Na23)*state.rho + screened_rates(k_p_Ne21_to_He4_F18_derived)*Y(Ne21)*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; jac.set(He4, H1, scratch); - scratch = -screened_rates(k_Ar36_He4_to_Ca40_approx)*Y(Ar36)*state.rho - screened_rates(k_Ca40_He4_to_Ti44_approx)*Y(Ca40)*state.rho - screened_rates(k_He4_C12_to_O16)*Y(C12)*state.rho - screened_rates(k_He4_Co55_to_p_Ni58)*Y(Co55)*state.rho - screened_rates(k_He4_Cr48_to_Fe52)*Y(Cr48)*state.rho - screened_rates(k_He4_Cr48_to_p_Mn51)*Y(Cr48)*state.rho - screened_rates(k_He4_F18_to_Na22)*Y(F18)*state.rho - screened_rates(k_He4_F18_to_p_Ne21)*Y(F18)*state.rho - screened_rates(k_He4_Fe52_to_Ni56)*Y(Fe52)*state.rho - screened_rates(k_He4_Fe52_to_p_Co55)*Y(Fe52)*state.rho - screened_rates(k_He4_Fe54_to_Ni58)*Y(Fe54)*state.rho - screened_rates(k_He4_Fe54_to_p_Co57_derived)*Y(Fe54)*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_Mg24_to_Si28)*Y(Mg24)*state.rho - screened_rates(k_He4_Mg24_to_p_Al27_derived)*Y(Mg24)*state.rho - screened_rates(k_He4_Mn51_to_Co55)*Y(Mn51)*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_N14_to_F18)*Y(N14)*state.rho - screened_rates(k_He4_Ne20_to_Mg24)*Y(Ne20)*state.rho - screened_rates(k_He4_Ne20_to_p_Na23_derived)*Y(Ne20)*state.rho - screened_rates(k_He4_O16_to_Ne20)*Y(O16)*state.rho - screened_rates(k_He4_Si28_to_S32)*Y(Si28)*state.rho - screened_rates(k_He4_Si28_to_p_P31_derived)*Y(Si28)*state.rho - screened_rates(k_S32_He4_to_Ar36_approx)*Y(S32)*state.rho - screened_rates(k_Ti44_He4_to_Cr48_approx)*Y(Ti44)*state.rho; + scratch = screened_rates(k_p_nse_Co55_to_He4_Fe52_derived)*Y(Co55)*state.rho + screened_rates(k_p_nse_Co57_to_He4_Fe54)*Y(Co57)*state.rho + screened_rates(k_p_nse_Fe54_to_He4_Mn51_derived)*Y(Fe54)*state.rho + screened_rates(k_p_nse_Mn51_to_He4_Cr48_derived)*Y(Mn51)*state.rho + screened_rates(k_p_nse_Ni58_to_He4_Co55_derived)*Y(Ni58)*state.rho; + jac.set(He4, P_nse, scratch); + + scratch = -screened_rates(k_Ar36_He4_to_Ca40_approx)*Y(Ar36)*state.rho - screened_rates(k_Ca40_He4_to_Ti44_approx)*Y(Ca40)*state.rho - screened_rates(k_He4_C12_to_O16)*Y(C12)*state.rho - screened_rates(k_He4_Co55_to_p_nse_Ni58)*Y(Co55)*state.rho - screened_rates(k_He4_Cr48_to_Fe52)*Y(Cr48)*state.rho - screened_rates(k_He4_Cr48_to_p_nse_Mn51)*Y(Cr48)*state.rho - screened_rates(k_He4_F18_to_Na22)*Y(F18)*state.rho - screened_rates(k_He4_F18_to_p_Ne21)*Y(F18)*state.rho - screened_rates(k_He4_Fe52_to_Ni56)*Y(Fe52)*state.rho - screened_rates(k_He4_Fe52_to_p_nse_Co55)*Y(Fe52)*state.rho - screened_rates(k_He4_Fe54_to_Ni58)*Y(Fe54)*state.rho - screened_rates(k_He4_Fe54_to_p_nse_Co57_derived)*Y(Fe54)*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_Mg24_to_Si28)*Y(Mg24)*state.rho - screened_rates(k_He4_Mg24_to_p_Al27_derived)*Y(Mg24)*state.rho - screened_rates(k_He4_Mn51_to_Co55)*Y(Mn51)*state.rho - screened_rates(k_He4_Mn51_to_p_nse_Fe54)*Y(Mn51)*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_Ne20_to_Mg24)*Y(Ne20)*state.rho - screened_rates(k_He4_Ne20_to_p_Na23_derived)*Y(Ne20)*state.rho - screened_rates(k_He4_O16_to_Ne20)*Y(O16)*state.rho - screened_rates(k_He4_Si28_to_S32)*Y(Si28)*state.rho - screened_rates(k_He4_Si28_to_p_P31_derived)*Y(Si28)*state.rho - screened_rates(k_S32_He4_to_Ar36_approx)*Y(S32)*state.rho - screened_rates(k_Ti44_He4_to_Cr48_approx)*Y(Ti44)*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 + 3.0*screened_rates(k_C12_to_He4_He4_He4_derived) - screened_rates(k_He4_C12_to_O16)*Y(He4)*state.rho; @@ -1688,28 +1707,28 @@ void jac_nuc(const burn_t& state, scratch = -screened_rates(k_Ti44_He4_to_Cr48_approx)*Y(He4)*state.rho + screened_rates(k_Ti44_to_Ca40_He4_approx); jac.set(He4, Ti44, scratch); - scratch = screened_rates(k_Cr48_to_Ti44_He4_approx) - screened_rates(k_He4_Cr48_to_Fe52)*Y(He4)*state.rho - screened_rates(k_He4_Cr48_to_p_Mn51)*Y(He4)*state.rho; + scratch = screened_rates(k_Cr48_to_Ti44_He4_approx) - screened_rates(k_He4_Cr48_to_Fe52)*Y(He4)*state.rho - screened_rates(k_He4_Cr48_to_p_nse_Mn51)*Y(He4)*state.rho; jac.set(He4, Cr48, scratch); - scratch = -screened_rates(k_He4_Mn51_to_Co55)*Y(He4)*state.rho - screened_rates(k_He4_Mn51_to_p_Fe54)*Y(He4)*state.rho + screened_rates(k_p_Mn51_to_He4_Cr48_derived)*Y(H1)*state.rho; + scratch = -screened_rates(k_He4_Mn51_to_Co55)*Y(He4)*state.rho - screened_rates(k_He4_Mn51_to_p_nse_Fe54)*Y(He4)*state.rho + screened_rates(k_p_nse_Mn51_to_He4_Cr48_derived)*Y(P_nse)*state.rho; jac.set(He4, Mn51, scratch); - scratch = screened_rates(k_Fe52_to_He4_Cr48_derived) - screened_rates(k_He4_Fe52_to_Ni56)*Y(He4)*state.rho - screened_rates(k_He4_Fe52_to_p_Co55)*Y(He4)*state.rho; + scratch = screened_rates(k_Fe52_to_He4_Cr48_derived) - screened_rates(k_He4_Fe52_to_Ni56)*Y(He4)*state.rho - screened_rates(k_He4_Fe52_to_p_nse_Co55)*Y(He4)*state.rho; jac.set(He4, Fe52, scratch); - scratch = -screened_rates(k_He4_Fe54_to_Ni58)*Y(He4)*state.rho - screened_rates(k_He4_Fe54_to_p_Co57_derived)*Y(He4)*state.rho + screened_rates(k_p_Fe54_to_He4_Mn51_derived)*Y(H1)*state.rho; + scratch = -screened_rates(k_He4_Fe54_to_Ni58)*Y(He4)*state.rho - screened_rates(k_He4_Fe54_to_p_nse_Co57_derived)*Y(He4)*state.rho + screened_rates(k_p_nse_Fe54_to_He4_Mn51_derived)*Y(P_nse)*state.rho; jac.set(He4, Fe54, scratch); - scratch = screened_rates(k_Co55_to_He4_Mn51_derived) - screened_rates(k_He4_Co55_to_p_Ni58)*Y(He4)*state.rho + screened_rates(k_p_Co55_to_He4_Fe52_derived)*Y(H1)*state.rho; + scratch = screened_rates(k_Co55_to_He4_Mn51_derived) - screened_rates(k_He4_Co55_to_p_nse_Ni58)*Y(He4)*state.rho + screened_rates(k_p_nse_Co55_to_He4_Fe52_derived)*Y(P_nse)*state.rho; jac.set(He4, Co55, scratch); - scratch = screened_rates(k_p_Co57_to_He4_Fe54)*Y(H1)*state.rho; + scratch = screened_rates(k_p_nse_Co57_to_He4_Fe54)*Y(P_nse)*state.rho; jac.set(He4, Co57, scratch); scratch = screened_rates(k_Ni56_to_He4_Fe52_derived); jac.set(He4, Ni56, scratch); - scratch = screened_rates(k_Ni58_to_He4_Fe54_derived) + screened_rates(k_p_Ni58_to_He4_Co55_derived)*Y(H1)*state.rho; + scratch = screened_rates(k_Ni58_to_He4_Fe54_derived) + screened_rates(k_p_nse_Ni58_to_He4_Co55_derived)*Y(P_nse)*state.rho; jac.set(He4, Ni58, scratch); scratch = -screened_rates(k_p_C12_to_N13)*Y(C12)*state.rho; @@ -2006,40 +2025,40 @@ void jac_nuc(const burn_t& state, scratch = screened_rates(k_Cr48_to_Ti44_He4_approx); jac.set(Ti44, Cr48, scratch); - scratch = screened_rates(k_p_Mn51_to_He4_Cr48_derived)*Y(Mn51)*state.rho; - jac.set(Cr48, H1, scratch); + scratch = screened_rates(k_p_nse_Mn51_to_He4_Cr48_derived)*Y(Mn51)*state.rho; + jac.set(Cr48, P_nse, scratch); - scratch = -screened_rates(k_He4_Cr48_to_Fe52)*Y(Cr48)*state.rho - screened_rates(k_He4_Cr48_to_p_Mn51)*Y(Cr48)*state.rho + screened_rates(k_Ti44_He4_to_Cr48_approx)*Y(Ti44)*state.rho; + scratch = -screened_rates(k_He4_Cr48_to_Fe52)*Y(Cr48)*state.rho - screened_rates(k_He4_Cr48_to_p_nse_Mn51)*Y(Cr48)*state.rho + screened_rates(k_Ti44_He4_to_Cr48_approx)*Y(Ti44)*state.rho; jac.set(Cr48, He4, scratch); scratch = screened_rates(k_Ti44_He4_to_Cr48_approx)*Y(He4)*state.rho; jac.set(Cr48, Ti44, scratch); - scratch = -screened_rates(k_Cr48_to_Ti44_He4_approx) - screened_rates(k_He4_Cr48_to_Fe52)*Y(He4)*state.rho - screened_rates(k_He4_Cr48_to_p_Mn51)*Y(He4)*state.rho; + scratch = -screened_rates(k_Cr48_to_Ti44_He4_approx) - screened_rates(k_He4_Cr48_to_Fe52)*Y(He4)*state.rho - screened_rates(k_He4_Cr48_to_p_nse_Mn51)*Y(He4)*state.rho; jac.set(Cr48, Cr48, scratch); - scratch = screened_rates(k_p_Mn51_to_He4_Cr48_derived)*Y(H1)*state.rho; + scratch = screened_rates(k_p_nse_Mn51_to_He4_Cr48_derived)*Y(P_nse)*state.rho; jac.set(Cr48, Mn51, scratch); scratch = screened_rates(k_Fe52_to_He4_Cr48_derived); jac.set(Cr48, Fe52, scratch); - scratch = screened_rates(k_p_Fe54_to_He4_Mn51_derived)*Y(Fe54)*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; - jac.set(Mn51, H1, scratch); + scratch = screened_rates(k_p_nse_Fe54_to_He4_Mn51_derived)*Y(Fe54)*state.rho - screened_rates(k_p_nse_Mn51_to_Fe52)*Y(Mn51)*state.rho - screened_rates(k_p_nse_Mn51_to_He4_Cr48_derived)*Y(Mn51)*state.rho; + jac.set(Mn51, P_nse, scratch); - scratch = screened_rates(k_He4_Cr48_to_p_Mn51)*Y(Cr48)*state.rho - screened_rates(k_He4_Mn51_to_Co55)*Y(Mn51)*state.rho - screened_rates(k_He4_Mn51_to_p_Fe54)*Y(Mn51)*state.rho; + scratch = screened_rates(k_He4_Cr48_to_p_nse_Mn51)*Y(Cr48)*state.rho - screened_rates(k_He4_Mn51_to_Co55)*Y(Mn51)*state.rho - screened_rates(k_He4_Mn51_to_p_nse_Fe54)*Y(Mn51)*state.rho; jac.set(Mn51, He4, scratch); - scratch = screened_rates(k_He4_Cr48_to_p_Mn51)*Y(He4)*state.rho; + scratch = screened_rates(k_He4_Cr48_to_p_nse_Mn51)*Y(He4)*state.rho; jac.set(Mn51, Cr48, scratch); - scratch = -screened_rates(k_He4_Mn51_to_Co55)*Y(He4)*state.rho - screened_rates(k_He4_Mn51_to_p_Fe54)*Y(He4)*state.rho - screened_rates(k_p_Mn51_to_Fe52)*Y(H1)*state.rho - screened_rates(k_p_Mn51_to_He4_Cr48_derived)*Y(H1)*state.rho; + scratch = -screened_rates(k_He4_Mn51_to_Co55)*Y(He4)*state.rho - screened_rates(k_He4_Mn51_to_p_nse_Fe54)*Y(He4)*state.rho - screened_rates(k_p_nse_Mn51_to_Fe52)*Y(P_nse)*state.rho - screened_rates(k_p_nse_Mn51_to_He4_Cr48_derived)*Y(P_nse)*state.rho; jac.set(Mn51, Mn51, scratch); - scratch = screened_rates(k_Fe52_to_p_Mn51_derived); + scratch = screened_rates(k_Fe52_to_p_nse_Mn51_derived); jac.set(Mn51, Fe52, scratch); - scratch = screened_rates(k_p_Fe54_to_He4_Mn51_derived)*Y(H1)*state.rho; + scratch = screened_rates(k_p_nse_Fe54_to_He4_Mn51_derived)*Y(P_nse)*state.rho; jac.set(Mn51, Fe54, scratch); scratch = screened_rates(k_Co55_to_He4_Mn51_derived); @@ -2048,25 +2067,25 @@ void jac_nuc(const burn_t& state, scratch = -2.0*screened_rates(k_Fe52_n_n_to_Fe54_approx)*Y(Fe52)*Y(N)*amrex::Math::powi<2>(state.rho); jac.set(Fe52, N, scratch); - scratch = screened_rates(k_p_Co55_to_He4_Fe52_derived)*Y(Co55)*state.rho + screened_rates(k_p_Mn51_to_Fe52)*Y(Mn51)*state.rho; - jac.set(Fe52, H1, scratch); + scratch = screened_rates(k_p_nse_Co55_to_He4_Fe52_derived)*Y(Co55)*state.rho + screened_rates(k_p_nse_Mn51_to_Fe52)*Y(Mn51)*state.rho; + jac.set(Fe52, P_nse, scratch); - scratch = screened_rates(k_He4_Cr48_to_Fe52)*Y(Cr48)*state.rho - screened_rates(k_He4_Fe52_to_Ni56)*Y(Fe52)*state.rho - screened_rates(k_He4_Fe52_to_p_Co55)*Y(Fe52)*state.rho; + scratch = screened_rates(k_He4_Cr48_to_Fe52)*Y(Cr48)*state.rho - screened_rates(k_He4_Fe52_to_Ni56)*Y(Fe52)*state.rho - screened_rates(k_He4_Fe52_to_p_nse_Co55)*Y(Fe52)*state.rho; jac.set(Fe52, He4, scratch); scratch = screened_rates(k_He4_Cr48_to_Fe52)*Y(He4)*state.rho; jac.set(Fe52, Cr48, scratch); - scratch = screened_rates(k_p_Mn51_to_Fe52)*Y(H1)*state.rho; + scratch = screened_rates(k_p_nse_Mn51_to_Fe52)*Y(P_nse)*state.rho; jac.set(Fe52, Mn51, scratch); - scratch = -screened_rates(k_Fe52_n_n_to_Fe54_approx)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho) - screened_rates(k_Fe52_to_He4_Cr48_derived) - screened_rates(k_Fe52_to_p_Mn51_derived) - screened_rates(k_He4_Fe52_to_Ni56)*Y(He4)*state.rho - screened_rates(k_He4_Fe52_to_p_Co55)*Y(He4)*state.rho; + scratch = -screened_rates(k_Fe52_n_n_to_Fe54_approx)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho) - screened_rates(k_Fe52_to_He4_Cr48_derived) - screened_rates(k_Fe52_to_p_nse_Mn51_derived) - screened_rates(k_He4_Fe52_to_Ni56)*Y(He4)*state.rho - screened_rates(k_He4_Fe52_to_p_nse_Co55)*Y(He4)*state.rho; jac.set(Fe52, Fe52, scratch); scratch = screened_rates(k_Fe54_to_Fe52_n_n_approx); jac.set(Fe52, Fe54, scratch); - scratch = screened_rates(k_p_Co55_to_He4_Fe52_derived)*Y(H1)*state.rho; + scratch = screened_rates(k_p_nse_Co55_to_He4_Fe52_derived)*Y(P_nse)*state.rho; jac.set(Fe52, Co55, scratch); scratch = screened_rates(k_Ni56_to_He4_Fe52_derived); @@ -2075,131 +2094,131 @@ void jac_nuc(const burn_t& state, scratch = 2.0*screened_rates(k_Fe52_n_n_to_Fe54_approx)*Y(Fe52)*Y(N)*amrex::Math::powi<2>(state.rho) - 2.0*screened_rates(k_Fe54_n_n_to_Fe56_approx)*Y(Fe54)*Y(N)*amrex::Math::powi<2>(state.rho); jac.set(Fe54, N, scratch); - scratch = screened_rates(k_p_Co57_to_He4_Fe54)*Y(Co57)*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; - jac.set(Fe54, H1, scratch); + scratch = screened_rates(k_p_nse_Co57_to_He4_Fe54)*Y(Co57)*state.rho - screened_rates(k_p_nse_Fe54_to_Co55)*Y(Fe54)*state.rho - screened_rates(k_p_nse_Fe54_to_He4_Mn51_derived)*Y(Fe54)*state.rho; + jac.set(Fe54, P_nse, scratch); - scratch = -screened_rates(k_He4_Fe54_to_Ni58)*Y(Fe54)*state.rho - screened_rates(k_He4_Fe54_to_p_Co57_derived)*Y(Fe54)*state.rho + screened_rates(k_He4_Mn51_to_p_Fe54)*Y(Mn51)*state.rho; + scratch = -screened_rates(k_He4_Fe54_to_Ni58)*Y(Fe54)*state.rho - screened_rates(k_He4_Fe54_to_p_nse_Co57_derived)*Y(Fe54)*state.rho + screened_rates(k_He4_Mn51_to_p_nse_Fe54)*Y(Mn51)*state.rho; jac.set(Fe54, He4, scratch); - scratch = screened_rates(k_He4_Mn51_to_p_Fe54)*Y(He4)*state.rho; + scratch = screened_rates(k_He4_Mn51_to_p_nse_Fe54)*Y(He4)*state.rho; jac.set(Fe54, Mn51, scratch); scratch = screened_rates(k_Fe52_n_n_to_Fe54_approx)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho); jac.set(Fe54, Fe52, scratch); - scratch = -screened_rates(k_Fe54_n_n_to_Fe56_approx)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho) - screened_rates(k_Fe54_to_Fe52_n_n_approx) - screened_rates(k_He4_Fe54_to_Ni58)*Y(He4)*state.rho - screened_rates(k_He4_Fe54_to_p_Co57_derived)*Y(He4)*state.rho - screened_rates(k_p_Fe54_to_Co55)*Y(H1)*state.rho - screened_rates(k_p_Fe54_to_He4_Mn51_derived)*Y(H1)*state.rho; + scratch = -screened_rates(k_Fe54_n_n_to_Fe56_approx)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho) - screened_rates(k_Fe54_to_Fe52_n_n_approx) - screened_rates(k_He4_Fe54_to_Ni58)*Y(He4)*state.rho - screened_rates(k_He4_Fe54_to_p_nse_Co57_derived)*Y(He4)*state.rho - screened_rates(k_p_nse_Fe54_to_Co55)*Y(P_nse)*state.rho - screened_rates(k_p_nse_Fe54_to_He4_Mn51_derived)*Y(P_nse)*state.rho; jac.set(Fe54, Fe54, scratch); scratch = screened_rates(k_Fe56_to_Fe54_n_n_approx); jac.set(Fe54, Fe56, scratch); - scratch = screened_rates(k_Co55_to_p_Fe54_derived); + scratch = screened_rates(k_Co55_to_p_nse_Fe54_derived); jac.set(Fe54, Co55, scratch); - scratch = screened_rates(k_p_Co57_to_He4_Fe54)*Y(H1)*state.rho; + scratch = screened_rates(k_p_nse_Co57_to_He4_Fe54)*Y(P_nse)*state.rho; jac.set(Fe54, Co57, scratch); scratch = screened_rates(k_Ni58_to_He4_Fe54_derived); jac.set(Fe54, Ni58, scratch); - scratch = 2.0*screened_rates(k_Fe54_n_n_to_Fe56_approx)*Y(Fe54)*Y(N)*amrex::Math::powi<2>(state.rho) + screened_rates(k_n_Co56_to_p_Fe56)*Y(Co56)*state.rho; + scratch = 2.0*screened_rates(k_Fe54_n_n_to_Fe56_approx)*Y(Fe54)*Y(N)*amrex::Math::powi<2>(state.rho) + screened_rates(k_n_Co56_to_p_nse_Fe56)*Y(Co56)*state.rho; jac.set(Fe56, N, scratch); - scratch = -screened_rates(k_p_Fe56_to_Co57)*Y(Fe56)*state.rho - screened_rates(k_p_Fe56_to_n_Co56_derived)*Y(Fe56)*state.rho; - jac.set(Fe56, H1, scratch); + scratch = -screened_rates(k_p_nse_Fe56_to_Co57)*Y(Fe56)*state.rho - screened_rates(k_p_nse_Fe56_to_n_Co56_derived)*Y(Fe56)*state.rho; + jac.set(Fe56, P_nse, scratch); scratch = screened_rates(k_Fe54_n_n_to_Fe56_approx)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho); jac.set(Fe56, Fe54, scratch); - scratch = -screened_rates(k_Fe56_to_Co56) - screened_rates(k_Fe56_to_Fe54_n_n_approx) - screened_rates(k_p_Fe56_to_Co57)*Y(H1)*state.rho - screened_rates(k_p_Fe56_to_n_Co56_derived)*Y(H1)*state.rho; + scratch = -screened_rates(k_Fe56_to_Co56) - screened_rates(k_Fe56_to_Fe54_n_n_approx) - screened_rates(k_p_nse_Fe56_to_Co57)*Y(P_nse)*state.rho - screened_rates(k_p_nse_Fe56_to_n_Co56_derived)*Y(P_nse)*state.rho; jac.set(Fe56, Fe56, scratch); - scratch = screened_rates(k_Co56_to_Fe56) + screened_rates(k_n_Co56_to_p_Fe56)*Y(N)*state.rho; + scratch = screened_rates(k_Co56_to_Fe56) + screened_rates(k_n_Co56_to_p_nse_Fe56)*Y(N)*state.rho; jac.set(Fe56, Co56, scratch); - scratch = screened_rates(k_Co57_to_p_Fe56_derived); + scratch = screened_rates(k_Co57_to_p_nse_Fe56_derived); jac.set(Fe56, Co57, scratch); scratch = -screened_rates(k_n_Co55_to_Co56)*Y(Co55)*state.rho; jac.set(Co55, N, scratch); - scratch = -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_Fe54_to_Co55)*Y(Fe54)*state.rho + screened_rates(k_p_Ni58_to_He4_Co55_derived)*Y(Ni58)*state.rho; - jac.set(Co55, H1, scratch); + scratch = -screened_rates(k_p_nse_Co55_to_He4_Fe52_derived)*Y(Co55)*state.rho - screened_rates(k_p_nse_Co55_to_Ni56)*Y(Co55)*state.rho + screened_rates(k_p_nse_Fe54_to_Co55)*Y(Fe54)*state.rho + screened_rates(k_p_nse_Ni58_to_He4_Co55_derived)*Y(Ni58)*state.rho; + jac.set(Co55, P_nse, scratch); - scratch = -screened_rates(k_He4_Co55_to_p_Ni58)*Y(Co55)*state.rho + screened_rates(k_He4_Fe52_to_p_Co55)*Y(Fe52)*state.rho + screened_rates(k_He4_Mn51_to_Co55)*Y(Mn51)*state.rho; + scratch = -screened_rates(k_He4_Co55_to_p_nse_Ni58)*Y(Co55)*state.rho + screened_rates(k_He4_Fe52_to_p_nse_Co55)*Y(Fe52)*state.rho + screened_rates(k_He4_Mn51_to_Co55)*Y(Mn51)*state.rho; jac.set(Co55, He4, scratch); scratch = screened_rates(k_He4_Mn51_to_Co55)*Y(He4)*state.rho; jac.set(Co55, Mn51, scratch); - scratch = screened_rates(k_He4_Fe52_to_p_Co55)*Y(He4)*state.rho; + scratch = screened_rates(k_He4_Fe52_to_p_nse_Co55)*Y(He4)*state.rho; jac.set(Co55, Fe52, scratch); - scratch = screened_rates(k_p_Fe54_to_Co55)*Y(H1)*state.rho; + scratch = screened_rates(k_p_nse_Fe54_to_Co55)*Y(P_nse)*state.rho; jac.set(Co55, Fe54, scratch); - scratch = -screened_rates(k_Co55_to_He4_Mn51_derived) - screened_rates(k_Co55_to_p_Fe54_derived) - screened_rates(k_He4_Co55_to_p_Ni58)*Y(He4)*state.rho - screened_rates(k_n_Co55_to_Co56)*Y(N)*state.rho - screened_rates(k_p_Co55_to_He4_Fe52_derived)*Y(H1)*state.rho - screened_rates(k_p_Co55_to_Ni56)*Y(H1)*state.rho; + scratch = -screened_rates(k_Co55_to_He4_Mn51_derived) - screened_rates(k_Co55_to_p_nse_Fe54_derived) - screened_rates(k_He4_Co55_to_p_nse_Ni58)*Y(He4)*state.rho - screened_rates(k_n_Co55_to_Co56)*Y(N)*state.rho - screened_rates(k_p_nse_Co55_to_He4_Fe52_derived)*Y(P_nse)*state.rho - screened_rates(k_p_nse_Co55_to_Ni56)*Y(P_nse)*state.rho; jac.set(Co55, Co55, scratch); scratch = screened_rates(k_Co56_to_n_Co55_derived); jac.set(Co55, Co56, scratch); - scratch = screened_rates(k_Ni56_to_p_Co55_derived); + scratch = screened_rates(k_Ni56_to_p_nse_Co55_derived); jac.set(Co55, Ni56, scratch); - scratch = screened_rates(k_p_Ni58_to_He4_Co55_derived)*Y(H1)*state.rho; + scratch = screened_rates(k_p_nse_Ni58_to_He4_Co55_derived)*Y(P_nse)*state.rho; jac.set(Co55, Ni58, scratch); - scratch = screened_rates(k_n_Co55_to_Co56)*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_Ni56_to_p_Co56)*Y(Ni56)*state.rho; + scratch = screened_rates(k_n_Co55_to_Co56)*Y(Co55)*state.rho - screened_rates(k_n_Co56_to_Co57)*Y(Co56)*state.rho - screened_rates(k_n_Co56_to_p_nse_Fe56)*Y(Co56)*state.rho + screened_rates(k_n_Ni56_to_p_nse_Co56)*Y(Ni56)*state.rho; jac.set(Co56, N, scratch); - scratch = -screened_rates(k_p_Co56_to_n_Ni56_derived)*Y(Co56)*state.rho + screened_rates(k_p_Fe56_to_n_Co56_derived)*Y(Fe56)*state.rho; - jac.set(Co56, H1, scratch); + scratch = -screened_rates(k_p_nse_Co56_to_n_Ni56_derived)*Y(Co56)*state.rho + screened_rates(k_p_nse_Fe56_to_n_Co56_derived)*Y(Fe56)*state.rho; + jac.set(Co56, P_nse, scratch); - scratch = screened_rates(k_Fe56_to_Co56) + screened_rates(k_p_Fe56_to_n_Co56_derived)*Y(H1)*state.rho; + scratch = screened_rates(k_Fe56_to_Co56) + screened_rates(k_p_nse_Fe56_to_n_Co56_derived)*Y(P_nse)*state.rho; jac.set(Co56, Fe56, scratch); scratch = screened_rates(k_n_Co55_to_Co56)*Y(N)*state.rho; jac.set(Co56, Co55, scratch); - scratch = -screened_rates(k_Co56_to_Fe56) - screened_rates(k_Co56_to_Ni56) - screened_rates(k_Co56_to_n_Co55_derived) - screened_rates(k_n_Co56_to_Co57)*Y(N)*state.rho - screened_rates(k_n_Co56_to_p_Fe56)*Y(N)*state.rho - screened_rates(k_p_Co56_to_n_Ni56_derived)*Y(H1)*state.rho; + scratch = -screened_rates(k_Co56_to_Fe56) - screened_rates(k_Co56_to_Ni56) - screened_rates(k_Co56_to_n_Co55_derived) - screened_rates(k_n_Co56_to_Co57)*Y(N)*state.rho - screened_rates(k_n_Co56_to_p_nse_Fe56)*Y(N)*state.rho - screened_rates(k_p_nse_Co56_to_n_Ni56_derived)*Y(P_nse)*state.rho; jac.set(Co56, Co56, scratch); scratch = screened_rates(k_Co57_to_n_Co56_derived); jac.set(Co56, Co57, scratch); - scratch = screened_rates(k_Ni56_to_Co56) + screened_rates(k_n_Ni56_to_p_Co56)*Y(N)*state.rho; + scratch = screened_rates(k_Ni56_to_Co56) + screened_rates(k_n_Ni56_to_p_nse_Co56)*Y(N)*state.rho; jac.set(Co56, Ni56, scratch); scratch = screened_rates(k_n_Co56_to_Co57)*Y(Co56)*state.rho; jac.set(Co57, N, scratch); - scratch = -screened_rates(k_p_Co57_to_He4_Fe54)*Y(Co57)*state.rho - screened_rates(k_p_Co57_to_Ni58)*Y(Co57)*state.rho + screened_rates(k_p_Fe56_to_Co57)*Y(Fe56)*state.rho; - jac.set(Co57, H1, scratch); + scratch = -screened_rates(k_p_nse_Co57_to_He4_Fe54)*Y(Co57)*state.rho - screened_rates(k_p_nse_Co57_to_Ni58)*Y(Co57)*state.rho + screened_rates(k_p_nse_Fe56_to_Co57)*Y(Fe56)*state.rho; + jac.set(Co57, P_nse, scratch); - scratch = screened_rates(k_He4_Fe54_to_p_Co57_derived)*Y(Fe54)*state.rho; + scratch = screened_rates(k_He4_Fe54_to_p_nse_Co57_derived)*Y(Fe54)*state.rho; jac.set(Co57, He4, scratch); - scratch = screened_rates(k_He4_Fe54_to_p_Co57_derived)*Y(He4)*state.rho; + scratch = screened_rates(k_He4_Fe54_to_p_nse_Co57_derived)*Y(He4)*state.rho; jac.set(Co57, Fe54, scratch); - scratch = screened_rates(k_p_Fe56_to_Co57)*Y(H1)*state.rho; + scratch = screened_rates(k_p_nse_Fe56_to_Co57)*Y(P_nse)*state.rho; jac.set(Co57, Fe56, scratch); scratch = screened_rates(k_n_Co56_to_Co57)*Y(N)*state.rho; jac.set(Co57, Co56, scratch); - scratch = -screened_rates(k_Co57_to_n_Co56_derived) - screened_rates(k_Co57_to_p_Fe56_derived) - screened_rates(k_p_Co57_to_He4_Fe54)*Y(H1)*state.rho - screened_rates(k_p_Co57_to_Ni58)*Y(H1)*state.rho; + scratch = -screened_rates(k_Co57_to_n_Co56_derived) - screened_rates(k_Co57_to_p_nse_Fe56_derived) - screened_rates(k_p_nse_Co57_to_He4_Fe54)*Y(P_nse)*state.rho - screened_rates(k_p_nse_Co57_to_Ni58)*Y(P_nse)*state.rho; jac.set(Co57, Co57, scratch); - scratch = screened_rates(k_Ni58_to_p_Co57_derived); + scratch = screened_rates(k_Ni58_to_p_nse_Co57_derived); jac.set(Co57, Ni58, scratch); - scratch = -2.0*screened_rates(k_Ni56_n_n_to_Ni58_approx)*Y(Ni56)*Y(N)*amrex::Math::powi<2>(state.rho) - screened_rates(k_n_Ni56_to_p_Co56)*Y(Ni56)*state.rho; + scratch = -2.0*screened_rates(k_Ni56_n_n_to_Ni58_approx)*Y(Ni56)*Y(N)*amrex::Math::powi<2>(state.rho) - screened_rates(k_n_Ni56_to_p_nse_Co56)*Y(Ni56)*state.rho; jac.set(Ni56, N, scratch); - scratch = screened_rates(k_p_Co55_to_Ni56)*Y(Co55)*state.rho + screened_rates(k_p_Co56_to_n_Ni56_derived)*Y(Co56)*state.rho; - jac.set(Ni56, H1, scratch); + scratch = screened_rates(k_p_nse_Co55_to_Ni56)*Y(Co55)*state.rho + screened_rates(k_p_nse_Co56_to_n_Ni56_derived)*Y(Co56)*state.rho; + jac.set(Ni56, P_nse, scratch); scratch = screened_rates(k_He4_Fe52_to_Ni56)*Y(Fe52)*state.rho; jac.set(Ni56, He4, scratch); @@ -2207,13 +2226,13 @@ void jac_nuc(const burn_t& state, scratch = screened_rates(k_He4_Fe52_to_Ni56)*Y(He4)*state.rho; jac.set(Ni56, Fe52, scratch); - scratch = screened_rates(k_p_Co55_to_Ni56)*Y(H1)*state.rho; + scratch = screened_rates(k_p_nse_Co55_to_Ni56)*Y(P_nse)*state.rho; jac.set(Ni56, Co55, scratch); - scratch = screened_rates(k_Co56_to_Ni56) + screened_rates(k_p_Co56_to_n_Ni56_derived)*Y(H1)*state.rho; + scratch = screened_rates(k_Co56_to_Ni56) + screened_rates(k_p_nse_Co56_to_n_Ni56_derived)*Y(P_nse)*state.rho; jac.set(Ni56, Co56, scratch); - scratch = -screened_rates(k_Ni56_n_n_to_Ni58_approx)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho) - screened_rates(k_Ni56_to_Co56) - screened_rates(k_Ni56_to_He4_Fe52_derived) - screened_rates(k_Ni56_to_p_Co55_derived) - screened_rates(k_n_Ni56_to_p_Co56)*Y(N)*state.rho; + scratch = -screened_rates(k_Ni56_n_n_to_Ni58_approx)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho) - screened_rates(k_Ni56_to_Co56) - screened_rates(k_Ni56_to_He4_Fe52_derived) - screened_rates(k_Ni56_to_p_nse_Co55_derived) - screened_rates(k_n_Ni56_to_p_nse_Co56)*Y(N)*state.rho; jac.set(Ni56, Ni56, scratch); scratch = screened_rates(k_Ni58_to_Ni56_n_n_approx); @@ -2222,25 +2241,25 @@ void jac_nuc(const burn_t& state, scratch = 2.0*screened_rates(k_Ni56_n_n_to_Ni58_approx)*Y(Ni56)*Y(N)*amrex::Math::powi<2>(state.rho); jac.set(Ni58, N, scratch); - scratch = screened_rates(k_p_Co57_to_Ni58)*Y(Co57)*state.rho - screened_rates(k_p_Ni58_to_He4_Co55_derived)*Y(Ni58)*state.rho; - jac.set(Ni58, H1, scratch); + scratch = screened_rates(k_p_nse_Co57_to_Ni58)*Y(Co57)*state.rho - screened_rates(k_p_nse_Ni58_to_He4_Co55_derived)*Y(Ni58)*state.rho; + jac.set(Ni58, P_nse, scratch); - scratch = screened_rates(k_He4_Co55_to_p_Ni58)*Y(Co55)*state.rho + screened_rates(k_He4_Fe54_to_Ni58)*Y(Fe54)*state.rho; + scratch = screened_rates(k_He4_Co55_to_p_nse_Ni58)*Y(Co55)*state.rho + screened_rates(k_He4_Fe54_to_Ni58)*Y(Fe54)*state.rho; jac.set(Ni58, He4, scratch); scratch = screened_rates(k_He4_Fe54_to_Ni58)*Y(He4)*state.rho; jac.set(Ni58, Fe54, scratch); - scratch = screened_rates(k_He4_Co55_to_p_Ni58)*Y(He4)*state.rho; + scratch = screened_rates(k_He4_Co55_to_p_nse_Ni58)*Y(He4)*state.rho; jac.set(Ni58, Co55, scratch); - scratch = screened_rates(k_p_Co57_to_Ni58)*Y(H1)*state.rho; + scratch = screened_rates(k_p_nse_Co57_to_Ni58)*Y(P_nse)*state.rho; jac.set(Ni58, Co57, scratch); scratch = screened_rates(k_Ni56_n_n_to_Ni58_approx)*amrex::Math::powi<2>(Y(N))*amrex::Math::powi<2>(state.rho); jac.set(Ni58, Ni56, scratch); - scratch = -screened_rates(k_Ni58_to_He4_Fe54_derived) - screened_rates(k_Ni58_to_Ni56_n_n_approx) - screened_rates(k_Ni58_to_p_Co57_derived) - screened_rates(k_p_Ni58_to_He4_Co55_derived)*Y(H1)*state.rho; + scratch = -screened_rates(k_Ni58_to_He4_Fe54_derived) - screened_rates(k_Ni58_to_Ni56_n_n_approx) - screened_rates(k_Ni58_to_p_nse_Co57_derived) - screened_rates(k_p_nse_Ni58_to_He4_Co55_derived)*Y(P_nse)*state.rho; jac.set(Ni58, Ni58, scratch); diff --git a/networks/he-burn/he-burn-31anp/he-burn-31anp.py b/networks/he-burn/he-burn-31anp/he-burn-31anp.py index b1a3118912..d8587c3147 100644 --- a/networks/he-burn/he-burn-31anp/he-burn-31anp.py +++ b/networks/he-burn/he-burn-31anp/he-burn-31anp.py @@ -117,6 +117,9 @@ net.make_nn_g_approx(intermediate_nuclei=["fe53", "fe55", "ni57"]) net.remove_nuclei(["fe53", "fe55", "ni57"]) +# make all rates with A >= 48 use NSE protons +net.make_nse_protons(48) + print(f"number of nuclei = {len(net.unique_nuclei)}") print(f"number of ReacLib rates = {len(net.reaclib_rates)}") print(f"number of tabular rates = {len(net.tabular_rates)}") diff --git a/networks/he-burn/he-burn-31anp/inputs.burn_cell.VODE b/networks/he-burn/he-burn-31anp/inputs.burn_cell.VODE index 98dc0cc869..29b7b866cc 100644 --- a/networks/he-burn/he-burn-31anp/inputs.burn_cell.VODE +++ b/networks/he-burn/he-burn-31anp/inputs.burn_cell.VODE @@ -54,3 +54,4 @@ unit_test.X27 = 0.0 unit_test.X28 = 0.0 unit_test.X29 = 0.0 unit_test.X30 = 0.0 +unit_test.X31 = 0.0 diff --git a/networks/he-burn/he-burn-31anp/partition_functions.H b/networks/he-burn/he-burn-31anp/partition_functions.H index d34ee02f2e..998d4b3d07 100644 --- a/networks/he-burn/he-burn-31anp/partition_functions.H +++ b/networks/he-burn/he-burn-31anp/partition_functions.H @@ -397,6 +397,7 @@ constexpr amrex::Real get_spin_state(const int inuc) { case N: case H1: + case P_nse: case N13: case P31: spin = 2; diff --git a/networks/he-burn/he-burn-31anp/pynucastro.net b/networks/he-burn/he-burn-31anp/pynucastro.net index 7773fdfd67..9db593b3f4 100644 --- a/networks/he-burn/he-burn-31anp/pynucastro.net +++ b/networks/he-burn/he-burn-31anp/pynucastro.net @@ -1,5 +1,6 @@ neutron n 1.0 0.0 hydrogen-1 H1 1.0 1.0 +proton-nse P_nse 1.0 1.0 helium-4 He4 4.0 2.0 carbon-12 C12 12.0 6.0 nitrogen-13 N13 13.0 7.0 diff --git a/networks/he-burn/he-burn-31anp/reaclib_rates.H b/networks/he-burn/he-burn-31anp/reaclib_rates.H index 780c60edd6..8762d5f5d4 100644 --- a/networks/he-burn/he-burn-31anp/reaclib_rates.H +++ b/networks/he-burn/he-burn-31anp/reaclib_rates.H @@ -764,9 +764,9 @@ void rate_He4_Cr48_to_Fe52(const tf_t& tfactors, amrex::Real& rate, amrex::Real& template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_p_Mn51_to_Fe52(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { +void rate_p_nse_Mn51_to_Fe52(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { - // Mn51 + p --> Fe52 + // Mn51 + p_nse --> Fe52 rate = 0.0; drate_dT = 0.0; @@ -860,9 +860,9 @@ void rate_He4_Fe52_to_Ni56(const tf_t& tfactors, amrex::Real& rate, amrex::Real& template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_p_Co55_to_Ni56(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { +void rate_p_nse_Co55_to_Ni56(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { - // Co55 + p --> Ni56 + // Co55 + p_nse --> Ni56 rate = 0.0; drate_dT = 0.0; @@ -1346,9 +1346,9 @@ void rate_p_P31_to_He4_Si28(const tf_t& tfactors, amrex::Real& rate, amrex::Real template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_He4_Cr48_to_p_Mn51(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { +void rate_He4_Cr48_to_p_nse_Mn51(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { - // Cr48 + He4 --> p + Mn51 + // Cr48 + He4 --> p_nse + Mn51 rate = 0.0; drate_dT = 0.0; @@ -1378,9 +1378,9 @@ void rate_He4_Cr48_to_p_Mn51(const tf_t& tfactors, amrex::Real& rate, amrex::Rea template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_He4_Fe52_to_p_Co55(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { +void rate_He4_Fe52_to_p_nse_Co55(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { - // Fe52 + He4 --> p + Co55 + // Fe52 + He4 --> p_nse + Co55 rate = 0.0; drate_dT = 0.0; @@ -1572,9 +1572,9 @@ void rate_C12_O16_to_Si28_modified(const tf_t& tfactors, amrex::Real& rate, amre template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_p_Fe54_to_Co55(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { +void rate_p_nse_Fe54_to_Co55(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { - // Fe54 + p --> Co55 + // Fe54 + p_nse --> Co55 rate = 0.0; drate_dT = 0.0; @@ -1636,9 +1636,9 @@ void rate_He4_Fe54_to_Ni58(const tf_t& tfactors, amrex::Real& rate, amrex::Real& template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_p_Fe56_to_Co57(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { +void rate_p_nse_Fe56_to_Co57(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { - // Fe56 + p --> Co57 + // Fe56 + p_nse --> Co57 rate = 0.0; drate_dT = 0.0; @@ -1732,9 +1732,9 @@ void rate_n_Co56_to_Co57(const tf_t& tfactors, amrex::Real& rate, amrex::Real& d template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_p_Co57_to_Ni58(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { +void rate_p_nse_Co57_to_Ni58(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { - // Co57 + p --> Ni58 + // Co57 + p_nse --> Ni58 rate = 0.0; drate_dT = 0.0; @@ -1764,9 +1764,9 @@ void rate_p_Co57_to_Ni58(const tf_t& tfactors, amrex::Real& rate, amrex::Real& d template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_He4_Mn51_to_p_Fe54(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { +void rate_He4_Mn51_to_p_nse_Fe54(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { - // Mn51 + He4 --> p + Fe54 + // Mn51 + He4 --> p_nse + Fe54 rate = 0.0; drate_dT = 0.0; @@ -1796,9 +1796,9 @@ void rate_He4_Mn51_to_p_Fe54(const tf_t& tfactors, amrex::Real& rate, amrex::Rea template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_He4_Co55_to_p_Ni58(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { +void rate_He4_Co55_to_p_nse_Ni58(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { - // Co55 + He4 --> p + Ni58 + // Co55 + He4 --> p_nse + Ni58 rate = 0.0; drate_dT = 0.0; @@ -1828,9 +1828,9 @@ void rate_He4_Co55_to_p_Ni58(const tf_t& tfactors, amrex::Real& rate, amrex::Rea template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_n_Co56_to_p_Fe56(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { +void rate_n_Co56_to_p_nse_Fe56(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { - // Co56 + n --> p + Fe56 + // Co56 + n --> p_nse + Fe56 rate = 0.0; drate_dT = 0.0; @@ -1860,9 +1860,9 @@ void rate_n_Co56_to_p_Fe56(const tf_t& tfactors, amrex::Real& rate, amrex::Real& template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_p_Co57_to_He4_Fe54(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { +void rate_p_nse_Co57_to_He4_Fe54(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { - // Co57 + p --> He4 + Fe54 + // Co57 + p_nse --> He4 + Fe54 rate = 0.0; drate_dT = 0.0; @@ -1892,9 +1892,9 @@ void rate_p_Co57_to_He4_Fe54(const tf_t& tfactors, amrex::Real& rate, amrex::Rea template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_n_Ni56_to_p_Co56(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { +void rate_n_Ni56_to_p_nse_Co56(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT) { - // Ni56 + n --> p + Co56 + // Ni56 + n --> p_nse + Co56 rate = 0.0; drate_dT = 0.0; @@ -3678,9 +3678,9 @@ void rate_S32_to_He4_Si28_derived(const tf_t& tfactors, amrex::Real& rate, amrex template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_Fe52_to_p_Mn51_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { +void rate_Fe52_to_p_nse_Mn51_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { - // Fe52 --> p + Mn51 + // Fe52 --> p_nse + Mn51 rate = 0.0; drate_dT = 0.0; @@ -3711,19 +3711,19 @@ void rate_Fe52_to_p_Mn51_derived(const tf_t& tfactors, amrex::Real& rate, amrex: // interpolating Mn51 partition function get_partition_function_cached(Mn51, tfactors, pf_cache, Mn51_pf, dMn51_pf_dT); - amrex::Real p_pf, dp_pf_dT; - // setting p partition function to 1.0 by default, independent of T - p_pf = 1.0_rt; - dp_pf_dT = 0.0_rt; + amrex::Real p_nse_pf, dp_nse_pf_dT; + // setting p_nse partition function to 1.0 by default, independent of T + p_nse_pf = 1.0_rt; + dp_nse_pf_dT = 0.0_rt; amrex::Real Fe52_pf, dFe52_pf_dT; // interpolating Fe52 partition function get_partition_function_cached(Fe52, tfactors, pf_cache, Fe52_pf, dFe52_pf_dT); - amrex::Real z_r = p_pf * Mn51_pf; + amrex::Real z_r = p_nse_pf * Mn51_pf; amrex::Real z_p = Fe52_pf; - amrex::Real dz_r_dT = Mn51_pf * dp_pf_dT + p_pf * dMn51_pf_dT; + amrex::Real dz_r_dT = Mn51_pf * dp_nse_pf_dT + p_nse_pf * dMn51_pf_dT; amrex::Real dz_p_dT = dFe52_pf_dT; amrex::Real dzterm_dT = (z_p * dz_r_dT - z_r * dz_p_dT) / (z_p * z_p); @@ -3849,9 +3849,9 @@ void rate_Co55_to_He4_Mn51_derived(const tf_t& tfactors, amrex::Real& rate, amre template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_Ni56_to_p_Co55_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { +void rate_Ni56_to_p_nse_Co55_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { - // Ni56 --> p + Co55 + // Ni56 --> p_nse + Co55 rate = 0.0; drate_dT = 0.0; @@ -3882,19 +3882,19 @@ void rate_Ni56_to_p_Co55_derived(const tf_t& tfactors, amrex::Real& rate, amrex: // interpolating Ni56 partition function get_partition_function_cached(Ni56, tfactors, pf_cache, Ni56_pf, dNi56_pf_dT); - amrex::Real p_pf, dp_pf_dT; - // setting p partition function to 1.0 by default, independent of T - p_pf = 1.0_rt; - dp_pf_dT = 0.0_rt; + amrex::Real p_nse_pf, dp_nse_pf_dT; + // setting p_nse partition function to 1.0 by default, independent of T + p_nse_pf = 1.0_rt; + dp_nse_pf_dT = 0.0_rt; amrex::Real Co55_pf, dCo55_pf_dT; // interpolating Co55 partition function get_partition_function_cached(Co55, tfactors, pf_cache, Co55_pf, dCo55_pf_dT); - amrex::Real z_r = p_pf * Co55_pf; + amrex::Real z_r = p_nse_pf * Co55_pf; amrex::Real z_p = Ni56_pf; - amrex::Real dz_r_dT = Co55_pf * dp_pf_dT + p_pf * dCo55_pf_dT; + amrex::Real dz_r_dT = Co55_pf * dp_nse_pf_dT + p_nse_pf * dCo55_pf_dT; amrex::Real dz_p_dT = dNi56_pf_dT; amrex::Real dzterm_dT = (z_p * dz_r_dT - z_r * dz_p_dT) / (z_p * z_p); @@ -4464,9 +4464,9 @@ void rate_He4_Si28_to_p_P31_derived(const tf_t& tfactors, amrex::Real& rate, amr template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_p_Mn51_to_He4_Cr48_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { +void rate_p_nse_Mn51_to_He4_Cr48_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { - // Mn51 + p --> He4 + Cr48 + // Mn51 + p_nse --> He4 + Cr48 rate = 0.0; drate_dT = 0.0; @@ -4493,10 +4493,10 @@ void rate_p_Mn51_to_He4_Cr48_derived(const tf_t& tfactors, amrex::Real& rate, am } - amrex::Real p_pf, dp_pf_dT; - // setting p partition function to 1.0 by default, independent of T - p_pf = 1.0_rt; - dp_pf_dT = 0.0_rt; + amrex::Real p_nse_pf, dp_nse_pf_dT; + // setting p_nse partition function to 1.0 by default, independent of T + p_nse_pf = 1.0_rt; + dp_nse_pf_dT = 0.0_rt; amrex::Real Cr48_pf, dCr48_pf_dT; // interpolating Cr48 partition function @@ -4512,10 +4512,10 @@ void rate_p_Mn51_to_He4_Cr48_derived(const tf_t& tfactors, amrex::Real& rate, am get_partition_function_cached(Mn51, tfactors, pf_cache, Mn51_pf, dMn51_pf_dT); amrex::Real z_r = He4_pf * Cr48_pf; - amrex::Real z_p = p_pf * Mn51_pf; + amrex::Real z_p = p_nse_pf * Mn51_pf; amrex::Real dz_r_dT = Cr48_pf * dHe4_pf_dT + He4_pf * dCr48_pf_dT; - amrex::Real dz_p_dT = Mn51_pf * dp_pf_dT + p_pf * dMn51_pf_dT; + amrex::Real dz_p_dT = Mn51_pf * dp_nse_pf_dT + p_nse_pf * dMn51_pf_dT; amrex::Real dzterm_dT = (z_p * dz_r_dT - z_r * dz_p_dT) / (z_p * z_p); @@ -4526,9 +4526,9 @@ void rate_p_Mn51_to_He4_Cr48_derived(const tf_t& tfactors, amrex::Real& rate, am template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_p_Co55_to_He4_Fe52_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { +void rate_p_nse_Co55_to_He4_Fe52_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { - // Co55 + p --> He4 + Fe52 + // Co55 + p_nse --> He4 + Fe52 rate = 0.0; drate_dT = 0.0; @@ -4555,10 +4555,10 @@ void rate_p_Co55_to_He4_Fe52_derived(const tf_t& tfactors, amrex::Real& rate, am } - amrex::Real p_pf, dp_pf_dT; - // setting p partition function to 1.0 by default, independent of T - p_pf = 1.0_rt; - dp_pf_dT = 0.0_rt; + amrex::Real p_nse_pf, dp_nse_pf_dT; + // setting p_nse partition function to 1.0 by default, independent of T + p_nse_pf = 1.0_rt; + dp_nse_pf_dT = 0.0_rt; amrex::Real He4_pf, dHe4_pf_dT; // setting He4 partition function to 1.0 by default, independent of T @@ -4574,10 +4574,10 @@ void rate_p_Co55_to_He4_Fe52_derived(const tf_t& tfactors, amrex::Real& rate, am get_partition_function_cached(Fe52, tfactors, pf_cache, Fe52_pf, dFe52_pf_dT); amrex::Real z_r = He4_pf * Fe52_pf; - amrex::Real z_p = p_pf * Co55_pf; + amrex::Real z_p = p_nse_pf * Co55_pf; amrex::Real dz_r_dT = Fe52_pf * dHe4_pf_dT + He4_pf * dFe52_pf_dT; - amrex::Real dz_p_dT = Co55_pf * dp_pf_dT + p_pf * dCo55_pf_dT; + amrex::Real dz_p_dT = Co55_pf * dp_nse_pf_dT + p_nse_pf * dCo55_pf_dT; amrex::Real dzterm_dT = (z_p * dz_r_dT - z_r * dz_p_dT) / (z_p * z_p); @@ -4588,9 +4588,9 @@ void rate_p_Co55_to_He4_Fe52_derived(const tf_t& tfactors, amrex::Real& rate, am template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_Co55_to_p_Fe54_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { +void rate_Co55_to_p_nse_Fe54_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { - // Co55 --> p + Fe54 + // Co55 --> p_nse + Fe54 rate = 0.0; drate_dT = 0.0; @@ -4617,10 +4617,10 @@ void rate_Co55_to_p_Fe54_derived(const tf_t& tfactors, amrex::Real& rate, amrex: } - amrex::Real p_pf, dp_pf_dT; - // setting p partition function to 1.0 by default, independent of T - p_pf = 1.0_rt; - dp_pf_dT = 0.0_rt; + amrex::Real p_nse_pf, dp_nse_pf_dT; + // setting p_nse partition function to 1.0 by default, independent of T + p_nse_pf = 1.0_rt; + dp_nse_pf_dT = 0.0_rt; amrex::Real Co55_pf, dCo55_pf_dT; // interpolating Co55 partition function @@ -4630,10 +4630,10 @@ void rate_Co55_to_p_Fe54_derived(const tf_t& tfactors, amrex::Real& rate, amrex: // interpolating Fe54 partition function get_partition_function_cached(Fe54, tfactors, pf_cache, Fe54_pf, dFe54_pf_dT); - amrex::Real z_r = p_pf * Fe54_pf; + amrex::Real z_r = p_nse_pf * Fe54_pf; amrex::Real z_p = Co55_pf; - amrex::Real dz_r_dT = Fe54_pf * dp_pf_dT + p_pf * dFe54_pf_dT; + amrex::Real dz_r_dT = Fe54_pf * dp_nse_pf_dT + p_nse_pf * dFe54_pf_dT; amrex::Real dz_p_dT = dCo55_pf_dT; amrex::Real dzterm_dT = (z_p * dz_r_dT - z_r * dz_p_dT) / (z_p * z_p); @@ -4759,9 +4759,9 @@ void rate_Co57_to_n_Co56_derived(const tf_t& tfactors, amrex::Real& rate, amrex: template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_Co57_to_p_Fe56_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { +void rate_Co57_to_p_nse_Fe56_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { - // Co57 --> p + Fe56 + // Co57 --> p_nse + Fe56 rate = 0.0; drate_dT = 0.0; @@ -4788,10 +4788,10 @@ void rate_Co57_to_p_Fe56_derived(const tf_t& tfactors, amrex::Real& rate, amrex: } - amrex::Real p_pf, dp_pf_dT; - // setting p partition function to 1.0 by default, independent of T - p_pf = 1.0_rt; - dp_pf_dT = 0.0_rt; + amrex::Real p_nse_pf, dp_nse_pf_dT; + // setting p_nse partition function to 1.0 by default, independent of T + p_nse_pf = 1.0_rt; + dp_nse_pf_dT = 0.0_rt; amrex::Real Fe56_pf, dFe56_pf_dT; // interpolating Fe56 partition function @@ -4801,10 +4801,10 @@ void rate_Co57_to_p_Fe56_derived(const tf_t& tfactors, amrex::Real& rate, amrex: // interpolating Co57 partition function get_partition_function_cached(Co57, tfactors, pf_cache, Co57_pf, dCo57_pf_dT); - amrex::Real z_r = p_pf * Fe56_pf; + amrex::Real z_r = p_nse_pf * Fe56_pf; amrex::Real z_p = Co57_pf; - amrex::Real dz_r_dT = Fe56_pf * dp_pf_dT + p_pf * dFe56_pf_dT; + amrex::Real dz_r_dT = Fe56_pf * dp_nse_pf_dT + p_nse_pf * dFe56_pf_dT; amrex::Real dz_p_dT = dCo57_pf_dT; amrex::Real dzterm_dT = (z_p * dz_r_dT - z_r * dz_p_dT) / (z_p * z_p); @@ -4816,9 +4816,9 @@ void rate_Co57_to_p_Fe56_derived(const tf_t& tfactors, amrex::Real& rate, amrex: template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_Ni58_to_p_Co57_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { +void rate_Ni58_to_p_nse_Co57_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { - // Ni58 --> p + Co57 + // Ni58 --> p_nse + Co57 rate = 0.0; drate_dT = 0.0; @@ -4845,10 +4845,10 @@ void rate_Ni58_to_p_Co57_derived(const tf_t& tfactors, amrex::Real& rate, amrex: } - amrex::Real p_pf, dp_pf_dT; - // setting p partition function to 1.0 by default, independent of T - p_pf = 1.0_rt; - dp_pf_dT = 0.0_rt; + amrex::Real p_nse_pf, dp_nse_pf_dT; + // setting p_nse partition function to 1.0 by default, independent of T + p_nse_pf = 1.0_rt; + dp_nse_pf_dT = 0.0_rt; amrex::Real Ni58_pf, dNi58_pf_dT; // interpolating Ni58 partition function @@ -4858,10 +4858,10 @@ void rate_Ni58_to_p_Co57_derived(const tf_t& tfactors, amrex::Real& rate, amrex: // interpolating Co57 partition function get_partition_function_cached(Co57, tfactors, pf_cache, Co57_pf, dCo57_pf_dT); - amrex::Real z_r = p_pf * Co57_pf; + amrex::Real z_r = p_nse_pf * Co57_pf; amrex::Real z_p = Ni58_pf; - amrex::Real dz_r_dT = Co57_pf * dp_pf_dT + p_pf * dCo57_pf_dT; + amrex::Real dz_r_dT = Co57_pf * dp_nse_pf_dT + p_nse_pf * dCo57_pf_dT; amrex::Real dz_p_dT = dNi58_pf_dT; amrex::Real dzterm_dT = (z_p * dz_r_dT - z_r * dz_p_dT) / (z_p * z_p); @@ -4930,9 +4930,9 @@ void rate_Ni58_to_He4_Fe54_derived(const tf_t& tfactors, amrex::Real& rate, amre template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_p_Fe54_to_He4_Mn51_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { +void rate_p_nse_Fe54_to_He4_Mn51_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { - // Fe54 + p --> He4 + Mn51 + // Fe54 + p_nse --> He4 + Mn51 rate = 0.0; drate_dT = 0.0; @@ -4959,10 +4959,10 @@ void rate_p_Fe54_to_He4_Mn51_derived(const tf_t& tfactors, amrex::Real& rate, am } - amrex::Real p_pf, dp_pf_dT; - // setting p partition function to 1.0 by default, independent of T - p_pf = 1.0_rt; - dp_pf_dT = 0.0_rt; + amrex::Real p_nse_pf, dp_nse_pf_dT; + // setting p_nse partition function to 1.0 by default, independent of T + p_nse_pf = 1.0_rt; + dp_nse_pf_dT = 0.0_rt; amrex::Real Mn51_pf, dMn51_pf_dT; // interpolating Mn51 partition function @@ -4978,10 +4978,10 @@ void rate_p_Fe54_to_He4_Mn51_derived(const tf_t& tfactors, amrex::Real& rate, am get_partition_function_cached(Fe54, tfactors, pf_cache, Fe54_pf, dFe54_pf_dT); amrex::Real z_r = He4_pf * Mn51_pf; - amrex::Real z_p = p_pf * Fe54_pf; + amrex::Real z_p = p_nse_pf * Fe54_pf; amrex::Real dz_r_dT = Mn51_pf * dHe4_pf_dT + He4_pf * dMn51_pf_dT; - amrex::Real dz_p_dT = Fe54_pf * dp_pf_dT + p_pf * dFe54_pf_dT; + amrex::Real dz_p_dT = Fe54_pf * dp_nse_pf_dT + p_nse_pf * dFe54_pf_dT; amrex::Real dzterm_dT = (z_p * dz_r_dT - z_r * dz_p_dT) / (z_p * z_p); @@ -4992,9 +4992,9 @@ void rate_p_Fe54_to_He4_Mn51_derived(const tf_t& tfactors, amrex::Real& rate, am template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_He4_Fe54_to_p_Co57_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { +void rate_He4_Fe54_to_p_nse_Co57_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { - // Fe54 + He4 --> p + Co57 + // Fe54 + He4 --> p_nse + Co57 rate = 0.0; drate_dT = 0.0; @@ -5021,10 +5021,10 @@ void rate_He4_Fe54_to_p_Co57_derived(const tf_t& tfactors, amrex::Real& rate, am } - amrex::Real p_pf, dp_pf_dT; - // setting p partition function to 1.0 by default, independent of T - p_pf = 1.0_rt; - dp_pf_dT = 0.0_rt; + amrex::Real p_nse_pf, dp_nse_pf_dT; + // setting p_nse partition function to 1.0 by default, independent of T + p_nse_pf = 1.0_rt; + dp_nse_pf_dT = 0.0_rt; amrex::Real He4_pf, dHe4_pf_dT; // setting He4 partition function to 1.0 by default, independent of T @@ -5039,10 +5039,10 @@ void rate_He4_Fe54_to_p_Co57_derived(const tf_t& tfactors, amrex::Real& rate, am // interpolating Fe54 partition function get_partition_function_cached(Fe54, tfactors, pf_cache, Fe54_pf, dFe54_pf_dT); - amrex::Real z_r = p_pf * Co57_pf; + amrex::Real z_r = p_nse_pf * Co57_pf; amrex::Real z_p = He4_pf * Fe54_pf; - amrex::Real dz_r_dT = Co57_pf * dp_pf_dT + p_pf * dCo57_pf_dT; + amrex::Real dz_r_dT = Co57_pf * dp_nse_pf_dT + p_nse_pf * dCo57_pf_dT; amrex::Real dz_p_dT = Fe54_pf * dHe4_pf_dT + He4_pf * dFe54_pf_dT; amrex::Real dzterm_dT = (z_p * dz_r_dT - z_r * dz_p_dT) / (z_p * z_p); @@ -5054,9 +5054,9 @@ void rate_He4_Fe54_to_p_Co57_derived(const tf_t& tfactors, amrex::Real& rate, am template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_p_Fe56_to_n_Co56_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { +void rate_p_nse_Fe56_to_n_Co56_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { - // Fe56 + p --> n + Co56 + // Fe56 + p_nse --> n + Co56 rate = 0.0; drate_dT = 0.0; @@ -5088,10 +5088,10 @@ void rate_p_Fe56_to_n_Co56_derived(const tf_t& tfactors, amrex::Real& rate, amre n_pf = 1.0_rt; dn_pf_dT = 0.0_rt; - amrex::Real p_pf, dp_pf_dT; - // setting p partition function to 1.0 by default, independent of T - p_pf = 1.0_rt; - dp_pf_dT = 0.0_rt; + amrex::Real p_nse_pf, dp_nse_pf_dT; + // setting p_nse partition function to 1.0 by default, independent of T + p_nse_pf = 1.0_rt; + dp_nse_pf_dT = 0.0_rt; amrex::Real Fe56_pf, dFe56_pf_dT; // interpolating Fe56 partition function @@ -5102,10 +5102,10 @@ void rate_p_Fe56_to_n_Co56_derived(const tf_t& tfactors, amrex::Real& rate, amre get_partition_function_cached(Co56, tfactors, pf_cache, Co56_pf, dCo56_pf_dT); amrex::Real z_r = n_pf * Co56_pf; - amrex::Real z_p = p_pf * Fe56_pf; + amrex::Real z_p = p_nse_pf * Fe56_pf; amrex::Real dz_r_dT = Co56_pf * dn_pf_dT + n_pf * dCo56_pf_dT; - amrex::Real dz_p_dT = Fe56_pf * dp_pf_dT + p_pf * dFe56_pf_dT; + amrex::Real dz_p_dT = Fe56_pf * dp_nse_pf_dT + p_nse_pf * dFe56_pf_dT; amrex::Real dzterm_dT = (z_p * dz_r_dT - z_r * dz_p_dT) / (z_p * z_p); @@ -5116,9 +5116,9 @@ void rate_p_Fe56_to_n_Co56_derived(const tf_t& tfactors, amrex::Real& rate, amre template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_p_Co56_to_n_Ni56_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { +void rate_p_nse_Co56_to_n_Ni56_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { - // Co56 + p --> n + Ni56 + // Co56 + p_nse --> n + Ni56 rate = 0.0; drate_dT = 0.0; @@ -5154,20 +5154,20 @@ void rate_p_Co56_to_n_Ni56_derived(const tf_t& tfactors, amrex::Real& rate, amre // interpolating Ni56 partition function get_partition_function_cached(Ni56, tfactors, pf_cache, Ni56_pf, dNi56_pf_dT); - amrex::Real p_pf, dp_pf_dT; - // setting p partition function to 1.0 by default, independent of T - p_pf = 1.0_rt; - dp_pf_dT = 0.0_rt; + amrex::Real p_nse_pf, dp_nse_pf_dT; + // setting p_nse partition function to 1.0 by default, independent of T + p_nse_pf = 1.0_rt; + dp_nse_pf_dT = 0.0_rt; amrex::Real Co56_pf, dCo56_pf_dT; // interpolating Co56 partition function get_partition_function_cached(Co56, tfactors, pf_cache, Co56_pf, dCo56_pf_dT); amrex::Real z_r = n_pf * Ni56_pf; - amrex::Real z_p = p_pf * Co56_pf; + amrex::Real z_p = p_nse_pf * Co56_pf; amrex::Real dz_r_dT = Ni56_pf * dn_pf_dT + n_pf * dNi56_pf_dT; - amrex::Real dz_p_dT = Co56_pf * dp_pf_dT + p_pf * dCo56_pf_dT; + amrex::Real dz_p_dT = Co56_pf * dp_nse_pf_dT + p_nse_pf * dCo56_pf_dT; amrex::Real dzterm_dT = (z_p * dz_r_dT - z_r * dz_p_dT) / (z_p * z_p); @@ -5178,9 +5178,9 @@ void rate_p_Co56_to_n_Ni56_derived(const tf_t& tfactors, amrex::Real& rate, amre template AMREX_GPU_HOST_DEVICE AMREX_INLINE -void rate_p_Ni58_to_He4_Co55_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { +void rate_p_nse_Ni58_to_He4_Co55_derived(const tf_t& tfactors, amrex::Real& rate, amrex::Real& drate_dT, [[maybe_unused]] part_fun::pf_cache_t& pf_cache) { - // Ni58 + p --> He4 + Co55 + // Ni58 + p_nse --> He4 + Co55 rate = 0.0; drate_dT = 0.0; @@ -5207,10 +5207,10 @@ void rate_p_Ni58_to_He4_Co55_derived(const tf_t& tfactors, amrex::Real& rate, am } - amrex::Real p_pf, dp_pf_dT; - // setting p partition function to 1.0 by default, independent of T - p_pf = 1.0_rt; - dp_pf_dT = 0.0_rt; + amrex::Real p_nse_pf, dp_nse_pf_dT; + // setting p_nse partition function to 1.0 by default, independent of T + p_nse_pf = 1.0_rt; + dp_nse_pf_dT = 0.0_rt; amrex::Real He4_pf, dHe4_pf_dT; // setting He4 partition function to 1.0 by default, independent of T @@ -5226,10 +5226,10 @@ void rate_p_Ni58_to_He4_Co55_derived(const tf_t& tfactors, amrex::Real& rate, am get_partition_function_cached(Ni58, tfactors, pf_cache, Ni58_pf, dNi58_pf_dT); amrex::Real z_r = He4_pf * Co55_pf; - amrex::Real z_p = p_pf * Ni58_pf; + amrex::Real z_p = p_nse_pf * Ni58_pf; amrex::Real dz_r_dT = Co55_pf * dHe4_pf_dT + He4_pf * dCo55_pf_dT; - amrex::Real dz_p_dT = Ni58_pf * dp_pf_dT + p_pf * dNi58_pf_dT; + amrex::Real dz_p_dT = Ni58_pf * dp_nse_pf_dT + p_nse_pf * dNi58_pf_dT; amrex::Real dzterm_dT = (z_p * dz_r_dT - z_r * dz_p_dT) / (z_p * z_p); @@ -6818,10 +6818,10 @@ fill_reaclib_rates(const tf_t& tfactors, T& rate_eval) rate_eval.dscreened_rates_dT(k_He4_Cr48_to_Fe52) = drate_dT; } - rate_p_Mn51_to_Fe52(tfactors, rate, drate_dT); - rate_eval.screened_rates(k_p_Mn51_to_Fe52) = rate; + rate_p_nse_Mn51_to_Fe52(tfactors, rate, drate_dT); + rate_eval.screened_rates(k_p_nse_Mn51_to_Fe52) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_p_Mn51_to_Fe52) = drate_dT; + rate_eval.dscreened_rates_dT(k_p_nse_Mn51_to_Fe52) = drate_dT; } rate_He4_Mn51_to_Co55(tfactors, rate, drate_dT); @@ -6836,10 +6836,10 @@ fill_reaclib_rates(const tf_t& tfactors, T& rate_eval) rate_eval.dscreened_rates_dT(k_He4_Fe52_to_Ni56) = drate_dT; } - rate_p_Co55_to_Ni56(tfactors, rate, drate_dT); - rate_eval.screened_rates(k_p_Co55_to_Ni56) = rate; + rate_p_nse_Co55_to_Ni56(tfactors, rate, drate_dT); + rate_eval.screened_rates(k_p_nse_Co55_to_Ni56) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_p_Co55_to_Ni56) = drate_dT; + rate_eval.dscreened_rates_dT(k_p_nse_Co55_to_Ni56) = drate_dT; } rate_C12_C12_to_p_Na23(tfactors, rate, drate_dT); @@ -6908,16 +6908,16 @@ fill_reaclib_rates(const tf_t& tfactors, T& rate_eval) rate_eval.dscreened_rates_dT(k_p_P31_to_He4_Si28) = drate_dT; } - rate_He4_Cr48_to_p_Mn51(tfactors, rate, drate_dT); - rate_eval.screened_rates(k_He4_Cr48_to_p_Mn51) = rate; + rate_He4_Cr48_to_p_nse_Mn51(tfactors, rate, drate_dT); + rate_eval.screened_rates(k_He4_Cr48_to_p_nse_Mn51) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_He4_Cr48_to_p_Mn51) = drate_dT; + rate_eval.dscreened_rates_dT(k_He4_Cr48_to_p_nse_Mn51) = drate_dT; } - rate_He4_Fe52_to_p_Co55(tfactors, rate, drate_dT); - rate_eval.screened_rates(k_He4_Fe52_to_p_Co55) = rate; + rate_He4_Fe52_to_p_nse_Co55(tfactors, rate, drate_dT); + rate_eval.screened_rates(k_He4_Fe52_to_p_nse_Co55) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_He4_Fe52_to_p_Co55) = drate_dT; + rate_eval.dscreened_rates_dT(k_He4_Fe52_to_p_nse_Co55) = drate_dT; } rate_He4_He4_He4_to_C12(tfactors, rate, drate_dT); @@ -6944,10 +6944,10 @@ fill_reaclib_rates(const tf_t& tfactors, T& rate_eval) rate_eval.dscreened_rates_dT(k_C12_O16_to_Si28_modified) = drate_dT; } - rate_p_Fe54_to_Co55(tfactors, rate, drate_dT); - rate_eval.screened_rates(k_p_Fe54_to_Co55) = rate; + rate_p_nse_Fe54_to_Co55(tfactors, rate, drate_dT); + rate_eval.screened_rates(k_p_nse_Fe54_to_Co55) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_p_Fe54_to_Co55) = drate_dT; + rate_eval.dscreened_rates_dT(k_p_nse_Fe54_to_Co55) = drate_dT; } rate_He4_Fe54_to_Ni58(tfactors, rate, drate_dT); @@ -6956,10 +6956,10 @@ fill_reaclib_rates(const tf_t& tfactors, T& rate_eval) rate_eval.dscreened_rates_dT(k_He4_Fe54_to_Ni58) = drate_dT; } - rate_p_Fe56_to_Co57(tfactors, rate, drate_dT); - rate_eval.screened_rates(k_p_Fe56_to_Co57) = rate; + rate_p_nse_Fe56_to_Co57(tfactors, rate, drate_dT); + rate_eval.screened_rates(k_p_nse_Fe56_to_Co57) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_p_Fe56_to_Co57) = drate_dT; + rate_eval.dscreened_rates_dT(k_p_nse_Fe56_to_Co57) = drate_dT; } rate_n_Co55_to_Co56(tfactors, rate, drate_dT); @@ -6974,40 +6974,40 @@ fill_reaclib_rates(const tf_t& tfactors, T& rate_eval) rate_eval.dscreened_rates_dT(k_n_Co56_to_Co57) = drate_dT; } - rate_p_Co57_to_Ni58(tfactors, rate, drate_dT); - rate_eval.screened_rates(k_p_Co57_to_Ni58) = rate; + rate_p_nse_Co57_to_Ni58(tfactors, rate, drate_dT); + rate_eval.screened_rates(k_p_nse_Co57_to_Ni58) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_p_Co57_to_Ni58) = drate_dT; + rate_eval.dscreened_rates_dT(k_p_nse_Co57_to_Ni58) = drate_dT; } - rate_He4_Mn51_to_p_Fe54(tfactors, rate, drate_dT); - rate_eval.screened_rates(k_He4_Mn51_to_p_Fe54) = rate; + rate_He4_Mn51_to_p_nse_Fe54(tfactors, rate, drate_dT); + rate_eval.screened_rates(k_He4_Mn51_to_p_nse_Fe54) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_He4_Mn51_to_p_Fe54) = drate_dT; + rate_eval.dscreened_rates_dT(k_He4_Mn51_to_p_nse_Fe54) = drate_dT; } - rate_He4_Co55_to_p_Ni58(tfactors, rate, drate_dT); - rate_eval.screened_rates(k_He4_Co55_to_p_Ni58) = rate; + rate_He4_Co55_to_p_nse_Ni58(tfactors, rate, drate_dT); + rate_eval.screened_rates(k_He4_Co55_to_p_nse_Ni58) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_He4_Co55_to_p_Ni58) = drate_dT; + rate_eval.dscreened_rates_dT(k_He4_Co55_to_p_nse_Ni58) = drate_dT; } - rate_n_Co56_to_p_Fe56(tfactors, rate, drate_dT); - rate_eval.screened_rates(k_n_Co56_to_p_Fe56) = rate; + rate_n_Co56_to_p_nse_Fe56(tfactors, rate, drate_dT); + rate_eval.screened_rates(k_n_Co56_to_p_nse_Fe56) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_n_Co56_to_p_Fe56) = drate_dT; + rate_eval.dscreened_rates_dT(k_n_Co56_to_p_nse_Fe56) = drate_dT; } - rate_p_Co57_to_He4_Fe54(tfactors, rate, drate_dT); - rate_eval.screened_rates(k_p_Co57_to_He4_Fe54) = rate; + rate_p_nse_Co57_to_He4_Fe54(tfactors, rate, drate_dT); + rate_eval.screened_rates(k_p_nse_Co57_to_He4_Fe54) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_p_Co57_to_He4_Fe54) = drate_dT; + rate_eval.dscreened_rates_dT(k_p_nse_Co57_to_He4_Fe54) = drate_dT; } - rate_n_Ni56_to_p_Co56(tfactors, rate, drate_dT); - rate_eval.screened_rates(k_n_Ni56_to_p_Co56) = rate; + rate_n_Ni56_to_p_nse_Co56(tfactors, rate, drate_dT); + rate_eval.screened_rates(k_n_Ni56_to_p_nse_Co56) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_n_Ni56_to_p_Co56) = drate_dT; + rate_eval.dscreened_rates_dT(k_n_Ni56_to_p_nse_Co56) = drate_dT; } rate_He4_S32_to_Ar36_removed(tfactors, rate, drate_dT); @@ -7190,10 +7190,10 @@ fill_reaclib_rates(const tf_t& tfactors, T& rate_eval) rate_eval.dscreened_rates_dT(k_S32_to_He4_Si28_derived) = drate_dT; } - rate_Fe52_to_p_Mn51_derived(tfactors, rate, drate_dT, pf_cache); - rate_eval.screened_rates(k_Fe52_to_p_Mn51_derived) = rate; + rate_Fe52_to_p_nse_Mn51_derived(tfactors, rate, drate_dT, pf_cache); + rate_eval.screened_rates(k_Fe52_to_p_nse_Mn51_derived) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_Fe52_to_p_Mn51_derived) = drate_dT; + rate_eval.dscreened_rates_dT(k_Fe52_to_p_nse_Mn51_derived) = drate_dT; } rate_Fe52_to_He4_Cr48_derived(tfactors, rate, drate_dT, pf_cache); @@ -7208,10 +7208,10 @@ fill_reaclib_rates(const tf_t& tfactors, T& rate_eval) rate_eval.dscreened_rates_dT(k_Co55_to_He4_Mn51_derived) = drate_dT; } - rate_Ni56_to_p_Co55_derived(tfactors, rate, drate_dT, pf_cache); - rate_eval.screened_rates(k_Ni56_to_p_Co55_derived) = rate; + rate_Ni56_to_p_nse_Co55_derived(tfactors, rate, drate_dT, pf_cache); + rate_eval.screened_rates(k_Ni56_to_p_nse_Co55_derived) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_Ni56_to_p_Co55_derived) = drate_dT; + rate_eval.dscreened_rates_dT(k_Ni56_to_p_nse_Co55_derived) = drate_dT; } rate_Ni56_to_He4_Fe52_derived(tfactors, rate, drate_dT, pf_cache); @@ -7256,22 +7256,22 @@ fill_reaclib_rates(const tf_t& tfactors, T& rate_eval) rate_eval.dscreened_rates_dT(k_He4_Si28_to_p_P31_derived) = drate_dT; } - rate_p_Mn51_to_He4_Cr48_derived(tfactors, rate, drate_dT, pf_cache); - rate_eval.screened_rates(k_p_Mn51_to_He4_Cr48_derived) = rate; + rate_p_nse_Mn51_to_He4_Cr48_derived(tfactors, rate, drate_dT, pf_cache); + rate_eval.screened_rates(k_p_nse_Mn51_to_He4_Cr48_derived) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_p_Mn51_to_He4_Cr48_derived) = drate_dT; + rate_eval.dscreened_rates_dT(k_p_nse_Mn51_to_He4_Cr48_derived) = drate_dT; } - rate_p_Co55_to_He4_Fe52_derived(tfactors, rate, drate_dT, pf_cache); - rate_eval.screened_rates(k_p_Co55_to_He4_Fe52_derived) = rate; + rate_p_nse_Co55_to_He4_Fe52_derived(tfactors, rate, drate_dT, pf_cache); + rate_eval.screened_rates(k_p_nse_Co55_to_He4_Fe52_derived) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_p_Co55_to_He4_Fe52_derived) = drate_dT; + rate_eval.dscreened_rates_dT(k_p_nse_Co55_to_He4_Fe52_derived) = drate_dT; } - rate_Co55_to_p_Fe54_derived(tfactors, rate, drate_dT, pf_cache); - rate_eval.screened_rates(k_Co55_to_p_Fe54_derived) = rate; + rate_Co55_to_p_nse_Fe54_derived(tfactors, rate, drate_dT, pf_cache); + rate_eval.screened_rates(k_Co55_to_p_nse_Fe54_derived) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_Co55_to_p_Fe54_derived) = drate_dT; + rate_eval.dscreened_rates_dT(k_Co55_to_p_nse_Fe54_derived) = drate_dT; } rate_Co56_to_n_Co55_derived(tfactors, rate, drate_dT, pf_cache); @@ -7286,16 +7286,16 @@ fill_reaclib_rates(const tf_t& tfactors, T& rate_eval) rate_eval.dscreened_rates_dT(k_Co57_to_n_Co56_derived) = drate_dT; } - rate_Co57_to_p_Fe56_derived(tfactors, rate, drate_dT, pf_cache); - rate_eval.screened_rates(k_Co57_to_p_Fe56_derived) = rate; + rate_Co57_to_p_nse_Fe56_derived(tfactors, rate, drate_dT, pf_cache); + rate_eval.screened_rates(k_Co57_to_p_nse_Fe56_derived) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_Co57_to_p_Fe56_derived) = drate_dT; + rate_eval.dscreened_rates_dT(k_Co57_to_p_nse_Fe56_derived) = drate_dT; } - rate_Ni58_to_p_Co57_derived(tfactors, rate, drate_dT, pf_cache); - rate_eval.screened_rates(k_Ni58_to_p_Co57_derived) = rate; + rate_Ni58_to_p_nse_Co57_derived(tfactors, rate, drate_dT, pf_cache); + rate_eval.screened_rates(k_Ni58_to_p_nse_Co57_derived) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_Ni58_to_p_Co57_derived) = drate_dT; + rate_eval.dscreened_rates_dT(k_Ni58_to_p_nse_Co57_derived) = drate_dT; } rate_Ni58_to_He4_Fe54_derived(tfactors, rate, drate_dT, pf_cache); @@ -7304,34 +7304,34 @@ fill_reaclib_rates(const tf_t& tfactors, T& rate_eval) rate_eval.dscreened_rates_dT(k_Ni58_to_He4_Fe54_derived) = drate_dT; } - rate_p_Fe54_to_He4_Mn51_derived(tfactors, rate, drate_dT, pf_cache); - rate_eval.screened_rates(k_p_Fe54_to_He4_Mn51_derived) = rate; + rate_p_nse_Fe54_to_He4_Mn51_derived(tfactors, rate, drate_dT, pf_cache); + rate_eval.screened_rates(k_p_nse_Fe54_to_He4_Mn51_derived) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_p_Fe54_to_He4_Mn51_derived) = drate_dT; + rate_eval.dscreened_rates_dT(k_p_nse_Fe54_to_He4_Mn51_derived) = drate_dT; } - rate_He4_Fe54_to_p_Co57_derived(tfactors, rate, drate_dT, pf_cache); - rate_eval.screened_rates(k_He4_Fe54_to_p_Co57_derived) = rate; + rate_He4_Fe54_to_p_nse_Co57_derived(tfactors, rate, drate_dT, pf_cache); + rate_eval.screened_rates(k_He4_Fe54_to_p_nse_Co57_derived) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_He4_Fe54_to_p_Co57_derived) = drate_dT; + rate_eval.dscreened_rates_dT(k_He4_Fe54_to_p_nse_Co57_derived) = drate_dT; } - rate_p_Fe56_to_n_Co56_derived(tfactors, rate, drate_dT, pf_cache); - rate_eval.screened_rates(k_p_Fe56_to_n_Co56_derived) = rate; + rate_p_nse_Fe56_to_n_Co56_derived(tfactors, rate, drate_dT, pf_cache); + rate_eval.screened_rates(k_p_nse_Fe56_to_n_Co56_derived) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_p_Fe56_to_n_Co56_derived) = drate_dT; + rate_eval.dscreened_rates_dT(k_p_nse_Fe56_to_n_Co56_derived) = drate_dT; } - rate_p_Co56_to_n_Ni56_derived(tfactors, rate, drate_dT, pf_cache); - rate_eval.screened_rates(k_p_Co56_to_n_Ni56_derived) = rate; + rate_p_nse_Co56_to_n_Ni56_derived(tfactors, rate, drate_dT, pf_cache); + rate_eval.screened_rates(k_p_nse_Co56_to_n_Ni56_derived) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_p_Co56_to_n_Ni56_derived) = drate_dT; + rate_eval.dscreened_rates_dT(k_p_nse_Co56_to_n_Ni56_derived) = drate_dT; } - rate_p_Ni58_to_He4_Co55_derived(tfactors, rate, drate_dT, pf_cache); - rate_eval.screened_rates(k_p_Ni58_to_He4_Co55_derived) = rate; + rate_p_nse_Ni58_to_He4_Co55_derived(tfactors, rate, drate_dT, pf_cache); + rate_eval.screened_rates(k_p_nse_Ni58_to_He4_Co55_derived) = rate; if constexpr (std::is_same_v) { - rate_eval.dscreened_rates_dT(k_p_Ni58_to_He4_Co55_derived) = drate_dT; + rate_eval.dscreened_rates_dT(k_p_nse_Ni58_to_He4_Co55_derived) = drate_dT; } rate_He4_S32_to_p_Cl35_derived_removed(tfactors, rate, drate_dT, pf_cache);