83 const REAL8 eccentricity,
96 const REAL8 eta = m1 * m2 / (
m *
m);
107 htilde = *htilde_out;
132 #if __GNUC__ >= 7 && !defined __INTEL_COMPILER
138 #if __GNUC__ >= 7 && !defined __INTEL_COMPILER
144 #if __GNUC__ >= 7 && !defined __INTEL_COMPILER
149 #if __GNUC__ >= 7 && !defined __INTEL_COMPILER
154 #if __GNUC__ >= 7 && !defined __INTEL_COMPILER
159 #if __GNUC__ >= 7 && !defined __INTEL_COMPILER
164 #if __GNUC__ >= 7 && !defined __INTEL_COMPILER
210 #if __GNUC__ >= 7 && !defined __INTEL_COMPILER
215 #if __GNUC__ >= 7 && !defined __INTEL_COMPILER
220 #if __GNUC__ >= 7 && !defined __INTEL_COMPILER
225 #if __GNUC__ >= 7 && !defined __INTEL_COMPILER
230 #if __GNUC__ >= 7 && !defined __INTEL_COMPILER
269 REAL8 v_ecc_ref = 0.0;
271 if( eccentricity > 0) {
272 v_ecc_ref = cbrt(piM*f_ecc);
282 REAL8 ref_phasing = 0.;
285 const REAL8 vref = cbrt(piM*f_ref);
286 const REAL8 logvref = log(vref);
287 const REAL8 v2ref = vref * vref;
288 const REAL8 v3ref = vref * v2ref;
289 const REAL8 v4ref = vref * v3ref;
290 const REAL8 v5ref = vref * v4ref;
291 const REAL8 v6ref = vref * v5ref;
292 const REAL8 v7ref = vref * v6ref;
293 const REAL8 v8ref = vref * v7ref;
294 const REAL8 v9ref = vref * v8ref;
295 const REAL8 v10ref = vref * v9ref;
296 const REAL8 v12ref = v2ref * v10ref;
297 const REAL8 v13ref = vref * v12ref;
298 const REAL8 v14ref = vref * v13ref;
299 const REAL8 v15ref = vref * v14ref;
300 ref_phasing += pfa7 * v7ref;
301 ref_phasing += (pfa6 + pfl6 * logvref) * v6ref;
302 ref_phasing += (pfa5 + pfl5 * logvref) * v5ref;
303 ref_phasing += pfa4 * v4ref;
304 ref_phasing += pfa3 * v3ref;
305 ref_phasing += pfa2 * v2ref;
306 ref_phasing += pfa1 * vref;
310 ref_phasing += pft15 * v15ref;
311 ref_phasing += pft14 * v14ref;
312 ref_phasing += pft13 * v13ref;
313 ref_phasing += pft12 * v12ref;
314 ref_phasing += pft10 * v10ref;
317 if( eccentricity > 0 ) {
321 ref_phasing /= v5ref;
324 #pragma omp parallel for
327 const REAL8 v = cbrt(piM*f);
328 const REAL8 logv = log(v);
329 const REAL8 v2 = v * v;
330 const REAL8 v3 = v * v2;
331 const REAL8 v4 = v * v3;
332 const REAL8 v5 = v * v4;
333 const REAL8 v6 = v * v5;
334 const REAL8 v7 = v * v6;
335 const REAL8 v8 = v * v7;
336 const REAL8 v9 = v * v8;
337 const REAL8 v10 = v * v9;
338 const REAL8 v12 = v2 * v10;
339 const REAL8 v13 = v * v12;
340 const REAL8 v14 = v * v13;
341 const REAL8 v15 = v * v14;
347 phasing += pfa7 * v7;
348 phasing += (pfa6 + pfl6 * logv) * v6;
349 phasing += (pfa5 + pfl5 * logv) * v5;
350 phasing += pfa4 * v4;
351 phasing += pfa3 * v3;
352 phasing += pfa2 * v2;
357 phasing += pft15 * v15;
358 phasing += pft14 * v14;
359 phasing += pft13 * v13;
360 phasing += pft12 * v12;
361 phasing += pft10 * v10;
364 if( eccentricity > 0 ) {
381 #if __GNUC__ >= 7 && !defined __INTEL_COMPILER
385 flux += (FTa6 + FTl6*logv) * v6;
386 dEnergy += dETa3 * v6;
387 #if __GNUC__ >= 7 && !defined __INTEL_COMPILER
392 #if __GNUC__ >= 7 && !defined __INTEL_COMPILER
397 dEnergy += dETa2 * v4;
398 #if __GNUC__ >= 7 && !defined __INTEL_COMPILER
403 #if __GNUC__ >= 7 && !defined __INTEL_COMPILER
408 dEnergy += dETa1 * v2;
409 #if __GNUC__ >= 7 && !defined __INTEL_COMPILER
419 dEnergy *= dETaN * v;
421 phasing += shft * f - 2.*phi_ref - ref_phasing;
422 amp = amp0 * sqrt(-dEnergy/flux) * v;
424 - amp * sin(phasing -
LAL_PI_4) * 1.0j;
427 *htilde_out = htilde;
460 const REAL8 eccentricity,
471 const REAL8 vISCO = 1. / sqrt(6.);
472 const REAL8 fISCO = vISCO * vISCO * vISCO / piM;
502 if ( (fEnd == 0.) && ( tideO == 0))
504 else if ( (fEnd == 0.) && ( tideO != 0 )) {
506 f_max = (fCONT > fISCO) ? fISCO : fCONT;
512 n = (size_t) (
f_max / deltaF + 1);
520 iStart = (
INT4) ceil(fStart / deltaF);
528 #pragma omp parallel for
529 for (
i = iStart;
i < n;
i++) {
530 freqs->
data[
i-iStart] =
i * deltaF;
537 f_ref, shft,
r, eccentricity,
p, &pfa);
541 *htilde_out = htilde;
int XLALSimInspiralSetQuadMonParamsFromLambdas(LALDict *LALparams)
if you do NOT provide a quadparam[1,2] term and you DO provide lamdba[1,2] then we calculate quad-mon...
REAL8 XLALSimInspiralContactFrequency(REAL8 m1_intr, REAL8 barlambda1, REAL8 m2_intr, REAL8 barlambda2)
This function estimates the radius for a binary of given masses and tidal deformability parameters.
static REAL8 UNUSED XLALSimInspiralPNFlux_5PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNEnergy_4PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_0PNCoeff(REAL8 eta)
Computes the flux PN Coefficients.
static REAL8 UNUSED XLALSimInspiralPNFlux_3PNCoeff(REAL8 UNUSED eta)
static REAL8 UNUSED XLALSimInspiralPNEnergy_2PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_4PNCoeff(REAL8 eta)
static REAL8 UNUSED eccentricityPhasing_F2(REAL8 v, REAL8 v0, REAL8 ecc, REAL8 eta, INT4 ecc_order)
Compute eccentric phase correction term using eccPNCeoffs[k][i][j].
static REAL8 UNUSED XLALSimInspiralPNEnergy_0PNCoeff(REAL8 eta)
Computes the PN Coefficients for using in the PN energy equation.
static REAL8 UNUSED XLALSimInspiralPNFlux_2PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_7PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNCoeff(REAL8 eta)
static REAL8 UNUSED XLALSimInspiralPNFlux_6PNLogCoeff(REAL8 UNUSED eta)
static void UNUSED XLALSimInspiralPNPhasing_F2(PNPhasingSeries *pfa, const REAL8 m1, const REAL8 m2, const REAL8 chi1L, const REAL8 chi2L, const REAL8 chi1sq, const REAL8 chi2sq, const REAL8 chi1dotchi2, LALDict *p)
static REAL8 UNUSED XLALSimInspiralPNEnergy_6PNCoeff(REAL8 eta)
COMPLEX16FrequencySeries * XLALCreateCOMPLEX16FrequencySeries(const CHAR *name, const LIGOTimeGPS *epoch, REAL8 f0, REAL8 deltaF, const LALUnit *sampleUnits, size_t length)
@ LAL_SIM_INSPIRAL_TIDAL_ORDER_5PN
@ LAL_SIM_INSPIRAL_TIDAL_ORDER_6PN
@ LAL_SIM_INSPIRAL_TIDAL_ORDER_ALL
@ LAL_SIM_INSPIRAL_TIDAL_ORDER_75PN
@ LAL_SIM_INSPIRAL_TIDAL_ORDER_7PN
@ LAL_SIM_INSPIRAL_TIDAL_ORDER_0PN
@ LAL_SIM_INSPIRAL_TIDAL_ORDER_65PN
int XLALSimInspiralTaylorF2Ecc(COMPLEX16FrequencySeries **htilde_out, const REAL8 phi_ref, const REAL8 deltaF, const REAL8 m1_SI, const REAL8 m2_SI, const REAL8 S1z, const REAL8 S2z, const REAL8 fStart, const REAL8 fEnd, const REAL8 f_ref, const REAL8 r, const REAL8 eccentricity, LALDict *p)
Computes the stationary phase approximation to the Fourier transform of a chirp waveform with eccentr...
int XLALSimInspiralTaylorF2CoreEcc(COMPLEX16FrequencySeries **htilde_out, const REAL8Sequence *freqs, const REAL8 phi_ref, const REAL8 m1_SI, const REAL8 m2_SI, const REAL8 f_ref, const REAL8 shft, const REAL8 r, const REAL8 eccentricity, LALDict *p, PNPhasingSeries *pfaP)
void XLALDestroyREAL8Sequence(REAL8Sequence *sequence)
REAL8Sequence * XLALCreateREAL8Sequence(size_t length)
const LALUnit lalStrainUnit
const LALUnit lalSecondUnit
LALUnit * XLALUnitMultiply(LALUnit *output, const LALUnit *unit1, const LALUnit *unit2)
#define XLAL_CHECK(assertion,...)
LIGOTimeGPS * XLALGPSAdd(LIGOTimeGPS *epoch, REAL8 dt)
REAL8 vlogv[PN_PHASING_SERIES_MAX_ORDER+1]
REAL8 v[PN_PHASING_SERIES_MAX_ORDER+1]