Commit 8144dc3e authored by Ian Harry's avatar Ian Harry
Browse files

First set of changes to add TaylorF2 templates to coh_PTF

Original: 60e3d68b4cf822a5148b3a29e423d6efd9944bcd
parent d4ab0d3d
......@@ -416,6 +416,12 @@ void coh_PTF_template (
FindChirpTmpltParams *params
);
void coh_PTF_template_PTF (
FindChirpTemplate *fcTmplt,
InspiralTemplate *InspTmplt,
FindChirpTmpltParams *params
);
void coh_PTF_normalize(
struct coh_PTF_params *params,
FindChirpTemplate *fcTmplt,
......
......@@ -360,6 +360,15 @@ int main( int argc, char **argv )
fcTmpltParams->PTFe2 = XLALCreateVectorSequence( 3, numPoints );
fcTmpltParams->fwdPlan = XLALCreateForwardREAL4FFTPlan( numPoints, 0 );
fcTmpltParams->deltaT = 1.0/params->sampleRate;
fcTmpltParams->xfacVec = XLALCreateVector(numPoints / 2 + 1 );
/* Set the values of xfacVec This is k^(-1/3) */
const REAL4 xfacExponent = -1.0/3.0;
REAL4 *xfac = NULL;
xfac = fcTmpltParams->xfacVec->data;
xfac[0] = 0;
for (ui = 1; ui < fcTmpltParams->xfacVec->length; ++ui)
xfac[ui] = pow( (REAL4) ui, xfacExponent );
for( ifoNumber = 0; ifoNumber < LAL_NUM_IFO; ifoNumber++ )
{
if ( params->haveTrig[ifoNumber] )
......
......@@ -73,8 +73,17 @@ RCSID( "$Id$" );
NRCSID(FINDCHIRPPTFTEMPLATEC, "$Id: FindChirpPTFTemplate.c,v 1.7 2008/06/26 19:05:07 dfazi Exp $");
void coh_PTF_template (
FindChirpTemplate *fcTmplt,
InspiralTemplate *InspTmplt,
FindChirpTmpltParams *params
)
{
coh_PTF_template_PTF (fcTmplt,InspTmplt,params);
}
void
coh_PTF_template (
coh_PTF_template_PTF (
FindChirpTemplate *fcTmplt,
InspiralTemplate *InspTmplt,
FindChirpTmpltParams *params
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment