From 8c0eb66e0bccd24de281ca590f8eebdf02c895fc Mon Sep 17 00:00:00 2001 From: Stephen Mardson McQuay Date: Mon, 1 Nov 2010 20:19:09 -0600 Subject: [PATCH] added a bunch of random point checks --- bin/parse_gmsh.py | 21 ++++++++++++++++++++- 1 file changed, 20 insertions(+), 1 deletion(-) 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