Reference values in tests, perturbation tests fixed

dG3D/benchmarks/EMT_Cylinder_Unit_Test/cyl_LinearEMT.py

@@ -87,8 +87,8 @@ magpot0_y_vac, magpot0_z_vac, Irms, freq, nTurnsCoil, coilLx, coilLy, coilLz, co
 # creation of ElasticField
 SMPfield = 10 # number of the field (physical number of SMP region)
 Vacfield = 20 # number of the field (physical number of Vacuum)
-SurfSMPStart = 6 # delta voltage surface
-SurfSMPEnd = 5 # ground voltage surface
+SurfSMPStart = 3 # delta voltage surface
+SurfSMPEnd = 4 # ground voltage surface
 
 space1 = 0 # function space (Lagrange=0)
 
@@ -146,7 +146,7 @@ mysolver.displacementBC("Volume",Vacfield,2,0.0)
 cyclicFunctionTemp1=cycleFunctionTime(0., tinitial,ftime,tinitial);
 mysolver.initialBC("Volume","Position",SMPfield,3,tinitial)
 mysolver.initialBC("Volume","Position",Vacfield,3,tinitial)
-mysolver.displacementBC("Face",7,3,cyclicFunctionTemp1)
+mysolver.displacementBC("Face",8,3,cyclicFunctionTemp1)
 #mysolver.displacementBC("Face",5,3,cyclicFunctionTemp1)
 #mysolver.displacementBC("Face",6,3,cyclicFunctionTemp1)
 #mysolver.displacementBC("Face",SurfSMPStart,3,cyclicFunctionTemp1)
@@ -165,12 +165,12 @@ mysolver.displacementBC("Face",SurfSMPEnd,4,cyclicFunctionvolt2)
 #magentic
 mysolver.curlDisplacementBC("Face",SurfSMPStart,5,0.0) # comp may also be 5
 mysolver.curlDisplacementBC("Face",SurfSMPEnd,5,0.0) # comp may also be 5
-mysolver.curlDisplacementBC("Face",3,5,0.0) # comp may also be 5
-mysolver.curlDisplacementBC("Face",4,5,0.0) # comp may also be 5
-mysolver.curlDisplacementBC("Face",7,5,0.0) # comp may also be 5
+mysolver.curlDisplacementBC("Face",5,5,0.0) # comp may also be 5
+mysolver.curlDisplacementBC("Face",6,5,0.0) # comp may also be 5
+mysolver.curlDisplacementBC("Face",8,5,0.0) # comp may also be 5
 #Gauging the edges using tree-cotree method
 PhysicalCurves = "" 		# input required as a string of comma separated ints
-PhysicalSurfaces = "3,4,5,6,7" # input required as a string of comma separated ints
+PhysicalSurfaces = "3,4,5,6,8" # input required as a string of comma separated ints
 PhysicalVolumes = "10,20" 	# input required as a string of comma separated ints
 OutputPhysical = 55 		# input required as a int
 mysolver.createTreeForBC(PhysicalCurves,PhysicalSurfaces,PhysicalVolumes,OutputPhysical)
@@ -206,9 +206,9 @@ mysolver.archivingForceOnPhysicalGroup("Face", 7, 5, 1);
 mysolver.solve()
 
 check = TestCheck()
-check.equal(2.200194e+09,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSMPStart, 3),1.e-5)
-check.equal(-1.232263e+09,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSMPStart, 4),1.e-5)
-check.equal(1.232263e+09,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSMPEnd, 4),1.e-5)
-check.equal(4838094703.964682,mysolver.getArchivedForceOnPhysicalGroup("Face", 7, 5),1.e-5)
+check.equal(0.000000e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSMPStart, 3),1.e-5)
+check.equal(-1.273695e+09,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSMPStart, 4),1.e-5)
+check.equal(1.273695e+09,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSMPEnd, 4),1.e-5)
+check.equal(1.387338e+08,mysolver.getArchivedForceOnPhysicalGroup("Face", 7, 5),1.e-5)
 
 

dG3D/benchmarks/linearEMT/LinearElecMagTherm.py

@@ -302,8 +302,8 @@ mysolver.archivingNodeDisplacement(SmpRefPoint,4,1);
 mysolver.solve()
 check = TestCheck()
 check.equal(0.000000e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 3),1.e-5)
-check.equal(-9.427355e-10,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 4),1.e1)
-check.equal(9.484542e-05,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 4),1.e1)
-check.equal(-7.110362e-08,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 4),1.e-1)
-check.equal(-4.387605e-01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 5),1.e-5)
+check.equal(-3.599847e-08,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 4),1.e1)
+check.equal(-3.318414e-05,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 4),1.e1)
+check.equal(1.0116289e-07,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 4),1.e-1)
+check.equal(0.1130046,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 5),1.e-5)
 

dG3D/benchmarks/linearEMT/smp.geo

@@ -48,9 +48,9 @@ BooleanFragments{ Volume{vAir(), vCoreE(), vCoil()}; Delete; }{} // This needs t
 
 // Adapting mesh size
 // characteristic lengths
-lc0 = Pi*Rint/3; // /10; // air
-lc1  = wcoil/2; /// /3; // coil 
-lc2  = R/3; ////6; // smp
+lc0 = Pi*Rint/5; // air
+lc1  = wcoil/3; // coil
+lc2  = R/2; // smp
 
 Characteristic Length { PointsOf{ Volume{vAir(0)}; } } = lc0;
 Characteristic Length { PointsOf{ Volume{vCoil()}; } } = lc1;

dG3D/benchmarks/linearEMT_Pert/LinearElecMagTherm_Pert.py

@@ -304,12 +304,12 @@ mysolver.solve()
 
 check = TestCheck()
 check.equal(0.000000e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 3),1.e-5)
-check.equal(-3.650220e-08,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 4),1.e-2)
+check.equal(-3.588918e-08,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 4),1.e-2)
 check.equal(-3.318414e-05,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 4),1.e-5)
 check.equal(1.011629e-07,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 4),1.e-5)
-check.equal(2.085980e+01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 5),1.e-5)
-check.equal(-8.683505e+01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 5),1.e-5)
-check.equal(-1.150760e+02,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 5),1.e-5)
+check.equal(2.429194e+01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 5),1.e-5)
+check.equal(-8.91641e+01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 5),1.e-5)
+check.equal(-1.15082e+02,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 5),1.e-5)
 check.equal(0.0,mysolver.getArchivedNodalValue(SmpRefPoint,0,mysolver.displacement),1.e-3)
 check.equal(0.0,mysolver.getArchivedNodalValue(SmpRefPoint,1,mysolver.displacement),1.e-3)
 check.equal(0.0,mysolver.getArchivedNodalValue(SmpRefPoint,2,mysolver.displacement),1.e-3)