From 5ef402982c685f80609b3489e9e9663be5ab1594 Mon Sep 17 00:00:00 2001
From: Stephen Mardson McQuay <stephen@mcquay.me>
Date: Fri, 20 May 2011 18:05:49 -0600
Subject: [PATCH] added a script to anaylze output from marylou runs

---
 gmsh/analyze_gmsh_output.py | 24 ++++++++++++++++++++++++
 1 file changed, 24 insertions(+)
 create mode 100644 gmsh/analyze_gmsh_output.py

diff --git a/gmsh/analyze_gmsh_output.py b/gmsh/analyze_gmsh_output.py
new file mode 100644
index 0000000..e02daae
--- /dev/null
+++ b/gmsh/analyze_gmsh_output.py
@@ -0,0 +1,24 @@
+import sys
+import numpy as np
+from interp.tools import rms
+
+ORDER = 0
+EP    = 1
+FINAL = 7
+EXACT = 8
+
+orders = [2,3,4,5]
+extras = [4,6,8,12, 16,20,32,48, 64,96,128,192, 256]
+
+params = (
+           (2,32),
+           (3,48),
+           (4,64),
+           (5,96),
+         )
+
+if __name__ == '__main__':
+  data = np.loadtxt(sys.argv[1])
+
+  for order, ep in params:
+    print rms([np.abs(i[FINAL] - i[EXACT]) for i in data if i[ORDER] == order and i[EP] == ep]),