Merge branch 'master' of gitlab.aei.uni-hannover.de:swetashah/xtalklisapropagation

testing/initial_test.py

+import matplotlib as mpl
+mpl.use('Agg')
+import matplotlib.pyplot as plt
+plt.rcParams['agg.path.chunksize'] = 10000
+
+import numpy as np
+from math import *
+import numpy.fft as FFT
+import random
+import scipy
+from numpy import transpose
+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]
+
+fig = plt.figure()
+ax = fig.add_subplot(111)
+#ax.plot(range(10))
+ax.loglog(freqmn, spXmn, ls='-', color='green')
+fig.savefig('/vagrant/swetas_precious_results.png')