diff --git a/bin/parse_gmsh.py b/bin/parse_gmsh.py index df096ec..3b6d2d3 100755 --- a/bin/parse_gmsh.py +++ b/bin/parse_gmsh.py @@ -6,6 +6,7 @@ import numpy as np from interp.grid.gmsh import gmsh_grid from interp.tools import exact_func, improved_answer +from interp import bootstrap if __name__ == '__main__': if len(sys.argv) != 2: @@ -13,7 +14,7 @@ if __name__ == '__main__': sys.exit(1) g = gmsh_grid(sys.argv[1]) -# g.dump_to_blender_files() + # g.dump_to_blender_files() X = np.array([0.2, 0.5, 0.0]) R = g.get_containing_simplex(X) @@ -63,3 +64,21 @@ if __name__ == '__main__': e = exact_func(X) print e print improved_answer(a, e) + sys.exit(0) + + results = {True:0, False:0} + for i in xrange(1000): + X = np.random.random((1,3))[0] + + try: + a = g.run_baker(X, order=2, extra_points = 10) + e = exact_func(X) + ia = improved_answer(a, e) + if not ia: + print a['final'] - e, a, e + results[ia] += 1 + except Exception,e: + print >>sys.stderr, e + print >>sys.stderr, X + + print results