Loading [MathJax]/extensions/TeX/AMSsymbols.js
LALPulsar 7.1.1.1-b246709
•All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Modules Pages
GeneralMeshTest.c File Reference
Pulsar Metric Routines » Header PtoleMetric.h

Tests and showcases the combination of a LAL metric function (of the user's specification) and Header TwoDMesh.h modules by producing a template mesh. More...

Detailed Description

Tests and showcases the combination of a LAL metric function (of the user's specification) and Header TwoDMesh.h modules by producing a template mesh.

Author
Jones, D. I., Owen, B. J.

Program <tt>GeneralMeshTest.c</tt>

Usage

GeneralMeshTest

Description

The -b option sets the beginning GPS time of integration to the option argument. (Default is \( 731265908 \) seconds, chosen to lie within the S2 run).

The -c option determins the center of the patch. This option is hardcoded to use equatorial coordinates and the argument should be given in hh:mm:ss:dd:mm:ss format. (Default is the center of the globular cluster 47 Tuc).

The -f option sets the maximum frequency of integration (in Hz) to the option argument. (The default value is 1000.)

The -l option determines, for a rectangular search region, the limits in right ascension and declination of the grid. The argument should be given in degrees as RA(min):RA(max):dec(min):dec(max). (The default is the octant of the sky defined by \( 0 < \textrm{RA} < 90 \) and \( 0< \textrm{dec} < 85 \) ; this avoids the coordinate singularity at the poles.) This option automatically overrides whatever is specified by the -r option.

The -m option sets the maximum mismatch of the mesh to the option argument. (Default is 0.02.)

The texttt{-p} option causes the coordinates of the nodes to be written to a file mesh.dat, for the benifit of users who don't have xmgrace installed. The format is one node per line, (RA, DEC), with the angles in degrees.

The -r option sets the radius (in arcminutes) of the circular sky patch. If you specify radius zero you will get a search over a rectangular region whose limits in RA and dec are specified by the -l option. (The default value is the radius of the globular cluster 47 Tuc).

The -t option sets the duration of integration in seconds. (The default is \( 39600 \) seconds \( = 11 \) hours, which is chosen because it is of the right size for S2 analyses).

Algorithm

Uses

lalDebugLevel LALDCreateVector()
LALCheckMemoryLeaks() LALDDestroyVector()
LALProjectMetric() XLALGetEarthTimes()
LALPtoleMetric() XLALInitBarycenter()
LALCreateTwoDMesh() LALDestroyTwoDMesh()
XLALGetEarthTimes
int XLALGetEarthTimes(const LIGOTimeGPS *tepoch, REAL8 *tMidnight, REAL8 *tAutumn)
This function takes a GPS time from tepoch and uses it to assign tAutumn and tMidnight,...
Definition: GetEarthTimes.c:76
XLALInitBarycenter
EphemerisData * XLALInitBarycenter(const CHAR *earthEphemerisFile, const CHAR *sunEphemerisFile)
XLAL interface to reading ephemeris files 'earth' and 'sun', and return ephemeris-data in old backwar...
Definition: LALInitBarycenter.c:206
lalDebugLevel
#define lalDebugLevel
LALPtoleMetric
void LALPtoleMetric(LALStatus *status, REAL8Vector *metric, PtoleMetricIn *input)
Computes metric components for a pulsar search in the `‘Ptolemaic’' approximation; both the Earth's s...
Definition: PtoleMetric.c:91
LALProjectMetric
void LALProjectMetric(LALStatus *stat, REAL8Vector *metric, BOOLEAN errors)
Project out the zeroth dimension of a metric.
Definition: PtoleMetric.c:732
LALCreateTwoDMesh
void LALCreateTwoDMesh(LALStatus *stat, TwoDMeshNode **mesh, TwoDMeshParamStruc *params)
Definition: TwoDMesh.c:106
LALDestroyTwoDMesh
void LALDestroyTwoDMesh(LALStatus *stat, TwoDMeshNode **mesh, UINT4 *nFree)
Definition: TwoDMesh.c:159
LALDCreateVector
void LALDCreateVector(LALStatus *, REAL8Vector **, UINT4)
LALDDestroyVector
void LALDDestroyVector(LALStatus *, REAL8Vector **)

Notes

For most regions of parameter space the three metric codes seem to agree well. However, for short (less than one day) runs, they are all capable of returning (unphysical) negative determinant metrics for points very close to the equator.

Definition in file GeneralMeshTest.c.

Go to the source code of this file.

Macros

Error Codes
#define GENERALMESHTESTC_EMEM   1
 
#define GENERALMESHTESTC_ERNG   2
 
#define GENERALMESHTESTC_EFIO   3
 
#define GENERALMESHTESTC_EOPT   4
 
#define GENERALMESHTESTC_EMET   5
 
#define GENERALMESHTESTC_MSGEMEM   "memory (de)allocation error"
 
#define GENERALMESHTESTC_MSGERNG   "value out of range"
 
#define GENERALMESHTESTC_MSGEFIO   "file I/O error"
 
#define GENERALMESHTESTC_MSGEOPT   "unknown command-line option"
 
#define GENERALMESHTESTC_MSGEMET   "determinant of projected metric negative"
 

Macro Definition Documentation

◆ GENERALMESHTESTC_EMEM

#define GENERALMESHTESTC_EMEM   1

Definition at line 100 of file GeneralMeshTest.c.

◆ GENERALMESHTESTC_ERNG

#define GENERALMESHTESTC_ERNG   2

Definition at line 101 of file GeneralMeshTest.c.

◆ GENERALMESHTESTC_EFIO

#define GENERALMESHTESTC_EFIO   3

Definition at line 102 of file GeneralMeshTest.c.

◆ GENERALMESHTESTC_EOPT

#define GENERALMESHTESTC_EOPT   4

Definition at line 103 of file GeneralMeshTest.c.

◆ GENERALMESHTESTC_EMET

#define GENERALMESHTESTC_EMET   5

Definition at line 104 of file GeneralMeshTest.c.

◆ GENERALMESHTESTC_MSGEMEM

#define GENERALMESHTESTC_MSGEMEM   "memory (de)allocation error"

Definition at line 106 of file GeneralMeshTest.c.

◆ GENERALMESHTESTC_MSGERNG

#define GENERALMESHTESTC_MSGERNG   "value out of range"

Definition at line 107 of file GeneralMeshTest.c.

◆ GENERALMESHTESTC_MSGEFIO

#define GENERALMESHTESTC_MSGEFIO   "file I/O error"

Definition at line 108 of file GeneralMeshTest.c.

◆ GENERALMESHTESTC_MSGEOPT

#define GENERALMESHTESTC_MSGEOPT   "unknown command-line option"

Definition at line 109 of file GeneralMeshTest.c.

◆ GENERALMESHTESTC_MSGEMET

#define GENERALMESHTESTC_MSGEMET   "determinant of projected metric negative"

Definition at line 110 of file GeneralMeshTest.c.