Commit 2f586195 authored by Marcus Haberland's avatar Marcus Haberland
Browse files

actual progress in our numerical computations and new insights

parent 6a8142fc
This diff is collapsed.
......@@ -148,7 +148,7 @@ class binary:
phi_S_det : The azimuth angle phi_S of the source in the detector frame
'''
arctan_term = np.arctan2((np.sqrt(3)*np.cos(self.theta_S_Ec)+np.sin(self.theta_S_Ec)*np.cos(omega_E*t-self.phi_S_Ec)),(2*np.sin(self.theta_S_Ec)*np.sin(omega_E*t-self.phi_S_Ec)))
return (omega_E*t - arctan_term) % (2*pi)
return (omega_E*t + arctan_term) % (2*pi)
def psi_S(self,t=0):
'''
......@@ -235,10 +235,9 @@ class binary:
raise ValueError
def h_ixh_j(self,i,j,diff=1e-6):
funcA = self.h_i(i,diff)
funcB = self.h_i(j,diff)
return integrate.quad(lambda t: funcA(t)*funcB(t),0,self.T_obs,limit=int(self.T_obs*self.f_GW*20),epsabs=1.49e-16)[0]
return integrate.quad(lambda t: funcA(t)*1e20*funcB(t)*1e20,0,self.T_obs,limit=int(self.T_obs*self.f_GW*20),epsrel=1.49e-2,epsabs=1.49e-2)[0]*1e-40
def h_i(self,i,diff=1e-6):
B = None
......@@ -276,7 +275,7 @@ class binary:
mat[i][j] += self.h_ixh_j(i, j,diff)
if j != i:
mat[j][i] = mat[i][j]
self.Fisher = mat #2*mat/S_n
self.Fisher = 2*mat/self.S_n() #2*mat/S_n
return self.Fisher
def rel_uncertainty(self,mode=1,diff=1e-6):
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment