introduce blast wave setup

Author: glesur

test/HD/SedovBlastWave/definitions.hpp

@@ -0,0 +1,4 @@
+#define     COMPONENTS      3
+#define     DIMENSIONS      3
+
+#define     GEOMETRY        CARTESIAN

test/HD/SedovBlastWave/idefix.ini

@@ -0,0 +1,25 @@
+[Grid]
+X1-grid    1  -0.5  256  u  0.5
+X2-grid    1  -0.5  256  u  0.5
+X3-grid    1  -0.5  256  u  0.5
+
+[TimeIntegrator]
+CFL         0.9
+tstop       5.0
+first_dt    1.e-6
+nstages     2
+
+[Hydro]
+solver    hllc
+gamma     1.666666666666666666
+
+[Boundary]
+X1-beg    periodic
+X1-end    periodic
+X2-beg    periodic
+X2-end    periodic
+X3-beg    periodic
+X3-end    periodic
+
+[Output]
+vtk    0.01

test/HD/SedovBlastWave/setup.cpp

@@ -0,0 +1,55 @@
+#include "idefix.hpp"
+#include "setup.hpp"
+
+
+// Default constructor
+
+
+// Initialisation routine. Can be used to allocate
+// Arrays or variables which are used later on
+Setup::Setup(Input &input, Grid &grid, DataBlock &data, Output &output) {
+
+}
+
+// This routine initialize the flow
+// Note that data is on the device.
+// One can therefore define locally
+// a datahost and sync it, if needed
+void Setup::InitFlow(DataBlock &data) {
+    // Create a host copy
+    DataBlockHost d(data);
+
+
+    for(int k = 0; k < d.np_tot[KDIR] ; k++) {
+        for(int j = 0; j < d.np_tot[JDIR] ; j++) {
+            for(int i = 0; i < d.np_tot[IDIR] ; i++) {
+              real x = d.x[IDIR](i);
+              real y = d.x[JDIR](j);
+              real z = d.x[KDIR](k);
+
+                // Sedov Blast Wave Following Stone+2018, 3.4.4
+                d.Vc(RHO,k,j,i) = 1.0;
+                d.Vc(VX1,k,j,i) = 0.0;
+                d.Vc(VX2,k,j,i) = 0.0;
+                d.Vc(VX3,k,j,i) = 0.0;
+
+                d.Vc(PRS,k,j,i) = 0.01;
+
+                real r=sqrt(x*x+y*y+z*z);
+                if(r<0.01) {
+                  d.Vc(PRS,k,j,i) = 1.6e5;
+                }
+
+
+            }
+        }
+    }
+
+    // Send it all, if needed
+    d.SyncToDevice();
+}
+
+// Analyse data to produce an output
+void MakeAnalysis(DataBlock & data) {
+
+}