make sure the blast wave works

Sanity check on the linearwavetest

Author: glesur

…SedovBlastWave/definitions_spherical.hpp …SedovBlastWave/definitions-spherical.hpptest/HD/SedovBlastWave/definitions_spherical.hpp renamed to test/HD/SedovBlastWave/definitions-spherical.hpp test/HD/SedovBlastWave/definitions_spherical.hpp renamed to test/HD/SedovBlastWave/definitions-spherical.hpp

@@ -1,7 +1,7 @@
 [Grid]
-X1-grid    1  -0.5  128  u  0.5
-X2-grid    1  -0.5  128  u  0.5
-X3-grid    1  -0.5  128  u  0.5
+X1-grid    1  -0.5  64  u  0.5
+X2-grid    1  -0.5  64  u  0.5
+X3-grid    1  -0.5  64  u  0.5
 
 [TimeIntegrator]
 CFL         0.9

test/HD/SedovBlastWave/idefix.ini

@@ -9,6 +9,17 @@
 from pytools.dump_io import readDump
 import numpy as np
 import matplotlib.pyplot as plt
+import argparse
+
+parser = argparse.ArgumentParser()
+parser.add_argument("-noplot",
+                    default=False,
+                    help="disable plotting",
+                    action="store_true")
+
+
+args, unknown=parser.parse_known_args()
+
 
 def getError(rep):
     try:
@@ -29,9 +40,10 @@ def getError(rep):
         Q2=Q2-np.mean(Q2)
         print(key+" error=%e"%((np.mean(np.abs(Q1-Q2)))*1e6))
         err=err + (np.mean(np.abs(Q1-Q2)))**2
-        plt.figure()
-        plt.contourf(V.data[key][:,:,16]-U.data[key][:,:,16])
-        plt.title(key)
+        if(not args.noplot):
+          plt.figure()
+          plt.contourf(V.data[key][:,:,16]-U.data[key][:,:,16])
+          plt.title(key)
 
 
     err = err/len(keylist)
@@ -41,4 +53,13 @@ def getError(rep):
 
 error = getError("..")
 print("Error=%e"%(error*1e6))
-plt.show()
+
+if(not args.noplot):
+  plt.show()
+
+if error*1e6<3e-2:
+  print("SUCCESS")
+  sys.exit(0)
+else:
+  print("Failed")
+  sys.exit(1)

test/MHD/LinearWaveTest/python/testidefix.py

@@ -1,7 +1,7 @@
 #!/bin/bash
 
 rep_2D_mpi_list="HD/MachReflection HD/ViscousFlowPastCylinder HD/FargoPlanet MHD/OrszagTang"
-rep_3D_mpi_list="HD/SedovBlastWave MHD/AmbipolarCshock3D MHD/AxisFluxTube MHD/LinearWaveTest MHD/FargoMHDSpherical MHD/OrszagTang3D"
+rep_3D_mpi_list="MHD/AmbipolarCshock3D MHD/AxisFluxTube MHD/LinearWaveTest MHD/FargoMHDSpherical MHD/OrszagTang3D"
 
 # refer to the parent dir of this file, wherever this is called from
 # a python equivalent is e.g.
@@ -102,5 +102,62 @@ for rep in $rep_3D_mpi_list; do
     cd $TEST_DIR
 done
 
+# Blast wave test
+rep="HD/SedovBlastWave"
+
+## Cartesian blast
+cp -R $TEST_DIR/$rep $TMP_DIR/$rep
+cd $TMP_DIR/$rep
+echo "***********************************************"
+echo "Configuring  $rep"
+echo "Using $TMP_DIR/$rep as working directory"
+echo "***********************************************"
+rm -f CMakeCache.txt
+cmake $IDEFIX_DIR -DIdefix_MPI=ON $options
+echo "***********************************************"
+echo "Making  $rep"
+echo "***********************************************"
+make clean; make -j 10
+
+ini="idefix.ini"
+echo "***********************************************"
+echo "Running  $rep with $ini"
+echo "***********************************************"
+mpirun -np 8 ./idefix -i $ini -dec 2 2 2 -nolog
+
+cd python
+echo "***********************************************"
+echo "Testing  $rep with $ini"
+echo "***********************************************"
+python3 testidefix.py -noplot -i ../$ini
+cd ..
+
+## Spherical blast
+echo "***********************************************"
+echo "Configuring  $rep in spherical geometry"
+echo "Using $TMP_DIR/$rep as working directory"
+echo "***********************************************"
+rm -f CMakeCache.txt
+cmake $IDEFIX_DIR -DIdefix_MPI=ON -DIdefix_DEFS=definitions-spherical.hpp $options
+echo "***********************************************"
+echo "Making  $rep"
+echo "***********************************************"
+make clean; make -j 10
+
+ini="idefix-spherical.ini"
+echo "***********************************************"
+echo "Running  $rep with $ini"
+echo "***********************************************"
+mpirun -np 8 ./idefix -i $ini -dec 2 2 2 -nolog
+
+cd python
+echo "***********************************************"
+echo "Testing  $rep with $ini"
+echo "***********************************************"
+python3 testidefix.py -noplot -i ../$ini
+
+cd $TEST_DIR
+
+
 echo "Cleaning temporary directory $TMP_DIR"
 rm -rf $TMP_DIR