LALSimIMREOBNQCTables.c File Reference

Prototypes
static UNUSED int	XLALSimIMRGetEOBCalibratedSpinNQC (EOBNonQCCoeffs *restrict coeffs, INT4 UNUSED l, INT4 UNUSED m, REAL8 eta, REAL8 a)
	More recent versions of the EOB models, such as EOBNRv2 and SEOBNRv1, utilise a non-quasicircular correction (NQC) to bring the peak of the EOB frequency into agreement with that of NR simulations. This file contains the functions used to calculate these NQC corrections, described in DCC document T1100433. The fits to NR peak amplitude, frequency, and their derivatives, are taken from Pan et al. PRD 84 124052 (2011), for EOBNRv2, and from Taracchini et al. PRD 86, 024011 (2012), for SEOBNRv1. More...
static UNUSED int	XLALSimIMRGetEOBCalibratedSpinNQCv2 (EOBNonQCCoeffs *restrict coeffs, INT4 UNUSED l, INT4 UNUSED m, REAL8 eta, REAL8 a)
	Function to interpolate known amplitude NQC coeffcients of spin terms for SEOBNRv2, namely a3s, a4 and a5. More...
static UNUSED int	XLALSimIMRGetEOBCalibratedSpinNQCv2chiAmax (EOBNonQCCoeffs *restrict coeffs, INT4 UNUSED l, INT4 UNUSED m, REAL8 eta, REAL8 a)
static UNUSED int	XLALSimIMRGetEOBCalibratedSpinNQCv2chiAmed (EOBNonQCCoeffs *restrict coeffs, INT4 UNUSED l, INT4 UNUSED m, REAL8 eta, REAL8 a)
static UNUSED int	XLALSimIMRGetEOBCalibratedSpinNQCv2chiAmin (EOBNonQCCoeffs *restrict coeffs, INT4 UNUSED l, INT4 UNUSED m, REAL8 eta, REAL8 a)
static UNUSED int	XLALSimIMRGetEOBCalibratedSpinNQC3D (EOBNonQCCoeffs *restrict coeffs, INT4 UNUSED l, INT4 UNUSED m, REAL8 m1, REAL8 m2, REAL8 a, REAL8 chiAin)
	Function to 3D-interpolate known amplitude NQC coeffcients of spin terms for SEOBNRv2, namely a3s, a4 and a5. More...

Go to the source code of this file.

Macros
#define	OLDNSNQC   0
#define	NEWNSNQC   1

Macro Definition Documentation

OLDNSNQC

#define OLDNSNQC   0

Definition at line 1 of file LALSimIMREOBNQCTables.c.

NEWNSNQC

#define NEWNSNQC   1

Definition at line 2 of file LALSimIMREOBNQCTables.c.

Function Documentation

XLALSimIMRGetEOBCalibratedSpinNQC()

static UNUSED int XLALSimIMRGetEOBCalibratedSpinNQC ( EOBNonQCCoeffs *restrict  coeffs, INT4 UNUSED  l, INT4 UNUSED  m, REAL8  eta, REAL8  a  )

static

More recent versions of the EOB models, such as EOBNRv2 and SEOBNRv1, utilise a non-quasicircular correction (NQC) to bring the peak of the EOB frequency into agreement with that of NR simulations. This file contains the functions used to calculate these NQC corrections, described in DCC document T1100433. The fits to NR peak amplitude, frequency, and their derivatives, are taken from Pan et al. PRD 84 124052 (2011), for EOBNRv2, and from Taracchini et al. PRD 86, 024011 (2012), for SEOBNRv1.

Author

Craig Robinson, Yi Pan, Andrea Taracchini Function to interpolate known amplitude NQC coeffcients of spin terms, namely a3s, a4 and a5. The a3s, a4 and a5 values were calculated for 11 mass ratios q=1,1.5,2,3,4,5,6,10,20,50 and 100, and 19 spin ( \(\chi\) defined in Taracchini et al. PRD 86, 024011 (2012)) values chi = -1, -0.9, -0.8, ......, 0.3, 0.4, 0.5, 0.55, 0.6, 0.65. The calculation was done by Andrea Taracchini using a C++ code of the UMaryland group. In principle, these numbers can be automatically calculated iteratively by the LAL code. However, since such calcualtion increase the cost of each waveform generation by about an order of magnitude, we prepare these numbers in advance reduce cost. These number can be verified by confirming that the peak amplitude and frequency agree well with the NR-fits predicted values, and to get exact NR-fits predicted values, corrections on these numbers are ~1%.

Definition at line 52 of file LALSimIMREOBNQCTables.c.

XLALSimIMRGetEOBCalibratedSpinNQCv2()

static UNUSED int XLALSimIMRGetEOBCalibratedSpinNQCv2 ( EOBNonQCCoeffs *restrict  coeffs, INT4 UNUSED  l, INT4 UNUSED  m, REAL8  eta, REAL8  a  )

static

Function to interpolate known amplitude NQC coeffcients of spin terms for SEOBNRv2, namely a3s, a4 and a5.

The a3s, a4 and a5 values were calculated for 11 mass ratios q=1,1.5,2,4,6,8,10,11,12,14 and 20, and 19 spin ( \(\chi\) defined in Taracchini et al. PRD 86, 024011 (2012)) values chi = -1, -0.9, -0.8, ......, 0.3, 0.4, 0.5, 0.55, 0.6, 0.65. The calculation was done by Andrea Taracchini using a C++ code of the UMaryland group. In principle, these numbers can be automatically calculated iteratively by the LAL code. However, since such calcualtion increase the cost of each waveform generation by about an order of magnitude, we prepare these numbers in advance reduce cost. These number can be verified by confirming that the peak amplitude and frequency agree well with the NR-fits predicted values, and to get exact NR-fits predicted values, corrections on these numbers are ~1%.

Definition at line 369 of file LALSimIMREOBNQCTables.c.

XLALSimIMRGetEOBCalibratedSpinNQCv2chiAmax()

static UNUSED int XLALSimIMRGetEOBCalibratedSpinNQCv2chiAmax ( EOBNonQCCoeffs *restrict  coeffs, INT4 UNUSED  l, INT4 UNUSED  m, REAL8  eta, REAL8  a  )

static

---

Below are three interpolation functions for max, med and min possible chiA values with given eta and chi values. They provide 2-D interpolation results for the 3-D interpolation function of SEOBNRv2

---

Definition at line 667 of file LALSimIMREOBNQCTables.c.

XLALSimIMRGetEOBCalibratedSpinNQCv2chiAmed()

static UNUSED int XLALSimIMRGetEOBCalibratedSpinNQCv2chiAmed ( EOBNonQCCoeffs *restrict  coeffs, INT4 UNUSED  l, INT4 UNUSED  m, REAL8  eta, REAL8  a  )

static

Definition at line 1248 of file LALSimIMREOBNQCTables.c.

XLALSimIMRGetEOBCalibratedSpinNQCv2chiAmin()

static UNUSED int XLALSimIMRGetEOBCalibratedSpinNQCv2chiAmin ( EOBNonQCCoeffs *restrict  coeffs, INT4 UNUSED  l, INT4 UNUSED  m, REAL8  eta, REAL8  a  )

static

Definition at line 1829 of file LALSimIMREOBNQCTables.c.

XLALSimIMRGetEOBCalibratedSpinNQC3D()

static UNUSED int XLALSimIMRGetEOBCalibratedSpinNQC3D ( EOBNonQCCoeffs *restrict  coeffs, INT4 UNUSED  l, INT4 UNUSED  m, REAL8  m1, REAL8  m2, REAL8  a, REAL8  chiAin  )

static

Function to 3D-interpolate known amplitude NQC coeffcients of spin terms for SEOBNRv2, namely a3s, a4 and a5.

The a3s, a4 and a5 values were calculated for 11 mass ratios q=1,1.5,2,4,6,8,10,11,12,14 and 20, 19 spin ( \(\chi\) defined in Taracchini et al. PRD 86, 024011 (2012)) values chi = -1, -0.9, -0.8, ......, 0.3, 0.4, 0.5, 0.55, 0.6, 0.65, and 3 chiA values, the max and min allowed chiA for given q and chi, as well as the mid value. The calculation was done by Andrea Taracchini using a C++ code of the UMaryland group. In principle, these numbers can be automatically calculated iteratively by the LAL code. However, since such calcualtion increase the cost of each waveform generation by about an order of magnitude, we prepare these numbers in advance reduce cost. These number can be verified by confirming that the peak amplitude and frequency agree well with the NR-fits predicted values, and to get exact NR-fits predicted values, corrections on these numbers are ~1%.

Definition at line 2424 of file LALSimIMREOBNQCTables.c.