NewGetAMCoeffsTest.c File Reference

Author:

John Whelan (based on code by Reinhard Prix) Test for LALNewGetAMCoeffs(): compare results to older, well-tested (but less efficient, and harder to understand) function LALComputeAM()

More...

#include <math.h>
#include <sys/times.h>
#include <lal/ComputeFstat.h>
#include <lal/LALBarycenter.h>
#include <lal/LALInitBarycenter.h>
#include <lal/AVFactories.h>

Include dependency graph for NewGetAMCoeffsTest.c:

Go to the source code of this file.

Defines
#define	RELERROR(x, y)   fabs( 2.0 * ((x) - (y)) / ( (x) + (y) ) )
#define	SQ(x)   ( (x) * (x) )
#define	MYMAX(a, b)   ( (a) > (b) ? (a) : (b) )
#define	ERROR(code, msg, statement)
#define	INFO(statement)
#define	SUB(func, statusptr)
Error codes
#define	NEWGETAMCOEFFSTEST_ENORM   0
#define	NEWGETAMCOEFFSTEST_ESUB   1
#define	NEWGETAMCOEFFSTEST_MSGENORM   "Normal exit"
#define	NEWGETAMCOEFFSTEST_MSGESUB   "Subroutine failed"
Functions
	NRCSID (NEWGETAMCOEFFSTEST,"$Id: NewGetAMCoeffsTest.c,v 1.6 2006/10/27 15:18:23 reinhard Exp $")
int	main (int argc, char *argv[])
	MAIN function of SideBandMCMC code Compute the posterior pdfs of the orbital and nuisance parameters of a binary signal in Fstat form.
Variables
int	lalDebugLevel
	DECLARE AND SET GLOBAL DEBUG LEVEL.

---

Detailed Description

Author:

John Whelan (based on code by Reinhard Prix) Test for LALNewGetAMCoeffs(): compare results to older, well-tested (but less efficient, and harder to understand) function LALComputeAM()

This routine currently compares LALComputeAM(), LALGetAMCoeffs() and LALNewGetAMCoeffs() to each other, by computing average and maximal relative differences between the respective a's and b's.

Detector and sky-location are picked at random each time, which allows a minimal Monte-Carlo validation by simply running this script repeatedly.

Definition in file NewGetAMCoeffsTest.c.

---

Define Documentation

#define NEWGETAMCOEFFSTEST_ENORM   0

Definition at line 47 of file NewGetAMCoeffsTest.c.

#define NEWGETAMCOEFFSTEST_ESUB   1

Definition at line 48 of file NewGetAMCoeffsTest.c.

#define NEWGETAMCOEFFSTEST_MSGENORM   "Normal exit"

Definition at line 50 of file NewGetAMCoeffsTest.c.

#define NEWGETAMCOEFFSTEST_MSGESUB   "Subroutine failed"

Definition at line 51 of file NewGetAMCoeffsTest.c.

#define RELERROR	(	x,
		y		)	fabs( 2.0 * ((x) - (y)) / ( (x) + (y) ) )

Definition at line 55 of file NewGetAMCoeffsTest.c.

#define SQ

(

x

)

( (x) * (x) )

Definition at line 56 of file NewGetAMCoeffsTest.c.

#define MYMAX	(	a,
		b		)	( (a) > (b) ? (a) : (b) )

Definition at line 57 of file NewGetAMCoeffsTest.c.

#define ERROR	(	code,
		msg,
		statement		)

Value:

do {                                                                   if ( lalDebugLevel & LALERROR )                                        LALPrintError( "Error[0] %d: program %s, file %s, line %d, %s\n"                    "        %s %s\n", (code), *argv, __FILE__,                     __LINE__, NEWGETAMCOEFFSTEST, statement ? statement :                     "", (msg) );                                      } while (0)

Definition at line 63 of file NewGetAMCoeffsTest.c.

#define INFO

(

statement

)

Value:

do {                                                                   if ( lalDebugLevel & LALINFO )                                         LALPrintError( "Info[0]: program %s, file %s, line %d, %s\n"                        "        %s\n", *argv, __FILE__, __LINE__,                      NEWGETAMCOEFFSTEST, (statement) );                     } while (0)

Definition at line 72 of file NewGetAMCoeffsTest.c.

#define SUB	(	func,
		statusptr		)

Value:

do {                                                                   if ( (func), (statusptr)->statusCode ) {                               ERROR( NEWGETAMCOEFFSTEST_ESUB, NEWGETAMCOEFFSTEST_MSGESUB,                         "Function call \"" #func "\" failed:" );                      return NEWGETAMCOEFFSTEST_ESUB;                                       }                                                                  } while (0)

Definition at line 80 of file NewGetAMCoeffsTest.c.

---

Function Documentation

NRCSID	(	NEWGETAMCOEFFSTEST	,
		"$Id: NewGetAMCoeffsTest.	c,
		v 1.6 2006/10/27 15:18:23 reinhard Exp $"
	)

int main	(	int	argc,
		char *	argv[]
	)

MAIN function of SideBandMCMC code Compute the posterior pdfs of the orbital and nuisance parameters of a binary signal in Fstat form.

..

TEST NO. 1

TEST NO. 2

TEST NO. 3

TEST NO. 4

Definition at line 100 of file NewGetAMCoeffsTest.c.

---

Variable Documentation

int lalDebugLevel

DECLARE AND SET GLOBAL DEBUG LEVEL.

see the section (currently 7.4.1) of the LSD on "Status-reporting objects" for a list of predefined debug levels

Definition at line 62 of file BCVSpinTemplates.c.

---