Reactivity support - code speedup

* Code is still WIP, but working much faster.

prasopes/reactivitytools.py

@@ -31,16 +31,23 @@ def pop_dial(augCanvas, drls, graph, labels):
     names, times, intensities = drl.get_daughterset(augCanvas.ds, drls)
     if len(names) < 2:
         return
-    params = augCanvas.ms['params']
+    params = augCanvas.ms['params'][1]
+    parlen = len(params)
     pressures = []
+    lastpos = 0
     for time in times:
         toavg = []
-        for param in params[1]:
-            if param[0] == time:
-                toavg.append(param[36])
+        for i in range(lastpos,parlen):
+            if float(params[i][0]) == time:
+                #toavg.append(param[36])
+                toavg.append((-(float(params[i][1])-24.09))/1000)
+                lastpos = i
+            elif float(params[i][0]) > time:
+                break
         if len(toavg) != 0:
             pressures.append([time, np.average(toavg)])
-    goodtimes = np.where([t in np.asarray(pressures).T[0] for t in times])
+    nptpressures = np.asarray(pressures).T[0]
+    goodtimes = np.where([t in nptpressures for t in times])
     for i in range(1,len(intensities)):
         relint = np.divide(intensities[i], np.clip(np.sum(
             intensities, 0), np.finfo(np.float32).eps, None),