Commit 761ef9b8 authored by NV KRISHNENDU's avatar NV KRISHNENDU

LALSimIMRPhenomP.c PhenomPCore quadparams modified

parent cff4c2ac
Pipeline #136046 passed with stages
in 161 minutes and 20 seconds
......@@ -639,8 +639,8 @@ static int PhenomPCore(
REAL8 m1_SI, m2_SI;
REAL8 lambda1_in = 0.0;
REAL8 lambda2_in = 0.0;
REAL8 quadparam1_in = 1.0;
REAL8 quadparam2_in = 1.0;
REAL8 quadparam1_in = 1. + XLALSimInspiralWaveformParamsLookupdQuadMon1(extraParams);
REAL8 quadparam2_in = 1. + XLALSimInspiralWaveformParamsLookupdQuadMon2(extraParams);
if (IMRPhenomP_version == IMRPhenomPv2NRTidal_V) {
int retcode;
......@@ -648,8 +648,6 @@ static int PhenomPCore(
XLAL_CHECK(retcode == XLAL_SUCCESS, XLAL_EFUNC, "Failed to set quadparams from Universal relation.\n");
lambda1_in = XLALSimInspiralWaveformParamsLookupTidalLambda1(extraParams);
lambda2_in = XLALSimInspiralWaveformParamsLookupTidalLambda2(extraParams);
quadparam1_in = 1. + XLALSimInspiralWaveformParamsLookupdQuadMon1(extraParams);
quadparam2_in = 1. + XLALSimInspiralWaveformParamsLookupdQuadMon2(extraParams);
}
REAL8 lambda1, lambda2;
......@@ -809,12 +807,13 @@ static int PhenomPCore(
extraParams=XLALCreateDict();
}
XLALSimInspiralWaveformParamsInsertPNSpinOrder(extraParams, LAL_SIM_INSPIRAL_SPIN_ORDER_35PN);
XLALSimInspiralWaveformParamsInsertdQuadMon1(extraParams, quadparam1-1.);
XLALSimInspiralWaveformParamsInsertdQuadMon2(extraParams, quadparam2-1.);
// Start making changes here: use XLALSimInspiralWaveformParamsInsertdQuadMon1() function
if (IMRPhenomP_version == IMRPhenomPv2NRTidal_V) {
XLALSimInspiralWaveformParamsInsertTidalLambda1(extraParams, lambda1);
XLALSimInspiralWaveformParamsInsertTidalLambda2(extraParams, lambda2);
XLALSimInspiralWaveformParamsInsertdQuadMon1(extraParams, quadparam1-1.);
XLALSimInspiralWaveformParamsInsertdQuadMon2(extraParams, quadparam2-1.);
XLALSimInspiralTaylorF2AlignedPhasing(&pn, m1, m2, chi1_l, chi2_l, extraParams);
} else {
XLALSimInspiralTaylorF2AlignedPhasing(&pn, m1, m2, chi1_l, chi2_l, extraParams);
......
Markdown is supported
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