checks_*: +exit_on_failure option

Author: JamesMcClung

src/include/checks_params.hxx

@@ -7,6 +7,7 @@ struct CheckParams
   double err_threshold = 1e-13;     // maximum acceptable error
   bool print_max_err_always = true; // always print error, even if acceptable
   bool dump_always = false; // always dump compared fields, even if acceptable
+  bool exit_on_failure = false; // exit the program if check fails
 
   bool enabled() { return check_interval > 0; }
 

src/libpsc/psc_checks/checks_impl.hxx

@@ -1,4 +1,3 @@
-
 #pragma once
 
 #include "fields.hxx"
@@ -102,7 +101,10 @@ public:
       MPI_Barrier(grid.comm());
     }
 
-    assert(max_err < err_threshold);
+    if (exit_on_failure && max_err >= err_threshold) {
+      LOG_ERROR("exiting due to failed continuity check\n");
+    }
+
     last_max_err = max_err;
   }
 
@@ -185,7 +187,10 @@ public:
       writer_.end_step();
     }
 
-    assert(max_err < err_threshold);
+    if (exit_on_failure && max_err >= err_threshold) {
+      LOG_ERROR("exiting due to failed gauss check\n");
+    }
+
     last_max_err = max_err;
   }