3 files changed, 75 insertions, 0 deletions

diff --git a/utils/__init__.py b/utils/__init__.py
new file mode 100644
index 0000000..e69de29
--- /dev/null
+++ b/utils/__init__.py
diff --git a/utils/iqr.py b/utils/iqr.py
new file mode 100755
index 0000000..bbecdfa
--- /dev/null
+++ b/utils/iqr.py
@@ -0,0 +1,32 @@
+import numpy as np
+from scipy.stats import scoreatpercentile
+import bisect
+
+def find_lt(a, x):
+        i = bisect.bisect_left(a,x)
+        if i:
+                return i - 1
+        else:
+                return None
+
+def find_gt(a, x):
+        i = bisect.bisect_right(a, x)
+        if i != len(a):
+                return i
+        else:
+                return None
+
+def apply_iqr(seq, extent = 1.5):
+        q1 = scoreatpercentile(seq, 25)
+        q3 = scoreatpercentile(seq, 75)
+        iqr = q3 - q1
+        start = 0
+        end = len(seq) - 1
+        l = find_lt(seq, q1 - extent*iqr)
+        if l is not None:
+                start = l + 1
+        r = find_gt(seq, q3 + extent*iqr)
+        if r is not None:
+                end = r - 1
+        seq = seq[start:end+1]
+        return (seq, q1 - extent*iqr, q3 + extent*iqr)
diff --git a/utils/machines.py b/utils/machines.py
new file mode 100644
index 0000000..37a3c2f
--- /dev/null
+++ b/utils/machines.py
@@ -0,0 +1,43 @@
+machines = {
+        'ludwig': {
+        'cpu': 2133.0,   # mhz
+        'sockets': 4,
+        'cores_per_socket': 6,
+        'L1': 256,       # kb
+        'nL1': 1,   # private
+        'L2': 3*1024,
+        'nL2': 2,   # shared by two L1 modules
+        'L3': 12*1024,
+        'nL3': 3,   # shared by three L2 modules
+        'nMem': 4   # shared by four L3 modules
+        },
+        'bonham': {
+        'cpu': 2666.0,   # mhz
+        'sockets': 2,
+        'cores_per_socket': 6,
+        'L1': 32,        # kb
+        'nL1': 1,   # private
+        'L2': 2*1024,
+        'nL2': 1,   # private
+        'L3': 12*1024,
+        'nL3': 6,   # shared by six L2 modules
+        'nMem': 1   # shared by one L3 module (per mem)
+        },
+        'ringo': {
+        'cpu': 1400.0,   # mhz
+        'sockets': 1,
+        'cores_per_socket': 4,
+        'L1': 32,        # kb (just d-cache)
+        'nL1': 1,   # private
+        'L2': 1024,
+        'nL2': 4,   # shared by four cpus
+        'L3': 0,    # n/a
+        'nL3': 0,   # n/a
+        'nMem': 1
+        }
+}
+
+def cycles_to_us(machine_name, ncycles):
+        mhz = machines[machine_name]['cpu']
+        us = ncycles / mhz
+        return us