Skip to content
Snippets Groups Projects
Commit 80b2c2b7 authored by Ludovic Noels's avatar Ludovic Noels
Browse files

reference for cm3007

parent 4b2c10fd
No related branches found
No related tags found
No related merge requests found
...@@ -296,7 +296,7 @@ check.equal(0.000000e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSm ...@@ -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.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.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(3.144224e-03,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 4),1.e-5)
check.equal(-6.785826e-04,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 5),1.e-5) check.equal(-6.775097e-04,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 5),1.e-5)
check.equal(3.720132e-01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 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(5.949466e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 5),1.e-5)
check.equal(0.0,mysolver.getArchivedNodalValue(SmpRefPoint,0,mysolver.displacement),1.e-5) check.equal(0.0,mysolver.getArchivedNodalValue(SmpRefPoint,0,mysolver.displacement),1.e-5)
......
...@@ -40,7 +40,7 @@ meshfile="cube.msh" # name of mesh file ...@@ -40,7 +40,7 @@ meshfile="cube.msh" # name of mesh file
# solver # solver
sol = 2 # Gmm=0 (default) Taucs=1 PETsc=2 sol = 2 # Gmm=0 (default) Taucs=1 PETsc=2
soltype = 1 # StaticLinear=0 (default) StaticNonLinear=1 soltype = 1 # StaticLinear=0 (default) StaticNonLinear=1
nstep = 400 # number of step (used only if soltype=1) nstep = 1000 # number of step (used only if soltype=1)
ftime =1. # Final time (used only if soltype=1) ftime =1. # Final time (used only if soltype=1)
tol=1.e-6 # relative tolerance for NR scheme (used only if soltype=1) tol=1.e-6 # relative tolerance for NR scheme (used only if soltype=1)
nstepArch=1 # Number of step between 2 archiving (used only if soltype=1) nstepArch=1 # Number of step between 2 archiving (used only if soltype=1)
...@@ -108,6 +108,6 @@ mysolver.archivingForceOnPhysicalGroup("Face", 5678, 2) ...@@ -108,6 +108,6 @@ mysolver.archivingForceOnPhysicalGroup("Face", 5678, 2)
mysolver.solve() mysolver.solve()
check = TestCheck() check = TestCheck()
check.equal(6.316031e-01,mysolver.getArchivedForceOnPhysicalGroup("Face", 5678, 2),1.e-3) check.equal(7.499560e-01,mysolver.getArchivedForceOnPhysicalGroup("Face", 5678, 2),1.e-3)
...@@ -295,7 +295,7 @@ check.equal(0.000000e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSm ...@@ -295,7 +295,7 @@ check.equal(0.000000e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSm
check.equal(-9.427355e-10,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 4),1.e1) 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(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(-7.110362e-08,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 4),1.e-1)
check.equal(3.691623e-01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 5),1.e-5) check.equal(2.749549e-01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 5),1.e-5)
check.equal(2.895048e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 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(-1.265934e+02,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 5),1.e-5)
check.equal(0.0,mysolver.getArchivedNodalValue(SmpRefPoint,0,mysolver.displacement),1.e-5) check.equal(0.0,mysolver.getArchivedNodalValue(SmpRefPoint,0,mysolver.displacement),1.e-5)
......
...@@ -293,7 +293,7 @@ check = TestCheck() ...@@ -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, 3),1.e-5)
check.equal(0.000000e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 4),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(-3.322070e-05,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 4),1.e-5)
check.equal(1.017297e-07,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 4),1.e-5) check.equal(1.017325e-07,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 4),1.e-5)
check.equal(2.235121e+01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 5),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(-6.599484e+01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 5),1.e-5)
check.equal(-1.150381e+02,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 5),1.e-5) check.equal(-1.150381e+02,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 5),1.e-5)
......
...@@ -287,7 +287,7 @@ mysolver.solve() ...@@ -287,7 +287,7 @@ mysolver.solve()
check = TestCheck() check = TestCheck()
check.equal(-6.996190e+11,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 3),1.e-5) check.equal(-6.996190e+11,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 3),1.e-5)
check.equal(0.000000e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 4),1.e-5) 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.765437e-08,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 4),1.e-2)
check.equal(1.786247e-08,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 4),1.e-2) check.equal(1.786247e-08,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 4),1.e-2)
check.equal(2.286854e-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(-6.912008e-01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 5),1.e-5)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment