diff --git a/src/ModelTree.cc b/src/ModelTree.cc
index dd652c7958201f966e8a43156cfdc38fbe5c2387..f6b1b9365a2195603f8860f8460a44421036e1c9 100644
--- a/src/ModelTree.cc
+++ b/src/ModelTree.cc
@@ -346,7 +346,7 @@ ModelTree::evaluateAndReduceJacobian(const eval_context_t &eval_context) const
               cerr << endl;
               exit(EXIT_FAILURE);
             }
-          if (lag == 0)
+          if ((isnan(val) || fabs(val) >= cutoff) && lag == 0)
             contemporaneous_jacobian[{ eq, var }] = val;
         }
     }