LALInspiralPhasing3.c

Go to the documentation of this file.

/*
*  Copyright (C) 2007 David Churches, B.S. Sathyaprakash
*
*  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
*  the Free Software Foundation; either version 2 of the License, or
*  (at your option) any later version.
*
*  This program is distributed in the hope that it will be useful,
*  but WITHOUT ANY WARRANTY; without even the implied warranty of
*  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
*  GNU General Public License for more details.
*
*  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
*  Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston,
*  MA  02111-1307  USA
*/

/*  <lalVerbatim file="LALInspiralPhasing3CV">
Author: Sathyaprakash, B. S.
$Id: LALInspiralPhasing3.c,v 1.4 2007/06/08 14:41:49 bema Exp $
</lalVerbatim>  */

/*  <lalLaTeX>

\subsection{Module \texttt{LALInspiralPhasing3.c}}

The code \texttt{LALInspiralPhasing3.c} calculates the phase the waveform 
from an inspiralling binary system as a function of time up to second post-Nowtonian 
order.

\subsubsection*{Prototypes}

\vspace{0.1in}

\input{LALInspiralPhasing3CP}

\index{\verb&LALInspiralPhasing3()&}
\begin{itemize}
\item {\tt phase:} Output, the phase of the wave at the current epoch.
\item {\tt td:} Input, the PN expansion coefficients of phase $\phi^t_k$ as a function 
of time (cf. Table \ref{table:flux}).
\item {\tt ak:} Input containing PN expansion coefficients.
\end{itemize}


\subsubsection*{Description}
The phase of the inspiral wave corresponding to the {\tt Approximant} {\tt TaylorT2} 
as in Equation~\ref{eq:InspiralWavePhase3}.


\subsubsection*{Algorithm}
None.


\subsubsection*{Uses}
None.

\subsubsection*{Notes}
None.

\vfill{\footnotesize\input{LALInspiralPhasing3CV}}

</lalLaTeX>  */

#include <lal/LALStdlib.h>
#include <lal/LALInspiral.h>

NRCSID (LALINSPIRALPHASING3C, "$Id: LALInspiralPhasing3.c,v 1.4 2007/06/08 14:41:49 bema Exp $");


/*  <lalVerbatim file="LALInspiralPhasing3CP"> */

void 
LALInspiralPhasing3_0PN (
   LALStatus  *status,
   REAL8      *phase,
   REAL8       td,
   expnCoeffs *ak
   ) 
{ /* </lalVerbatim>  */

  REAL8 theta5;

  INITSTATUS (status, "LALInspiralPhasing3_0PN", LALINSPIRALPHASING3C);
  ATTATCHSTATUSPTR(status);
  ASSERT(ak, status, LALINSPIRALH_ENULL, LALINSPIRALH_MSGENULL);

  theta5 = pow(td,-0.625);
  *phase = (ak->ptaN/theta5);
  DETATCHSTATUSPTR(status);
  RETURN(status);
}

/*  <lalVerbatim file="LALInspiralPhasing3CP"> */

void 
LALInspiralPhasing3_2PN (
   LALStatus  *status,
   REAL8      *phase,
   REAL8       td,
   expnCoeffs *ak
   ) 
{ /* </lalVerbatim>  */

  REAL8 theta,theta2,theta5;

  INITSTATUS (status, "LALInspiralPhasing3_2PN", LALINSPIRALPHASING3C);
  ATTATCHSTATUSPTR(status);
  ASSERT(ak, status, LALINSPIRALH_ENULL, LALINSPIRALH_MSGENULL);

  theta = pow(td,-0.125);
  theta2 = theta*theta;
  theta5 = theta2*theta2*theta;

  *phase = (ak->ptaN/theta5) * (1. 
         + ak->pta2*theta2);
  DETATCHSTATUSPTR(status);
  RETURN(status);
}

/*  <lalVerbatim file="LALInspiralPhasing3CP"> */

void 
LALInspiralPhasing3_3PN (
   LALStatus  *status,
   REAL8      *phase,
   REAL8       td,
   expnCoeffs *ak
   ) 
{ /* </lalVerbatim>  */

  REAL8 theta,theta2,theta3,theta5;

  INITSTATUS (status, "LALInspiralPhasing3_3PN", LALINSPIRALPHASING3C);
  ATTATCHSTATUSPTR(status);
  ASSERT(ak, status, LALINSPIRALH_ENULL, LALINSPIRALH_MSGENULL);

  theta = pow(td,-0.125);
  theta2 = theta*theta;
  theta3 = theta2*theta;
  theta5 = theta2*theta3;

  *phase = (ak->ptaN/theta5) * (1. 
         + ak->pta2*theta2 
         + ak->pta3*theta3);
  DETATCHSTATUSPTR(status);
  RETURN(status);
}

/*  <lalVerbatim file="LALInspiralPhasing3CP"> */

void 
LALInspiralPhasing3_4PN (
   LALStatus  *status,
   REAL8      *phase,
   REAL8       td,
   expnCoeffs *ak
   ) 
{ /* </lalVerbatim>  */

  REAL8 theta,theta2,theta3,theta4,theta5;

  INITSTATUS (status, "LALInspiralPhasing3_4PN", LALINSPIRALPHASING3C);
  ATTATCHSTATUSPTR(status);
  ASSERT(ak, status, LALINSPIRALH_ENULL, LALINSPIRALH_MSGENULL);

  theta = pow(td,-0.125);
  theta2 = theta*theta;
  theta3 = theta2*theta;
  theta4 = theta3*theta;
  theta5 = theta4*theta;

  *phase = (ak->ptaN/theta5) * (1. 
         + ak->pta2*theta2 
         + ak->pta3*theta3
         + ak->pta4*theta4);
  DETATCHSTATUSPTR(status);
  RETURN(status);
}

/*  <lalVerbatim file="LALInspiralPhasing3CP"> */

void 
LALInspiralPhasing3_5PN (
   LALStatus  *status,
   REAL8      *phase,
   REAL8       td,
   expnCoeffs *ak
   ) 
{ /* </lalVerbatim>  */

  REAL8 theta,theta2,theta3,theta4,theta5;

  INITSTATUS (status, "LALInspiralPhasing3_5PN", LALINSPIRALPHASING3C);
  ATTATCHSTATUSPTR(status);
  ASSERT(ak, status, LALINSPIRALH_ENULL, LALINSPIRALH_MSGENULL);

  theta = pow(td,-0.125);
  theta2 = theta*theta;
  theta3 = theta2*theta;
  theta4 = theta3*theta;
  theta5 = theta4*theta;

  *phase = (ak->ptaN/theta5) * (1. 
         + ak->pta2*theta2 
         + ak->pta3*theta3
         + ak->pta4*theta4 
         + ak->pta5 * log(td/ak->tn) * theta5);
  DETATCHSTATUSPTR(status);
  RETURN(status);
}

/*  <lalVerbatim file="LALInspiralPhasing3CP"> */

void 
LALInspiralPhasing3_6PN (
   LALStatus  *status,
   REAL8      *phase,
   REAL8       td,
   expnCoeffs *ak
   ) 
{ /* </lalVerbatim>  */

  REAL8 theta,theta2,theta3,theta4,theta5,theta6;

  INITSTATUS (status, "LALInspiralPhasing3_6PN", LALINSPIRALPHASING3C);
  ATTATCHSTATUSPTR(status);
  ASSERT(ak, status, LALINSPIRALH_ENULL, LALINSPIRALH_MSGENULL);

  theta = pow(td,-0.125);
  theta2 = theta*theta;
  theta3 = theta2*theta;
  theta4 = theta3*theta;
  theta5 = theta4*theta;
  theta6 = theta5*theta;

  *phase = (ak->ptaN/theta5) * (1. 
         + ak->pta2*theta2 
         + ak->pta3*theta3
         + ak->pta4*theta4 
         + ak->pta5*log(td/ak->tn)*theta5
         +(ak->ptl6*log(td/256.) + ak->pta6)*theta6);
  DETATCHSTATUSPTR(status);
  RETURN(status);
}

/*  <lalVerbatim file="LALInspiralPhasing3CP"> */

void 
LALInspiralPhasing3_7PN (
   LALStatus  *status,
   REAL8      *phase,
   REAL8       td,
   expnCoeffs *ak
   ) 
{ /* </lalVerbatim>  */

  REAL8 theta,theta2,theta3,theta4,theta5,theta6,theta7;

  INITSTATUS (status, "LALInspiralPhasing3_7PN", LALINSPIRALPHASING3C);
  ATTATCHSTATUSPTR(status);
  ASSERT(ak, status, LALINSPIRALH_ENULL, LALINSPIRALH_MSGENULL);

  theta = pow(td,-0.125);
  theta2 = theta*theta;
  theta3 = theta2*theta;
  theta4 = theta3*theta;
  theta5 = theta4*theta;
  theta6 = theta5*theta;
  theta7 = theta6*theta;

  *phase = (ak->ptaN/theta5) * (1. 
         + ak->pta2*theta2 
         + ak->pta3*theta3
         + ak->pta4*theta4 
         + ak->pta5*log(td/ak->tn)*theta5
         +(ak->ptl6*log(td/256.) + ak->pta6)*theta6
         + ak->pta7*theta7);
  DETATCHSTATUSPTR(status);
  RETURN(status);
}

---