trimmed the test script and added a calltree functionality

testing/initial_test.py

@@ -13,29 +13,34 @@ from decimal import Decimal
 
 from synthlisa import *
 
-c = 2.99792458e8
-
-patches = 256
-OneYear = 2**25 
-NYear = 1.0
-stime = 4.0
-samples = int(NYear*OneYear/stime)
-
-lisa = EccentricInclined(2.5e9/2.99792458e8,0.0,0.,-1,0)
-
-Sl  = 0 # no laser noise since it would be canceled with first-generation TDI
-
-Sop = ((10.0e-12)**2)*(4.0*math.pi*math.pi/c**2) ## 4.4e-38
-#Spm = 0
-# Sop = 0
-Spm = 3.0e-15**2/(4.0*math.pi*math.pi*c**2)      ## 2.5e-48 
-Sop = 2*((10.0e-12)**2)*(4.0*math.pi*math.pi/c**2) ## 4.4e-38
-
-noiseTDI =  TDInoise(lisa, stime,Spm, stime,Sop, stime,0.0)
-[tn, Xmn] = np.transpose(getobs(samples,stime,[noiseTDI.t, noiseTDI.Xm]))
-specXmn = spect(Xmn,stime,patches=256)
-freqmn, spXmn = specXmn[1:,0], specXmn[1:,1]
-
+from pycallgraph import PyCallGraph
+from pycallgraph.output import GraphvizOutput
+
+with PyCallGraph(output=GraphvizOutput()):
+    
+    c = 2.99792458e8
+    
+    patches = 256
+    OneYear = 2**25 
+    NYear = 0.25
+    stime = 124.0
+    samples = int(NYear*OneYear/stime)
+    
+    lisa = EccentricInclined(2.5e9/2.99792458e8,0.0,0.,-1,0)
+    
+    Sl  = 0 # no laser noise since it would be canceled with first-generation TDI
+    
+    Sop = ((10.0e-12)**2)*(4.0*math.pi*math.pi/c**2) ## 4.4e-38
+    #Spm = 0
+    # Sop = 0
+    Spm = 3.0e-15**2/(4.0*math.pi*math.pi*c**2)      ## 2.5e-48 
+    Sop = 2*((10.0e-12)**2)*(4.0*math.pi*math.pi/c**2) ## 4.4e-38
+    
+    noiseTDI =  TDInoise(lisa, stime,Spm, stime,Sop, stime,0.0)
+    [tn, Xmn] = np.transpose(getobs(samples,stime,[noiseTDI.t, noiseTDI.Xm]))
+    specXmn = spect(Xmn,stime,patches=256)
+    freqmn, spXmn = specXmn[1:,0], specXmn[1:,1]
+    
 fig = plt.figure()
 ax = fig.add_subplot(111)
 #ax.plot(range(10))