Added tetrahedron 3D python test

diff --git a/CoreFlows/examples/Python/SinglePhase/SinglePhase_3DSphericalExplosion_TETRAHEDRON.py b/CoreFlows/examples/Python/SinglePhase/SinglePhase_3DSphericalExplosion_TETRAHEDRON.py
new file mode 100755 (executable)
index 0000000..b33e039
--- /dev/null
+++ b/CoreFlows/examples/Python/SinglePhase/SinglePhase_3DSphericalExplosion_TETRAHEDRON.py
@@ -0,0 +1,111 @@
+#!/usr/bin/env python3
+# -*-coding:utf-8 -*
+
+import CoreFlows as cf
+import cdmath
+
+
+def SinglePhase_3DSphericalExplosion_TETRAHEDRON():
+
+       inputfile="../resources/meshTetrahedron10.med";
+       my_mesh=cdmath.Mesh(inputfile);
+       spaceDim=3
+       
+       # Initial field data
+       nVar=2+spaceDim;
+       radius=0.1;
+       Center=cdmath.Vector(spaceDim);#default value is (0,0,0)
+       Vout = cdmath.Vector(nVar)
+       Vin  = cdmath.Vector(nVar)
+       
+       Pmin=100e5
+       Pmax=155e5
+       InitialTemperature = 563
+
+       Vin[0]=Pmax
+       Vin[1]=0;
+       Vin[2]=0;
+       Vin[3]=0
+       Vin[4]=InitialTemperature
+       Vout[0]=Pmin;
+       Vout[1]=0;
+       Vout[2]=0;
+       Vout[3]=0
+       Vout[4]=InitialTemperature
+
+       myProblem = cf.SinglePhase(cf.Liquid,cf.around155bars600K,spaceDim);
+
+       # Initial field creation
+       print ("Setting mesh and initial data" ) ;
+       myProblem.setInitialFieldSphericalStepFunction( my_mesh, Vin, Vout, radius, Center);
+
+       # set the boundary conditions
+       wallVelocityX=0;
+       wallVelocityY=0;
+       wallVelocityZ=0;
+       wallTemperature=InitialTemperature;
+
+       myProblem.setWallBoundaryCondition("BoundaryFaces", wallTemperature, wallVelocityX, wallVelocityY, wallVelocityZ);
+
+       # set the numerical method
+       myProblem.setNumericalScheme(cf.upwind, cf.Explicit);
+    
+       # name file save
+       fileName = "3DSphericalExplosion_TETRAHEDRON";
+
+       # parameters calculation
+       MaxNbOfTimeStep = 3 ;
+       freqSave = 1;
+       cfl = 0.3;
+       maxTime = 5;
+       precision = 1e-6;
+
+       myProblem.setCFL(cfl);
+       myProblem.setPrecision(precision);
+       myProblem.setMaxNbOfTimeStep(MaxNbOfTimeStep);
+       myProblem.setTimeMax(maxTime);
+       myProblem.setFreqSave(freqSave);
+       myProblem.setFileName(fileName);
+       myProblem.setNewtonSolver(precision,20);
+       myProblem.saveConservativeField(True);
+
+
+       # evolution
+       myProblem.initialize();
+
+       masseInitiale = abs( myProblem.getDensityField().integral()[0] )
+       initialTotalEnergy=abs(myProblem.getTotalEnergyField().integral()[0] )
+
+       ok = myProblem.run();
+       if (ok):
+               print( "Simulation python " + fileName + " is successful !" );
+               pass
+       else:
+               print( "Simulation python " + fileName + "  failed ! " );
+               pass
+
+       print( "------------ End of calculation !!! -----------" );
+
+       # Control of the results
+       PressureField=myProblem.getPressureField()
+       print( "La pression est bornée par le maximum de pression initiale, à 10% près :\n (PressureField.max() - InitialPressure.max())/InitialPressure.max()= ", (PressureField.max() - Pmax)/Pmax )
+       assert PressureField.max() < Pmax
+
+       temperatureField=myProblem.getTemperatureField()
+       print("La température est constante à 0.1 % près : erreur relative = ", max(abs(temperatureField.max() - InitialTemperature),abs(temperatureField.min() - InitialTemperature))/InitialTemperature )
+       assert abs(temperatureField.max() - InitialTemperature)/InitialTemperature < 0.001
+       assert abs(temperatureField.min() - InitialTemperature)/InitialTemperature < 0.001
+
+       densityField=myProblem.getDensityField()
+       print("La masse totale est conservée au zero machine près : erreur relative = ", abs( (abs(densityField.integral()[0]) - masseInitiale)/masseInitiale ) )
+       assert abs( (abs(densityField.integral()[0]) - masseInitiale)/masseInitiale ) < 1e-14
+       
+       totalEnergyField=myProblem.getTotalEnergyField()
+       print("L'énergie totale est conservée au zero machine près : erreur relative = ", abs( (abs(totalEnergyField.integral()[0]) - initialTotalEnergy)/initialTotalEnergy ) )
+       assert abs( (abs(totalEnergyField.integral()[0]) - initialTotalEnergy)/initialTotalEnergy ) < 1e-13
+
+       myProblem.terminate();
+       return ok
+
+if __name__ == """__main__""":
+    SinglePhase_3DSphericalExplosion_TETRAHEDRON()