minor: notes and moved a file for the resolution study

--HG--
rename : gmsh/resolution_3D.py => plots/resolution/resolution_3D.py

plots/resolution/resolution.3D.gmsh.out

@@ -1,3 +1,6 @@
+# generated by running resolution_3D.py, and then by calculating x by dividing
+# the length of the domain by the cube root of number of points.
+
 0.0454730588061  0.000270108596188 8.53500404066e-05 1.41791618738e-05 4.9979436523e-06
 0.00997060697355 2.0540042226e-06 1.69775361716e-07 1.01412678674e-08 8.0860277895e-10
 0.00461278449044 2.18695319382e-07 1.52441703341e-08 3.24157996636e-10 2.18935788738e-11

plots/resolution/resolution_3D.py

@@ -0,0 +1,27 @@
+import sys
+import sqlite3
+import numpy as np
+
+from interp.tools import rms
+
+con = sqlite3.connect(sys.argv[1])
+cur = con.cursor()
+
+orders = np.array(cur.execute('select distinct ord from results').fetchall())[:,0]
+eps    = np.array(cur.execute('select distinct ep  from results').fetchall())[:,0]
+res    = np.array(cur.execute('select distinct res from results').fetchall())[:,0]
+
+params = (
+           (2,64),
+           (3,64),
+           (4,64),
+           (5,64),
+         )
+for r in res:
+  for order, ep in params:
+    # print np.array([np.abs(i[FINAL] - i[EXACT]) for i in data if i[ORDER] == order and i[EP] == ep])),
+    cur_results = cur.execute('select abs(exact - final) from results where res = ? and ep = ? and ord = ?',
+                             (int(r), int(ep),int(order)))
+    cur_ = np.array(cur_results.fetchall())
+    print rms(np.array(cur_)),
+  print