switch how we deal with the zeroing.

Now we make it the job of the network
AMReX-Astro · Jul 31, 2024 · 14e234b · 14e234b
1 parent 391ffcb
commit 14e234b
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 0 deletions.
diff --git a/integration/integrator_rhs_sdc.H b/integration/integrator_rhs_sdc.H
@@ -105,6 +105,7 @@ void jac (const amrex::Real time, BurnT& state, T& int_state, MatrixType& pd)
 
 
     if (state.T <= EOSData::mintemp || state.T >= integrator_rp::MAX_TEMP) {
+        pd.zero();
         return;
     }
 

diff --git a/integration/integrator_rhs_strang.H b/integration/integrator_rhs_strang.H
@@ -119,6 +119,7 @@ void jac ([[maybe_unused]] const amrex::Real time, BurnT& state, T& int_state, M
     // bounds. Otherwise set the Jacobian to zero and return.
 
     if (state.T <= EOSData::mintemp || state.T >= integrator_rp::MAX_TEMP) {
+        pd.zero();
         return;
     }
-Original file line number
+Diff line change
@@ Expand Up @@
         if (state.T <= EOSData::mintemp || state.T >= integrator_rp::MAX_TEMP) {
+            pd.zero();
             return;
         }
@@ Expand Down @@