Return False in 2x regression tests when analyze() detects issue

Closes #174

Author: felker

src/task_list/task_list.hpp

@@ -155,7 +155,7 @@ class TimeIntegratorTaskList : public TaskList {
 
 
 //----------------------------------------------------------------------------------------
-// 64-bit integers with "1" in different bit positions used to ID  each hydro task.
+// 64-bit integers with "1" in different bit positions used to ID each hydro task.
 
 namespace HydroIntegratorTaskNames {
   const uint64_t NONE=0;

tst/regression/scripts/tests/curvilinear/blast_sph.py

@@ -33,5 +33,6 @@ def analyze():
     # check blast is spherical
     if data[0][3] > 1.0:
         print("Distortion of blast wave in spherical coords too large", data[0][3])
+        return False
 
     return True

tst/regression/scripts/tests/symmetry/hydro_linwave_aligned.py

@@ -98,6 +98,8 @@ def analyze():
         # (Differences in d_max can exceed 2e-15, but other differences are 2e-16 at most)
         atol = 5e-15
         rtol = 1e-8
+
+        # Useful optional diagnostics for determining if differences are meaningful in FP:
         # print(np.allclose(results_1D, results_2D, atol=atol, rtol=rtol))
         # print("numpy tolerance = {}".format(atol + 1e-10*abs(results_2D)))
         # print(np.allclose(results_2D, results_3D, atol=atol, rtol=rtol))
@@ -113,6 +115,6 @@ def analyze():
             print(results_2D)
             print(results_3D)
             print("Exhibit differences that are not close to round-off")
-            # return False
+            return False
 
     return True