From 34e32fa6bd04f5f498b5f32b7888590c8276205f Mon Sep 17 00:00:00 2001
From: ChunHuangPhy <chunhuang@mails.ccnu.edu.cn>
Date: Mon, 18 Nov 2024 08:40:21 -0600
Subject: [PATCH] add chirp massto front line

---
 InferenceWorkflow/Likelihood.py | 3 +--
 1 file changed, 1 insertion(+), 2 deletions(-)

diff --git a/InferenceWorkflow/Likelihood.py b/InferenceWorkflow/Likelihood.py
index 50939f3..aaca249 100644
--- a/InferenceWorkflow/Likelihood.py
+++ b/InferenceWorkflow/Likelihood.py
@@ -60,19 +60,18 @@ def TidalLikihood_kernel(eps_total,pres_total,x,d1):
     """
     kernelGW,chrip = x
     MRT = []
+    chrip_mass = chrip.resample(1)
     if d1 ==0 :
         likelihood = -1e101
     else:
         d1 = 10**(d1)
         if   all(x<y for x,y in zip(eps_total[:], eps_total[1:])) and all(x<y for x, y in zip(pres_total[:], pres_total[1:])):
             MRT = main.OutputMRTpoint(d1,eps_total,pres_total).T
-            chrip_mass = chrip.resample(1)
             M1 = TOV_solver.m1_from_mc_m2(chrip_mass, MRT[1][0])
             Tidal_line = main.OutputMRT('',eps_total,pres_total).T
         if len(MRT[0]) == False or len(Tidal_line[0]) == False:
             likelihood = -1e101
         else:
-            chrip_mass = chrip.resample(1)
             MTspline = interpolate.interp1d(Tidal_line[1]/Msun,Tidal_line[2])
             point = np.array([[chrip_mass[0][0]], [MRT[1][0] / M1[0][0]],[ MTspline(M1)[0][0]], [MRT[2][0]]])
             likelihood = np.log(kernelGW.evaluate(point))