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

new ref values

parent 1d2a38ec
No related branches found
No related tags found
No related merge requests found
......@@ -124,4 +124,4 @@ mysolver.archivingNodeDisplacement(86,0,1)
mysolver.solve()
check = TestCheck()
check.equal(1.060988e-08,mysolver.getEnergy(3),5.e-3)
check.equal(1.066476e-08,mysolver.getEnergy(3),5.e-3)
......@@ -86,7 +86,7 @@ mysolver = nonLinearMechSolver(numDomain1) # Solver associated with numSolver1
geofile="cube.geo"
meshfile= "cube.msh" # name of mesh file
mysolver.createModel(geofile,meshfile,3,1)
mysolver.loadModel(meshfile) # add mesh
#mysolver.loadModel(meshfile) # add mesh
mysolver.addDomain(field1) # add domain
mysolver.addMaterialLaw(BulkLaw1) # add material law
mysolver.Solver(2) # Library solving: Gmm=0 (default) Taucs=1 PETsc=2
......
......@@ -86,7 +86,7 @@ mysolver = nonLinearMechSolver(numDomain1) # Solver associated with numSolver1
geofile="twoHole.geo"
meshfile= "twoHole.msh" # name of mesh file
mysolver.createModel(geofile,meshfile,3,1)
mysolver.loadModel(meshfile) # add mesh
#mysolver.loadModel(meshfile) # add mesh
mysolver.addDomain(field1) # add domain
mysolver.addMaterialLaw(BulkLaw1) # add material law
mysolver.Solver(2) # Library solving: Gmm=0 (default) Taucs=1 PETsc=2
......@@ -133,7 +133,7 @@ mysolver.archivingNodeDisplacement(19,1, nstepArchForce)
mysolver.solve()
check = TestCheck()
check.equal(-6.777304e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", 84, 1),1.e-6)
check.equal(-6.779454e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", 84, 1),1.e-6)
......
......@@ -117,6 +117,6 @@ mysolver.archivingForceOnPhysicalGroup("Face", 5678, 2)
mysolver.solve()
check = TestCheck()
check.equal(9.612241e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", 5678, 2),1.e-6)
check.equal(9.612804e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", 5678, 2),1.e-6)
......@@ -102,6 +102,6 @@ check.equal(-2.122983e-06,mysolver.getArchivedNodalValue(8,0,mysolver.displaceme
check.equal(2.124192e-06,mysolver.getArchivedNodalValue(8,1,mysolver.displacement),1.e-5)
check.equal(4.519716e-07,mysolver.getArchivedNodalValue(8,2,mysolver.displacement),1.e-5)
check.equal(2.272999e+03,mysolver.getArchivedNodalValue(8,3,mysolver.displacement),1.e-5)
check.equal(-1.099895e-08,mysolver.getArchivedForceOnPhysicalGroup("Face", 1234, 2),1.e-4)
check.equal(-1.100211e-08,mysolver.getArchivedForceOnPhysicalGroup("Face", 1234, 2),1.e-4)
......@@ -114,4 +114,4 @@ mysolver.archivingForceOnPhysicalGroup('Face',1,1)
mysolver.solve()
check = TestCheck()
check.equal(-1.194313e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", 1, 1),1.e-4)
check.equal(-1.195570e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", 1, 1),1.e-4)
This diff is collapsed.
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment