From e1cf89185182d9273abfeb1c52d7eb692a80a89b Mon Sep 17 00:00:00 2001
From: Michael Zingale <michael.zingale@stonybrook.edu>
Date: Sun, 17 Dec 2023 15:24:33 -0500
Subject: [PATCH] this works with ROCm > 5.4.0

---
 integration/integrator_rhs_simplified_sdc.H | 4 ++++
 integration/integrator_rhs_strang.H         | 6 ++++--
 2 files changed, 8 insertions(+), 2 deletions(-)

diff --git a/integration/integrator_rhs_simplified_sdc.H b/integration/integrator_rhs_simplified_sdc.H
index e44c8f8e4c..c7fba04aa5 100644
--- 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
index 2018c76b10..d67a403056 100644
--- 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.