31 #ifndef _OPTIMIZEDLALSIMIMRSPINEOBFACTORIZEDFLUX_C
32 #define _OPTIMIZEDLALSIMIMRSPINEOBFACTORIZEDFLUX_C
35 #include <lal/LALSimInspiral.h>
36 #include <lal/LALSimIMR.h>
65 SpinAlignedEOBversion);
91 const UINT4 SpinAlignedEOBversion
115 omegaSq = omega * omega;
128 REAL8 rcrossp[3], rcrosspMag, s1dotL, s2dotL;
129 REAL8 chiS, chiA, tplspin;
138 sqrt (rcrossp[0] * rcrossp[0] + rcrossp[1] * rcrossp[1] +
139 rcrossp[2] * rcrossp[2]);
141 rcrossp[0] /= rcrosspMag;
142 rcrossp[1] /= rcrosspMag;
143 rcrossp[2] /= rcrosspMag;
145 s1dotL = values->
data[6] * rcrossp[0] + values->
data[7] * rcrossp[1]
146 + values->
data[8] * rcrossp[2];
147 s2dotL = values->
data[9] * rcrossp[0] + values->
data[10] * rcrossp[1]
148 + values->
data[11] * rcrossp[2];
150 chiS = 0.5 * (s1dotL + s2dotL);
151 chiA = 0.5 * (s1dotL - s2dotL);
155 switch (SpinAlignedEOBversion)
176 (
"XLAL Error - %s: Unknown SEOBNR version!\nAt present only v1 and v2 are available.\n",
207 for (
l = 2;
l <= lMax;
l++)
209 for (
m = 1;
m <=
l;
m++)
211 INT4 use_optimized_v2 = 1;
213 (&hLM, values, v,
H,
l,
m, ak, use_optimized_v2,
219 if (
l == 2 &&
m == 2)
243 (
REAL8) (
m *
m) * omegaSq * (creal (hLM) * creal (hLM) +
244 cimag (hLM) * cimag (hLM));
static UNUSED int XLALSimIMREOBNonQCCorrection(COMPLEX16 *restrict nqc, REAL8Vector *restrict values, const REAL8 omega, EOBNonQCCoeffs *restrict coeffs)
This function calculates the non-quasicircular correction to apply to the waveform.
static REAL8 XLALInspiralSpinFactorizedFluxOptimized(REAL8Vector *values, EOBNonQCCoeffs *nqcCoeffs, const REAL8 omega, SpinEOBParams *ak, const REAL8 H, const INT4 lMax, const UINT4 SpinAlignedEOBversion)
This function calculates the spin factorized-resummed GW energy flux for given dynamical variables.
static REAL8 XLALSimIMRSpinAlignedEOBNonKeplerCoeffOptimized(const REAL8 values[], SpinEOBParams *funcParams)
Function to calculate the non-Keplerian coefficient for the spin-aligned EOB model.
void XLALDestroyREAL8Vector(REAL8Vector *vector)
#define XLAL_ERROR_REAL8(...)
int XLALPrintError(const char *fmt,...) _LAL_GCC_PRINTF_FORMAT_(1
The coefficients which are used in calculating the non-quasicircular correction to the EOBNRv2 model.
FacWaveformCoeffs * hCoeffs
Parameters for the spinning EOB model.