this works with ROCm > 5.4.0

AMReX-Astro · Dec 17, 2023 · e1cf891 · e1cf891
1 parent 42fcc6e
commit e1cf891
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 2 deletions.
diff --git a/integration/integrator_rhs_simplified_sdc.H b/integration/integrator_rhs_simplified_sdc.H
@@ -42,6 +42,7 @@ void rhs(const Real time, BurnT& state, T& int_state, RArray1D& ydot,
 
     // make sure that the temperature is valid
 
+#if 0
     if (state.T <= EOSData::mintemp || state.T >= MAX_TEMP) {
 
         for (int n = 1; n <= INT_NEQS; ++n) {
@@ -51,6 +52,7 @@ void rhs(const Real time, BurnT& state, T& int_state, RArray1D& ydot,
         return;
 
     }
+#endif
 
     // call the specific network to get the RHS
 
@@ -104,6 +106,7 @@ void jac (const Real time, BurnT& state, T& int_state, MatrixType& pd)
     int_to_burn(time, int_state, state);
 
 
+#if 0
     if (state.T <= EOSData::mintemp || state.T >= MAX_TEMP) {
 
         for (int j = 1; j <= INT_NEQS; ++j) {
@@ -115,6 +118,7 @@ void jac (const Real time, BurnT& state, T& int_state, MatrixType& pd)
         return;
 
     }
+#endif
 
     // Call the specific network routine to get the Jacobian.
 

diff --git a/integration/integrator_rhs_strang.H b/integration/integrator_rhs_strang.H
@@ -50,7 +50,7 @@ void rhs (const Real time, BurnT& state, T& int_state, RArray1D& ydot, [[maybe_u
 
     // Only do the burn if the incoming temperature is within the temperature
     // bounds. Otherwise set the RHS to zero and return.
-
+#if 0
     if (state.T <= EOSData::mintemp || state.T >= MAX_TEMP) {
 
         for (int n = 1; n <= INT_NEQS; ++n) {
@@ -60,6 +60,7 @@ void rhs (const Real time, BurnT& state, T& int_state, RArray1D& ydot, [[maybe_u
         return;
 
     }
+#endif
 
     state.time = time;
 
@@ -119,7 +120,7 @@ void jac ([[maybe_unused]] const Real time, BurnT& state, T& int_state, MatrixTy
 
     // Only do the burn if the incoming temperature is within the temperature
     // bounds. Otherwise set the Jacobian to zero and return.
-
+#if 0
     if (state.T <= EOSData::mintemp || state.T >= MAX_TEMP) {
 
         for (int j = 1; j <= INT_NEQS; ++j) {
@@ -131,6 +132,7 @@ void jac ([[maybe_unused]] const Real time, BurnT& state, T& int_state, MatrixTy
         return;
 
     }
+#endif
 
     // Call the specific network routine to get the Jacobian.