Changeset 13451 for fact/tools/pyscripts/sandbox/kraehenb/CalFitsPerformanceWeave.py

Timestamp:

04/26/12 10:28:19 (13 years ago)

Author:

kraehenb

Message:

New comments [all], cleaned up unused code parts [CalFitsPerformanceWeave].

File:

• fact/tools/pyscripts/sandbox/kraehenb/CalFitsPerformanceWeave.py (modified) (4 diffs)

fact/tools/pyscripts/sandbox/kraehenb/CalFitsPerformanceWeave.py

@@ -12 +12 @@
 
 import numpy as np
-import itertools
 from scipy import weave
 from scipy.weave import converters
@@ -51 +50 @@
 caltest.datafile.PrintColumns()
 
-def offsum(x,y):
-    return np.sum(npcalevent[x,y+5:y+15])
-vecoffsum = np.vectorize(offsum) #specifying otypes=[np.float/double] does not improve speed
-
-#AND WE HAVE A WINNER: 35 Hz with scipy.weave!
+#AND WE HAVE A WINNER: 43 Hz with scipy.weave!
 while caltest.GetCalEvent():
     c_integrals = np.zeros(40000, np.float64) #Allocate the memory for about 40000 singles
@@ -83 +78 @@
     }
 //    std::cout << std::endl << singles << std::endl;
-    
     """
     weave.inline(cppcode, ['c_integrals', 'npcalevent'], type_converters=converters.blitz)
+    
 #    print "CalFitsPerformanceWeave"
 #    for i in range(10):
@@ -92 +87 @@
 #        if x:
 #            print x,
-    
-#    npcalevent = npcalevent.reshape((caltest.npix, caltest.nroi))
-#    npthr = npcalevent>2.5 #true if above threshold
-#    npthr = npthr[:,1:-14] * np.logical_not(npthr[:,:-15]) #only true for the zero crossings, shape (1400,285) [smaller due to necessary integration range]
-##    print [(x,y) for x,y in zip(npthr.nonzero()[0],npthr.nonzero()[1])] #print the coordinates, range 0-1399,0-284
-#    
-#    #Various versions to get the integral, all with approximately the same miserable speed (3 Hz), except for the last one (vectorized function) with ~4.3 Hz
-#    #Missing: add deadtime after an integration, remove start & end of the array
-##    integrals = [np.sum(npcalevent[x,y+5:y+15]) for x,y in np.transpose(npthr.nonzero())]
-##    integrals = [np.sum(npcalevent[x,y+5:y+15]) for x,y in zip(npthr.nonzero()[0],npthr.nonzero()[1])]
-##    integrals = [np.sum(npcalevent[x,y+5:y+15]) for x,y in itertools.izip(npthr.nonzero()[0],npthr.nonzero()[1])]
-##    integrals = map((lambda index_a,index_b: np.sum(npcalevent[index_a,index_b+5:index_b+15])),npthr.nonzero()[0],npthr.nonzero()[1])
-#    integrals = vecoffsum(npthr.nonzero()[0],npthr.nonzero()[1])
-#    print len(integrals)
-    
 #    print caltest.event_id, caltest.event_triggertype, caltest.event_caldata[10]
 #    pass