Commit 6afd29ee authored by NV KRISHNENDU's avatar NV KRISHNENDU
Browse files

rearranging terms in lalsimulation/lib/LALSimIMRPhenomX_internals.c

parent 015ffead
Pipeline #269469 failed with stages
in 35 minutes and 8 seconds
......@@ -1673,7 +1673,8 @@ int IMRPhenomXGetPhaseCoefficients(
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.0*powers_of_lalpi.two)/3 - 520.0*eta*powers_of_lalpi.two + 0.5*(chi1L - chi2L)*((26015.0*(qm_def1 - qm_def2))/28.0 - (44255/21.)*eta*(qm_def1 - qm_def2)- 240.0*eta*eta*(qm_def1 - qm_def2))) + 0.25*(chi1L - chi2L)*(chi1L - chi2L)*(-(1344475/2016.0) + (267815*eta)/252.0 - 240.0 * eta * eta + (26015.0*(qm_def1 + qm_def2))/56.0 - (44255/42.0)*eta*(qm_def1 + qm_def2)-120.0*eta*eta*(qm_def1 + qm_def2)) + 0.25*(chi1L + chi2L)*(chi1L + chi2L)*(-(1344475/2016.0) + (829705*eta)/504. + (3415*eta*eta)/9. + (26015 * ((qm_def1 + qm_def2))/56.0) -(44255/42.) * eta * (qm_def1 + qm_def2) - 120.0 * eta * eta * (qm_def1 + qm_def2)) + delta * (1188.57 * (chi1L - chi2L) + 0.25 * (chi1L - chi2L)* (chi1L - chi2L) * ((26015 * ((qm_def1 - qm_def2))/56.0) - (1495/12.) * eta * (qm_def1 - qm_def2)) + 0.25 * (chi1L + chi2L) * (chi1L + chi2L) * ((26015.0 * ((qm_def1 - qm_def2))/56.0) - (1495/12.) * eta * (qm_def1 - qm_def2)) + 0.25 * (chi1L - chi2L) * (chi1L + chi2L)*(-(1344475/1008.0) + ((745.0*eta)/18.0) + (26015.0 * (qm_def1 + qm_def2))/28. - (1495/6) * eta * (qm_def1 + qm_def2)))) * powers_of_lalpi.two;
pPhase->phi6 += (0.5 * (chi1L + chi2L) * (2377.14 - 520*eta *LAL_PI) + 0.25 * (chi1L - chi2L)* (chi1L - chi2L)* (-(1344475/2016) + (267815 * eta)/252 - 240 * eta * eta + (26015*qm_def1)/28 - (44255*eta*qm_def1)/21 - 240 * eta * eta * qm_def1) + 0.25 * (chi1L + chi2L)* (chi1L + chi2L)* (-(1344475/2016) + (829705*eta)/504 + (3415 *eta* eta)/9 + (26015*qm_def1)/28 - (44255*eta*qm_def1)/21 - 240*eta* eta* qm_def1) + delta * (0. + 1188.57 * (chi1L - chi2L) + 0.25 * (chi1L - chi2L)* (chi1L + chi2L) *(-(1344475/1008) + (745 * eta)/18 + (26015 * qm_def1)/14 - (1495 * eta * qm_def1)/3)))* 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;
......@@ -1682,10 +1683,10 @@ int IMRPhenomXGetPhaseCoefficients(
pPhase->phi6 += ( (5*(20*chi1L2L*eta*(11763 + 12488*eta) + 7*chi2L2*(-15103*(-1 + delta) + 2*(-21683 + 6580*delta)*eta - 9808*eta2) -
7*chi1L2*(-15103*(1 + delta) + 2*(21683 + 6580*delta)*eta + 9808*eta2)))/4032. ) * powers_of_lalpi.two;
3.0 PN, Log Term
3.0 PN, Log Term */
pPhase->phi6L = (-6848/63.) * powers_of_lalpi.two;
3.5PN, Non-Spinning */
/* 3.5PN, Non-Spinning */
pPhase->phi7 = ( (5*(15419335 + 168*(75703 - 29618*eta)*eta)*LAL_PI)/254016. ) * powers_of_lalpi.seven_thirds;
/* 3.5PN, Spin-Orbit */
......
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