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

used different exact func

Browse Source
This commit is contained in:
Stephen McQuay 2011-02-20 19:23:57 -07:00
parent da95afb14d
commit e8b15c14e2
1 changed files with 6 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
bin/parse_gmsh3D.py
View File

@ -6,7 +6,7 @@ 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, log from interp.tools import improved_answer, log
from interp.grid import grid from interp.grid import grid
@ -34,15 +34,17 @@ if __name__ == '__main__':
print "s", S print "s", S
results = {True:0, False:0} results = {True:0, False:0}
# import pdb; pdb.set_trace()
for i in xrange(1000): for i in xrange(1000):
X = np.random.random((1,3))[0] X = np.random.random((1,3))[0]
a = g.run_baker(X, order = 4, extra_points = 32) a = g.run_baker(X, order = 3, extra_points = 32)
e = exact_func_3D(X) e = exact(X)
results[improved_answer(a, e)] += 1 results[improved_answer(a, e)] += 1
print results print results
g.normalize_q() # g.normalize_q()
# print g.toxml() # print g.toxml()
igg = grid(g.verts, g.q) igg = grid(g.verts, g.q)