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

reference

parent 94a7d413
No related branches found
No related tags found
No related merge requests found
...@@ -281,7 +281,7 @@ check.equal(6.859167e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", 102, 0 ...@@ -281,7 +281,7 @@ check.equal(6.859167e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", 102, 0
data21 = csv.reader(open('force102comp0_part2.csv'), delimiter=';') data21 = csv.reader(open('force102comp0_part2.csv'), delimiter=';')
force2 = list(data21) force2 = list(data21)
check.equal(1.235179e+03,float(force2[-1][1]),1e-5) check.equal(1.234095e+03,float(force2[-1][1]),1e-5)
#mysolver.disableResetRestart() #only battery: the second solve mimics a restart #mysolver.disableResetRestart() #only battery: the second solve mimics a restart
......
import os,csv import os,csv
import sys import sys
def checkEqual(ref, cur, tol): def checkEqual(part,ref, cur, tol):
if(abs(ref-cur)>abs(tol*ref)): if(abs(ref-cur)>abs(tol*ref)):
print("Error : reference value ",ref," current value ",cur," relative error ", (ref-cur)/ref," tolerance ",tol) print("Error Part ", part,": reference value ",ref," current value ",cur," relative error ", (ref-cur)/ref," tolerance ",tol)
if sys.version_info[0] < 3: if sys.version_info[0] < 3:
...@@ -13,7 +13,7 @@ else: ...@@ -13,7 +13,7 @@ else:
data1 = csv.reader(open('force102comp0_part2.csv'), delimiter=';') data1 = csv.reader(open('force102comp0_part2.csv'), delimiter=';')
force = list(data1) force = list(data1)
checkEqual(1235.178839997943,float(force[-1][1]),1e-5) checkEqual(0,1235.178839997943,float(force[-1][1]),1e-5)
if sys.version_info[0] < 3: if sys.version_info[0] < 3:
os.system('mpiexec -np 4 python fullPlateHoleRes.py') os.system('mpiexec -np 4 python fullPlateHoleRes.py')
...@@ -22,5 +22,5 @@ else: ...@@ -22,5 +22,5 @@ else:
data21 = csv.reader(open('force102comp0_part2.csv'), delimiter=';') data21 = csv.reader(open('force102comp0_part2.csv'), delimiter=';')
force2 = list(data21) force2 = list(data21)
checkEqual(1235.178839997943,float(force2[-1][1]),1e-5) checkEqual(1,1.234095e+03,float(force2[-1][1]),1e-5)
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment