diff --git a/bin/parse_gmsh.py b/bin/parse_gmsh.py
index 633cece..5704b8f 100755
--- a/bin/parse_gmsh.py
+++ b/bin/parse_gmsh.py
@@ -5,6 +5,7 @@ import sys
 import numpy as np
 
 from interp.grid.gmsh import gmsh_grid
+from interp.tools import exact_func, improved_answer
 
 if __name__ == '__main__':
   if len(sys.argv) != 2:
@@ -14,7 +15,45 @@ if __name__ == '__main__':
   g = gmsh_grid(sys.argv[1])
   g.dump_to_blender_files()
 
-  X = np.array([0.2, 0.2, 0.0])
+  X = np.array([0.2, 0.5, 0.0])
   R = g.get_containing_simplex(X)
+  print R
+  R, S = g.get_simplex_and_nearest_points(X, 10)
+  print "r", R
+  print "s", S
 
-  print g
+  a = g.run_baker(X)
+  print a
+  e = exact_func(X)
+  print e
+  print improved_answer(a, e)
+
+  a = g.run_baker(X, extra_points = 5)
+  print a
+  e = exact_func(X)
+  print e
+  print improved_answer(a, e)
+
+  a = g.run_baker(X, extra_points = 7)
+  print a
+  e = exact_func(X)
+  print e
+  print improved_answer(a, e)
+
+  a = g.run_baker(X, order=3)
+  print a
+  e = exact_func(X)
+  print e
+  print improved_answer(a, e)
+
+  a = g.run_baker(X, order=3, extra_points = 5)
+  print a
+  e = exact_func(X)
+  print e
+  print improved_answer(a, e)
+
+  a = g.run_baker(X, order=3, extra_points = 7)
+  print a
+  e = exact_func(X)
+  print e
+  print improved_answer(a, e)