sanity check at every time step.

LLNL · May 9, 2024 · f846f7e · f846f7e
1 parent 80b5334
commit f846f7e
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/sketches/usns.cpp b/sketches/usns.cpp
@@ -279,6 +279,12 @@ class NavierSolver
 
       time += dt;
    }
+
+   void SanityCheck()
+   {
+      assert(!isnan(u->Min()) && !isnan(u->Max()));
+      assert(!isnan(p->Min()) && !isnan(p->Max()));
+   }
 };
 
 int main(int argc, char *argv[])
@@ -404,6 +410,7 @@ int main(int argc, char *argv[])
       }
 
       navier.Step(time, step);
+      navier.SanityCheck();
    }
 
    // 17. Free the used memory.