Loading [MathJax]/extensions/TeX/AMSsymbols.js
LALInspiral 5.0.3.1-ea7c608
All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Modules Pages
LALInspiralTest.c File Reference

Test routine for codes that generate inspiral waveform from non-spinning black hole binaries. More...

Prototypes

void printf_timeseries (int n, float *sig, double delta, double t0)
 
int main (void)
 

Detailed Description

Test routine for codes that generate inspiral waveform from non-spinning black hole binaries.

Author
Sathyaprakash, B. S.

Time domain signals are returned when the approximant is one of TaylorT1, TaylorT2, TaylorT3, PadeT1, EOB, SpinTaylorT3 and frequency domain signals are returned when the approximant is one of TaylorF1, TaylorF2, BCV. This code checks every available approximant at every order and reports whether or not there was any problem with the generation codes.

To generate a waveform first set the InspiralTemplate structure (see below for an example). Next, to measure the length of the array required call the function

void LALInspiralWaveLength(LALStatus *status, UINT4 *n, InspiralTemplate params)

The length will be returned in n. Finally, call the function

void LALInspiralWave(LALStatus *status, REAL4Vector *signalvec, InspiralTemplate *params)

to generate the wave, which will be returned in signal.

Example values of the parameters that can be set (with options in brackets) is:

params.OmegaS = 0.; (Unknown 3PN parameter in energy; shown to be 0 by DJS)
params.Theta = 0.; (Unknown 3PN parameter in flux; arbitrarily set to 0)
params.ieta=1; (1 for comparable masses model, 0 for test mass model)
params.mass1=1.4; (masses of the component stars in solar masses)
params.mass2=1.4;
params.startTime=0.0; (defined so that the instantaneous GW frequency
is params.fLower at params.startTime)
params.startPhase=0.0; (0 to LAL_PI_2)
params.fLower=40.0; (in Hz)
params.fCutoff=1000.0; (in Hz)
params.tSampling=4000.; (in Hz; should be larger than 2 fCutoff or 2 flso,
whichever is smaller)
params.signalAmplitude=1.0;
params.nStartPad=0; (number of leading zero bins)
params.nEndPad=0; (number of trailing zero bins)
params.approximant=TaylorF2; (TaylorT1, PadeT1=ODE solver;
TaylorT2=implicit phasing formula solved in quadrature;
TaylorT3=explicit time-domain phasing;
TaylorF1=stationary phase approx. using ODEs;
TaylorF2=usual stationary phase approx.;
EOB=effective-one-body approach)
params.order=twoPN; (also newtonian, onePN, oneAndHalfPN, twoPN,
twoAndHalfPN, threePN, threeAndHalfPN)
params.massChoice=m1Andm2; (also t0t2, t0t3, t0t4, totalMassAndEta,totalMassAndMu)
params.psi0 = 132250.; (parameter required to generate BCV detection templates)
params.psi3 = -1014.2; (parameter required to generate BCV detection templates)
params.alpha = 0.528; (amplitude correction used in BCV templates)
params.fFinal = 868.7; (frequency at which the BCV template is terminated)
static INT4 test(REAL4, REAL4, REAL4, REAL4, REAL4, REAL4, REAL4, REAL4)
@ totalMassAndMu
total mass and reduced mass
Definition: LALInspiral.h:182
@ totalMassAndEta
total mass and symmetric mass ratio
Definition: LALInspiral.h:180
@ m1Andm2
component masses
Definition: LALInspiral.h:179
EOB
PadeT1
TaylorF1
TaylorT3
TaylorF2
BCV
TaylorT1
TaylorT2

Definition in file LALInspiralTest.c.

Go to the source code of this file.

Function Documentation

◆ printf_timeseries()

void printf_timeseries ( int  n,
float *  sig,
double  delta,
double  t0 
)

Definition at line 89 of file LALInspiralTest.c.

◆ main()

int main ( void  )

Definition at line 105 of file LALInspiralTest.c.