LALSimBlackHoleRingdown.h   LALSimBlackHoleRingdown.h 
skipping to change at line 24 skipping to change at line 24
* You should have received a copy of the GNU General Public License * You should have received a copy of the GNU General Public License
* along with with program; see the file COPYING. If not, write to the * along with with program; see the file COPYING. If not, write to the
* Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, * Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston,
* MA 02111-1307 USA * MA 02111-1307 USA
*/ */
#ifndef _LALSIMBLACKHOLERINGDOWN_H #ifndef _LALSIMBLACKHOLERINGDOWN_H
#define _LALSIMBLACKHOLERINGDOWN_H #define _LALSIMBLACKHOLERINGDOWN_H
#include <lal/LALDatatypes.h> #include <lal/LALDatatypes.h>
#include <lal/LALSimInspiral.h>
#include <gsl/gsl_linalg.h>
#include <gsl/gsl_interp.h>
#include <gsl/gsl_spline.h>
#if defined(__cplusplus) #if defined(__cplusplus)
extern "C" { extern "C" {
#elif 0 #elif 0
} /* so that editors will match preceding brace */ } /* so that editors will match preceding brace */
#endif #endif
/* LOW-LEVEL ROUTINES (USE LEAVER'S CONVENSIONS G = c = 2M = 1) */ /* LOW-LEVEL ROUTINES (USE LEAVER'S CONVENSIONS G = c = 2M = 1) */
/** /**
skipping to change at line 157 skipping to change at line 161
double deltaT, /**< sampling interval (s) */ double deltaT, /**< sampling interval (s) */
double mass, /**< black hole mass (kg) */ double mass, /**< black hole mass (kg) */
double dimensionless_spin, /**< black hole dimensionless spin p arameter */ double dimensionless_spin, /**< black hole dimensionless spin p arameter */
double fractional_mass_loss, /**< fraction of mass radiated in th is mode */ double fractional_mass_loss, /**< fraction of mass radiated in th is mode */
double distance, /**< distance to source (m) */ double distance, /**< distance to source (m) */
double inclination, /**< inclination of source's spin ax is (rad) */ double inclination, /**< inclination of source's spin ax is (rad) */
int l, /**< polar mode number */ int l, /**< polar mode number */
int m /**< azimuthal mode number */ int m /**< azimuthal mode number */
); );
/**
* Computes the final mass and spin of the black hole resulting from merger
.
* They are given by fittings of NR simulations results. Specifically,
* for EOBNR, Table I of Buonanno et al. PRD76, 104049;
* for EOBNRv2 and EOBNRv2HM, Eqs. 29a and 29b of Pan et al. PRD84, 124052;
* for SEOBNRv1, Eq. 8 of Tichy and Marronetti PRD78, 081501 and
* Eqs. 1 and 3 of Barausse and Rezzolla ApJ704, L40.
*/
INT4 XLALSimIMREOBFinalMassSpin(
REAL8 *finalMass, /**<< OUTPUT, the final mass (scaled by origin
al total mass) */
REAL8 *finalSpin, /**<< OUTPUT, the final spin (scaled by final
mass) */
const REAL8 mass1, /**<< The mass of the 1st component of the sy
stem */
const REAL8 mass2, /**<< The mass of the 2nd component of the sy
stem */
const REAL8 spin1[3], /**<< The spin of the 1st object; only needed
for spin waveforms */
const REAL8 spin2[3], /**<< The spin of the 2nd object; only needed
for spin waveforms */
Approximant approximant /**<< The waveform approximant being used */
);
/**
* This function generates the quasinormal mode frequencies for a black
* hole ringdown. At present, this function works for the 22, 21, 33, 44
* and 55 modes, and includes 8 overtones. The final frequencies are
* computed by interpolating the data found on the webpage of
* Vitor Cardoso, http://centra.ist.utl.pt/~vitor/?page=ringdown
* In this page, frequecy data are given for positive final spins only.
* For a negative final spin chi<0 case, the (l,m) mode frequency is given
by
* the (l,-m) mode frequency of the positive final spin -chi case.
*/
INT4 XLALSimIMREOBGenerateQNMFreqV2(
COMPLEX16Vector *modefreqs, /**<< OUTPUT, complex freqs of overtones (sca
led by total mass) */
const REAL8 mass1, /**<< The mass of the 1st component (in Solar
masses) */
const REAL8 mass2, /**<< The mass of the 2nd component (in Solar
masses) */
const REAL8 spin1[3], /**<< The spin of the 1st object; only needed
for spin waveforms */
const REAL8 spin2[3], /**<< The spin of the 2nd object; only needed
for spin waveforms */
UINT4 l, /**<< The l value of the mode in question */
UINT4 m, /**<< The m value of the mode in question */
UINT4 nmodes, /**<< The number of overtones that should be
included (max 8) */
Approximant approximant/**<< The waveform approximant being used */
);
#if 0 #if 0
{ /* so that editors will match succeeding brace */ { /* so that editors will match succeeding brace */
#elif defined(__cplusplus) #elif defined(__cplusplus)
} }
#endif #endif
#endif /* _LALSIMBLACKHOLERINGDOWN_H */ #endif /* _LALSIMBLACKHOLERINGDOWN_H */
 End of changes. 2 change blocks. 
0 lines changed or deleted 58 lines changed or added


 LALSimIMR.h   LALSimIMR.h 
skipping to change at line 24 skipping to change at line 24
* You should have received a copy of the GNU General Public License * You should have received a copy of the GNU General Public License
* along with with program; see the file COPYING. If not, write to the * along with with program; see the file COPYING. If not, write to the
* Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, * Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston,
* MA 02111-1307 USA * MA 02111-1307 USA
*/ */
#ifndef _LALSIMIMR_H #ifndef _LALSIMIMR_H
#define _LALSIMIMR_H #define _LALSIMIMR_H
#include <lal/LALDatatypes.h> #include <lal/LALDatatypes.h>
#include <lal/LALSimInspiralWaveformFlags.h>
#if defined(__cplusplus) #if defined(__cplusplus)
extern "C" { extern "C" {
#elif 0 #elif 0
} /* so that editors will match preceding brace */ } /* so that editors will match preceding brace */
#endif #endif
/** /**
* The number of e-folds of ringdown which should be attached for * The number of e-folds of ringdown which should be attached for
* EOBNR models * EOBNR models
*/ */
#define EOB_RD_EFOLDS 10.0 #define EOB_RD_EFOLDS 10.0
/** /**
* Constant which comes up in some of the EOB models. Its value is * Constant which comes up in some of the EOB models. Its value is
* (94/3 -41/32*pi*pi) * (94/3 -41/32*pi*pi)
*/ */
#define ninty4by3etc 18.687902694437592603 #define ninty4by3etc 18.687902694437592603
/** /**
* Enumerator for choosing the reference frame associated with
* PSpinInspiralRD waveforms.
*/
typedef enum {
View,
TotalJ,
OrbitalL,
} InputAxis;
/**
* Driver routine to compute the non-spinning, inspiral-merger-ringdown * Driver routine to compute the non-spinning, inspiral-merger-ringdown
* phenomenological waveform IMRPhenomA in the frequency domain. * phenomenological waveform IMRPhenomA in the frequency domain.
* *
* Reference: * Reference:
* - Waveform: Eq.(4.13) and (4.16) of http://arxiv.org/pdf/0710.2335 * - Waveform: Eq.(4.13) and (4.16) of http://arxiv.org/pdf/0710.2335
* - Coefficients: Eq.(4.18) of http://arxiv.org/pdf/0710.2335 and * - Coefficients: Eq.(4.18) of http://arxiv.org/pdf/0710.2335 and
* Table I of http://arxiv.org/pdf/0712.0343 * Table I of http://arxiv.org/pdf/0712.0343
* *
* All input parameters should be SI units. * All input parameters should be SI units.
*/ */
int XLALSimIMRPhenomAGenerateFD( int XLALSimIMRPhenomAGenerateFD(
COMPLEX16FrequencySeries **htilde, /**< FD waveform */ COMPLEX16FrequencySeries **htilde, /**< FD waveform */
const REAL8 phi0, /**< initial phase */ const REAL8 phiPeak, /**< orbital phase at peak (rad) */
const REAL8 deltaF, /**< sampling interval */ const REAL8 deltaF, /**< sampling interval (Hz) */
const REAL8 m1_SI, /**< mass of companion 1 (kg) const REAL8 m1_SI, /**< mass of companion 1 (kg) */
*/ const REAL8 m2_SI, /**< mass of companion 2 (kg) */
const REAL8 m2_SI, /**< mass of companion 2 (kg) const REAL8 f_min, /**< starting GW frequency (Hz) */
*/ const REAL8 f_max, /**< end frequency; 0 defaults to ri
const REAL8 f_min, /**< start frequency */ ngdown cutoff freq */
const REAL8 f_max, /**< end frequency */ const REAL8 distance /**< distance of source (m) */
const REAL8 distance /**< distance of source (m) */
); );
/** /**
* Driver routine to compute the non-spinning, inspiral-merger-ringdown * Driver routine to compute the non-spinning, inspiral-merger-ringdown
* phenomenological waveform IMRPhenomA in the time domain. * phenomenological waveform IMRPhenomA in the time domain.
* *
* Reference: * Reference:
* - Waveform: Eq.(4.13) and (4.16) of http://arxiv.org/pdf/0710.2335 * - Waveform: Eq.(4.13) and (4.16) of http://arxiv.org/pdf/0710.2335
* - Coefficients: Eq.(4.18) of http://arxiv.org/pdf/0710.2335 and * - Coefficients: Eq.(4.18) of http://arxiv.org/pdf/0710.2335 and
* Table I of http://arxiv.org/pdf/0712.0343 * Table I of http://arxiv.org/pdf/0712.0343
* *
* All input parameters should be in SI units. Angles should be in radians. * All input parameters should be in SI units. Angles should be in radians.
*/ */
int XLALSimIMRPhenomAGenerateTD( int XLALSimIMRPhenomAGenerateTD(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarization waveform */
REAL8TimeSeries **hcross, /**< x-polarization waveform */ REAL8TimeSeries **hcross, /**< x-polarization waveform */
const REAL8 phiPeak, /**< phase at peak */ const REAL8 phiPeak, /**< orbital phase at peak (rad) */
const REAL8 deltaT, /**< sampling interval */ const REAL8 deltaT, /**< sampling interval (s) */
const REAL8 m1_SI, /**< mass of companion 1 (kg) */ const REAL8 m1_SI, /**< mass of companion 1 (kg) */
const REAL8 m2_SI, /**< mass of companion 2 (kg) */ const REAL8 m2_SI, /**< mass of companion 2 (kg) */
const REAL8 f_min, /**< start frequency */ const REAL8 f_min, /**< starting GW frequency (Hz) */
const REAL8 f_max, /**< end frequency */ const REAL8 f_max, /**< end frequency; 0 defaults to ringdown cu
const REAL8 distance, /**< distance of source (m) */ toff freq */
const REAL8 inclination /**< inclination of source */ const REAL8 distance, /**< distance of source (m) */
const REAL8 inclination /**< inclination of source (rad) */
); );
/** /**
* Compute the dimensionless, spin-aligned parameter chi as used in the * Compute the dimensionless, spin-aligned parameter chi as used in the
* IMRPhenomB waveform. This is different from chi in SpinTaylorRedSpin! * IMRPhenomB waveform. This is different from chi in SpinTaylorRedSpin!
* Reference: http://arxiv.org/pdf/0909.2867, paragraph 3. * Reference: http://arxiv.org/pdf/0909.2867, paragraph 3.
*/ */
double XLALSimIMRPhenomBComputeChi( double XLALSimIMRPhenomBComputeChi(
const REAL8 m1, /**< mass of companion 1 */ const REAL8 m1, /**< mass of companion 1 */
const REAL8 m2, /**< mass of companion 2 */ const REAL8 m2, /**< mass of companion 2 */
skipping to change at line 123 skipping to change at line 114
* phenomenological waveform IMRPhenomB in the frequency domain. * phenomenological waveform IMRPhenomB in the frequency domain.
* *
* Reference: http://arxiv.org/pdf/0909.2867 * Reference: http://arxiv.org/pdf/0909.2867
* - Waveform: Eq.(1) * - Waveform: Eq.(1)
* - Coefficients: Eq.(2) and Table I * - Coefficients: Eq.(2) and Table I
* *
* All input parameters should be in SI units. Angles should be in radians. * All input parameters should be in SI units. Angles should be in radians.
*/ */
int XLALSimIMRPhenomBGenerateFD( int XLALSimIMRPhenomBGenerateFD(
COMPLEX16FrequencySeries **htilde, /**< FD waveform */ COMPLEX16FrequencySeries **htilde, /**< FD waveform */
const REAL8 phi0, /**< initial phase */ const REAL8 phiPeak, /**< orbital phase at peak (rad) */
const REAL8 deltaF, /**< sampling interval */ const REAL8 deltaF, /**< sampling interval (Hz) */
const REAL8 m1_SI, /**< mass of companion 1 (kg) const REAL8 m1_SI, /**< mass of companion 1 (kg) */
*/ const REAL8 m2_SI, /**< mass of companion 2 (kg) */
const REAL8 m2_SI, /**< mass of companion 2 (kg) const REAL8 chi, /**< mass-weighted aligned-spin para
*/ meter */
const REAL8 chi, /**< mass-weighted aligned-spi const REAL8 f_min, /**< starting GW frequency (Hz) */
n parameter */ const REAL8 f_max, /**< end frequency; 0 defaults to ri
const REAL8 f_min, /**< start frequency */ ngdown cutoff freq */
const REAL8 f_max, /**< end frequency */ const REAL8 distance /**< distance of source (m) */
const REAL8 distance /**< distance of source (m) */
); );
/** /**
* Driver routine to compute the spin-aligned, inspiral-merger-ringdown * Driver routine to compute the spin-aligned, inspiral-merger-ringdown
* phenomenological waveform IMRPhenomB in the time domain. * phenomenological waveform IMRPhenomB in the time domain.
* *
* Reference: http://arxiv.org/pdf/0909.2867 * Reference: http://arxiv.org/pdf/0909.2867
* - Waveform: Eq.(1) * - Waveform: Eq.(1)
* - Coefficients: Eq.(2) and Table I * - Coefficients: Eq.(2) and Table I
* *
* All input parameters should be in SI units. Angles should be in radians. * All input parameters should be in SI units. Angles should be in radians.
*/ */
int XLALSimIMRPhenomBGenerateTD( int XLALSimIMRPhenomBGenerateTD(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarization waveform */
REAL8TimeSeries **hcross, /**< x-polarization waveform */ REAL8TimeSeries **hcross, /**< x-polarization waveform */
const REAL8 phiPeak, /**< phase at peak */ const REAL8 phiPeak, /**< orbital phase at peak (rad) */
const REAL8 deltaT, /**< sampling interval */ const REAL8 deltaT, /**< sampling interval (s) */
const REAL8 m1_SI, /**< mass of companion 1 (kg) */ const REAL8 m1_SI, /**< mass of companion 1 (kg) */
const REAL8 m2_SI, /**< mass of companion 2 (kg) */ const REAL8 m2_SI, /**< mass of companion 2 (kg) */
const REAL8 chi, /**< mass-weighted aligned-spin paramet const REAL8 chi, /**< mass-weighted aligned-spin parameter */
er */ const REAL8 f_min, /**< starting GW frequency (Hz) */
const REAL8 f_min, /**< start frequency */ const REAL8 f_max, /**< end frequency; 0 defaults to ringdown cu
const REAL8 f_max, /**< end frequency */ toff freq */
const REAL8 distance, /**< distance of source (m) */ const REAL8 distance, /**< distance of source (m) */
const REAL8 inclination /**< inclination of source */ const REAL8 inclination /**< inclination of source (rad) */
); );
/** /**
* This function generates the plus and cross polarizations for the dominan t * This function generates the plus and cross polarizations for the dominan t
* (2,2) mode of the EOBNRv2 approximant. This model is defined in Pan et a l, * (2,2) mode of the EOBNRv2 approximant. This model is defined in Pan et a l,
* arXiv:1106.1021v1 [gr-qc]. * arXiv:1106.1021v1 [gr-qc].
*/ */
int XLALSimIMREOBNRv2DominantMode( int XLALSimIMREOBNRv2DominantMode(
REAL8TimeSeries **hplus, /**<< The +-polarization waveform (return ed) */ REAL8TimeSeries **hplus, /**<< The +-polarization waveform (return ed) */
REAL8TimeSeries **hcross, /**<< The x-polarization waveform (return ed) */ REAL8TimeSeries **hcross, /**<< The x-polarization waveform (return ed) */
skipping to change at line 246 skipping to change at line 237
REAL8 f_min, /**< start frequency */ REAL8 f_min, /**< start frequency */
REAL8 r, /**< distance of source */ REAL8 r, /**< distance of source */
REAL8 iota, /**< inclination of source (rad) */ REAL8 iota, /**< inclination of source (rad) */
REAL8 s1x, /**< x-component of dimensionless spin for object 1 */ REAL8 s1x, /**< x-component of dimensionless spin for object 1 */
REAL8 s1y, /**< y-component of dimensionless spin for object 1 */ REAL8 s1y, /**< y-component of dimensionless spin for object 1 */
REAL8 s1z, /**< z-component of dimensionless spin for object 1 */ REAL8 s1z, /**< z-component of dimensionless spin for object 1 */
REAL8 s2x, /**< x-component of dimensionless spin for object 2 */ REAL8 s2x, /**< x-component of dimensionless spin for object 2 */
REAL8 s2y, /**< y-component of dimensionless spin for object 2 */ REAL8 s2y, /**< y-component of dimensionless spin for object 2 */
REAL8 s2z, /**< z-component of dimensionless spin for object 2 */ REAL8 s2z, /**< z-component of dimensionless spin for object 2 */
int phaseO, /**< twice post-Newtonian phase order */ int phaseO, /**< twice post-Newtonian phase order */
InputAxis axisChoice, /**< Choice of axis for input spin p arams */ LALSimInspiralFrameAxis axisChoice, /**< Choice of axis for input spin params */
int inspiralOnly /**< 0 generate RD, 1 generate inspiralOnly */ int inspiralOnly /**< 0 generate RD, 1 generate inspiralOnly */
); );
#if 0 #if 0
{ /* so that editors will match succeeding brace */ { /* so that editors will match succeeding brace */
#elif defined(__cplusplus) #elif defined(__cplusplus)
} }
#endif #endif
#endif /* _LALSIMIMR_H */ #endif /* _LALSIMIMR_H */
 End of changes. 7 change blocks. 
49 lines changed or deleted 39 lines changed or added


 LALSimInspiral.h   LALSimInspiral.h 
/* /*
* Copyright (C) 2008 J. Creighton, S. Fairhurst, B. Krishnan, L. Santamari a, E. Ochsner * Copyright (C) 2008 J. Creighton, S. Fairhurst, B. Krishnan, L. Santamari a, E. Ochsner, C. Pankow
* *
* This program is free software; you can redistribute it and/or modify * This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by * it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or * the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version. * (at your option) any later version.
* *
* This program is distributed in the hope that it will be useful, * This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of * but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details. * GNU General Public License for more details.
* *
* You should have received a copy of the GNU General Public License * You should have received a copy of the GNU General Public License
* along with with program; see the file COPYING. If not, write to the * along with with program; see the file COPYING. If not, write to the
* Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, * Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston,
* MA 02111-1307 USA * MA 02111-1307 USA
*/ */
#ifndef _LALSIMINSPIRAL_H #ifndef _LALSIMINSPIRAL_H
#define _LALSIMINSPIRAL_H #define _LALSIMINSPIRAL_H
/* remove SWIG interface directives */
#if !defined(SWIG) && !defined(SWIGLAL_STRUCT)
#define SWIGLAL_STRUCT(...)
#endif
#include <lal/LALDatatypes.h> #include <lal/LALDatatypes.h>
#include <lal/LALSimIMR.h> #include <lal/LALSimIMR.h>
#include <lal/LALSimInspiralWaveformFlags.h>
#include <lal/LALSimInspiralTestGRParams.h>
#include <lal/TimeSeries.h>
#include <gsl/gsl_matrix.h>
#if defined(__cplusplus) #if defined(__cplusplus)
extern "C" { extern "C" {
#elif 0 #elif 0
} /* so that editors will match preceding brace */ } /* so that editors will match preceding brace */
#endif #endif
#define LAL_PN_MODE_L_MAX 3 #define LAL_PN_MODE_L_MAX 3
/** Enum that specifies the PN approximant to be used in computing the wave form. /** Enum that specifies the PN approximant to be used in computing the wave form.
skipping to change at line 91 skipping to change at line 90
EOBNRv2, /**< UNDOCUMENTED */ EOBNRv2, /**< UNDOCUMENTED */
EOBNRv2HM, /**< UNDOCUMENTED */ EOBNRv2HM, /**< UNDOCUMENTED */
SEOBNRv1, /**< Spin-aligned EOBNR model */ SEOBNRv1, /**< Spin-aligned EOBNR model */
IMRPhenomA, /**< Time domain (non-spinning) inspiral-merger-ring down waveforms generated from the inverse FFT of IMRPhenomFA */ IMRPhenomA, /**< Time domain (non-spinning) inspiral-merger-ring down waveforms generated from the inverse FFT of IMRPhenomFA */
IMRPhenomB, /**< Time domain (non-precessing spins) inspiral-mer ger-ringdown waveforms generated from the inverse FFT of IMRPhenomFB */ IMRPhenomB, /**< Time domain (non-precessing spins) inspiral-mer ger-ringdown waveforms generated from the inverse FFT of IMRPhenomFB */
IMRPhenomFA, /**< Frequency domain (non-spinning) inspira l-merger-ringdown templates of Ajith et al [Ajith_2007kx] with phenomenolog ical coefficients defined in the Table I of [Ajith_2007xh]*/ IMRPhenomFA, /**< Frequency domain (non-spinning) inspira l-merger-ringdown templates of Ajith et al [Ajith_2007kx] with phenomenolog ical coefficients defined in the Table I of [Ajith_2007xh]*/
IMRPhenomFB, /**< Frequency domain (non-precessing spins) inspiral-merger-ringdown templates of Ajith et al [Ajith_2009bn] */ IMRPhenomFB, /**< Frequency domain (non-precessing spins) inspiral-merger-ringdown templates of Ajith et al [Ajith_2009bn] */
TaylorEt, /**< UNDOCUMENTED */ TaylorEt, /**< UNDOCUMENTED */
TaylorT4, /**< UNDOCUMENTED */ TaylorT4, /**< UNDOCUMENTED */
TaylorN, /**< UNDOCUMENTED */ TaylorN, /**< UNDOCUMENTED */
NumApproximants /**< UNDOCUMENTED */ NumApproximants /**< Number of elements in enum, useful for checking bounds */
} Approximant; } Approximant;
/** Enumeration to specify which interaction will be used in the waveform /** Enum of possible values to use for post-Newtonian order. */
* generation. Their combination also can be used by the bitwise or.
**/
typedef enum { typedef enum {
LAL_SIM_INSPIRAL_INTERACTION_NONE = 0, /**< No spin, tidal or other LAL_PNORDER_NEWTONIAN, /**< Newtonain (leading) order */
interactions */ LAL_PNORDER_HALF, /**< 0.5PN <==> O(v) */
LAL_SIM_INSPIRAL_INTERACTION_SPIN_ORBIT_15PN = 1, /**< Leading order LAL_PNORDER_ONE, /**< 1PN <==> O(v^2) */
spin-orbit interaction */ LAL_PNORDER_ONE_POINT_FIVE, /**< 1.5PN <==> O(v^3) */
LAL_SIM_INSPIRAL_INTERACTION_SPIN_SPIN_2PN = 1 << 1, /**< Spin-spin LAL_PNORDER_TWO, /**< 2PN <==> O(v^4) */
interaction */ LAL_PNORDER_TWO_POINT_FIVE, /**< 2.5PN <==> O(v^5) */
LAL_SIM_INSPIRAL_INTERACTION_SPIN_SPIN_SELF_2PN = 1 << 2, /**< LAL_PNORDER_THREE, /**< 3PN <==> O(v^6) */
Spin-spin-self interaction */ LAL_PNORDER_THREE_POINT_FIVE, /**< 3.5PN <==> O(v^7) */
LAL_SIM_INSPIRAL_INTERACTION_QUAD_MONO_2PN = 1 << 3, /**< Quadru LAL_PNORDER_PSEUDO_FOUR, /**< pseudo-4PN tuning coefficients included
pole-monopole interaction */ , true 4PN terms currently unknown */
LAL_SIM_INSPIRAL_INTERACTION_SPIN_ORBIT_25PN = 1 << 4, /**< Nex LAL_PNORDER_NUM_ORDER /**< Number of elements in enum, use
t-to-leading-order spin-orbit interaction */ ful for checking bounds */
LAL_SIM_INSPIRAL_INTERACTION_SPIN_ORBIT_3PN = 1 << 5, /**< Spin-spi } LALPNOrder;
n interaction */
LAL_SIM_INSPIRAL_INTERACTION_TIDAL_5PN = 1 << 6, /**< Leading-order
tidal interaction */
LAL_SIM_INSPIRAL_INTERACTION_TIDAL_6PN = 1 << 7, /**< Next-to-leadin
g-order tidal interaction */
LAL_SIM_INSPIRAL_INTERACTION_ALL_SPIN = (1 << 6) - 1, /**< all spin
interactions, no tidal interactions */
LAL_SIM_INSPIRAL_INTERACTION_ALL = (1 << 8) - 1 /**< all spin and ti
dal interactions */
} LALSimInspiralInteraction;
/** Enumeration to specify the tapering method to apply to the waveform */ /** Enumeration to specify the tapering method to apply to the waveform */
typedef enum typedef enum
{ {
LAL_SIM_INSPIRAL_TAPER_NONE, /**< No tapering */ LAL_SIM_INSPIRAL_TAPER_NONE, /**< No tapering */
LAL_SIM_INSPIRAL_TAPER_START, /**< Taper the start of the waveform */ LAL_SIM_INSPIRAL_TAPER_START, /**< Taper the start of the waveform */
LAL_SIM_INSPIRAL_TAPER_END, /**< Taper the end of the waveform * / LAL_SIM_INSPIRAL_TAPER_END, /**< Taper the end of the waveform * /
LAL_SIM_INSPIRAL_TAPER_STARTEND, /**< Taper the start and the end of the waveform */ LAL_SIM_INSPIRAL_TAPER_STARTEND, /**< Taper the start and the end of the waveform */
LAL_SIM_INSPIRAL_TAPER_NUM_OPTS /**< UNDOCUMENTED */ LAL_SIM_INSPIRAL_TAPER_NUM_OPTS /**< Number of elements in enum, use ful for checking bounds */
} LALSimInspiralApplyTaper; } LALSimInspiralApplyTaper;
/** /**
* Tapers a REAL4 inspiral waveform in the time domain. * Tapers a REAL4 inspiral waveform in the time domain.
*/ */
int XLALSimInspiralREAL4WaveTaper( int XLALSimInspiralREAL4WaveTaper(
REAL4Vector *signalvec, /**< pointer to wave form vector */ REAL4Vector *signalvec, /**< pointer to wave form vector */
LALSimInspiralApplyTaper bookends /**< taper type enum erator */ LALSimInspiralApplyTaper bookends /**< taper type enum erator */
); );
/** /**
* Tapers a REAL8 inspiral waveform in the time domain. * Tapers a REAL8 inspiral waveform in the time domain.
*/ */
int XLALSimInspiralREAL8WaveTaper( int XLALSimInspiralREAL8WaveTaper(
REAL8Vector *signalvec, /**< pointer to wave form vector */ REAL8Vector *signalvec, /**< pointer to wave form vector */
LALSimInspiralApplyTaper bookends /**< taper type enum erator */ LALSimInspiralApplyTaper bookends /**< taper type enum erator */
); );
/*
* Structyure to carry a collectio of spherical harmonic modes in COMPLEX16
* time series. Contains convenience getter and setter functions, as well a
s
* a convienence "maximum l mode" function. Implemented as a singly forward
* linked list.
*/
typedef struct tagSphHarmTimeSeries SphHarmTimeSeries;
/*
* Create a SphHarmTimeSeries. If appended is not NULL, this will prepend a
new
* structure to the list by duplicating the mode inmode, mode numbers l, an
d m,
* and then set the next pointer to the appended structure.
*/
SphHarmTimeSeries* XLALSphHarmTimeSeriesAddMode(
SphHarmTimeSeries *appended, /**< List structure to prepend
to */
const COMPLEX16TimeSeries* inmode, /**< mode series to cont
ain */
UINT4 l, /**< major mode number */
INT4 m /**< minor mode number */
);
/*
* Destroy a SphHarmTimeSeries. Note that this will destroy any
* COMPLEX16TimeSeries which it has references to.
*/
void XLALDestroySphHarmTimeSeries( SphHarmTimeSeries* ts );
/*
* Destroy a SphHarmTimeSeries. Note that this will destroy any
* COMPLEX16TimeSeries which it has references to.
*/
UINT4 XLALSphHarmTimeSeriesGetMaxL( SphHarmTimeSeries* ts );
/*
* Get the mode-decomposed time series corresponding to l,m.
*/
COMPLEX16TimeSeries* XLALSphHarmTimeSeriesGetMode(
SphHarmTimeSeries *ts,
UINT4 l,
INT4 m
);
/** /**
* Computes h(2,2) mode of spherical harmonic decomposition of * Computes h(2,2) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform. * the post-Newtonian inspiral waveform.
* *
* Implements Equation (79) of: * Implements Equation (79) of:
* Lawrence E. Kidder, \"Using Full Information When Computing Modes of * Lawrence E. Kidder, \"Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular * Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit\", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc]. * Orbit\", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/ */
COMPLEX16 XLALSimInspiralPNMode22( COMPLEX16TimeSeries *XLALSimInspiralPNMode22(
REAL8 v, /**< post-Newtonian parameter */ REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8 phi, /**< orbital phase */ REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter */ REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 */ REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 */ REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source */ REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */ int O /**< twice post-Newtonian order */
); );
/** /**
* Computes h(2,1) mode of spherical harmonic decomposition of * Computes h(2,1) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform. * the post-Newtonian inspiral waveform.
* *
* Implements Equation (80) of: * Implements Equation (80) of:
* Lawrence E. Kidder, \"Using Full Information When Computing Modes of * Lawrence E. Kidder, \"Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular * Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit\", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc]. * Orbit\", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/ */
COMPLEX16 XLALSimInspiralPNMode21( COMPLEX16TimeSeries *XLALSimInspiralPNMode21(
REAL8 v, /**< post-Newtonian parameter */ REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8 phi, /**< orbital phase */ REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter */ REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 */ REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 */ REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source */ REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */ int O /**< twice post-Newtonian order */
); );
/**
* Computes h(2,0) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* Implements Equation (81) of:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode20(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/** /**
* Computes h(3,3) mode of spherical harmonic decomposition of * Computes h(3,3) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform. * the post-Newtonian inspiral waveform.
* *
* Implements Equation (82) of: * Implements Equation (82) of:
* Lawrence E. Kidder, \"Using Full Information When Computing Modes of * Lawrence E. Kidder, \"Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular * Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit\", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc]. * Orbit\", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/ */
COMPLEX16 XLALSimInspiralPNMode33( COMPLEX16TimeSeries *XLALSimInspiralPNMode33(
REAL8 v, /**< post-Newtonian parameter */ REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8 phi, /**< orbital phase */ REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter */ REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 */ REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 */ REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source */ REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */ int O /**< twice post-Newtonian order */
); );
/** /**
* Computes h(3,2) mode of spherical harmonic decomposition of * Computes h(3,2) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform. * the post-Newtonian inspiral waveform.
* *
* Implements Equation (83) of: * Implements Equation (83) of:
* Lawrence E. Kidder, \"Using Full Information When Computing Modes of * Lawrence E. Kidder, \"Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular * Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit\", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc]. * Orbit\", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/ */
COMPLEX16 XLALSimInspiralPNMode32( COMPLEX16TimeSeries *XLALSimInspiralPNMode32(
REAL8 v, /**< post-Newtonian parameter */ REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8 phi, /**< orbital phase */ REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter */ REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 */ REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 */ REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source */ REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */ int O /**< twice post-Newtonian order */
); );
/** /**
* Computes h(3,1) mode of spherical harmonic decomposition of * Computes h(3,1) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform. * the post-Newtonian inspiral waveform.
* *
* Implements Equation (84) of: * Implements Equation (84) of:
* Lawrence E. Kidder, \"Using Full Information When Computing Modes of * Lawrence E. Kidder, \"Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular * Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit\", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc]. * Orbit\", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/ */
COMPLEX16 XLALSimInspiralPNMode31( COMPLEX16TimeSeries *XLALSimInspiralPNMode31(
REAL8 v, /**< post-Newtonian parameter */ REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8 phi, /**< orbital phase */ REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter */ REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 */ REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 */ REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source */ REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */ int O /**< twice post-Newtonian order */
); );
/**
* Computes h(3,0) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* Implements Equation (85) of:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode30(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/**
* Computes h(4,4) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* Implements Equation (86) of:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode44(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/**
* Computes h(4,3) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* Implements Equation (87) of:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode43(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/**
* Computes h(4,2) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* Implements Equation (88) of:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode42(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/**
* Computes h(4,1) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* Implements Equation (89) of:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode41(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/**
* Computes h(4,0) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* Implements Equation (90) of:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode40(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/**
* Computes h(5,5) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* Implements Equation (91) of:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode55(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/**
* Computes h(5,4) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* Implements Equation (92) of:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode54(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/**
* Computes h(5,3) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* Implements Equation (93) of:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode53(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/**
* Computes h(5,2) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* Implements Equation (94) of:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode52(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/**
* Computes h(5,1) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* Implements Equation (95) of:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode51(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/**
* Computes h(5,0) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* THIS MODE IS ZERO TO THE ORDER CONSIDERED IN:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode50(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/**
* Computes h(6,6) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* Implements Equation (96) of:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode66(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/**
* Computes h(6,5) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* Implements Equation (97) of:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode65(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/**
* Computes h(6,4) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* Implements Equation (98) of:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode64(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/**
* Computes h(6,3) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* Implements Equation (99) of:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode63(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/**
* Computes h(6,2) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* Implements Equation (100) of:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode62(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/**
* Computes h(6,1) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* Implements Equation (101) of:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode61(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/**
* Computes h(6,0) mode of spherical harmonic decomposition of
* the post-Newtonian inspiral waveform.
*
* THIS MODE IS ZERO TO THE ORDER CONSIDERED IN:
* Lawrence E. Kidder, "Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/
COMPLEX16TimeSeries *XLALSimInspiralPNMode60(
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8 v0, /**< tail gauge parameter (default=1) */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */
int O /**< twice post-Newtonian order */
);
/** /**
* Multiplies a mode h(l,m) by a spin-2 weighted spherical harmonic * Multiplies a mode h(l,m) by a spin-2 weighted spherical harmonic
* to obtain hplus - i hcross, which is added to the time series. * to obtain hplus - i hcross, which is added to the time series.
* *
* Implements the sum of a single term of Eq. (11) of: * Implements the sum of a single term of Eq. (11) of:
* Lawrence E. Kidder, \"Using Full Information When Computing Modes of * Lawrence E. Kidder, \"Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular * Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit\", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc]. * Orbit\", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
* *
skipping to change at line 267 skipping to change at line 684
* the post-Newtonian inspiral waveform. * the post-Newtonian inspiral waveform.
* *
* See Eqns. (79)-(116) of: * See Eqns. (79)-(116) of:
* Lawrence E. Kidder, \"Using Full Information When Computing Modes of * Lawrence E. Kidder, \"Using Full Information When Computing Modes of
* Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular * Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular
* Orbit\", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc]. * Orbit\", Physical Review D 77, 044016 (2008), arXiv:0710.0614v1 [gr-qc].
*/ */
COMPLEX16TimeSeries *XLALCreateSimInspiralPNModeCOMPLEX16TimeSeries( COMPLEX16TimeSeries *XLALCreateSimInspiralPNModeCOMPLEX16TimeSeries(
REAL8TimeSeries *v, /**< post-Newtonian parameter */ REAL8TimeSeries *v, /**< post-Newtonian parameter */
REAL8TimeSeries *phi, /**< orbital phase */ REAL8TimeSeries *phi, /**< orbital phase */
REAL8 v0, /**< tail-term gauge choice (default = REAL8 v0, /**< tail gauge parameter (default = 1
1) */ ) */
REAL8 m1, /**< mass of companion 1 */ REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 */ REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source */ REAL8 r, /**< distance of source (m) */
int O, /**< twice post-Newtonain order */ int O, /**< twice post-Newtonain order */
int l, /**< mode number l */ int l, /**< mode number l */
int m /**< mode number m */ int m /**< mode number m */
); );
/** /**
* Given time series for a binary's orbital dynamical variables, * Given time series for a binary's orbital dynamical variables,
* construct the waveform polarizations h+ and hx as a sum of * construct the waveform polarizations h+ and hx as a sum of
* -2 spin-weighted spherical harmonic modes, h_lm. * -2 spin-weighted spherical harmonic modes, h_lm.
* NB: Valid only for non-precessing systems! * NB: Valid only for non-precessing systems!
skipping to change at line 374 skipping to change at line 791
REAL8TimeSeries *E1y, /**< orbital plane basis vector y comp. */ REAL8TimeSeries *E1y, /**< orbital plane basis vector y comp. */
REAL8TimeSeries *E1z, /**< orbital plane basis vector z comp. */ REAL8TimeSeries *E1z, /**< orbital plane basis vector z comp. */
REAL8 m1, /**< mass of companion 1 (kg) */ REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */ REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r, /**< distance of source (m) */ REAL8 r, /**< distance of source (m) */
REAL8 v0, /**< tail-term gauge choice (default = 1) */ REAL8 v0, /**< tail-term gauge choice (default = 1) */
INT4 ampO /**< twice amp. post-Newtonian order */ INT4 ampO /**< twice amp. post-Newtonian order */
); );
/** /**
* Compute the physical template family "Q" vectors for a spinning, precess
ing
* binary when provided time series of all the dynamical quantities.
* These vectors always supplied to dominant order.
*
* Based on Pan, Buonanno, Chan and Vallisneri PRD69 104017, (see also thes
es
* of Diego Fazi and Ian Harry)
*
* NOTE: The vectors MUST be given in the so-called radiation frame where
* Z is the direction of propagation, X is the principal '+' axis and Y = Z
x X
*/
int XLALSimInspiralPrecessingPTFQWaveforms(
REAL8TimeSeries **Q1, /**< PTF-Q1 waveform [returned] */
REAL8TimeSeries **Q2, /**< PTF-Q2 waveform [returned] */
REAL8TimeSeries **Q3, /**< PTF-Q2 waveform [returned] */
REAL8TimeSeries **Q4, /**< PTF-Q2 waveform [returned] */
REAL8TimeSeries **Q5, /**< PTF-Q2 waveform [returned] */
REAL8TimeSeries *V, /**< post-Newtonian parameter */
REAL8TimeSeries *Phi, /**< orbital phase */
REAL8TimeSeries *S1x, /**< Spin1 vector x component */
REAL8TimeSeries *S1y, /**< Spin1 vector y component */
REAL8TimeSeries *S1z, /**< Spin1 vector z component */
REAL8TimeSeries *S2x, /**< Spin2 vector x component */
REAL8TimeSeries *S2y, /**< Spin2 vector y component */
REAL8TimeSeries *S2z, /**< Spin2 vector z component */
REAL8TimeSeries *LNhatx, /**< unit orbital ang. mom. x comp. */
REAL8TimeSeries *LNhaty, /**< unit orbital ang. mom. y comp. */
REAL8TimeSeries *LNhatz, /**< unit orbital ang. mom. z comp. */
REAL8TimeSeries *E1x, /**< orbital plane basis vector x comp. */
REAL8TimeSeries *E1y, /**< orbital plane basis vector y comp. */
REAL8TimeSeries *E1z, /**< orbital plane basis vector z comp. */
REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 r /**< distance of source (m) */
);
/**
* Compute the length of an inspiral waveform assuming the Taylor dEnergy a nd Flux equations * Compute the length of an inspiral waveform assuming the Taylor dEnergy a nd Flux equations
*/ */
REAL8 REAL8
XLALSimInspiralTaylorLength( XLALSimInspiralTaylorLength(
REAL8 deltaT, /**< sampling interval */ REAL8 deltaT, /**< sampling interval */
REAL8 m1, /**< mass of companion 1 */ REAL8 m1, /**< mass of companion 1 */
REAL8 m2, /**< mass of companion 2 */ REAL8 m2, /**< mass of companion 2 */
REAL8 f_min, /**< start frequency */ REAL8 f_min, /**< start frequency */
int O /**< twice post-Newtonian order */ int O /**< twice post-Newtonian order */
); );
skipping to change at line 406 skipping to change at line 860
/** /**
* Checks whether the given approximant is implemented in lalsimulation's X LALSimInspiralChooseFDWaveform(). * Checks whether the given approximant is implemented in lalsimulation's X LALSimInspiralChooseFDWaveform().
* *
* returns 1 if the approximant is implemented, 0 otherwise. * returns 1 if the approximant is implemented, 0 otherwise.
*/ */
int XLALSimInspiralImplementedFDApproximants( int XLALSimInspiralImplementedFDApproximants(
Approximant approximant /**< post-Newtonian approximant for use in wave form production */ Approximant approximant /**< post-Newtonian approximant for use in wave form production */
); );
/** /**
* XLAL function to determine approximant from a string. The string need n
ot
* match exactly, only contain a member of the Approximant enum.
*/
int XLALGetApproximantFromString(const CHAR *inString);
/**
* XLAL function to determine string from approximant enum.
* This function needs to be updated when new approximants are added.
*/
char* XLALGetStringFromApproximant(Approximant approximant);
/**
* XLAL function to determine PN order from a string. The string need not
* match exactly, only contain a member of the LALPNOrder enum.
*/
int XLALGetOrderFromString(const CHAR *inString);
/**
* XLAL function to determine tapering flag from a string. The string must
* match exactly with a member of the LALSimInspiralApplyTaper enum.
*/
int XLALGetTaperFromString(const CHAR *inString);
/**
* XLAL function to determine LALSimInspiralInteraction from a string.
*
* TODO: return the bit sum if the string is a concatenation of several
* interaction terms. Also make names match cases of enum.
*/
int XLALGetInteractionFromString(const CHAR *inString);
/** XLAL function to determine axis choice flag from a string */
int XLALGetFrameAxisFromString(const CHAR *inString);
/**
* XLAL function to determine adaptive integration flag from a string. Ret
urns
* 1 if string contains 'fixedStep', otherwise returns 0 to signal
* adaptive integration should be used.
*/
int XLALGetAdaptiveIntFromString(const CHAR *inString);
/**
* XLAL function to determine inspiral-only flag from a string. Returns
* 1 if string contains 'inspiralOnly', otherwise returns 0 to signal
* full inspiral-merger-ringdown waveform should be generated.
*/
int XLALGetInspiralOnlyFromString(const CHAR *inString);
/**
* DEPRECATED: USE XLALSimInspiralChooseTDWaveform() INSTEAD * DEPRECATED: USE XLALSimInspiralChooseTDWaveform() INSTEAD
* *
* Chooses between different approximants when requesting a waveform to be generated * Chooses between different approximants when requesting a waveform to be generated
* For spinning waveforms, all known spin effects up to given PN order are included * For spinning waveforms, all known spin effects up to given PN order are included
* *
* The parameters passed must be in SI units. * The parameters passed must be in SI units.
*/ */
int XLALSimInspiralChooseWaveform( int XLALSimInspiralChooseWaveform(
REAL8TimeSeries **hplus, /**< +-polarization wavefor m */ REAL8TimeSeries **hplus, /**< +-polarization wavefor m */
REAL8TimeSeries **hcross, /**< x-polarization wavefor m */ REAL8TimeSeries **hcross, /**< x-polarization wavefor m */
REAL8 phi0, /**< peak phase */ REAL8 phiRef, /**< reference orbital phas
REAL8 deltaT, /**< sampling interval */ e (rad) */
REAL8 m1, /**< mass of companion 1 */ REAL8 deltaT, /**< sampling interval (s)
REAL8 m2, /**< mass of companion 2 */ */
REAL8 m1, /**< mass of companion 1 (k
g) */
REAL8 m2, /**< mass of companion 2 (k
g) */
REAL8 S1x, /**< x-component of the dim ensionless spin of object 1 */ REAL8 S1x, /**< x-component of the dim ensionless spin of object 1 */
REAL8 S1y, /**< y-component of the dim ensionless spin of object 1 */ REAL8 S1y, /**< y-component of the dim ensionless spin of object 1 */
REAL8 S1z, /**< z-component of the dim ensionless spin of object 1 */ REAL8 S1z, /**< z-component of the dim ensionless spin of object 1 */
REAL8 S2x, /**< x-component of the dim ensionless spin of object 2 */ REAL8 S2x, /**< x-component of the dim ensionless spin of object 2 */
REAL8 S2y, /**< y-component of the dim ensionless spin of object 2 */ REAL8 S2y, /**< y-component of the dim ensionless spin of object 2 */
REAL8 S2z, /**< z-component of the dim ensionless spin of object 2 */ REAL8 S2z, /**< z-component of the dim ensionless spin of object 2 */
REAL8 f_min, /**< start frequency */ REAL8 f_min, /**< starting GW frequency
REAL8 r, /**< distance of source */ (Hz) */
REAL8 f_ref, /**< reference GW frequency
(Hz) */
REAL8 r, /**< distance of source (m)
*/
REAL8 i, /**< inclination of source (rad) */ REAL8 i, /**< inclination of source (rad) */
REAL8 lambda1, /**< (tidal deformability o REAL8 lambda1, /**< (tidal deformability o
f mass 1) / (total mass)^5 (dimensionless) */ f mass 1) / m1^5 (dimensionless) */
REAL8 lambda2, /**< (tidal deformability o REAL8 lambda2, /**< (tidal deformability o
f mass 2) / (total mass)^5 (dimensionless) */ f mass 2) / m2^5 (dimensionless) */
LALSimInspiralInteraction interactionFlags, /**< flag to control spin a LALSimInspiralWaveformFlags *waveFlags, /**< Set of flags to contro
nd tidal effects */ l special behavior of some waveform families. Pass in NULL (or None in pyth
on) for default flags */
LALSimInspiralTestGRParam *nonGRparams, /**< Linked list of non-GR p
arameters. Pass in NULL (or None in python) for standard GR waveforms */
int amplitudeO, /**< twice post-Newtonian a mplitude order */ int amplitudeO, /**< twice post-Newtonian a mplitude order */
int phaseO, /**< twice post-Newtonian o rder */ int phaseO, /**< twice post-Newtonian o rder */
Approximant approximant /**< post-Newtonian approxi mant to use for waveform production */ Approximant approximant /**< post-Newtonian approxi mant to use for waveform production */
); );
/** /**
* Chooses between different approximants when requesting a waveform to be generated * Chooses between different approximants when requesting a waveform to be generated
* For spinning waveforms, all known spin effects up to given PN order are included * For spinning waveforms, all known spin effects up to given PN order are included
* Returns the waveform in the time domain. * Returns the waveform in the time domain.
* *
* The parameters passed must be in SI units. * The parameters passed must be in SI units.
*/ */
int XLALSimInspiralChooseTDWaveform( int XLALSimInspiralChooseTDWaveform(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarization waveform */
REAL8TimeSeries **hcross, /**< x-polarization waveform */ REAL8TimeSeries **hcross, /**< x-polarization waveform */
REAL8 phi0, /**< reference phase */ REAL8 phiRef, /**< reference orbital phase (rad) */
REAL8 deltaT, /**< sampling interval */ REAL8 deltaT, /**< sampling interval (s) */
REAL8 m1, /**< mass of companion 1 */ REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 */ REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 s1x, /**< x-component of the dimensionless spin of object 1 */ REAL8 s1x, /**< x-component of the dimensionless spin of object 1 */
REAL8 s1y, /**< y-component of the dimensionless spin of object 1 */ REAL8 s1y, /**< y-component of the dimensionless spin of object 1 */
REAL8 s1z, /**< z-component of the dimensionless spin of object 1 */ REAL8 s1z, /**< z-component of the dimensionless spin of object 1 */
REAL8 s2x, /**< x-component of the dimensionless spin of object 2 */ REAL8 s2x, /**< x-component of the dimensionless spin of object 2 */
REAL8 s2y, /**< y-component of the dimensionless spin of object 2 */ REAL8 s2y, /**< y-component of the dimensionless spin of object 2 */
REAL8 s2z, /**< z-component of the dimensionless spin of object 2 */ REAL8 s2z, /**< z-component of the dimensionless spin of object 2 */
REAL8 f_min, /**< start frequency */ REAL8 f_min, /**< starting GW frequency (Hz) */
REAL8 r, /**< distance of source */ REAL8 f_ref, /**< reference GW frequency (Hz) */
REAL8 r, /**< distance of source (m) */
REAL8 i, /**< inclination of source (rad) */ REAL8 i, /**< inclination of source (rad) */
REAL8 lambda1, /**< (tidal deformability of mass 1) / (tot REAL8 lambda1, /**< (tidal deformability of mass 1) / m1^5
al mass)^5 (dimensionless) */ (dimensionless) */
REAL8 lambda2, /**< (tidal deformability of mass 2) / (tot REAL8 lambda2, /**< (tidal deformability of mass 2) / m2^5
al mass)^5 (dimensionless) */ (dimensionless) */
LALSimInspiralInteraction interactionFlags, /**< flag to control spin a LALSimInspiralWaveformFlags *waveFlags, /**< Set of flags to control sp
nd tidal effects */ ecial behavior of some waveform families. Pass in NULL (or None in python)
for default flags */
LALSimInspiralTestGRParam *nonGRparams, /**< Linked list of non-GR para
meters. Pass in NULL (or None in python) for standard GR waveforms */
int amplitudeO, /**< twice post-Newtonian amplitude order * / int amplitudeO, /**< twice post-Newtonian amplitude order * /
int phaseO, /**< twice post-Newtonian phase order */ int phaseO, /**< twice post-Newtonian phase order */
Approximant approximant /**< post-Newtonian approximant to use for waveform production */ Approximant approximant /**< post-Newtonian approximant to use for waveform production */
); );
/** /**
* Chooses between different approximants when requesting a waveform to be generated * Chooses between different approximants when requesting a waveform to be generated
* For spinning waveforms, all known spin effects up to given PN order are included * For spinning waveforms, all known spin effects up to given PN order are included
* Returns the waveform in the frequency domain. * Returns the waveform in the frequency domain.
* *
* The parameters passed must be in SI units. * The parameters passed must be in SI units.
*/ */
int XLALSimInspiralChooseFDWaveform( int XLALSimInspiralChooseFDWaveform(
COMPLEX16FrequencySeries **htilde, /**< FD waveform */ COMPLEX16FrequencySeries **htilde, /**< FD waveform */
REAL8 phi0, /**< peak phase */ REAL8 phiRef, /**< reference orbital phas
REAL8 deltaF, /**< sampling interval */ e (rad) */
REAL8 m1, /**< mass of companion 1 */ REAL8 deltaF, /**< sampling interval (Hz)
REAL8 m2, /**< mass of companion 2 */ */
REAL8 m1, /**< mass of companion 1 (k
g) */
REAL8 m2, /**< mass of companion 2 (k
g) */
REAL8 S1x, /**< x-component of the dim ensionless spin of object 1 */ REAL8 S1x, /**< x-component of the dim ensionless spin of object 1 */
REAL8 S1y, /**< y-component of the dim ensionless spin of object 1 */ REAL8 S1y, /**< y-component of the dim ensionless spin of object 1 */
REAL8 S1z, /**< z-component of the dim ensionless spin of object 1 */ REAL8 S1z, /**< z-component of the dim ensionless spin of object 1 */
REAL8 S2x, /**< x-component of the dim ensionless spin of object 2 */ REAL8 S2x, /**< x-component of the dim ensionless spin of object 2 */
REAL8 S2y, /**< y-component of the dim ensionless spin of object 2 */ REAL8 S2y, /**< y-component of the dim ensionless spin of object 2 */
REAL8 S2z, /**< z-component of the dim ensionless spin of object 2 */ REAL8 S2z, /**< z-component of the dim ensionless spin of object 2 */
REAL8 f_min, /**< start frequency */ REAL8 f_min, /**< starting GW frequency
REAL8 f_max, /**< end frequency */ (Hz) */
REAL8 r, /**< distance of source */ REAL8 f_max, /**< ending GW frequency (H
z) */
REAL8 r, /**< distance of source (m)
*/
REAL8 i, /**< inclination of source (rad) */ REAL8 i, /**< inclination of source (rad) */
REAL8 lambda1, /**< (tidal deformability o REAL8 lambda1, /**< (tidal deformability o
f mass 1) / (total mass)^5 (dimensionless) */ f mass 1) / m1^5 (dimensionless) */
REAL8 lambda2, /**< (tidal deformability o REAL8 lambda2, /**< (tidal deformability o
f mass 2) / (total mass)^5 (dimensionless) */ f mass 2) / m2^5 (dimensionless) */
LALSimInspiralInteraction interactionFlags, /**< flag to control spin a LALSimInspiralWaveformFlags *waveFlags, /**< Set of flags to contro
nd tidal effects */ l special behavior of some waveform families. Pass in NULL (or None in pyth
on) for default flags */
LALSimInspiralTestGRParam *nonGRparams, /**< Linked list of non-GR p
arameters. Pass in NULL (or None in python) for standard GR waveforms */
int amplitudeO, /**< twice post-Newtonian a mplitude order */ int amplitudeO, /**< twice post-Newtonian a mplitude order */
int phaseO, /**< twice post-Newtonian o rder */ int phaseO, /**< twice post-Newtonian o rder */
Approximant approximant /**< post-Newtonian approxi mant to use for waveform production */ Approximant approximant /**< post-Newtonian approxi mant to use for waveform production */
); );
/**
* Interface to compute -2 spin-weighted spherical harmonic modes for a bin
ary
* inspiral of any available amplitude and phase PN order.
* The phasing is computed with any of the TaylorT1, T2, T3, T4 methods.
*
* FIXME: Interface will be changed to return a collection of modes.
*/
COMPLEX16TimeSeries *XLALSimInspiralChooseTDModes(
REAL8 phiRef, /**< reference orbital phas
e (rad) */
REAL8 deltaT, /**< sampling interval (s)
*/
REAL8 m1, /**< mass of companion 1 (k
g) */
REAL8 m2, /**< mass of companion 2 (k
g) */
REAL8 f_min, /**< starting GW frequency
(Hz) */
REAL8 f_ref, /**< reference GW frequency
(Hz) */
REAL8 r, /**< distance of source (m)
*/
REAL8 lambda1, /**< (tidal deformability o
f mass 1) / m1^5 (dimensionless) */
REAL8 lambda2, /**< (tidal deformability o
f mass 2) / m2^5 (dimensionless) */
LALSimInspiralWaveformFlags *waveFlags, /**< Set of flags to contro
l special behavior of some waveform families. Pass in NULL (or None in pyth
on) for default flags */
LALSimInspiralTestGRParam *nonGRparams, /**< Linked list of non-GR p
arameters. Pass in NULL (or None in python) for standard GR waveforms */
int amplitudeO, /**< twice post-Newtonian a
mplitude order */
int phaseO, /**< twice post-Newtonian o
rder */
int l, /**< l index of mode - repl
ace with a struct of several integer pairs */
int m, /**< m index of mode - ditt
o */
Approximant approximant /**< post-Newtonian approxi
mant to use for waveform production */
);
/* TaylorT4 functions */ /* TaylorT4 functions */
/** /**
* Evolves a post-Newtonian orbit using the Taylor T4 method. * Evolves a post-Newtonian orbit using the Taylor T4 method.
* *
* See: * See:
* Michael Boyle, Duncan A. Brown, Lawrence E. Kidder, Abdul H. Mroue, * Michael Boyle, Duncan A. Brown, Lawrence E. Kidder, Abdul H. Mroue,
* Harald P. Pfeiffer, Mark A. Scheel, Gregory B. Cook, and Saul A. Teukolsk y * Harald P. Pfeiffer, Mark A. Scheel, Gregory B. Cook, and Saul A. Teukolsk y
* \"High-accuracy comparison of numerical relativity simulations with * \"High-accuracy comparison of numerical relativity simulations with
* post-Newtonian expansions\" * post-Newtonian expansions\"
* <a href="http://arxiv.org/abs/0710.0158v2">arXiv:0710.0158v2</a>. * <a href="http://arxiv.org/abs/0710.0158v2">arXiv:0710.0158v2</a>.
*/ */
int XLALSimInspiralTaylorT4PNEvolveOrbit( int XLALSimInspiralTaylorT4PNEvolveOrbit(
REAL8TimeSeries **V, /**< post-Newtonian parameter [return REAL8TimeSeries **v, /**< post-Newton
ed] */ ian parameter [returned] */
REAL8TimeSeries **phi, /**< orbital phase [returned] */ REAL8TimeSeries **phi, /**< orbital pha
REAL8 phic, /**< coalescence phase */ se [returned] */
REAL8 deltaT, /**< sampling interval */ REAL8 phiRef, /**< reference o
REAL8 m1, /**< mass of companion 1 */ rbital phase (rad) */
REAL8 m2, /**< mass of companion 2 */ REAL8 deltaT, /**< sampling in
REAL8 f_min, /**< start frequency */ terval (s) */
int O /**< twice post-Newtonian order */ REAL8 m1, /**< mass of com
panion 1 (kg) */
REAL8 m2, /**< mass of com
panion 2 (kg) */
REAL8 f_min, /**< start frequ
ency (Hz) */
REAL8 fRef, /**< reference f
requency (Hz) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int O /**< twice post-
Newtonian order */
); );
/** /**
* Driver routine to compute the post-Newtonian inspiral waveform. * Driver routine to compute the post-Newtonian inspiral waveform.
* *
* This routine allows the user to specify different pN orders * This routine allows the user to specify different pN orders
* for phasing calcuation vs. amplitude calculations. * for phasing calcuation vs. amplitude calculations.
*/ */
int XLALSimInspiralTaylorT4PNGenerator( int XLALSimInspiralTaylorT4PNGenerator(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarizat
REAL8TimeSeries **hcross, /**< x-polarization waveform */ ion waveform */
REAL8 phic, /**< coalescence phase */ REAL8TimeSeries **hcross, /**< x-polarizat
REAL8 v0, /**< tail-term gauge choice (defau ion waveform */
lt = 1) */ REAL8 phiRef, /**< reference o
REAL8 deltaT, /**< sampling interval */ rbital phase (rad) */
REAL8 m1, /**< mass of companion 1 */ REAL8 v0, /**< tail-term g
REAL8 m2, /**< mass of companion 2 */ auge choice (default = 1) */
REAL8 f_min, /**< start frequency */ REAL8 deltaT, /**< sampling in
REAL8 r, /**< distance of source */ terval (s) */
REAL8 i, /**< inclination of source (rad) * REAL8 m1, /**< mass of com
/ panion 1 (kg) */
int amplitudeO, /**< twice post-Newtonian amplitud REAL8 m2, /**< mass of com
e order */ panion 2 (kg) */
int phaseO /**< twice post-Newtonian phase or REAL8 f_min, /**< starting GW
der */ frequency (Hz) */
REAL8 fRef, /**< reference G
W frequency (Hz) */
REAL8 r, /**< distance of
source (m) */
REAL8 i, /**< inclination
of source (rad) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int amplitudeO, /**< twice post-
Newtonian amplitude order */
int phaseO /**< twice post-
Newtonian phase order */
);
/**
* Driver routine to compute the -2 spin-weighted spherical harmonic mode
* using TaylorT4 phasing.
*/
COMPLEX16TimeSeries *XLALSimInspiralTaylorT4PNModes(
REAL8 phiRef, /**< reference o
rbital phase (rad) */
REAL8 v0, /**< tail-term g
auge choice (default = 1) */
REAL8 deltaT, /**< sampling in
terval (s) */
REAL8 m1, /**< mass of com
panion 1 (kg) */
REAL8 m2, /**< mass of com
panion 2 (kg) */
REAL8 f_min, /**< starting GW
frequency (Hz) */
REAL8 fRef, /**< reference G
W frequency (Hz) */
REAL8 r, /**< distance of
source (m) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int amplitudeO, /**< twice post-
Newtonian amplitude order */
int phaseO, /**< twice post-
Newtonian phase order */
int l, /**< l index of
mode */
int m /**< m index of
mode */
); );
/** /**
* Driver routine to compute the post-Newtonian inspiral waveform. * Driver routine to compute the post-Newtonian inspiral waveform.
* *
* This routine uses the same pN order for phasing and amplitude * This routine uses the same pN order for phasing and amplitude
* (unless the order is -1 in which case the highest available * (unless the order is -1 in which case the highest available
* order is used for both of these -- which might not be the same). * order is used for both of these -- which might not be the same).
* *
* Constant log term in amplitude set to 1. This is a gauge choice. * Constant log term in amplitude set to 1. This is a gauge choice.
*/ */
int XLALSimInspiralTaylorT4PN( int XLALSimInspiralTaylorT4PN(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarizat
REAL8TimeSeries **hcross, /**< x-polarization waveform */ ion waveform */
REAL8 phic, /**< coalescence phase */ REAL8TimeSeries **hcross, /**< x-polarizat
REAL8 deltaT, /**< sampling interval */ ion waveform */
REAL8 m1, /**< mass of companion 1 */ REAL8 phiRef, /**< reference o
REAL8 m2, /**< mass of companion 2 */ rbital phase (rad) */
REAL8 f_min, /**< start frequency */ REAL8 deltaT, /**< sampling in
REAL8 r, /**< distance of source */ terval (Hz) */
REAL8 i, /**< inclination of source (rad) * REAL8 m1, /**< mass of com
/ panion 1 (kg) */
int O /**< twice post-Newtonian order */ REAL8 m2, /**< mass of com
panion 2 (kg) */
REAL8 f_min, /**< start frequ
ency (Hz) */
REAL8 fRef, /**< reference f
requency (Hz) */
REAL8 r, /**< distance of
source (m) */
REAL8 i, /**< inclination
of source (rad) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int O /**< twice post-
Newtonian order */
); );
/** /**
* Driver routine to compute the restricted post-Newtonian inspiral wavefor m. * Driver routine to compute the restricted post-Newtonian inspiral wavefor m.
* *
* This routine computes the phasing to the specified order, but * This routine computes the phasing to the specified order, but
* only computes the amplitudes to the Newtonian (quadrupole) order. * only computes the amplitudes to the Newtonian (quadrupole) order.
* *
* Constant log term in amplitude set to 1. This is a gauge choice. * Constant log term in amplitude set to 1. This is a gauge choice.
*/ */
int XLALSimInspiralTaylorT4PNRestricted( int XLALSimInspiralTaylorT4PNRestricted(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarizat
REAL8TimeSeries **hcross, /**< x-polarization waveform */ ion waveform */
REAL8 phic, /**< coalescence phase */ REAL8TimeSeries **hcross, /**< x-polarizat
REAL8 deltaT, /**< sampling interval */ ion waveform */
REAL8 m1, /**< mass of companion 1 */ REAL8 phiRef, /**< reference o
REAL8 m2, /**< mass of companion 2 */ rbital phase (rad) */
REAL8 f_min, /**< start frequency */ REAL8 deltaT, /**< sampling in
REAL8 r, /**< distance of source */ terval (s) */
REAL8 i, /**< inclination of source (rad) * REAL8 m1, /**< mass of com
/ panion 1 (kg) */
int O /**< twice post-Newtonian phase or REAL8 m2, /**< mass of com
der */ panion 2 (kg) */
REAL8 f_min, /**< start frequ
ency (Hz) */
REAL8 fRef, /**< reference f
requency (Hz) */
REAL8 r, /**< distance of
source (m) */
REAL8 i, /**< inclination
of source (rad) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int O /**< twice post-
Newtonian phase order */
); );
/* TaylorT3 functions */ /* TaylorT3 functions */
/** /**
* Evolves a post-Newtonian orbit using the Taylor T3 method. * Evolves a post-Newtonian orbit using the Taylor T3 method.
*/ */
int XLALSimInspiralTaylorT3PNEvolveOrbit( int XLALSimInspiralTaylorT3PNEvolveOrbit(
REAL8TimeSeries **V, /**< post-Newtonian parameter [return REAL8TimeSeries **V, /**< post-Newton
ed] */ ian parameter [returned] */
REAL8TimeSeries **phi, /**< orbital phase [returned] */ REAL8TimeSeries **phi, /**< orbital pha
REAL8 phic, /**< coalescence phase */ se [returned] */
REAL8 deltaT, /**< sampling interval */ REAL8 phiRef, /**< reference o
REAL8 m1, /**< mass of companion 1 */ rbital phase (rad) */
REAL8 m2, /**< mass of companion 2 */ REAL8 deltaT, /**< sampling in
REAL8 f_min, /**< start frequency */ terval (s) */
int O /**< twice post-Newtonian order */ REAL8 m1, /**< mass of com
panion 1 (kg) */
REAL8 m2, /**< mass of com
panion 2 (kg) */
REAL8 f_min, /**< starting GW
frequency (Hz) */
REAL8 fRef, /**< reference G
W frequency (Hz) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int O /**< twice post-
Newtonian order */
); );
/** /**
* Driver routine to compute the post-Newtonian inspiral waveform. * Driver routine to compute the post-Newtonian inspiral waveform.
* *
* This routine allows the user to specify different pN orders * This routine allows the user to specify different pN orders
* for phasing calcuation vs. amplitude calculations. * for phasing calcuation vs. amplitude calculations.
*/ */
int XLALSimInspiralTaylorT3PNGenerator( int XLALSimInspiralTaylorT3PNGenerator(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarizat
REAL8TimeSeries **hcross, /**< x-polarization waveform */ ion waveform */
REAL8 phic, /**< coalescence phase */ REAL8TimeSeries **hcross, /**< x-polarizat
REAL8 v0, /**< tail-term gauge choice (defau ion waveform */
lt = 1) */ REAL8 phiRef, /**< reference o
REAL8 deltaT, /**< sampling interval */ rbital phase (rad) */
REAL8 m1, /**< mass of companion 1 */ REAL8 v0, /**< tail-term g
REAL8 m2, /**< mass of companion 2 */ auge choice (default = 1) */
REAL8 f_min, /**< start frequency */ REAL8 deltaT, /**< sampling in
REAL8 r, /**< distance of source */ terval (s) */
REAL8 i, /**< inclination of source (rad) * REAL8 m1, /**< mass of com
/ panion 1 (kg) */
int amplitudeO, /**< twice post-Newtonian amplitud REAL8 m2, /**< mass of com
e order */ panion 2 (kg) */
int phaseO /**< twice post-Newtonian phase or REAL8 f_min, /**< starting GW
der */ frequency (Hz) */
REAL8 fRef, /**< reference G
W frequency (Hz) */
REAL8 r, /**< distance of
source (m) */
REAL8 i, /**< inclination
of source (rad) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int amplitudeO, /**< twice post-
Newtonian amplitude order */
int phaseO /**< twice post-
Newtonian phase order */
);
/**
* Driver routine to compute the -2 spin-weighted spherical harmonic mode
* using TaylorT3 phasing.
*/
COMPLEX16TimeSeries *XLALSimInspiralTaylorT3PNModes(
REAL8 phiRef, /**< reference o
rbital phase (rad) */
REAL8 v0, /**< tail-term g
auge choice (default = 1) */
REAL8 deltaT, /**< sampling in
terval (s) */
REAL8 m1, /**< mass of com
panion 1 (kg) */
REAL8 m2, /**< mass of com
panion 2 (kg) */
REAL8 f_min, /**< starting GW
frequency (Hz) */
REAL8 fRef, /**< reference G
W frequency (Hz) */
REAL8 r, /**< distance of
source (m) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int amplitudeO, /**< twice post-
Newtonian amplitude order */
int phaseO, /**< twice post-
Newtonian phase order */
int l, /**< l index of
mode */
int m /**< m index of
mode */
); );
/** /**
* Driver routine to compute the post-Newtonian inspiral waveform. * Driver routine to compute the post-Newtonian inspiral waveform.
* *
* This routine uses the same pN order for phasing and amplitude * This routine uses the same pN order for phasing and amplitude
* (unless the order is -1 in which case the highest available * (unless the order is -1 in which case the highest available
* order is used for both of these -- which might not be the same). * order is used for both of these -- which might not be the same).
* *
* Constant log term in amplitude set to 1. This is a gauge choice. * Constant log term in amplitude set to 1. This is a gauge choice.
*/ */
int XLALSimInspiralTaylorT3PN( int XLALSimInspiralTaylorT3PN(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarizat
REAL8TimeSeries **hcross, /**< x-polarization waveform */ ion waveform */
REAL8 phic, /**< coalescence phase */ REAL8TimeSeries **hcross, /**< x-polarizat
REAL8 deltaT, /**< sampling interval */ ion waveform */
REAL8 m1, /**< mass of companion 1 */ REAL8 phiRef, /**< reference o
REAL8 m2, /**< mass of companion 2 */ rbital phase (rad) */
REAL8 f_min, /**< start frequency */ REAL8 deltaT, /**< sampling in
REAL8 r, /**< distance of source */ terval (s) */
REAL8 i, /**< inclination of source (rad) * REAL8 m1, /**< mass of com
/ panion 1 (kg) */
int O /**< twice post-Newtonian order */ REAL8 m2, /**< mass of com
panion 2 (kg) */
REAL8 f_min, /**< starting GW
frequency (Hz) */
REAL8 fRef, /**< reference G
W frequency (Hz) */
REAL8 r, /**< distance of
source (m) */
REAL8 i, /**< inclination
of source (rad) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int O /**< twice post-
Newtonian order */
); );
/** /**
* Driver routine to compute the restricted post-Newtonian inspiral wavefor m. * Driver routine to compute the restricted post-Newtonian inspiral wavefor m.
* *
* This routine computes the phasing to the specified order, but * This routine computes the phasing to the specified order, but
* only computes the amplitudes to the Newtonian (quadrupole) order. * only computes the amplitudes to the Newtonian (quadrupole) order.
* *
* Constant log term in amplitude set to 1. This is a gauge choice. * Constant log term in amplitude set to 1. This is a gauge choice.
*/ */
int XLALSimInspiralTaylorT3PNRestricted( int XLALSimInspiralTaylorT3PNRestricted(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarizat
REAL8TimeSeries **hcross, /**< x-polarization waveform */ ion waveform */
REAL8 phic, /**< coalescence phase */ REAL8TimeSeries **hcross, /**< x-polarizat
REAL8 deltaT, /**< sampling interval */ ion waveform */
REAL8 m1, /**< mass of companion 1 */ REAL8 phiRef, /**< reference o
REAL8 m2, /**< mass of companion 2 */ rbital phase (rad) */
REAL8 f_min, /**< start frequency */ REAL8 deltaT, /**< sampling in
REAL8 r, /**< distance of source */ terval (s) */
REAL8 i, /**< inclination of source (rad) * REAL8 m1, /**< mass of com
/ panion 1 (kg) */
int O /**< twice post-Newtonian phase or REAL8 m2, /**< mass of com
der */ panion 2 (kg) */
REAL8 f_min, /**< starting GW
frequency (Hz) */
REAL8 fRef, /**< reference G
W frequency (Hz) */
REAL8 r, /**< distance of
source (m)*/
REAL8 i, /**< inclination
of source (rad) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int O /**< twice post-
Newtonian phase order */
); );
/* TaylorT2 functions */ /* TaylorT2 functions */
/** /**
* Evolves a post-Newtonian orbit using the Taylor T2 method. * Evolves a post-Newtonian orbit using the Taylor T2 method.
*/ */
int XLALSimInspiralTaylorT2PNEvolveOrbit( int XLALSimInspiralTaylorT2PNEvolveOrbit(
REAL8TimeSeries **V, /**< post-Newtonian parameter [return REAL8TimeSeries **V, /**< post-Newton
ed] */ ian parameter [returned] */
REAL8TimeSeries **phi, /**< orbital phase [returned] */ REAL8TimeSeries **phi, /**< orbital pha
REAL8 phic, /**< coalescence phase */ se [returned] */
REAL8 deltaT, /**< sampling interval */ REAL8 phiRef, /**< reference o
REAL8 m1, /**< mass of companion 1 */ rbital phase (rad) */
REAL8 m2, /**< mass of companion 2 */ REAL8 deltaT, /**< sampling in
REAL8 f_min, /**< start frequency */ terval (s) */
int O /**< twice post-Newtonian order */ REAL8 m1, /**< mass of com
panion 1 (kg) */
REAL8 m2, /**< mass of com
panion 2 (kg) */
REAL8 f_min, /**< starting GW
frequency (Hz) */
REAL8 fRef, /**< reference G
W frequency (Hz) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int O /**< twice post-
Newtonian order */
); );
/** /**
* Driver routine to compute the post-Newtonian inspiral waveform. * Driver routine to compute the post-Newtonian inspiral waveform.
* *
* This routine allows the user to specify different pN orders * This routine allows the user to specify different pN orders
* for phasing calcuation vs. amplitude calculations. * for phasing calcuation vs. amplitude calculations.
*/ */
int XLALSimInspiralTaylorT2PNGenerator( int XLALSimInspiralTaylorT2PNGenerator(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarizat
REAL8TimeSeries **hcross, /**< x-polarization waveform */ ion waveform */
REAL8 phic, /**< coalescence phase */ REAL8TimeSeries **hcross, /**< x-polarizat
REAL8 v0, /**< tail-term gauge choice (defau ion waveform */
lt = 1) */ REAL8 phiRef, /**< reference o
REAL8 deltaT, /**< sampling interval */ rbital phase (rad) */
REAL8 m1, /**< mass of companion 1 */ REAL8 v0, /**< tail-term g
REAL8 m2, /**< mass of companion 2 */ auge choice (default = 1) */
REAL8 f_min, /**< start frequency */ REAL8 deltaT, /**< sampling in
REAL8 r, /**< distance of source */ terval (s) */
REAL8 i, /**< inclination of source (rad) * REAL8 m1, /**< mass of com
/ panion 1 (kg) */
int amplitudeO, /**< twice post-Newtonian amplitud REAL8 m2, /**< mass of com
e order */ panion 2 (kg) */
int phaseO /**< twice post-Newtonian phase or REAL8 f_min, /**< starting GW
der */ frequency (Hz) */
REAL8 fRef, /**< reference G
W frequency (Hz) */
REAL8 r, /**< distance of
source (m) */
REAL8 i, /**< inclination
of source (rad) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int amplitudeO, /**< twice post-
Newtonian amplitude order */
int phaseO /**< twice post-
Newtonian phase order */
);
/**
* Driver routine to compute the -2 spin-weighted spherical harmonic mode
* using TaylorT2 phasing.
*/
COMPLEX16TimeSeries *XLALSimInspiralTaylorT2PNModes(
REAL8 phiRef, /**< reference o
rbital phase (rad) */
REAL8 v0, /**< tail-term g
auge choice (default = 1) */
REAL8 deltaT, /**< sampling in
terval (s) */
REAL8 m1, /**< mass of com
panion 1 (kg) */
REAL8 m2, /**< mass of com
panion 2 (kg) */
REAL8 f_min, /**< starting GW
frequency (Hz) */
REAL8 fRef, /**< reference G
W frequency (Hz) */
REAL8 r, /**< distance of
source (m) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int amplitudeO, /**< twice post-
Newtonian amplitude order */
int phaseO, /**< twice post-
Newtonian phase order */
int l, /**< l index of
mode */
int m /**< m index of
mode */
); );
/** /**
* Driver routine to compute the post-Newtonian inspiral waveform. * Driver routine to compute the post-Newtonian inspiral waveform.
* *
* This routine uses the same pN order for phasing and amplitude * This routine uses the same pN order for phasing and amplitude
* (unless the order is -1 in which case the highest available * (unless the order is -1 in which case the highest available
* order is used for both of these -- which might not be the same). * order is used for both of these -- which might not be the same).
* *
* Constant log term in amplitude set to 1. This is a gauge choice. * Constant log term in amplitude set to 1. This is a gauge choice.
*/ */
int XLALSimInspiralTaylorT2PN( int XLALSimInspiralTaylorT2PN(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarizat
REAL8TimeSeries **hcross, /**< x-polarization waveform */ ion waveform */
REAL8 phic, /**< coalescence phase */ REAL8TimeSeries **hcross, /**< x-polarizat
REAL8 deltaT, /**< sampling interval */ ion waveform */
REAL8 m1, /**< mass of companion 1 */ REAL8 phiRef, /**< reference o
REAL8 m2, /**< mass of companion 2 */ rbital phase (rad) */
REAL8 f_min, /**< start frequency */ REAL8 deltaT, /**< sampling in
REAL8 r, /**< distance of source */ terval (s) */
REAL8 i, /**< inclination of source (rad) * REAL8 m1, /**< mass of com
/ panion 1 (kg) */
int O /**< twice post-Newtonian order */ REAL8 m2, /**< mass of com
panion 2 (kg) */
REAL8 f_min, /**< starting GW
frequency (Hz)*/
REAL8 fRef, /**< reference G
W frequency (Hz)*/
REAL8 r, /**< distance of
source (m) */
REAL8 i, /**< inclination
of source (rad) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int O /**< twice post-
Newtonian order */
); );
/** /**
* Driver routine to compute the restricted post-Newtonian inspiral wavefor m. * Driver routine to compute the restricted post-Newtonian inspiral wavefor m.
* *
* This routine computes the phasing to the specified order, but * This routine computes the phasing to the specified order, but
* only computes the amplitudes to the Newtonian (quadrupole) order. * only computes the amplitudes to the Newtonian (quadrupole) order.
* *
* Constant log term in amplitude set to 1. This is a gauge choice. * Constant log term in amplitude set to 1. This is a gauge choice.
*/ */
int XLALSimInspiralTaylorT2PNRestricted( int XLALSimInspiralTaylorT2PNRestricted(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarizat
REAL8TimeSeries **hcross, /**< x-polarization waveform */ ion waveform */
REAL8 phic, /**< coalescence phase */ REAL8TimeSeries **hcross, /**< x-polarizat
REAL8 deltaT, /**< sampling interval */ ion waveform */
REAL8 m1, /**< mass of companion 1 */ REAL8 phiRef, /**< reference o
REAL8 m2, /**< mass of companion 2 */ rbital phase (rad) */
REAL8 f_min, /**< start frequency */ REAL8 deltaT, /**< sampling in
REAL8 r, /**< distance of source */ terval (s) */
REAL8 i, /**< inclination of source (rad) * REAL8 m1, /**< mass of com
/ panion 1 (kg) */
int O /**< twice post-Newtonian phase or REAL8 m2, /**< mass of com
der */ panion 2 (kg) */
REAL8 f_min, /**< starting GW
frequency (Hz) */
REAL8 fRef, /**< reference G
W frequency (Hz) */
REAL8 r, /**< distance of
source (m) */
REAL8 i, /**< inclination
of source (rad) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int O /**< twice post-
Newtonian phase order */
); );
/* TaylorT1 functions */ /* TaylorT1 functions */
/** /**
* Evolves a post-Newtonian orbit using the Taylor T1 method. * Evolves a post-Newtonian orbit using the Taylor T1 method.
*/ */
int XLALSimInspiralTaylorT1PNEvolveOrbit( int XLALSimInspiralTaylorT1PNEvolveOrbit(
REAL8TimeSeries **V, /**< post-Newtonian parameter [return REAL8TimeSeries **V, /**< post-Newton
ed] */ ian parameter [returned] */
REAL8TimeSeries **phi, /**< orbital phase [returned] */ REAL8TimeSeries **phi, /**< orbital pha
REAL8 phic, /**< coalescence phase */ se [returned] */
REAL8 deltaT, /**< sampling interval */ REAL8 phiRef, /**< reference o
REAL8 m1, /**< mass of companion 1 */ rbital phase (rad) */
REAL8 m2, /**< mass of companion 2 */ REAL8 deltaT, /**< sampling in
REAL8 f_min, /**< start frequency */ terval (s) */
int O /**< twice post-Newtonian order */ REAL8 m1, /**< mass of com
panion 1 (kg) */
REAL8 m2, /**< mass of com
panion 2 (kg) */
REAL8 f_min, /**< start frequ
ency (Hz) */
REAL8 fRef, /**< reference f
requency (Hz) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int O /**< twice post-
Newtonian order */
); );
/** /**
* Driver routine to compute the post-Newtonian inspiral waveform. * Driver routine to compute the post-Newtonian inspiral waveform.
* *
* This routine allows the user to specify different pN orders * This routine allows the user to specify different pN orders
* for phasing calcuation vs. amplitude calculations. * for phasing calcuation vs. amplitude calculations.
*/ */
int XLALSimInspiralTaylorT1PNGenerator( int XLALSimInspiralTaylorT1PNGenerator(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarizat
REAL8TimeSeries **hcross, /**< x-polarization waveform */ ion waveform */
REAL8 phic, /**< coalescence phase */ REAL8TimeSeries **hcross, /**< x-polarizat
REAL8 v0, /**< tail-term gauge choice (defau ion waveform */
lt = 1) */ REAL8 phiRef, /**< reference o
REAL8 deltaT, /**< sampling interval */ rbital phase (rad) */
REAL8 m1, /**< mass of companion 1 */ REAL8 v0, /**< tail-term g
REAL8 m2, /**< mass of companion 2 */ auge choice (default = 1) */
REAL8 f_min, /**< start frequency */ REAL8 deltaT, /**< sampling in
REAL8 r, /**< distance of source */ terval (s) */
REAL8 i, /**< inclination of source (rad) * REAL8 m1, /**< mass of com
/ panion 1 (kg) */
int amplitudeO, /**< twice post-Newtonian amplitud REAL8 m2, /**< mass of com
e order */ panion 2 (kg) */
int phaseO /**< twice post-Newtonian phase or REAL8 f_min, /**< starting GW
der */ frequency (Hz) */
REAL8 fRef, /**< reference G
W frequency (Hz) */
REAL8 r, /**< distance of
source (m) */
REAL8 i, /**< inclination
of source (rad) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int amplitudeO, /**< twice post-
Newtonian amplitude order */
int phaseO /**< twice post-
Newtonian phase order */
);
/**
* Driver routine to compute the -2 spin-weighted spherical harmonic mode
* using TaylorT1 phasing.
*/
COMPLEX16TimeSeries *XLALSimInspiralTaylorT1PNModes(
REAL8 phiRef, /**< reference o
rbital phase (rad) */
REAL8 v0, /**< tail-term g
auge choice (default = 1) */
REAL8 deltaT, /**< sampling in
terval (s) */
REAL8 m1, /**< mass of com
panion 1 (kg) */
REAL8 m2, /**< mass of com
panion 2 (kg) */
REAL8 f_min, /**< starting GW
frequency (Hz) */
REAL8 fRef, /**< reference G
W frequency (Hz) */
REAL8 r, /**< distance of
source (m) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int amplitudeO, /**< twice post-
Newtonian amplitude order */
int phaseO, /**< twice post-
Newtonian phase order */
int l, /**< l index of
mode */
int m /**< m index of
mode */
); );
/** /**
* Driver routine to compute the post-Newtonian inspiral waveform. * Driver routine to compute the post-Newtonian inspiral waveform.
* *
* This routine uses the same pN order for phasing and amplitude * This routine uses the same pN order for phasing and amplitude
* (unless the order is -1 in which case the highest available * (unless the order is -1 in which case the highest available
* order is used for both of these -- which might not be the same). * order is used for both of these -- which might not be the same).
* *
* Constant log term in amplitude set to 1. This is a gauge choice. * Constant log term in amplitude set to 1. This is a gauge choice.
*/ */
int XLALSimInspiralTaylorT1PN( int XLALSimInspiralTaylorT1PN(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarizat
REAL8TimeSeries **hcross, /**< x-polarization waveform */ ion waveform */
REAL8 phic, /**< coalescence phase */ REAL8TimeSeries **hcross, /**< x-polarizat
REAL8 deltaT, /**< sampling interval */ ion waveform */
REAL8 m1, /**< mass of companion 1 */ REAL8 phiRef, /**< reference o
REAL8 m2, /**< mass of companion 2 */ rbital phase (rad) */
REAL8 f_min, /**< start frequency */ REAL8 deltaT, /**< sampling in
REAL8 r, /**< distance of source */ terval (s) */
REAL8 i, /**< inclination of source (rad) * REAL8 m1, /**< mass of com
/ panion 1 (kg) */
int O /**< twice post-Newtonian order */ REAL8 m2, /**< mass of com
panion 2 (kg) */
REAL8 f_min, /**< start frequ
ency (Hz) */
REAL8 fRef, /**< reference f
requency (Hz) */
REAL8 r, /**< distance of
source (m) */
REAL8 i, /**< inclination
of source (rad) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int O /**< twice post-
Newtonian order */
); );
/** /**
* Driver routine to compute the restricted post-Newtonian inspiral wavefor m. * Driver routine to compute the restricted post-Newtonian inspiral wavefor m.
* *
* This routine computes the phasing to the specified order, but * This routine computes the phasing to the specified order, but
* only computes the amplitudes to the Newtonian (quadrupole) order. * only computes the amplitudes to the Newtonian (quadrupole) order.
* *
* Constant log term in amplitude set to 1. This is a gauge choice. * Constant log term in amplitude set to 1. This is a gauge choice.
*/ */
int XLALSimInspiralTaylorT1PNRestricted( int XLALSimInspiralTaylorT1PNRestricted(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarizat
REAL8TimeSeries **hcross, /**< x-polarization waveform */ ion waveform */
REAL8 phic, /**< coalescence phase */ REAL8TimeSeries **hcross, /**< x-polarizat
REAL8 deltaT, /**< sampling interval */ ion waveform */
REAL8 m1, /**< mass of companion 1 */ REAL8 phiRef, /**< reference o
REAL8 m2, /**< mass of companion 2 */ rbital phase (rad) */
REAL8 f_min, /**< start frequency */ REAL8 deltaT, /**< sampling in
REAL8 r, /**< distance of source */ terval (s) */
REAL8 i, /**< inclination of source (rad) * REAL8 m1, /**< mass of com
/ panion 1 (kg) */
int O /**< twice post-Newtonian phase or REAL8 m2, /**< mass of com
der */ panion 2 (kg) */
REAL8 f_min, /**< starting GW
frequency (Hz) */
REAL8 fRef, /**< reference G
W frequency (Hz) */
REAL8 r, /**< distance of
source (m)*/
REAL8 i, /**< inclination
of source (rad) */
REAL8 lambda1, /**< (tidal defo
rmability of body 1)/(mass of body 1)^5 */
REAL8 lambda2, /**< (tidal defo
rmability of body 2)/(mass of body 2)^5 */
LALSimInspiralInteraction interactionFlags, /**< flag to con
trol spin and tidal effects */
int O /**< twice post-
Newtonian phase order */
); );
/* TaylorEt functions */ /* TaylorEt functions */
/** /**
* Evolves a post-Newtonian orbit using the Taylor T1 method. * Evolves a post-Newtonian orbit using the Taylor T1 method.
*/ */
int XLALSimInspiralTaylorEtPNEvolveOrbit( int XLALSimInspiralTaylorEtPNEvolveOrbit(
REAL8TimeSeries **V, /**< post-Newtonian parameter [return ed] */ REAL8TimeSeries **V, /**< post-Newtonian parameter [return ed] */
REAL8TimeSeries **phi, /**< orbital phase [returned] */ REAL8TimeSeries **phi, /**< orbital phase [returned] */
REAL8 phic, /**< coalescence phase */ REAL8 phic, /**< orbital phase at end */
REAL8 deltaT, /**< sampling interval */ REAL8 deltaT, /**< sampling interval */
REAL8 m1, /**< mass of companion 1 */ REAL8 m1, /**< mass of companion 1 */
REAL8 m2, /**< mass of companion 2 */ REAL8 m2, /**< mass of companion 2 */
REAL8 f_min, /**< start frequency */ REAL8 f_min, /**< start frequency */
int O /**< twice post-Newtonian order */ int O /**< twice post-Newtonian order */
); );
/** /**
* Driver routine to compute the post-Newtonian inspiral waveform. * Driver routine to compute the post-Newtonian inspiral waveform.
* *
* This routine allows the user to specify different pN orders * This routine allows the user to specify different pN orders
* for phasing calcuation vs. amplitude calculations. * for phasing calcuation vs. amplitude calculations.
*/ */
int XLALSimInspiralTaylorEtPNGenerator( int XLALSimInspiralTaylorEtPNGenerator(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarization waveform */
REAL8TimeSeries **hcross, /**< x-polarization waveform */ REAL8TimeSeries **hcross, /**< x-polarization waveform */
REAL8 phic, /**< coalescence phase */ REAL8 phic, /**< orbital phase at end */
REAL8 x0, /**< tail-term gauge choice thing REAL8 x0, /**< tail-term gauge choice (if yo
(if you don't know, just set it to zero) */ u don't know, just set it to zero) */
REAL8 deltaT, /**< sampling interval */ REAL8 deltaT, /**< sampling interval */
REAL8 m1, /**< mass of companion 1 */ REAL8 m1, /**< mass of companion 1 */
REAL8 m2, /**< mass of companion 2 */ REAL8 m2, /**< mass of companion 2 */
REAL8 f_min, /**< start frequency */ REAL8 f_min, /**< start frequency */
REAL8 r, /**< distance of source */ REAL8 r, /**< distance of source */
REAL8 i, /**< inclination of source (rad) * / REAL8 i, /**< inclination of source (rad) * /
int amplitudeO, /**< twice post-Newtonian amplitud e order */ int amplitudeO, /**< twice post-Newtonian amplitud e order */
int phaseO /**< twice post-Newtonian phase or der */ int phaseO /**< twice post-Newtonian phase or der */
); );
skipping to change at line 875 skipping to change at line 1561
* *
* This routine uses the same pN order for phasing and amplitude * This routine uses the same pN order for phasing and amplitude
* (unless the order is -1 in which case the highest available * (unless the order is -1 in which case the highest available
* order is used for both of these -- which might not be the same). * order is used for both of these -- which might not be the same).
* *
* Log terms in amplitudes are ignored. This is a gauge choice. * Log terms in amplitudes are ignored. This is a gauge choice.
*/ */
int XLALSimInspiralTaylorEtPN( int XLALSimInspiralTaylorEtPN(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarization waveform */
REAL8TimeSeries **hcross, /**< x-polarization waveform */ REAL8TimeSeries **hcross, /**< x-polarization waveform */
REAL8 phic, /**< coalescence phase */ REAL8 phic, /**< orbital phase at end */
REAL8 deltaT, /**< sampling interval */ REAL8 deltaT, /**< sampling interval */
REAL8 m1, /**< mass of companion 1 */ REAL8 m1, /**< mass of companion 1 */
REAL8 m2, /**< mass of companion 2 */ REAL8 m2, /**< mass of companion 2 */
REAL8 f_min, /**< start frequency */ REAL8 f_min, /**< start frequency */
REAL8 r, /**< distance of source */ REAL8 r, /**< distance of source */
REAL8 i, /**< inclination of source (rad) * / REAL8 i, /**< inclination of source (rad) * /
int O /**< twice post-Newtonian order */ int O /**< twice post-Newtonian order */
); );
/** /**
* Driver routine to compute the restricted post-Newtonian inspiral wavefor m. * Driver routine to compute the restricted post-Newtonian inspiral wavefor m.
* *
* This routine computes the phasing to the specified order, but * This routine computes the phasing to the specified order, but
* only computes the amplitudes to the Newtonian (quadrupole) order. * only computes the amplitudes to the Newtonian (quadrupole) order.
* *
* Log terms in amplitudes are ignored. This is a gauge choice. * Log terms in amplitudes are ignored. This is a gauge choice.
*/ */
int XLALSimInspiralTaylorEtPNRestricted( int XLALSimInspiralTaylorEtPNRestricted(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarization waveform */
REAL8TimeSeries **hcross, /**< x-polarization waveform */ REAL8TimeSeries **hcross, /**< x-polarization waveform */
REAL8 phic, /**< coalescence phase */ REAL8 phic, /**< orbital phase at end */
REAL8 deltaT, /**< sampling interval */ REAL8 deltaT, /**< sampling interval */
REAL8 m1, /**< mass of companion 1 */ REAL8 m1, /**< mass of companion 1 */
REAL8 m2, /**< mass of companion 2 */ REAL8 m2, /**< mass of companion 2 */
REAL8 f_min, /**< start frequency */ REAL8 f_min, /**< start frequency */
REAL8 r, /**< distance of source */ REAL8 r, /**< distance of source */
REAL8 i, /**< inclination of source (rad) * / REAL8 i, /**< inclination of source (rad) * /
int O /**< twice post-Newtonian phase or der */ int O /**< twice post-Newtonian phase or der */
); );
/** /**
* Computes the stationary phase approximation to the Fourier transform of
* a chirp waveform with phase given by Eq.\eqref{eq_InspiralFourierPhase_f
2}
* and amplitude given by expanding \f$1/\sqrt{\dot{F}}\f$. If the PN order
is
* set to -1, then the highest implemented order is used.
*/
int XLALSimInspiralTaylorF2(
COMPLEX16FrequencySeries **htilde, /**< FD waveform */
const REAL8 phic, /**< orbital coalescence ph
ase (rad) */
const REAL8 deltaF, /**< frequency resolution *
/
const REAL8 m1_SI, /**< mass of companion 1 (k
g) */
const REAL8 m2_SI, /**< mass of companion 2 (k
g) */
const REAL8 fStart, /**< start GW frequency (Hz
) */
const REAL8 r, /**< distance of source (m)
*/
const REAL8 lambda1, /**< (tidal deformation of
body 1)/(mass of body 1)^5 */
const REAL8 lambda2, /**< (tidal deformation of
body 2)/(mass of body 2)^5 */
const LALSimInspiralInteraction interactionFlags, /**< flag
to control spin and tidal effects */
const INT4 phaseO, /**< twice PN phase order *
/
const INT4 amplitudeO /**< twice PN amplitude ord
er */
);
/**
* Functions for generic spinning waveforms. * Functions for generic spinning waveforms.
* Reproduce and extend old SpinTaylor(Frameless) and SQTPN waveforms * Reproduce and extend old SpinTaylor(Frameless) and SQTPN waveforms
*/ */
/** /**
* This function evolves the orbital equations for a precessing binary usin g * This function evolves the orbital equations for a precessing binary usin g
* the \"TaylorT4\" approximant for solving the orbital dynamics * the \"TaylorT4\" approximant for solving the orbital dynamics
* (see arXiv:0907.0700 for a review of the various PN approximants). * (see arXiv:0907.0700 for a review of the various PN approximants).
* *
* It returns time series of the \"orbital velocity\", orbital phase, * It returns time series of the \"orbital velocity\", orbital phase,
skipping to change at line 949 skipping to change at line 1656
REAL8TimeSeries **S1z, /**< -- z component [returned]*/ REAL8TimeSeries **S1z, /**< -- z component [returned]*/
REAL8TimeSeries **S2x, /**< Spin2 vector x component [returned]*/ REAL8TimeSeries **S2x, /**< Spin2 vector x component [returned]*/
REAL8TimeSeries **S2y, /**< -- y component [returned]*/ REAL8TimeSeries **S2y, /**< -- y component [returned]*/
REAL8TimeSeries **S2z, /**< -- z component [returned]*/ REAL8TimeSeries **S2z, /**< -- z component [returned]*/
REAL8TimeSeries **LNhatx, /**< unit orbital ang. mom. x [returned]*/ REAL8TimeSeries **LNhatx, /**< unit orbital ang. mom. x [returned]*/
REAL8TimeSeries **LNhaty, /**< -- y component [returned]*/ REAL8TimeSeries **LNhaty, /**< -- y component [returned]*/
REAL8TimeSeries **LNhatz, /**< -- z component [returned]*/ REAL8TimeSeries **LNhatz, /**< -- z component [returned]*/
REAL8TimeSeries **E1x, /**< orb. plane basis vector x[returned]*/ REAL8TimeSeries **E1x, /**< orb. plane basis vector x[returned]*/
REAL8TimeSeries **E1y, /**< -- y component [returned]*/ REAL8TimeSeries **E1y, /**< -- y component [returned]*/
REAL8TimeSeries **E1z, /**< -- z component [returned]*/ REAL8TimeSeries **E1z, /**< -- z component [returned]*/
REAL8 phi_end, /**< orbital phase at last sample */
REAL8 deltaT, /**< sampling interval (s) */ REAL8 deltaT, /**< sampling interval (s) */
REAL8 m1, /**< mass of companion 1 (kg) */ REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */ REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 fStart, /**< start frequency */ REAL8 fStart, /**< starting GW frequency */
REAL8 fEnd, /**< ending GW frequency, fEnd=0 means int
egrate as far forward as possible */
REAL8 s1x, /**< initial value of S1x */ REAL8 s1x, /**< initial value of S1x */
REAL8 s1y, /**< initial value of S1y */ REAL8 s1y, /**< initial value of S1y */
REAL8 s1z, /**< initial value of S1z */ REAL8 s1z, /**< initial value of S1z */
REAL8 s2x, /**< initial value of S2x */ REAL8 s2x, /**< initial value of S2x */
REAL8 s2y, /**< initial value of S2y */ REAL8 s2y, /**< initial value of S2y */
REAL8 s2z, /**< initial value of S2z */ REAL8 s2z, /**< initial value of S2z */
REAL8 lnhatx, /**< initial value of LNhatx */ REAL8 lnhatx, /**< initial value of LNhatx */
REAL8 lnhaty, /**< initial value of LNhaty */ REAL8 lnhaty, /**< initial value of LNhaty */
REAL8 lnhatz, /**< initial value of LNhatz */ REAL8 lnhatz, /**< initial value of LNhatz */
REAL8 e1x, /**< initial value of E1x */ REAL8 e1x, /**< initial value of E1x */
skipping to change at line 977 skipping to change at line 1684
REAL8 lambda2, /**< (tidal deformability of mass 2) / (total mass)^5 (dimensionless) */ REAL8 lambda2, /**< (tidal deformability of mass 2) / (total mass)^5 (dimensionless) */
LALSimInspiralInteraction interactionFlags, /**< flag to control spin and tidal effects */ LALSimInspiralInteraction interactionFlags, /**< flag to control spin and tidal effects */
INT4 phaseO /**< twice post-Newtonian order */ INT4 phaseO /**< twice post-Newtonian order */
); );
/** /**
* Driver routine to compute a precessing post-Newtonian inspiral waveform * Driver routine to compute a precessing post-Newtonian inspiral waveform
* with phasing computed from energy balance using the so-called \"T4\" met hod. * with phasing computed from energy balance using the so-called \"T4\" met hod.
* *
* This routine allows the user to specify different pN orders * This routine allows the user to specify different pN orders
* for phasing calcuation vs. amplitude calculations. * for the phasing and amplitude of the waveform.
*
* The reference frequency fRef is used as follows:
* 1) if fRef = 0: The initial values of s1, s2, lnhat and e1 will be the
* values at frequency fStart. The orbital phase of the last sample is s
et
* to phiRef (i.e. phiRef is the "coalescence phase", roughly speaking).
* THIS IS THE DEFAULT BEHAVIOR CONSISTENT WITH OTHER APPROXIMANTS
*
* 2) If fRef = fStart: The initial values of s1, s2, lnhat and e1 will be
the
* values at frequency fStart. phiRef is used to set the orbital phase
* of the first sample at fStart.
*
* 3) If fRef > fStart: The initial values of s1, s2, lnhat and e1 will be
the
* values at frequency fRef. phiRef is used to set the orbital phase at
fRef.
* The code will integrate forwards and backwards from fRef and stitch t
he
* two together to create a complete waveform. This allows one to specif
y
* the orientation of the binary in-band (or at any arbitrary point).
* Otherwise, the user can only directly control the initial orientation
.
*
* 4) fRef < 0 or fRef >= Schwarz. ISCO are forbidden and the code will abo
rt.
*/ */
int XLALSimInspiralSpinTaylorT4( int XLALSimInspiralSpinTaylorT4(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ REAL8TimeSeries **hplus, /**< +-polarization waveform */
REAL8TimeSeries **hcross, /**< x-polarization waveform */ REAL8TimeSeries **hcross, /**< x-polarization waveform */
REAL8 phi_end, /**< GW phase of final sample (rad REAL8 phiRef, /**< orbital phase at reference pt.
) */ */
REAL8 v0, /**< tail gauge term (default = 1) REAL8 v0, /**< tail gauge term (default = 1) *
*/ /
REAL8 deltaT, /**< sampling interval (s) */ REAL8 deltaT, /**< sampling interval (s) */
REAL8 m1, /**< mass of companion 1 (kg) */ REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */ REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 fStart, /**< start GW frequency (Hz) */ REAL8 fStart, /**< start GW frequency (Hz) */
REAL8 r, /**< distance of source (m) */ REAL8 fRef, /**< reference GW frequency (Hz) */
REAL8 s1x, /**< initial value of S1x */ REAL8 r, /**< distance of source (m) */
REAL8 s1y, /**< initial value of S1y */ REAL8 s1x, /**< initial value of S1x */
REAL8 s1z, /**< initial value of S1z */ REAL8 s1y, /**< initial value of S1y */
REAL8 s2x, /**< initial value of S2x */ REAL8 s1z, /**< initial value of S1z */
REAL8 s2y, /**< initial value of S2y */ REAL8 s2x, /**< initial value of S2x */
REAL8 s2z, /**< initial value of S2z */ REAL8 s2y, /**< initial value of S2y */
REAL8 lnhatx, /**< initial value of LNhatx */ REAL8 s2z, /**< initial value of S2z */
REAL8 lnhaty, /**< initial value of LNhaty */ REAL8 lnhatx, /**< initial value of LNhatx */
REAL8 lnhatz, /**< initial value of LNhatz */ REAL8 lnhaty, /**< initial value of LNhaty */
REAL8 e1x, /**< initial value of E1x */ REAL8 lnhatz, /**< initial value of LNhatz */
REAL8 e1y, /**< initial value of E1y */ REAL8 e1x, /**< initial value of E1x */
REAL8 e1z, /**< initial value of E1z */ REAL8 e1y, /**< initial value of E1y */
REAL8 lambda1, /**< (tidal deformability of mass REAL8 e1z, /**< initial value of E1z */
1) / (total mass)^5 (dimensionless) */ REAL8 lambda1, /**< (tidal deformability of mass 1)
REAL8 lambda2, /**< (tidal deformability of mass / (total mass)^5 (dimensionless) */
2) / (total mass)^5 (dimensionless) */ REAL8 lambda2, /**< (tidal deformability of mass 2)
LALSimInspiralInteraction interactionFlags, /**< flag to con / (total mass)^5 (dimensionless) */
trol spin and tidal effects */ LALSimInspiralInteraction interactionFlags, /**< flag to control spi
int phaseO, /**< twice PN phase order */ n and tidal effects */
int amplitudeO /**< twice PN amplitude order */ int phaseO, /**< twice PN phase order */
); int amplitudeO /**< twice PN amplitude order */
);
/** /**
* Driver routine to compute a precessing post-Newtonian inspiral waveform * Driver routine to compute the physical template family "Q" vectors using
* with phasing computed from energy balance using the so-called \"T4\" met * the \"T4\" method. Note that PTF describes single spin systems
hod.
* *
* This routine assumes leading-order amplitude dependence (restricted wave * This routine requires leading-order amplitude dependence
form) * but allows the user to specify the phase PN order
* but allows hte user to specify the phase PN order
*/ */
int XLALSimInspiralRestrictedSpinTaylorT4(
REAL8TimeSeries **hplus, /**< +-polarization waveform */ int XLALSimInspiralSpinTaylorT4PTFQVecs(
REAL8TimeSeries **hcross, /**< x-polarization waveform */ REAL8TimeSeries **Q1, /**< Q1 output vector */
REAL8 phi_end, /**< GW phase of final sample (rad REAL8TimeSeries **Q2, /**< Q2 output vector */
) */ REAL8TimeSeries **Q3, /**< Q3 output vector */
REAL8 v0, /**< tail gauge term (default = 1) REAL8TimeSeries **Q4, /**< Q4 output vector */
*/ REAL8TimeSeries **Q5, /**< Q5 output vector */
REAL8 deltaT, /**< sampling interval (s) */ REAL8 deltaT, /**< sampling interval (s) */
REAL8 m1, /**< mass of companion 1 (kg) */ REAL8 m1, /**< mass of companion 1 (kg) */
REAL8 m2, /**< mass of companion 2 (kg) */ REAL8 m2, /**< mass of companion 2 (kg) */
REAL8 fStart, /**< start GW frequency (Hz) */ REAL8 chi1, /**< spin magnitude (|S1|) */
REAL8 r, /**< distance of source (m) */ REAL8 kappa, /**< L . S (1 if they are aligned)
REAL8 s1x, /**< initial value of S1x */ */
REAL8 s1y, /**< initial value of S1y */ REAL8 fStart, /**< start GW frequency (Hz) */
REAL8 s1z, /**< initial value of S1z */ REAL8 lambda1, /**< (tidal deformability of mass 1
REAL8 s2x, /**< initial value of S2x */ ) / (total mass)^5 (dimensionless) */
REAL8 s2y, /**< initial value of S2y */ REAL8 lambda2, /**< (tidal deformability of mass 2
REAL8 s2z, /**< initial value of S2z */ ) / (total mass)^5 (dimensionless) */
REAL8 lnhatx, /**< initial value of LNhatx */ LALSimInspiralInteraction interactionFlags, /**< flag to control sp
REAL8 lnhaty, /**< initial value of LNhaty */ in and tidal effects */
REAL8 lnhatz, /**< initial value of LNhatz */ int phaseO /**< twice PN phase order */
REAL8 e1x, /**< initial value of E1x */ );
REAL8 e1y, /**< initial value of E1y */
REAL8 e1z, /**< initial value of E1z */
REAL8 lambda1, /**< (tidal deformability of mass
1) / (total mass)^5 (dimensionless) */
REAL8 lambda2, /**< (tidal deformability of mass
2) / (total mass)^5 (dimensionless) */
LALSimInspiralInteraction interactionFlags, /**< flag to control
spin and tidal effects */
int phaseO /**< twice PN phase order */
);
/** /**
* Function to specify the desired orientation of a precessing binary in te rms * Function to specify the desired orientation of a precessing binary in te rms
* of several angles and then compute the vector components in the so-calle d * of several angles and then compute the vector components in the so-calle d
* \"radiation frame\" (with the z-axis along the direction of propagation) as * \"radiation frame\" (with the z-axis along the direction of propagation) as
* needed for initial conditions for the SpinTaylorT4 waveform routines. * needed to specify binary configuration for ChooseTDWaveform.
* *
* Input: * Input:
* thetaJN, phiJN are angles describing the desired orientation of the * thetaJN is the inclination between total angular momentum (J) and th
* total angular momentum (J) relative to direction of propagation (N) e
* theta1, phi1, theta2, phi2 are angles describing the desired orienta * direction of propagation (N)
tion * theta1 and theta2 are the inclinations of S1 and S2
* of spin 1 and 2 relative to the Newtonian orbital angular momentum (L_N) * measured from the Newtonian orbital angular momentum (L_N)
* m1, m2, f0 are the component masses and initial GW frequency, * phi12 is the difference in azimuthal angles of S1 and S2.
* they are needed to compute the magnitude of L_N, and thus J * chi1, chi2 are the dimensionless spin magnitudes ( \f$0 \le chi1,2 \
* chi1, chi2 are the dimensionless spin magnitudes ( 0 <= chi1,2 <= 1) le 1\f$)
, * phiJL is the azimuthal angle of L_N on its cone about J.
* they are needed to compute the magnitude of S1 and S2, and thus J * m1, m2, f_ref are the component masses and reference GW frequency,
* they are needed to compute the magnitude of L_N, and thus J.
* *
* Output: * Output:
* x, y, z components of LNhat (unit vector along orbital angular momen tum), * incl - inclination angle of L_N relative to N
* x, y, z components of E1 (unit vector in the initial orbital plane) * x, y, z components of E1 (unit vector in the initial orbital plane)
* x, y, z components S1 and S2 (unit spin vectors times their * x, y, z components S1 and S2 (unit spin vectors times their
* dimensionless spin magnitudes - i.e. they have unit magnitude for * dimensionless spin magnitudes - i.e. they have unit magnitude fo
* extremal BHs and smaller magnitude for slower spins) r
* extremal BHs and smaller magnitude for slower spins).
* *
* NOTE: Here the \"total\" angular momentum is computed as * NOTE: Here the \"total\" angular momentum is computed as
* J = L_N + S1 + S2 * J = L_N + S1 + S2
* where L_N is the Newtonian orbital angular momentum. In fact, there are * where L_N is the Newtonian orbital angular momentum. In fact, there are
* PN corrections to L which contribute to J that are NOT ACCOUNTED FOR * PN corrections to L which contribute to J that are NOT ACCOUNTED FOR
* in this function. This is done so the function does not need to know abo ut * in this function. This is done so the function does not need to know abo ut
* the PN order of the system and to avoid subtleties with spin-orbit * the PN order of the system and to avoid subtleties with spin-orbit
* contributions to L. Also, it is believed that the difference in Jhat * contributions to L. Also, it is believed that the difference in Jhat
* with or without these PN corrections to L is quite small. * with or without these PN corrections to L is quite small.
*
* NOTE: fRef = 0 is not a valid choice. If you will pass fRef=0 into
* ChooseWaveform, then here pass in f_min, the starting GW frequency
*
* The various rotations in this transformation are described in more detai
l
* in a Mathematica notebook available here:
* https://www.lsc-group.phys.uwm.edu/ligovirgo/cbcnote/Waveforms/Transform
PrecessingInitialConditions
*/ */
int XLALSimInspiralTransformPrecessingInitialConditions( int XLALSimInspiralTransformPrecessingInitialConditions(
REAL8 *LNhatx, /**< LNhat x component (returned) */ REAL8 *incl, /**< Inclination angle of L_N (returned) */
REAL8 *LNhaty, /**< LNhat y component (returned) */
REAL8 *LNhatz, /**< LNhat z component (returned) */
REAL8 *E1x, /**< E1 x component (returned) */
REAL8 *E1y, /**< E1 y component (returned) */
REAL8 *E1z, /**< E1 z component (returned) */
REAL8 *S1x, /**< S1 x component (returned) */ REAL8 *S1x, /**< S1 x component (returned) */
REAL8 *S1y, /**< S1 y component (returned) */ REAL8 *S1y, /**< S1 y component (returned) */
REAL8 *S1z, /**< S1 z component (returned) */ REAL8 *S1z, /**< S1 z component (returned) */
REAL8 *S2x, /**< S2 x component (returned) */ REAL8 *S2x, /**< S2 x component (returned) */
REAL8 *S2y, /**< S2 y component (returned) */ REAL8 *S2y, /**< S2 y component (returned) */
REAL8 *S2z, /**< S2 z component (returned) */ REAL8 *S2z, /**< S2 z component (returned) */
REAL8 thetaJN, /**< zenith angle between J and N */ REAL8 thetaJN, /**< zenith angle between J and N (rad) */
REAL8 phiJN, /**< azimuth angle between J and N */ REAL8 phiJL, /**< azimuthal angle of L_N on its cone abou
REAL8 theta1, /**< zenith angle between S1 and LNhat */ t J (rad) */
REAL8 phi1, /**< azimuth angle between S1 and LNhat */ REAL8 theta1, /**< zenith angle between S1 and LNhat (rad)
REAL8 theta2, /**< zenith angle between S2 and LNhat */ */
REAL8 phi2, /**< azimuth angle between S2 and LNhat */ REAL8 theta2, /**< zenith angle between S2 and LNhat (rad)
*/
REAL8 phi12, /**< difference in azimuthal angle btwn S1,
S2 (rad) */
REAL8 chi1, /**< dimensionless spin of body 1 */
REAL8 chi2, /**< dimensionless spin of body 2 */
REAL8 m1, /**< mass of body 1 (kg) */ REAL8 m1, /**< mass of body 1 (kg) */
REAL8 m2, /**< mass of body 2 (kg) */ REAL8 m2, /**< mass of body 2 (kg) */
REAL8 f0, /**< initial GW frequency (Hz) */ REAL8 fRef /**< reference GW frequency (Hz) */
REAL8 chi1, /**< dimensionless spin of body 1 */
REAL8 chi2 /**< dimensionless spin of body 2 */
); );
/** /**
* Driver routine to compute a non-precessing post-Newtonian inspiral wavef orm * Driver routine to compute a non-precessing post-Newtonian inspiral wavef orm
* in the frequency domain, described in http://arxiv.org/abs/1107.1267. * in the frequency domain, described in http://arxiv.org/abs/1107.1267.
* *
* The chi parameter should be determined from * The chi parameter should be determined from
* XLALSimInspiralTaylorF2ReducedSpinComputeChi. * XLALSimInspiralTaylorF2ReducedSpinComputeChi.
* *
* A note from Evan Ochsner on differences with respect to TaylorF2: * A note from Evan Ochsner on differences with respect to TaylorF2:
skipping to change at line 1135 skipping to change at line 1855
* In particular, the F2ReducedSpin will always have only the 2nd harmonic, but * In particular, the F2ReducedSpin will always have only the 2nd harmonic, but
* F2 will have multiple harmonics starting at ampO = 0.5PN. Even if you we re * F2 will have multiple harmonics starting at ampO = 0.5PN. Even if you we re
* to compare just the 2nd harmonic, you would have a difference starting a t * to compare just the 2nd harmonic, you would have a difference starting a t
* 1PN ampO, because the F2 has a 1PN TD amp. correction to the 2nd harmoni c * 1PN ampO, because the F2 has a 1PN TD amp. correction to the 2nd harmoni c
* (alpha/beta_2^(2)(t)) which will not be accounted for by the F2ReducedSp in. * (alpha/beta_2^(2)(t)) which will not be accounted for by the F2ReducedSp in.
* So, the two should agree when ampO=0, but will be different in any other * So, the two should agree when ampO=0, but will be different in any other
* case. * case.
*/ */
int XLALSimInspiralTaylorF2ReducedSpin( int XLALSimInspiralTaylorF2ReducedSpin(
COMPLEX16FrequencySeries **htilde, /**< FD waveform */ COMPLEX16FrequencySeries **htilde, /**< FD waveform */
const REAL8 phiStart, /**< initial GW phase (rad) const REAL8 phic, /**< orbital coalescence phase (rad
*/ ) */
const REAL8 deltaF, /**< frequency resolution * const REAL8 deltaF, /**< frequency resolution (Hz) */
/ const REAL8 m1_SI, /**< mass of companion 1 (kg) */
const REAL8 m1_SI, /**< mass of companion 1 (k const REAL8 m2_SI, /**< mass of companion 2 (kg) */
g) */ const REAL8 chi, /**< dimensionless aligned-spin par
const REAL8 m2_SI, /**< mass of companion 2 (k am */
g) */ const REAL8 fStart, /**< start GW frequency (Hz) */
const REAL8 chi, /**< dimensionless aligned- const REAL8 r, /**< distance of source (m) */
spin param */ const INT4 phaseO, /**< twice PN phase order */
const REAL8 fStart, /**< start GW frequency (Hz const INT4 ampO /**< twice PN amplitude order */
) */
const REAL8 r, /**< distance of source (m)
*/
const UINT4 phaseO, /**< twice PN phase order *
/
const UINT4 ampO /**< twice PN amplitude ord
er */
); );
/** /**
* Generate the \"reduced-spin templates\" proposed in http://arxiv.org/abs /1107.1267 * Generate the \"reduced-spin templates\" proposed in http://arxiv.org/abs /1107.1267
* Add the tidal phase terms from http://arxiv.org/abs/1101.1673 (Eqs. 3.9, 3.10) * Add the tidal phase terms from http://arxiv.org/abs/1101.1673 (Eqs. 3.9, 3.10)
* The chi parameter should be determined from XLALSimInspiralTaylorF2Reduc edSpinComputeChi. * The chi parameter should be determined from XLALSimInspiralTaylorF2Reduc edSpinComputeChi.
*/ */
int XLALSimInspiralTaylorF2ReducedSpinTidal( int XLALSimInspiralTaylorF2ReducedSpinTidal(
COMPLEX16FrequencySeries **htilde, /**< FD waveform */ COMPLEX16FrequencySeries **htilde, /**< FD waveform */
const REAL8 phiStart, /**< initial GW phase (rad) const REAL8 phic, /**< orbital coalescence phase (rad
*/ ) */
const REAL8 deltaF, /**< frequency resolution * const REAL8 deltaF, /**< frequency resolution (Hz) */
/ const REAL8 m1_SI, /**< mass of companion 1 (kg) */
const REAL8 m1_SI, /**< mass of companion 1 (k const REAL8 m2_SI, /**< mass of companion 2 (kg) */
g) */ const REAL8 chi, /**< dimensionless aligned-spin par
const REAL8 m2_SI, /**< mass of companion 2 (k am */
g) */ const REAL8 lam1, /**< dimensionless deformability of
const REAL8 chi, /**< dimensionless aligned- 1 */
spin param */ const REAL8 lam2, /**< dimensionless deformability of
const REAL8 lam1, /**< dimensionless deformab 2 */
ility of 1 */ const REAL8 fStart, /**< start GW frequency (Hz) */
const REAL8 lam2, /**< dimensionless deformab const REAL8 r, /**< distance of source (m) */
ility of 2 */ const INT4 phaseO, /**< twice PN phase order */
const REAL8 fStart, /**< start GW frequency (Hz const INT4 ampO /**< twice PN amplitude order */
) */
const REAL8 r, /**< distance of source (m)
*/
const UINT4 phaseO, /**< twice PN phase order *
/
const UINT4 ampO /**< twice PN amplitude ord
er */
); );
/** /**
* Compute the chirp time of the \"reduced-spin\" templates, described in * Compute the chirp time of the \"reduced-spin\" templates, described in
* http://arxiv.org/abs/1107.1267. * http://arxiv.org/abs/1107.1267.
*/ */
REAL8 XLALSimInspiralTaylorF2ReducedSpinChirpTime( REAL8 XLALSimInspiralTaylorF2ReducedSpinChirpTime(
const REAL8 fStart, /**< start GW frequency (Hz) */ const REAL8 fStart, /**< start GW frequency (Hz) */
const REAL8 m1_SI, /**< mass of companion 1 (kg) */ const REAL8 m1_SI, /**< mass of companion 1 (kg) */
const REAL8 m2_SI, /**< mass of companion 2 (kg) */ const REAL8 m2_SI, /**< mass of companion 2 (kg) */
const REAL8 chi, /**< dimensionless aligned-spin param * / const REAL8 chi, /**< dimensionless aligned-spin param * /
const UINT4 O /**< twice PN phase order */ const INT4 O /**< twice PN phase order */
); );
/** /**
* Compute the dimensionless, spin-aligned parameter chi as used in the * Compute the dimensionless, spin-aligned parameter chi as used in the
* TaylorF2RedSpin waveform. This is different from chi in IMRPhenomB! * TaylorF2RedSpin waveform. This is different from chi in IMRPhenomB!
* Reference: http://arxiv.org/abs/1107.1267, paragraph 3. * Reference: http://arxiv.org/abs/1107.1267, paragraph 3.
*/ */
REAL8 XLALSimInspiralTaylorF2ReducedSpinComputeChi( REAL8 XLALSimInspiralTaylorF2ReducedSpinComputeChi(
const REAL8 m1, /**< mass of companion 1 */ const REAL8 m1, /**< mass of companion 1 */
const REAL8 m2, /**< mass of companion 2 */ const REAL8 m2, /**< mass of companion 2 */
skipping to change at line 1207 skipping to change at line 1927
REAL8 *gamma11, /**< template metric coeff. 11 in mChirp-eta-chi */ REAL8 *gamma11, /**< template metric coeff. 11 in mChirp-eta-chi */
REAL8 *gamma12, /**< template metric coeff. 12/21 in mChirp-eta-chi */ REAL8 *gamma12, /**< template metric coeff. 12/21 in mChirp-eta-chi */
REAL8 *gamma22, /**< template metric coeff. 22 in mChirp-eta-chi */ REAL8 *gamma22, /**< template metric coeff. 22 in mChirp-eta-chi */
const REAL8 mc, /**< chirp mass (in solar mass) */ const REAL8 mc, /**< chirp mass (in solar mass) */
const REAL8 eta, /**< symmetric mass ratio */ const REAL8 eta, /**< symmetric mass ratio */
const REAL8 chi, /**< reduced-spin parameter */ const REAL8 chi, /**< reduced-spin parameter */
const REAL8 fLow, /**< low-frequency cutoff (Hz) */ const REAL8 fLow, /**< low-frequency cutoff (Hz) */
const REAL8FrequencySeries *Sh const REAL8FrequencySeries *Sh
); );
/**
* Compute the Fisher information matrix of "reduced-spin" PN templates in
* theta0, theta3, theta3s, t0, phi0 parameter space, for an SNR=1/sqrt(2)
signal.
*/
gsl_matrix *XLALSimInspiralTaylorF2RedSpinFisherMatrixChirpTimes(
const REAL8 theta0, /**< dimensionless parameter related to the chi
rp time by theta0 = 2 pi fLow tau0 */
const REAL8 theta3, /**< dimensionless parameter related to the chi
rp time by theta3 = -2 pi fLow tau3 */
const REAL8 theta3s, /**< dimensionless parameter related to the chi
rp time by theta3s = 2 pi fLow tau3s */
const REAL8 fLow, /**< low-frequency cutoff (Hz) */
const REAL8 df, /**< frequency resolution of the noise moment v
ectors (Hz) */
REAL8Vector *momI_0, /**< noise moments: momI_0(f) = \int_f0^f (f'/
f0)^{(0-17)/3} df' */
REAL8Vector *momI_2, /**< noise moments: momI_2(f) = \int_f0^f (f'/
f0)^{(2-17)/3} df' */
REAL8Vector *momI_3, /**< noise moments: momI_3(f) = \int_f0^f (f'/
f0)^{(3-17)/3} df' */
REAL8Vector *momI_4, /**< noise moments: momI_4(f) = \int_f0^f (f'/
f0)^{(4-17)/3} df' */
REAL8Vector *momI_5, /**< noise moments: momI_5(f) = \int_f0^f (f'/
f0)^{(5-17)/3} df' */
REAL8Vector *momI_6, /**< noise moments: momI_6(f) = \int_f0^f (f'/
f0)^{(6-17)/3} df' */
REAL8Vector *momI_7, /**< noise moments: momI_7(f) = \int_f0^f (f'/
f0)^{(7-17)/3} df' */
REAL8Vector *momI_8, /**< noise moments: momI_8(f) = \int_f0^f (f'/
f0)^{(8-17)/3} df' */
REAL8Vector *momI_9, /**< noise moments: momI_9(f) = \int_f0^f (f'/
f0)^{(9-17)/3} df' */
REAL8Vector *momI_10, /**< noise moments: momI_10(f) = \int_f0^f (f'
/f0)^{(10-17)/3} df' */
REAL8Vector *momI_11, /**< noise moments: momI_11(f) = \int_f0^f (f'
/f0)^{(11-17)/3} df' */
REAL8Vector *momI_12, /**< noise moments: momI_12(f) = \int_f0^f (f'
/f0)^{(12-17)/3} df' */
REAL8Vector *momI_13, /**< noise moments: momI_13(f) = \int_f0^f (f'
/f0)^{(13-17)/3} df' */
REAL8Vector *momI_14, /**< noise moments: momI_14(f) = \int_f0^f (f'
/f0)^{(14-17)/3} df' */
REAL8Vector *momI_15, /**< noise moments: momI_15(f) = \int_f0^f (f'
/f0)^{(15-17)/3} df' */
REAL8Vector *momI_16, /**< noise moments: momI_16(f) = \int_f0^f (f'
/f0)^{(16-17)/3} df' */
REAL8Vector *momJ_5, /**< noise moments: momJ_5(f) = \int_f0^f (f'/
f0)^{(5-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_6, /**< noise moments: momJ_6(f) = \int_f0^f (f'/
f0)^{(6-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_7, /**< noise moments: momJ_7(f) = \int_f0^f (f'/
f0)^{(7-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_8, /**< noise moments: momJ_8(f) = \int_f0^f (f'/
f0)^{(8-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_9, /**< noise moments: momJ_9(f) = \int_f0^f (f'/
f0)^{(9-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_10, /**< noise moments: momJ_10(f) = \int_f0^f (f'
/f0)^{(10-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_11, /**< noise moments: momJ_11(f) = \int_f0^f (f'
/f0)^{(11-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_12, /**< noise moments: momJ_12(f) = \int_f0^f (f'
/f0)^{(12-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_13, /**< noise moments: momJ_13(f) = \int_f0^f (f'
/f0)^{(13-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_14, /**< noise moments: momJ_14(f) = \int_f0^f (f'
/f0)^{(14-17)/3} log(f'/f0) df' */
REAL8Vector *momK_10, /**< noise moments: momK_14(f) = \int_f0^f (f'
/f0)^{(14-17)/3} log(f'/f0) log(f'/f0) df' */
REAL8Vector *momK_11, /**< noise moments: momK_15(f) = \int_f0^f (f'
/f0)^{(15-17)/3} log(f'/f0) log(f'/f0) df' */
REAL8Vector *momK_12 /**< noise moments: momK_16(f) = \int_f0^f (f'
/f0)^{(16-17)/3} log(f'/f0) log(f'/f0) df' */
);
/**
* Compute the template-space metric of "reduced-spin" PN templates in
* theta0, theta3, theta3s parameter space.
*/
int XLALSimInspiralTaylorF2RedSpinMetricChirpTimes(
REAL8 *gamma00, /**< template metric coeff. 00 in theta0-theta3
-theta3s*/
REAL8 *gamma01, /**< template metric coeff. 01/10 in theta0-the
ta3-theta3s */
REAL8 *gamma02, /**< template metric coeff. 02/20 in theta0-the
ta3-theta3s */
REAL8 *gamma11, /**< template metric coeff. 11 in theta0-theta3
-theta3s */
REAL8 *gamma12, /**< template metric coeff. 12/21 in theta0-the
ta3-theta3s */
REAL8 *gamma22, /**< template metric coeff. 22 in theta0-theta3
-theta3s */
const REAL8 theta0, /**< dimensionless parameter related to the chi
rp time by theta0 = 2 pi fLow tau0 */
const REAL8 theta3, /**< dimensionless parameter related to the chi
rp time by theta3 = -2 pi fLow tau3 */
const REAL8 theta3s, /**< dimensionless parameter related to the chi
rp time by theta3s = 2 pi fLow tau3s */
const REAL8 fLow, /**< low-frequency cutoff (Hz) */
const REAL8 df, /**< frequency resolution of the noise moment v
ectors (Hz) */
REAL8Vector *momI_0, /**< noise moments: momI_0(f) = \int_f0^f (f'/
f0)^{(0-17)/3} df' */
REAL8Vector *momI_2, /**< noise moments: momI_2(f) = \int_f0^f (f'/
f0)^{(2-17)/3} df' */
REAL8Vector *momI_3, /**< noise moments: momI_3(f) = \int_f0^f (f'/
f0)^{(3-17)/3} df' */
REAL8Vector *momI_4, /**< noise moments: momI_4(f) = \int_f0^f (f'/
f0)^{(4-17)/3} df' */
REAL8Vector *momI_5, /**< noise moments: momI_5(f) = \int_f0^f (f'/
f0)^{(5-17)/3} df' */
REAL8Vector *momI_6, /**< noise moments: momI_6(f) = \int_f0^f (f'/
f0)^{(6-17)/3} df' */
REAL8Vector *momI_7, /**< noise moments: momI_7(f) = \int_f0^f (f'/
f0)^{(7-17)/3} df' */
REAL8Vector *momI_8, /**< noise moments: momI_8(f) = \int_f0^f (f'/
f0)^{(8-17)/3} df' */
REAL8Vector *momI_9, /**< noise moments: momI_9(f) = \int_f0^f (f'/
f0)^{(9-17)/3} df' */
REAL8Vector *momI_10, /**< noise moments: momI_10(f) = \int_f0^f (f'
/f0)^{(10-17)/3} df' */
REAL8Vector *momI_11, /**< noise moments: momI_11(f) = \int_f0^f (f'
/f0)^{(11-17)/3} df' */
REAL8Vector *momI_12, /**< noise moments: momI_12(f) = \int_f0^f (f'
/f0)^{(12-17)/3} df' */
REAL8Vector *momI_13, /**< noise moments: momI_13(f) = \int_f0^f (f'
/f0)^{(13-17)/3} df' */
REAL8Vector *momI_14, /**< noise moments: momI_14(f) = \int_f0^f (f'
/f0)^{(14-17)/3} df' */
REAL8Vector *momI_15, /**< noise moments: momI_15(f) = \int_f0^f (f'
/f0)^{(15-17)/3} df' */
REAL8Vector *momI_16, /**< noise moments: momI_16(f) = \int_f0^f (f'
/f0)^{(16-17)/3} df' */
REAL8Vector *momJ_5, /**< noise moments: momJ_5(f) = \int_f0^f (f'/
f0)^{(5-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_6, /**< noise moments: momJ_6(f) = \int_f0^f (f'/
f0)^{(6-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_7, /**< noise moments: momJ_7(f) = \int_f0^f (f'/
f0)^{(7-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_8, /**< noise moments: momJ_8(f) = \int_f0^f (f'/
f0)^{(8-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_9, /**< noise moments: momJ_9(f) = \int_f0^f (f'/
f0)^{(9-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_10, /**< noise moments: momJ_10(f) = \int_f0^f (f'
/f0)^{(10-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_11, /**< noise moments: momJ_11(f) = \int_f0^f (f'
/f0)^{(11-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_12, /**< noise moments: momJ_12(f) = \int_f0^f (f'
/f0)^{(12-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_13, /**< noise moments: momJ_13(f) = \int_f0^f (f'
/f0)^{(13-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_14, /**< noise moments: momJ_14(f) = \int_f0^f (f'
/f0)^{(14-17)/3} log(f'/f0) df' */
REAL8Vector *momK_10, /**< noise moments: momK_14(f) = \int_f0^f (f'
/f0)^{(14-17)/3} log(f'/f0) log(f'/f0) df' */
REAL8Vector *momK_11, /**< noise moments: momK_15(f) = \int_f0^f (f'
/f0)^{(15-17)/3} log(f'/f0) log(f'/f0) df' */
REAL8Vector *momK_12 /**< noise moments: momK_16(f) = \int_f0^f (f'
/f0)^{(16-17)/3} log(f'/f0) log(f'/f0) df' */
);
int XLALSimInspiralTaylorF2RedSpinComputeNoiseMoments(
REAL8Vector *momI_0, /**< noise moments: momI_0(f) = \int_f0^f (f'/f
0)^{(0-17)/3} df' */
REAL8Vector *momI_2, /**< noise moments: momI_2(f) = \int_f0^f (f'/f
0)^{(2-17)/3} df' */
REAL8Vector *momI_3, /**< noise moments: momI_3(f) = \int_f0^f (f'/f
0)^{(3-17)/3} df' */
REAL8Vector *momI_4, /**< noise moments: momI_4(f) = \int_f0^f (f'/f
0)^{(4-17)/3} df' */
REAL8Vector *momI_5, /**< noise moments: momI_5(f) = \int_f0^f (f'/f
0)^{(5-17)/3} df' */
REAL8Vector *momI_6, /**< noise moments: momI_6(f) = \int_f0^f (f'/f
0)^{(6-17)/3} df' */
REAL8Vector *momI_7, /**< noise moments: momI_7(f) = \int_f0^f (f'/f
0)^{(7-17)/3} df' */
REAL8Vector *momI_8, /**< noise moments: momI_8(f) = \int_f0^f (f'/f
0)^{(8-17)/3} df' */
REAL8Vector *momI_9, /**< noise moments: momI_9(f) = \int_f0^f (f'/f
0)^{(9-17)/3} df' */
REAL8Vector *momI_10, /**< noise moments: momI_10(f) = \int_f0^f (f'/
f0)^{(10-17)/3} df' */
REAL8Vector *momI_11, /**< noise moments: momI_11(f) = \int_f0^f (f'/
f0)^{(11-17)/3} df' */
REAL8Vector *momI_12, /**< noise moments: momI_12(f) = \int_f0^f (f'/
f0)^{(12-17)/3} df' */
REAL8Vector *momI_13, /**< noise moments: momI_13(f) = \int_f0^f (f'/
f0)^{(13-17)/3} df' */
REAL8Vector *momI_14, /**< noise moments: momI_14(f) = \int_f0^f (f'/
f0)^{(14-17)/3} df' */
REAL8Vector *momI_15, /**< noise moments: momI_15(f) = \int_f0^f (f'/
f0)^{(15-17)/3} df' */
REAL8Vector *momI_16, /**< noise moments: momI_16(f) = \int_f0^f (f'/
f0)^{(16-17)/3} df' */
REAL8Vector *momJ_5, /**< noise moments: momJ_5(f) = \int_f0^f (f'/f
0)^{(5-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_6, /**< noise moments: momJ_6(f) = \int_f0^f (f'/f
0)^{(6-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_7, /**< noise moments: momJ_7(f) = \int_f0^f (f'/f
0)^{(7-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_8, /**< noise moments: momJ_8(f) = \int_f0^f (f'/f
0)^{(8-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_9, /**< noise moments: momJ_9(f) = \int_f0^f (f'/f
0)^{(9-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_10, /**< noise moments: momJ_10(f) = \int_f0^f (f'/
f0)^{(10-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_11, /**< noise moments: momJ_11(f) = \int_f0^f (f'/
f0)^{(11-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_12, /**< noise moments: momJ_12(f) = \int_f0^f (f'/
f0)^{(12-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_13, /**< noise moments: momJ_13(f) = \int_f0^f (f'/
f0)^{(13-17)/3} log(f'/f0) df' */
REAL8Vector *momJ_14, /**< noise moments: momJ_14(f) = \int_f0^f (f'/
f0)^{(14-17)/3} log(f'/f0) df' */
REAL8Vector *momK_10, /**< noise moments: momK_10(f) = \int_f0^f (f'/
f0)^{(10-17)/3} log(f'/f0) log(f'/f0) df' */
REAL8Vector *momK_11, /**< noise moments: momK_11(f) = \int_f0^f (f'/
f0)^{(11-17)/3} log(f'/f0) log(f'/f0) df' */
REAL8Vector *momK_12, /**< noise moments: momK_12(f) = \int_f0^f (f'/
f0)^{(12-17)/3} log(f'/f0) log(f'/f0) df' */
REAL8Vector *Sh, /**< one sided PSD of the detector noise: Sh(f
) for f = [fLow, fNyq] */
REAL8 fLow, /**< low frequency cutoff (Hz) */
REAL8 df
);
/* compute theta0, theta3, theta3s from mc, eta, chi */
void XLALSimInspiralTaylorF2RedSpinChirpTimesFromMchirpEtaChi(
double *theta0, /**< dimensionless parameter related to the chirp time
by theta0 = 2 pi fLow tau0 */
double *theta3, /**< dimensionless parameter related to the chirp time
by theta3 = -2 pi fLow tau3 */
double *theta3s,/**< dimensionless parameter related to the chirp time
by theta3s = 2 pi fLow tau3s */
double mc, /**< chirp mass (M_sun) */
double eta, /**< symmetric mass ratio */
double chi, /**< reduced-spin parameter */
double fLow /**< low-frequency cutoff (Hz) */
);
/* compute mc, eta, chi from theta0, theta3, theta3s */
void XLALSimInspiralTaylorF2RedSpinMchirpEtaChiFromChirpTimes(
double *mc, /**< chirp mass (M_sun) */
double *eta, /**< symmetric mass ratio */
double *chi, /**< reduced-spin parameter */
double theta0, /**< dimensionless parameter related to the chirp time
by theta0 = 2 pi fLow tau0 */
double theta3, /**< dimensionless parameter related to the chirp time
by theta3 = -2 pi fLow tau3 */
double theta3s, /**< dimensionless parameter related to the chirp time
by theta3s = 2 pi fLow tau3s */
double fLow /**< low-frequency cutoff (Hz) */
);
#if 0 #if 0
{ /* so that editors will match succeeding brace */ { /* so that editors will match succeeding brace */
#elif defined(__cplusplus) #elif defined(__cplusplus)
} }
#endif #endif
#endif /* _LALSIMINSPIRAL_H */ #endif /* _LALSIMINSPIRAL_H */
 End of changes. 66 change blocks. 
462 lines changed or deleted 1739 lines changed or added


 LALSimNoise.h   LALSimNoise.h 
skipping to change at line 33 skipping to change at line 33
#include <stddef.h> #include <stddef.h>
#include <lal/LALDatatypes.h> #include <lal/LALDatatypes.h>
#include <gsl/gsl_rng.h> #include <gsl/gsl_rng.h>
#if defined(__cplusplus) #if defined(__cplusplus)
extern "C" { extern "C" {
#elif 0 #elif 0
} /* so that editors will match preceding brace */ } /* so that editors will match preceding brace */
#endif #endif
#ifdef SWIG // SWIG interface directives
SWIGLAL(FUNCTION_POINTER(XLALSimNoisePSDiLIGOSRD, XLALSimNoisePSDiLIGOModel
,
XLALSimNoisePSDeLIGOModel, XLALSimNoisePSDVirgo, XLALSimNoisePSDGEO,
XLALSimNoisePSDTAMA, XLALSimNoisePSDaLIGONoSRMLowPower,
XLALSimNoisePSDaLIGONoSRMHighPower,
XLALSimNoisePSDaLIGOZeroDetLowPower,
XLALSimNoisePSDaLIGOZeroDetHighPower, XLALSimNoisePSDaLIGONSNSOpt,
XLALSimNoisePSDaLIGOBHBH20Deg, XLALSimNoisePSDaLIGOHighFrequency,
XLALSimNoisePSDKAGRA, XLALSimNoisePSDAdvVirgo));
#endif
/* /*
* *
* *
* PSD GENERATION FUNCTIONS * PSD GENERATION FUNCTIONS
* *
* *
*/ */
/* /*
* *
 End of changes. 1 change blocks. 
0 lines changed or deleted 12 lines changed or added


 LALSimulation.h   LALSimulation.h 
skipping to change at line 32 skipping to change at line 32
#include <lal/LALDatatypes.h> #include <lal/LALDatatypes.h>
#include <lal/LALDetectors.h> #include <lal/LALDetectors.h>
#if defined(__cplusplus) #if defined(__cplusplus)
extern "C" { extern "C" {
#elif 0 #elif 0
} /* so that editors will match preceding brace */ } /* so that editors will match preceding brace */
#endif #endif
const LALDetector *XLALInstrumentNameToLALDetector( const LALDetector *XLALDetectorPrefixToLALDetector(
const char *string const char *string
); );
/* FIXME: compatibility wrapper. remove when not needed */
const LALDetector *XLALInstrumentNameToLALDetector(const char *string);
#ifdef SWIG
%newobject XLALSimDetectorStrainREAL8TimeSeries;
#endif
REAL8TimeSeries *XLALSimDetectorStrainREAL8TimeSeries( REAL8TimeSeries *XLALSimDetectorStrainREAL8TimeSeries(
const REAL8TimeSeries *hplus, const REAL8TimeSeries *hplus,
const REAL8TimeSeries *hcross, const REAL8TimeSeries *hcross,
REAL8 right_ascension, REAL8 right_ascension,
REAL8 declination, REAL8 declination,
REAL8 psi, REAL8 psi,
LALDetector *detector LALDetector *detector
); );
int XLALSimAddInjectionREAL8TimeSeries( int XLALSimAddInjectionREAL8TimeSeries(
REAL8TimeSeries *target, REAL8TimeSeries *target,
REAL8TimeSeries *h, REAL8TimeSeries *h,
const COMPLEX16FrequencySeries *response const COMPLEX16FrequencySeries *response
); );
int XLALSimAddInjectionREAL4TimeSeries(
REAL4TimeSeries *target,
REAL4TimeSeries *h,
const COMPLEX8FrequencySeries *response
);
#if 0 #if 0
{ /* so that editors will match succeeding brace */ { /* so that editors will match succeeding brace */
#elif defined(__cplusplus) #elif defined(__cplusplus)
} }
#endif #endif
#endif /* _LALSIMULATION_H */ #endif /* _LALSIMULATION_H */
 End of changes. 4 change blocks. 
4 lines changed or deleted 8 lines changed or added


 LALSimulationConfig.h   LALSimulationConfig.h 
/* src/LALSimulationConfig.h. Generated from LALSimulationConfig.h.in by c onfigure. */ /* src/LALSimulationConfig.h. Generated from LALSimulationConfig.h.in by c onfigure. */
/* only include this file if LALSimulation's config.h has not been included */ /* only include this file if LALSimulation's config.h has not been included */
#ifndef _LALSIMULATIONCONFIG_H #ifndef _LALSIMULATIONCONFIG_H
#define _LALSIMULATIONCONFIG_H #define _LALSIMULATIONCONFIG_H
/* LALFrame Version */ /* LALFrame Version */
#define LALSIMULATION_VERSION "0.3.0" #define LALSIMULATION_VERSION "0.6.0"
/* LALFrame Version Major Number */ /* LALFrame Version Major Number */
#define LALSIMULATION_VERSION_MAJOR 0 #define LALSIMULATION_VERSION_MAJOR 0
/* LALFrame Version Minor Number */ /* LALFrame Version Minor Number */
#define LALSIMULATION_VERSION_MINOR 3 #define LALSIMULATION_VERSION_MINOR 6
/* LALFrame Version Micro Number */ /* LALFrame Version Micro Number */
#define LALSIMULATION_VERSION_MICRO 0 #define LALSIMULATION_VERSION_MICRO 0
/* LALFrame Version Devel Number */ /* LALFrame Version Devel Number */
#define LALSIMULATION_VERSION_DEVEL 0 #define LALSIMULATION_VERSION_DEVEL 0
/* LALFrame Configure Arguments */ /* LALFrame Configure Arguments */
#define LALSIMULATION_CONFIGURE_ARGS " '--enable-shared' '--prefix=/home/an drey/upstream-tracker/testing/lalsimulation/0.3.0' '--exec-prefix=/home/and rey/upstream-tracker/testing/lalsimulation/0.3.0' '--sysconfdir=/home/andre y/upstream-tracker/testing/lalsimulation/0.3.0' '--datadir=/home/andrey/ups tream-tracker/testing/lalsimulation/0.3.0' 'CFLAGS=-w -fpermissive' 'CXXFLA GS=-w -fpermissive'" #define LALSIMULATION_CONFIGURE_ARGS " '--enable-shared' '--prefix=/home/ut /testing/lalsimulation/0.6.0' '--exec-prefix=/home/ut/testing/lalsimulation /0.6.0' '--sysconfdir=/home/ut/testing/lalsimulation/0.6.0' '--datadir=/hom e/ut/testing/lalsimulation/0.6.0/share' 'CFLAGS=-w -fpermissive' 'CXXFLAGS= -w -fpermissive'"
/* LALFrame Configure Date */ /* LALFrame Configure Date */
#define LALSIMULATION_CONFIGURE_DATE "2012-05-12T12:31:46+0400" #define LALSIMULATION_CONFIGURE_DATE "2012-12-28T13:09:00+0400"
#endif /* _LALSIMULATIONCONFIG_H */ #endif /* _LALSIMULATIONCONFIG_H */
 End of changes. 4 change blocks. 
4 lines changed or deleted 4 lines changed or added


 LALSimulationVCSInfo.h   LALSimulationVCSInfo.h 
/* /*
* LALSimulationVCSInfo.h.in - LALSimulation VCS Information Header * LALSimulationVCSInfo.h - LALSimulation VCS Information Header
* *
* This program is free software; you can redistribute it and/or modify * This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by * it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or (at * the Free Software Foundation; either version 2 of the License, or (at
* your option) any later version. * your option) any later version.
* *
* This program is distributed in the hope that it will be useful, but * This program is distributed in the hope that it will be useful, but
* WITHOUT ANY WARRANTY; without even the implied warranty of * WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* General Public License for more details. * General Public License for more details.
skipping to change at line 26 skipping to change at line 26
* Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, * Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston,
* MA 02111-1307 USA * MA 02111-1307 USA
* *
* Copyright (C) 2009-2011 Adam Mercer * Copyright (C) 2009-2011 Adam Mercer
*/ */
#ifndef _LALSIMULATIONVCSINFO_H #ifndef _LALSIMULATIONVCSINFO_H
#define _LALSIMULATIONVCSINFO_H #define _LALSIMULATIONVCSINFO_H
#include <lal/LALVCSInfo.h> #include <lal/LALVCSInfo.h>
#include <lal/LALSimulationConfig.h> #include <lal/LALSimulationLibVCSInfo.h>
#ifdef __cplusplus #ifdef __cplusplus
extern "C" { extern "C" {
#endif #endif
/* vcs information defines */ /* vcs information defines */
#define LALSIMULATION_VCS_NAME "LALSimulation" #define LALSIMULATION_VCS_NAME "LALSimulation"
#define LALSIMULATION_VCS_ID "f93dc7b9184083729ad83b6e60b202dc57471948" #define LALSIMULATION_VCS_ID "6fb2f15038f3af147ffc456268b5d60f2b7d2c03"
#define LALSIMULATION_VCS_DATE "2012-05-11 16:55:41 +0000" #define LALSIMULATION_VCS_DATE "2012-12-11 17:06:25 +0000"
#define LALSIMULATION_VCS_BRANCH "None" #define LALSIMULATION_VCS_BRANCH "None"
#define LALSIMULATION_VCS_TAG "lalsimulation-v0.3.0" #define LALSIMULATION_VCS_TAG "lalsimulation-v0.6.0"
#define LALSIMULATION_VCS_AUTHOR "Adam Mercer <adam.mercer@ligo.org>" #define LALSIMULATION_VCS_AUTHOR "Adam Mercer <adam.mercer@ligo.org>"
#define LALSIMULATION_VCS_COMMITTER "Adam Mercer <adam.mercer@ligo.org>" #define LALSIMULATION_VCS_COMMITTER "Adam Mercer <adam.mercer@ligo.org>"
#define LALSIMULATION_VCS_STATUS "CLEAN: All modifications committed" #define LALSIMULATION_VCS_STATUS "CLEAN: All modifications committed"
/* vcs information defines - identable*/ /* vcs information defines - identable*/
#define LALSIMULATION_VCS_IDENT_ID "$" "LALSimulationId: f93dc7b9184083729a #define LALSIMULATION_VCS_IDENT_ID "$" "LALSimulationId: 6fb2f15038f3af147f
d83b6e60b202dc57471948 " "$" fc456268b5d60f2b7d2c03 " "$"
#define LALSIMULATION_VCS_IDENT_DATE "$" "LALSimulationDate: 2012-05-11 16: #define LALSIMULATION_VCS_IDENT_DATE "$" "LALSimulationDate: 2012-12-11 17:
55:41 +0000 " "$" 06:25 +0000 " "$"
#define LALSIMULATION_VCS_IDENT_BRANCH "$" "LALSimulationBranch: None " "$" #define LALSIMULATION_VCS_IDENT_BRANCH "$" "LALSimulationBranch: None " "$"
#define LALSIMULATION_VCS_IDENT_TAG "$" "LALSimulationTag: lalsimulation-v0 .3.0 " "$" #define LALSIMULATION_VCS_IDENT_TAG "$" "LALSimulationTag: lalsimulation-v0 .6.0 " "$"
#define LALSIMULATION_VCS_IDENT_AUTHOR "$" "LALSimulationAuthor: Adam Merce r <adam.mercer@ligo.org> " "$" #define LALSIMULATION_VCS_IDENT_AUTHOR "$" "LALSimulationAuthor: Adam Merce r <adam.mercer@ligo.org> " "$"
#define LALSIMULATION_VCS_IDENT_COMMITTER "$" "LALSimulationCommitter: Adam Mercer <adam.mercer@ligo.org> " "$" #define LALSIMULATION_VCS_IDENT_COMMITTER "$" "LALSimulationCommitter: Adam Mercer <adam.mercer@ligo.org> " "$"
#define LALSIMULATION_VCS_IDENT_STATUS "$" "LALSimulationStatus: CLEAN: All modifications committed " "$" #define LALSIMULATION_VCS_IDENT_STATUS "$" "LALSimulationStatus: CLEAN: All modifications committed " "$"
/* vcs information structures */ /* vcs information structures */
extern struct tagLALVCSInfo lalSimulationVCSInfo;
static const struct tagLALVCSInfo lalSimulationHeaderVCSInfo = { \ static const struct tagLALVCSInfo lalSimulationHeaderVCSInfo = { \
LALSIMULATION_VCS_NAME, \ LALSIMULATION_VCS_NAME, \
LALSIMULATION_VERSION, \ LALSIMULATION_VERSION, \
LALSIMULATION_VCS_ID, \ LALSIMULATION_VCS_ID, \
LALSIMULATION_VCS_DATE, \ LALSIMULATION_VCS_DATE, \
LALSIMULATION_VCS_BRANCH, \ LALSIMULATION_VCS_BRANCH, \
LALSIMULATION_VCS_TAG, \ LALSIMULATION_VCS_TAG, \
LALSIMULATION_VCS_AUTHOR, \ LALSIMULATION_VCS_AUTHOR, \
LALSIMULATION_VCS_COMMITTER, \ LALSIMULATION_VCS_COMMITTER, \
LALSIMULATION_VCS_STATUS \ LALSIMULATION_VCS_STATUS \
 End of changes. 7 change blocks. 
11 lines changed or deleted 10 lines changed or added

This html diff was produced by rfcdiff 1.41. The latest version is available from http://tools.ietf.org/tools/rfcdiff/