From 17b22c173d774ac01f63c997e2a1ad0e4f873bc3 Mon Sep 17 00:00:00 2001 From: Stephen Mardson McQuay Date: Mon, 1 Nov 2010 15:51:40 -0600 Subject: [PATCH] added a 3d gmsh test script --- bin/parse_gmsh3D.py | 33 +++++++++++++++++++++++++++++++++ 1 file changed, 33 insertions(+) create mode 100755 bin/parse_gmsh3D.py diff --git a/bin/parse_gmsh3D.py b/bin/parse_gmsh3D.py new file mode 100755 index 0000000..90cfb11 --- /dev/null +++ b/bin/parse_gmsh3D.py @@ -0,0 +1,33 @@ +#!/usr/bin/env python + +import sys + +import numpy as np + +import interp.bootstrap +from interp.grid.gmsh import gmsh_grid3D +from interp.tools import exact_func_3D, improved_answer + +if __name__ == '__main__': + if len(sys.argv) != 2: + print >> sys.stderr, "usage: %s " % sys.argv[0] + sys.exit(1) + + g = gmsh_grid3D(sys.argv[1]) + + X = np.array([0.2, 0.5, 0.01]) + R = g.get_containing_simplex(X) + print R + R, S = g.get_simplex_and_nearest_points(X, 10) + print "r", R + print "s", S + + + for i in xrange(10): + X = np.array([0.2,0.5,0.01])#np.random.random((1,3)) + + a = g.run_baker(X) + print a + e = exact_func_3D(X) + print e + print improved_answer(a, e)