sm/smbinterp
1
0
Fork 0
Code Issues Pull Requests Activity

still proving out the 3D gmsh class, added correct funcs (3D) to the driver

Browse Source
This commit is contained in:
Stephen McQuay 2010-11-01 20:20:05 -06:00
parent 8c0eb66e0b
commit 8dff2747fd
1 changed files with 11 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
bin/parse_gmsh3D.py
View File

@ -5,8 +5,9 @@ import sys
import numpy as np import numpy as np
import interp.bootstrap import interp.bootstrap
from interp.grid.gmsh import gmsh_grid3D from interp.grid.gmsh import gmsh_grid3D
from interp.tools import exact_func_3D, improved_answer from interp.tools import exact_func_3D, improved_answer, log
if __name__ == '__main__': if __name__ == '__main__':
if len(sys.argv) != 2: if len(sys.argv) != 2:
@ -22,12 +23,13 @@ if __name__ == '__main__':
print "r", R print "r", R
print "s", S print "s", S
# results = {True:0, False:0}
# for i in xrange(10):
# X = np.random.random((1,3))[0]
for i in xrange(10): # a = g.run_baker(X, order = 5, extra_points = 25)
X = np.array([0.2,0.5,0.01])#np.random.random((1,3)) # e = exact_func_3D(X)
# ia = improved_answer(a, e)
# results[ia] += 1
a = g.run_baker(X) # print results
print a
e = exact_func_3D(X)
print e
print improved_answer(a, e)