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

reference values

parent 6071fe3b
Branches
Tags ver3.2
No related merge requests found
...@@ -194,4 +194,4 @@ mysolver.archivingIPOnPhysicalGroup("Face", 700, IPField.LOCAL_POROSITY,IPField. ...@@ -194,4 +194,4 @@ mysolver.archivingIPOnPhysicalGroup("Face", 700, IPField.LOCAL_POROSITY,IPField.
mysolver.solve() mysolver.solve()
check = TestCheck() check = TestCheck()
check.equal(-8.920918e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-4) check.equal(-8.848617e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-4)
...@@ -189,4 +189,4 @@ mysolver.archivingIPOnPhysicalGroup("Face", 700, IPField.LOCAL_POROSITY,IPField. ...@@ -189,4 +189,4 @@ mysolver.archivingIPOnPhysicalGroup("Face", 700, IPField.LOCAL_POROSITY,IPField.
mysolver.solve() mysolver.solve()
check = TestCheck() check = TestCheck()
check.equal(-2.027542e+03,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-6) check.equal(-2.496462e+03,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-6)
...@@ -198,7 +198,7 @@ check = TestCheck() ...@@ -198,7 +198,7 @@ check = TestCheck()
import csv import csv
data = csv.reader(open('IPVolume11val_LOCAL_POROSITYMax.csv'), delimiter=';') data = csv.reader(open('IPVolume11val_LOCAL_POROSITYMax.csv'), delimiter=';')
porosity = list(data) porosity = list(data)
check.equal(4.828137e-01,float(porosity[-1][1]),1e-6) check.equal(4.597634e-01,float(porosity[-1][1]),1e-6)
...@@ -156,7 +156,7 @@ mysolver.solve() ...@@ -156,7 +156,7 @@ mysolver.solve()
# Test # Test
# =========== # ===========
check = TestCheck() check = TestCheck()
check.equal(9.528820e+02,mysolver.getArchivedForceOnPhysicalGroup("Face", 31, 0),1.e-6) check.equal(9.741993e+02,mysolver.getArchivedForceOnPhysicalGroup("Face", 31, 0),1.e-6)
import csv import csv
data = csv.reader(open('IPVolume29val_LOCAL_POROSITYMean.csv'), delimiter=';') data = csv.reader(open('IPVolume29val_LOCAL_POROSITYMean.csv'), delimiter=';')
......
...@@ -159,4 +159,4 @@ mysolver.archivingIPOnPhysicalGroup("Face", 700, IPField.LOCAL_POROSITY,IPField. ...@@ -159,4 +159,4 @@ mysolver.archivingIPOnPhysicalGroup("Face", 700, IPField.LOCAL_POROSITY,IPField.
mysolver.solve() mysolver.solve()
check = TestCheck() check = TestCheck()
check.equal(-4.332014e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-6) check.equal(-4.303072e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-6)
...@@ -149,7 +149,7 @@ import csv ...@@ -149,7 +149,7 @@ import csv
data = csv.reader(open('IPVolume29val_LOCAL_POROSITYMean.csv'), delimiter=';') data = csv.reader(open('IPVolume29val_LOCAL_POROSITYMean.csv'), delimiter=';')
porosity = list(data) porosity = list(data)
check.equal(3.406171e-01,float(porosity[-1][1]),1e-5) check.equal(3.407290e-01,float(porosity[-1][1]),1e-5)
data = csv.reader(open('IPVolume29val_LIGAMENT_RATIOMean.csv'), delimiter=';') data = csv.reader(open('IPVolume29val_LIGAMENT_RATIOMean.csv'), delimiter=';')
porosity = list(data) porosity = list(data)
......
...@@ -108,6 +108,6 @@ check.equal(9.200879e+02,mysolver.getHomogenizedTangent(0,1,0,1),1.e-4) ...@@ -108,6 +108,6 @@ check.equal(9.200879e+02,mysolver.getHomogenizedTangent(0,1,0,1),1.e-4)
check.equal(1.720444e+02,mysolver.getHomogenizedSecondTangent(0,0,0,0,0,0),1.e-4) check.equal(1.720444e+02,mysolver.getHomogenizedSecondTangent(0,0,0,0,0,0),1.e-4)
check.equal(8.453953e+01,mysolver.getHomogenizedSecondTangent(1,0,1,1,0,1),1.e-4) check.equal(8.453953e+01,mysolver.getHomogenizedSecondTangent(1,0,1,1,0,1),1.e-4)
check.equal(78.2909010986,mysolver.getHomogenizedSecondTangent(1,0,0,1,0,0),1.e-4) check.equal(7.866244e+01,mysolver.getHomogenizedSecondTangent(1,0,0,1,0,0),1.e-4)
check.equal(28.491995815,mysolver.getHomogenizedSecondTangent(0,1,0,1,0,0),1.e-4) check.equal(28.491995815,mysolver.getHomogenizedSecondTangent(0,1,0,1,0,0),1.e-4)
...@@ -248,7 +248,7 @@ mysolver.solve() ...@@ -248,7 +248,7 @@ mysolver.solve()
check = TestCheck() check = TestCheck()
check.equal(-6.476339e+04,mysolver.getArchivedForceOnPhysicalGroup("Face", 85, 3),1.e-5) check.equal(-6.476339e+04,mysolver.getArchivedForceOnPhysicalGroup("Face", 85, 3),1.e-5)
check.equal(-2.097632e+05,mysolver.getArchivedForceOnPhysicalGroup("Face", 85, 4),1.e-5) check.equal(-2.097693e+05,mysolver.getArchivedForceOnPhysicalGroup("Face", 85, 4),1.e-5)
check.equal(3.355705e-03,mysolver.getArchivedNodalValue(88,3,mysolver.displacement),1.e-5) check.equal(3.355705e-03,mysolver.getArchivedNodalValue(88,3,mysolver.displacement),1.e-5)
check.equal(-3.355705e-02,mysolver.getArchivedNodalValue(88,4,mysolver.displacement),1.e-5) check.equal(-3.355705e-02,mysolver.getArchivedNodalValue(88,4,mysolver.displacement),1.e-5)
check.equal(3.276348e-03,mysolver.getArchivedNodalValue(89,3,mysolver.displacement),1.e-5) check.equal(3.276348e-03,mysolver.getArchivedNodalValue(89,3,mysolver.displacement),1.e-5)
......
...@@ -226,7 +226,7 @@ check.equal(-2.462150e-05,mysolver.getArchivedForceOnPhysicalGroup("Face", 247, ...@@ -226,7 +226,7 @@ check.equal(-2.462150e-05,mysolver.getArchivedForceOnPhysicalGroup("Face", 247,
check.equal(3.584923e-05,mysolver.getArchivedForceOnPhysicalGroup("Face", 247, 2),1.e-6) check.equal(3.584923e-05,mysolver.getArchivedForceOnPhysicalGroup("Face", 247, 2),1.e-6)
check.equal(4.282709e+04,mysolver.getArchivedForceOnPhysicalGroup("Face", 247, 3),1.e-6) check.equal(4.282709e+04,mysolver.getArchivedForceOnPhysicalGroup("Face", 247, 3),1.e-6)
check.equal(4.429112e+05,mysolver.getArchivedForceOnPhysicalGroup("Face", 247, 4),1.e-6) check.equal(4.429112e+05,mysolver.getArchivedForceOnPhysicalGroup("Face", 247, 4),1.e-6)
check.equal(6.422382e-09,mysolver.getArchivedNodalValue(93,0,mysolver.displacement),1.e-6) check.equal(6.028935e-09,mysolver.getArchivedNodalValue(93,0,mysolver.displacement),1.e-1)
check.equal(3.570000e+02,mysolver.getArchivedNodalValue(97,3,mysolver.displacement),1.e-6) check.equal(3.570000e+02,mysolver.getArchivedNodalValue(97,3,mysolver.displacement),1.e-6)
check.equal(5.800000e-02,mysolver.getArchivedNodalValue(97,4,mysolver.displacement),1.e-6) check.equal(5.800000e-02,mysolver.getArchivedNodalValue(97,4,mysolver.displacement),1.e-6)
check.equal(3.585401e+02,mysolver.getArchivedNodalValue(96,3,mysolver.displacement),1.e-6) check.equal(3.585401e+02,mysolver.getArchivedNodalValue(96,3,mysolver.displacement),1.e-6)
......
...@@ -263,12 +263,12 @@ mysolver.createRestartBySteps(300); ...@@ -263,12 +263,12 @@ mysolver.createRestartBySteps(300);
mysolver.solve() mysolver.solve()
check = TestCheck() check = TestCheck()
check.equal(6.859186e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", 102, 0),5.e-3) check.equal(6.859188e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", 102, 0),5.e-3)
mysolver.disableResetRestart() #only battery: the second solve mimics a restart mysolver.disableResetRestart() #only battery: the second solve mimics a restart
mysolver.createRestartBySteps(300); mysolver.createRestartBySteps(300);
mysolver.solve() mysolver.solve()
check = TestCheck() check = TestCheck()
check.equal(6.859186e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", 102, 0),5.e-3) check.equal(6.859188e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", 102, 0),5.e-3)
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment