more syncing

AMReX-Astro · Feb 6, 2025 · 60bde11 · 60bde11
1 parent 27b71fc
commit 60bde11
Show file tree

Hide file tree

Showing 3 changed files with 6 additions and 4 deletions.
diff --git a/unit_test/burn_cell_sdc/burn_cell.H b/unit_test/burn_cell_sdc/burn_cell.H
@@ -20,7 +20,7 @@ void burn_cell_c()
     // Make sure user set all the mass fractions to values in the interval [0, 1]
     for (int n = 1; n <= NumSpec; ++n) {
 
-        massfractions[n-1] = get_xn(n);
+        massfractions[n-1] = get_xn(n, unit_test_rp::uniform_xn);
 
         if (massfractions[n-1] < 0 || massfractions[n-1] > 1) {
             amrex::Error("mass fraction for " + short_spec_names_cxx[n-1] + " not initialized in the interval [0,1]!");

diff --git a/unit_test/eos_cell/eos_cell.H b/unit_test/eos_cell/eos_cell.H
@@ -24,7 +24,7 @@ void eos_cell_c()
     // Make sure user set all the mass fractions to values in the interval [0, 1]
     for (int n = 1; n <= NumSpec; ++n) {
 
-        massfractions[n-1] = get_xn(n);
+        massfractions[n-1] = get_xn(n, unit_test_rp::uniform_xn);
 
         if (massfractions[n-1] < 0 || massfractions[n-1] > 1) {
             amrex::Error("mass fraction for " + short_spec_names_cxx[n-1] + " not initialized in the interval [0,1]!");

diff --git a/unit_test/jac_cell/jac_cell.H b/unit_test/jac_cell/jac_cell.H
@@ -22,7 +22,7 @@ void jac_cell_c()
     // Make sure user set all the mass fractions to values in the interval [0, 1]
     for (int n = 1; n <= NumSpec; ++n) {
 
-        massfractions[n-1] = get_xn(n);
+        massfractions[n-1] = get_xn(n, unit_test_rp::uniform_xn);
 
         if (massfractions[n-1] < 0 || massfractions[n-1] > 1) {
             amrex::Error("mass fraction for " + short_spec_names_cxx[n-1] + " not initialized in the interval [0,1]!");
@@ -38,7 +38,9 @@ void jac_cell_c()
         burn_state.xn[n] = massfractions[n];
     }
 
-    normalize_abundances_burn(burn_state);
+    if (! unit_test_rp::skip_initial_normalization) {
+        normalize_abundances_burn(burn_state);
+    }
 
     // get the energy