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/