21 #include <lal/LALStdlib.h>
22 #include <lal/LALConstants.h>
23 #include <lal/SkyCoordinates.h>
25 #define LAL_ALPHAGAL (3.366032942)
26 #define LAL_DELTAGAL (0.473477302)
27 #define LAL_LGAL (0.576)
156 REAL8 sinB, cosB, sinL, cosL;
157 REAL8 sinD, sinA, cosA;
178 sinD = cosB*cosDGal*sinL + sinB*sinDGal;
180 cosA = sinB*cosDGal - cosB*sinL*sinDGal;
184 output->latitude = asin( sinD );
204 REAL8 sinD, cosD, sinA, cosA;
205 REAL8 sinB, sinL, cosL;
226 sinB = cosD*cosDGal*cosA + sinD*sinDGal;
227 sinL = sinD*cosDGal - cosD*cosA*sinDGal;
232 output->latitude = asin( sinB );
251 REAL8 sinB, cosB, sinL, cosL;
252 REAL8 sinD, sinA, cosA;
272 sinD = cosB*sinL*sinE + sinB*cosE;
273 sinA = cosB*sinL*cosE - sinB*sinE;
278 output->latitude = asin( sinD );
279 output->longitude = atan2( sinA, cosA );
282 if (
output->longitude < 0.0 )
296 REAL8 sinD, cosD, sinA, cosA;
297 REAL8 sinB, sinL, cosL;
317 sinB = sinD*cosE - cosD*sinA*sinE;
318 sinL = cosD*sinA*cosE + sinD*sinE;
323 output->latitude = asin( sinB );
324 output->longitude = atan2( sinL, cosL );
327 if (
output->longitude < 0.0 )
#define ABORT(statusptr, code, mesg)
#define ASSERT(assertion, statusptr, code, mesg)
#define INITSTATUS(statusptr)
#define RETURN(statusptr)
void LALEquatorialToGalactic(LALStatus *stat, SkyPosition *output, SkyPosition *input)
void LALEquatorialToEcliptic(LALStatus *stat, SkyPosition *output, SkyPosition *input)
void LALEclipticToEquatorial(LALStatus *stat, SkyPosition *output, SkyPosition *input)
void LALGalacticToEquatorial(LALStatus *stat, SkyPosition *output, SkyPosition *input)
#define LAL_IEARTH
Earth inclination (2000), radians.
#define LAL_TWOPI
2*pi is circumference of a circle divided by its radius
double REAL8
Double precision real floating-point number (8 bytes).
#define SKYCOORDINATESH_ENUL
Unexpected null pointer in arguments.
#define SKYCOORDINATESH_ESYS
Wrong coordinate system in input.
@ COORDINATESYSTEM_GALACTIC
The galactic coordinate system.
@ COORDINATESYSTEM_ECLIPTIC
The ecliptic coordinate system.
@ COORDINATESYSTEM_EQUATORIAL
The sky-fixed equatorial coordinate system.
LAL status structure, see The LALStatus structure for more details.
This structure stores the two spherical coordinates of a sky position; ie a generic latitude and long...
REAL8 longitude
The longitudinal coordinate (in radians), as defined above.
REAL8 latitude
The latitudinal coordinate (in radians), as defined above.
CoordinateSystem system
The coordinate system in which latitude/longitude are expressed.
void output(int gps_sec, int output_type)