update reference values

dG3D/benchmarks/GursonPlasticInstabilities_PathFollowing/Plane_notch.py

@@ -190,4 +190,4 @@ mysolver.archivingIPOnPhysicalGroup("Face", 700, IPField.LOCAL_POROSITY,IPField.
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-1.957777e+03,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-6)
+check.equal(-1.957901e+03,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-6)

dG3D/benchmarks/GursonThomasonMultipleNonlocalVar_PathFollowing/Plane_notch.py

@@ -190,4 +190,4 @@ mysolver.archivingIPOnPhysicalGroup("Face", 700, IPField.LOCAL_POROSITY,IPField.
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-5.458367e+01,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-6)
+check.equal(-5.458283e+01,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-6)

dG3D/benchmarks/GursonThomason_cube_I1J2J3_Transition/cube.py

@@ -71,7 +71,7 @@ nstepArchIP=1		# Number of step between 2 archiving
 nstepArchForce = 1	# Number of step between 2 force archiving
 nstepArchEnergy = nstepArchForce # Number of step between 2 energy computation and archiving
 
-MaxIter = 12		# Maximum number of iterations
+MaxIter = 16		# Maximum number of iterations
 StepIncrease = 2	# Number of successfull timestep before reincreasing it
 StepReducFactor = 2 	# Timestep reduction factor
 NumberReduction = 100	# Maximum number of timespep reduction: max reduction = pow(StepReducFactor,this)

dG3D/benchmarks/GursonThomason_cube_Transition/cube.py

@@ -68,7 +68,7 @@ nstepArchIP=1		# Number of step between 2 archiving
 nstepArchForce = 1	# Number of step between 2 force archiving
 nstepArchEnergy = nstepArchForce # Number of step between 2 energy computation and archiving
 
-MaxIter = 12		# Maximum number of iterations
+MaxIter = 16		# Maximum number of iterations
 StepIncrease = 2	# Number of successfull timestep before reincreasing it
 StepReducFactor = 2 	# Timestep reduction factor
 NumberReduction = 100	# Maximum number of timespep reduction: max reduction = pow(StepReducFactor,this)

dG3D/benchmarks/GursonThomason_cube_withLodeEnhanced/cube.py

@@ -197,7 +197,7 @@ check = TestCheck()
 import csv
 data = csv.reader(open('IPVolume11val_LOCAL_POROSITYMax.csv'), delimiter=';')
 porosity = list(data)
-check.equal(4.420505e-01,float(porosity[-1][1]),1e-6)
+check.equal(4.415553e-01,float(porosity[-1][1]),1e-6)
 
 
 

dG3D/benchmarks/Inno_EMT_test/InnoElecMagTherm.py

@@ -296,7 +296,7 @@ check.equal(0.000000e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSm
 #check.equal(3.399153e-04,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 4),1.e-5)
 check.equal(-3.484000e-03,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 4),1.e-2)
 check.equal(3.144224e-03,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 4),1.e-5)
-check.equal(-6.802441e-04,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 5),1.e-5)
+check.equal(-6.785826e-04,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 5),1.e-5)
 check.equal(3.720132e-01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 5),1.e-5)
 check.equal(5.949466e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 5),1.e-5)
 check.equal(0.0,mysolver.getArchivedNodalValue(SmpRefPoint,0,mysolver.displacement),1.e-5)

dG3D/benchmarks/defoDefoContactSlaveEdge2DQuad/squares.py

@@ -59,7 +59,7 @@ mysolver.displacementBC("Face",nfield1,2,0.)
 mysolver.displacementBC("Face",nfield2,2,0.)
 #compression along z
 mysolver.displacementBC("Edge",12,0,0.)
-mysolver.displacementBC("Edge",12,2,0.)
+mysolver.displacementBC("Edge",12,1,0.)
 d1=-0.00025
 cyclicFunction1=cycleFunctionTime(0.,0.,ftime/4., d1/2., ftime/2., 0., 3.*ftime/4., d1/2., ftime, d1);
 #cyclicFunction1=cycleFunctionTime(0.,0.,ftime, d1);
@@ -93,6 +93,6 @@ mysolver.archivingForceOnPhysicalGroup("Edge", 156, 1)
 mysolver.solve()
 
 check = TestCheck()
-check.equal(2.451079e+07,mysolver.getArchivedForceOnPhysicalGroup("Edge", 12, 1),1.e-6)
+check.equal(2.450488e+07,mysolver.getArchivedForceOnPhysicalGroup("Edge", 12, 1),1.e-6)
 
 

dG3D/benchmarks/defoDefoContactSlaveEdge2DQuadPert/squares.py

@@ -95,6 +95,6 @@ mysolver.archivingForceOnPhysicalGroup("Edge", 156, 1)
 mysolver.solve()
 
 check = TestCheck()
-check.equal(2.451079e+07,mysolver.getArchivedForceOnPhysicalGroup("Edge", 12, 1),1.e-6)
+check.equal(2.450488e+07,mysolver.getArchivedForceOnPhysicalGroup("Edge", 12, 1),1.e-6)
 
 

dG3D/benchmarks/linearEMT/LinearElecMagTherm.py

@@ -290,13 +290,12 @@ mysolver.archivingNodeDisplacement(SmpRefPoint,4,1);
 #mysolver.archivingNodeIP(SmpRef, IPField.ELECTRICALFLUX_X,IPField.MEAN_VALUE,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(3.662529e-01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 5),1.e-5)
+check.equal(3.691623e-01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 5),1.e-5)
 check.equal(2.895048e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 5),1.e-5)
 check.equal(-1.265934e+02,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 5),1.e-5)
 check.equal(0.0,mysolver.getArchivedNodalValue(SmpRefPoint,0,mysolver.displacement),1.e-5)

dG3D/benchmarks/linearEMT_Pert/LinearElecMagTherm_Pert.py

@@ -293,7 +293,7 @@ check = TestCheck()
 check.equal(0.000000e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 3),1.e-5)
 check.equal(0.000000e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 4),1.e-5)
 check.equal(-3.322070e-05,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 4),1.e-5)
-check.equal(1.017284e-07,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 4),1.e-5)
+check.equal(1.017297e-07,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 4),1.e-5)
 check.equal(2.235121e+01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 5),1.e-5)
 check.equal(-6.599484e+01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 5),1.e-5)
 check.equal(-1.150381e+02,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 5),1.e-5)

dG3D/benchmarks/linearElectroMagnetic/LinearElecMag.py

@@ -289,7 +289,7 @@ check.equal(-6.996190e+11,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfS
 check.equal(0.000000e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 4),1.e-5)
 check.equal(-1.789120e-08,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 4),1.e-2)
 check.equal(1.786247e-08,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 4),1.e-2)
-check.equal(2.291965e-01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 5),1.e-5)
+check.equal(2.286854e-01 ,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 5),1.e-5)
 check.equal(-6.912008e-01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 5),1.e-5)
 check.equal(-3.522141e+02,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 5),1.e-5)
 check.equal(0.0,mysolver.getArchivedNodalValue(SmpRefPoint,0,mysolver.displacement),1.e-5)
@@ -297,4 +297,3 @@ check.equal(0.0,mysolver.getArchivedNodalValue(SmpRefPoint,1,mysolver.displaceme
 check.equal(0.0,mysolver.getArchivedNodalValue(SmpRefPoint,2,mysolver.displacement),1.e-5)
 check.equal(2.980000e+02,mysolver.getArchivedNodalValue(SmpRefPoint,3,mysolver.displacement),1.e-5)
 check.equal(-6.173236e-04,mysolver.getArchivedNodalValue(SmpRefPoint,4,mysolver.displacement),1.e-5)
-

dG3D/benchmarks/multiscale_2order_shearlayer_pert/shearlayerPert.py

@@ -149,5 +149,5 @@ mysolver.solve()
 
 #test check
 check = TestCheck()
-check.equal(4.649837e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge",2,1),1.e-5)
+check.equal(4.649913e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge",2,1),1.e-5)
 

dG3D/benchmarks/nanoIndentation/ni.py

@@ -166,7 +166,7 @@ mysolver.archivingRigidContactForce(22, 2, 1)
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-2.811155e-04,mysolver.getArchivedNodalValue(1001,2,mysolver.displacement),1.e-4)
+check.equal(-2.816262e-04,mysolver.getArchivedNodalValue(1001,2,mysolver.displacement),1.e-4)
 
 
 

dG3D/benchmarks/shiftMultiSystem/run.py

@@ -13,7 +13,7 @@ else:
 
 data1 = csv.reader(open('force102comp0_part2.csv'), delimiter=';')
 force = list(data1)
-checkEqual(1235.178136860775,float(force[-1][1]),1e-5)
+checkEqual(1.236021e+03,float(force[-1][1]),1e-5)
 
 if sys.version_info[0] < 3:
   os.system('mpiexec -np 4 python fullPlateHoleRes.py')
@@ -22,5 +22,5 @@ else:
 
 data21 = csv.reader(open('force102comp0_part2.csv'), delimiter=';')
 force2 = list(data21)
-checkEqual(1.235168e+03,float(force2[-1][1]),1e-5)
+checkEqual(1.236021e+03,float(force2[-1][1]),1e-5)
 

dG3D/benchmarks/torchUniaxial/macro_cpp.py

@@ -55,7 +55,7 @@ macrogeofile="squares.geo"
 # solver
 sol = 2  # Gmm=0 (default) Taucs=1 PETsc=2
 soltype =1 # StaticLinear=0 (default) StaticNonLinear=1
-nstep = 50  # number of step (used only if soltype=1)
+nstep = 50000  # number of step (used only if soltype=1)
 ftime =1.   # Final time (used only if soltype=1)
 tol=1.e-4  # relative tolerance for NR scheme (used only if soltype=1)
 nstepArch=1 # Number of step between 2 archiving (used only if soltype=1)
@@ -82,7 +82,7 @@ mysolver.addMaterialLaw(macromat1)
 #mysolver.addMaterialLaw(pertmacromat1)
 mysolver.Scheme(soltype)
 mysolver.Solver(sol)
-mysolver.snlData(nstep,ftime,tol)
+mysolver.snlData(nstep,ftime,tol,tol/100.)
 
 # boundary condition
 mysolver.displacementBC("Face",nfield1,2,0.0)