Commit d132237f authored by NV Krishnendu's avatar NV Krishnendu
Browse files

rewriting the 3pn ss so term

parent a900445f
Pipeline #419604 failed
......@@ -1672,7 +1672,28 @@ int IMRPhenomXGetPhaseCoefficients(
/* 3.0 PN, Non-Spinning */
pPhase->phi6 = ( 11583231236531/4.69421568e9 - (5*eta*(3147553127 + 588*eta*(-45633 + 102260*eta)))/3.048192e6 - (6848*LAL_GAMMA)/21. -
(640*powers_of_lalpi.two)/3. + (2255*eta*powers_of_lalpi.two)/12. - (13696*log(2))/21. - (6848*powers_of_lalpi.log)/63. ) * powers_of_lalpi.two;
pPhase->phi6 += (0.5 * (chi1L + chi2L) * ((2270 *LAL_PI)/3 - 520 * eta * LAL_PI + 0.5 * (chi1L - chi2L) * ((26015 * (qm_def1 - qm_def2))/28 - (44255/21) * eta * (qm_def1 - qm_def2) - 240.0 * eta * eta * (qm_def1 - qm_def2))) + 0.25 * (chi1L - chi2L) * (chi1L - chi2L) * (-(1344475/2016) + (267815 * eta)/252 - 240 * eta * eta + (26015 * (qm_def1 + qm_def2))/56 - (44255/42) * eta * (qm_def1 + qm_def2) - 120 * eta * eta * (qm_def1 + qm_def2)) + 0.25 * (chi1L + chi2L) *(chi1L + chi2L) * (-(1344475/2016) + (829705 * eta)/504 + (3415 * eta* eta)/9 + (26015 * (qm_def1 + qm_def2))/56 - (44255/42) * eta * (qm_def1 + qm_def2) - 120 * eta * eta * (qm_def1 + qm_def2)) + delta * (1188.57 * (chi1L - chi2L) + 0.25 * (chi1L - chi2L) * (chi1L - chi2L) * ((26015 * (qm_def1 - qm_def2))/56 - (1495/12) * eta * (qm_def1 - qm_def2)) + 0.25 * (chi1L + chi2L) * (chi1L + chi2L) * ((26015 * (qm_def1 - qm_def2))/56 - (1495/12) * eta * (qm_def1 - qm_def2)) + 0.25 * (chi1L - chi2L) * (chi1L + chi2L) * (-(1344475/1008) + (745 * eta)/18 + (26015 * (qm_def1 + qm_def2))/28 - (1495/6) * eta * (qm_def1 + qm_def2))))* powers_of_lalpi.two;
/* 3.0 PN, Spin-Orbit 3.0 PN, Spin-Spin */
/*the terms below has decimal numbers and changing it on 1 july 22. New terms below. */
/* pPhase->phi6 += (0.5 * (chi1L + chi2L) * ((2270 *LAL_PI)/3 - 520 * eta * LAL_PI + 0.5 * (chi1L - chi2L) * ((26015 * (qm_def1 - qm_def2))/28 - (44255/21) * eta * (qm_def1 - qm_def2) - 240.0 * eta * eta * (qm_def1 - qm_def2))) + 0.25 * (chi1L - chi2L) * (chi1L - chi2L) * (-(1344475/2016) + (267815 * eta)/252 - 240 * eta * eta + (26015 * (qm_def1 + qm_def2))/56 - (44255/42) * eta * (qm_def1 + qm_def2) - 120 * eta * eta * (qm_def1 + qm_def2)) + 0.25 * (chi1L + chi2L) *(chi1L + chi2L) * (-(1344475/2016) + (829705 * eta)/504 + (3415 * eta* eta)/9 + (26015 * (qm_def1 + qm_def2))/56 - (44255/42) * eta * (qm_def1 + qm_def2) - 120 * eta * eta * (qm_def1 + qm_def2)) + delta * (1188.57 * (chi1L - chi2L) + 0.25 * (chi1L - chi2L) * (chi1L - chi2L) * ((26015 * (qm_def1 - qm_def2))/56 - (1495/12) * eta * (qm_def1 - qm_def2)) + 0.25 * (chi1L + chi2L) * (chi1L + chi2L) * ((26015 * (qm_def1 - qm_def2))/56 - (1495/12) * eta * (qm_def1 - qm_def2)) + 0.25 * (chi1L - chi2L) * (chi1L + chi2L) * (-(1344475/1008) + (745 * eta)/18 + (26015 * (qm_def1 + qm_def2))/28 - (1495/6) * eta * (qm_def1 + qm_def2))))* powers_of_lalpi.two; */
pPhase->phi6 += (chi1L*chi2L*((32675*eta)/112 + (5575*eta*eta)/18) +
chi1L*chi1L*(-1344475/4032 - (1344475*delta)/4032 + (1365335*eta)/2016 +
(745*delta*eta)/72 + (1255*eta*eta)/36 + (26015*(qm_def1 - qm_def2))/112 +
(26015*delta*(qm_def1 - qm_def2))/112 - (44255*eta*(qm_def1 - qm_def2))/84 -
(1495*delta*eta*(qm_def1 - qm_def2))/24 - 60*eta*eta*(qm_def1 - qm_def2) +
(26015*(qm_def1 + qm_def2))/112 + (26015*delta*(qm_def1 + qm_def2))/112 -
(44255*eta*(qm_def1 + qm_def2))/84 - (1495*delta*eta*(qm_def1 + qm_def2))/24 -
60*eta*eta*(qm_def1 + qm_def2)) + chi2L*chi2L*(-1344475/4032 + (1344475*delta)/4032 +
(1365335*eta)/2016 - (745*delta*eta)/72 + (1255*eta*eta)/36 -
(26015*(qm_def1 - qm_def2))/112 + (26015*delta*(qm_def1 - qm_def2))/112 +
(44255*eta*(qm_def1 - qm_def2))/84 - (1495*delta*eta*(qm_def1 - qm_def2))/24 +
60*eta*eta*(qm_def1 - qm_def2) + (26015*(qm_def1 + qm_def2))/112 -
(26015*delta*(qm_def1 + qm_def2))/112 - (44255*eta*(qm_def1 + qm_def2))/84 +
(1495*delta*eta*(qm_def1 + qm_def2))/24 - 60*eta*eta*(qm_def1 + qm_def2)) +
chi2L*((1135*(LAL_PI))/3 - (1135*delta*(LAL_PI))/3 - 260*eta*(LAL_PI)) +
chi1L*((1135*(LAL_PI))/3 + (1135*delta*(LAL_PI))/3 - 260*eta*(LAL_PI)))* powers_of_lalpi.two;
/* 3.0 PN, Spin-Orbit
pPhase->phi6 += ( (5*(227*(chi1L + chi2L + chi1L*delta - chi2L*delta) - 156*(chi1L + chi2L)*eta)*LAL_PI)/3. ) * powers_of_lalpi.two;
......
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