Skip to content
Snippets Groups Projects
Commit 84e33fa7 authored by Jonathan Thompson's avatar Jonathan Thompson
Browse files

removing debug output and fixing memory leak

parent 44210e7e
No related branches found
No related tags found
1 merge request!2122Add IMRPhenomXO4a
......@@ -436,7 +436,7 @@ IMRPhenomX_UsefulPowers powers_of_lalpi;
const double eta = pWF->eta;
const double MfRD = pWF->fRING;
// const double MfRD = XLALSimPhenomUtilsHztoMf(fRD, M);
const double theta = pWF->theta_antisymmetric * LAL_PI/180;
const double theta = pWF->theta_antisymmetric * LAL_PI/180.0;
const double Chi = pPrec->chi_singleSpin_antisymmetric;
/* coefficients for phenemenological fit of amplitude ratio in PN-NR transition*/
......@@ -447,10 +447,6 @@ IMRPhenomX_UsefulPowers powers_of_lalpi;
double b = b0 + b1*eta + b2*theta + b3*eta*theta;
FILE *fp5 = fopen("debugfile.txt","w");
fprintf(fp5,"%f",m2);
fclose(fp5);
const double vRD = cbrt (LAL_PI * MfRD * (2.0 / 2.0) );
const double kappaPNRD = (21 * vRD*vRD * (1 + delta) * Chi * sin(theta))/( 2 * (42 + vRD * vRD * (55 * eta - 107)) + vRD * vRD * vRD *(84 * LAL_PI - 28 * (1 + delta - eta) * Chi * cos ( theta )));
const double kappaRD = kappaPNRD * (1 + b * vRD * vRD * vRD * vRD * vRD);
......@@ -581,6 +577,8 @@ IMRPhenomX_UsefulPowers powers_of_lalpi;
*phi_A0 = pPrec-> alpha_offset;
*phi_B0 = alpha_MfT - phi_MfT/2 + *A0 * MfT + *phi_A0;
LALFree(alphaParams);
return XLAL_SUCCESS;
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment