| 1 | #ifndef PALMACDEF
|
|---|
| 2 | #define PALMACDEF
|
|---|
| 3 |
|
|---|
| 4 | /*
|
|---|
| 5 | *+
|
|---|
| 6 | * Name:
|
|---|
| 7 | * palmac.h
|
|---|
| 8 |
|
|---|
| 9 | * Purpose:
|
|---|
| 10 | * Macros used by the PAL library
|
|---|
| 11 |
|
|---|
| 12 | * Language:
|
|---|
| 13 | * Starlink ANSI C
|
|---|
| 14 |
|
|---|
| 15 | * Type of Module:
|
|---|
| 16 | * Include file
|
|---|
| 17 |
|
|---|
| 18 | * Description:
|
|---|
| 19 | * A collection of useful macros provided and used by the PAL library
|
|---|
| 20 |
|
|---|
| 21 | * Authors:
|
|---|
| 22 | * TIMJ: Tim Jenness (JAC, Hawaii)
|
|---|
| 23 | * DSB: David Berry (JAC, Hawaii)
|
|---|
| 24 | * {enter_new_authors_here}
|
|---|
| 25 |
|
|---|
| 26 | * Notes:
|
|---|
| 27 | *
|
|---|
| 28 |
|
|---|
| 29 | * History:
|
|---|
| 30 | * 2012-02-08 (TIMJ):
|
|---|
| 31 | * Initial version.
|
|---|
| 32 | * Adapted with permission from the Fortran SLALIB library.
|
|---|
| 33 | * 2012-04-13 (DSB):
|
|---|
| 34 | * Added PAL__DR2H and PAL__DR2S
|
|---|
| 35 | * {enter_further_changes_here}
|
|---|
| 36 |
|
|---|
| 37 | * Copyright:
|
|---|
| 38 | * Copyright (C) 2012 Science and Technology Facilities Council.
|
|---|
| 39 | * All Rights Reserved.
|
|---|
| 40 |
|
|---|
| 41 | * Licence:
|
|---|
| 42 | * This program is free software: you can redistribute it and/or
|
|---|
| 43 | * modify it under the terms of the GNU Lesser General Public
|
|---|
| 44 | * License as published by the Free Software Foundation, either
|
|---|
| 45 | * version 3 of the License, or (at your option) any later
|
|---|
| 46 | * version.
|
|---|
| 47 | *
|
|---|
| 48 | * This program is distributed in the hope that it will be useful,
|
|---|
| 49 | * but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|---|
| 50 | * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|---|
| 51 | * GNU Lesser General Public License for more details.
|
|---|
| 52 | *
|
|---|
| 53 | * You should have received a copy of the GNU Lesser General
|
|---|
| 54 | * License along with this program. If not, see
|
|---|
| 55 | * <http://www.gnu.org/licenses/>.
|
|---|
| 56 |
|
|---|
| 57 | * Bugs:
|
|---|
| 58 | * {note_any_bugs_here}
|
|---|
| 59 | *-
|
|---|
| 60 | */
|
|---|
| 61 |
|
|---|
| 62 | /* Pi */
|
|---|
| 63 | static const double PAL__DPI = 3.1415926535897932384626433832795028841971693993751;
|
|---|
| 64 |
|
|---|
| 65 | /* 2Pi */
|
|---|
| 66 | static const double PAL__D2PI = 6.2831853071795864769252867665590057683943387987502;
|
|---|
| 67 |
|
|---|
| 68 | /* pi/2: 90 degrees in radians */
|
|---|
| 69 | static const double PAL__DPIBY2 = 1.5707963267948966192313216916397514420985846996876;
|
|---|
| 70 |
|
|---|
| 71 | /* pi/180: degrees to radians */
|
|---|
| 72 | static const double PAL__DD2R = 0.017453292519943295769236907684886127134428718885417;
|
|---|
| 73 |
|
|---|
| 74 | /* Radians to arcseconds */
|
|---|
| 75 | static const double PAL__DR2AS = 2.0626480624709635515647335733077861319665970087963e5;
|
|---|
| 76 |
|
|---|
| 77 | /* Arcseconds to radians */
|
|---|
| 78 | static const double PAL__DAS2R = 4.8481368110953599358991410235794797595635330237270e-6;
|
|---|
| 79 |
|
|---|
| 80 | /* Radians to degrees */
|
|---|
| 81 | static const double PAL__DR2D = 57.295779513082320876798154814105170332405472466564;
|
|---|
| 82 |
|
|---|
| 83 | /* Hours to radians */
|
|---|
| 84 | static const double PAL__DH2R = 0.26179938779914943653855361527329190701643078328126;
|
|---|
| 85 |
|
|---|
| 86 | /* Radians to hours */
|
|---|
| 87 | static const double PAL__DR2H = 3.8197186342054880584532103209403446888270314977709;
|
|---|
| 88 |
|
|---|
| 89 | /* Radians to seconds of time */
|
|---|
| 90 | static const double PAL__DR2S = 1.3750987083139757010431557155385240879777313391975e4;
|
|---|
| 91 |
|
|---|
| 92 | /* Seconds of time to radians */
|
|---|
| 93 | static const double PAL__DS2R = 7.272205216643039903848712e-5;
|
|---|
| 94 |
|
|---|
| 95 | /* Start of SLA modified Julian date epoch */
|
|---|
| 96 | static const double PAL__MJD0 = 2400000.5;
|
|---|
| 97 |
|
|---|
| 98 | /* Light time for 1 AU (sec) */
|
|---|
| 99 | static const double PAL__CR = 499.004782;
|
|---|
| 100 |
|
|---|
| 101 | /* Seconds per day */
|
|---|
| 102 | static const double PAL__SPD = 86400.0;
|
|---|
| 103 |
|
|---|
| 104 | /* Km per sec to AU per tropical century
|
|---|
| 105 | = 86400 * 36524.2198782 / 149597870 */
|
|---|
| 106 | static const double PAL__VF = 21.095;
|
|---|
| 107 |
|
|---|
| 108 | /* Radians per year to arcsec per century. This needs to be a macro since it
|
|---|
| 109 | is an expression including other constants. */
|
|---|
| 110 | #define PAL__PMF (100.0*60.0*60.0*360.0/PAL__D2PI);
|
|---|
| 111 |
|
|---|
| 112 | /* Mean sidereal rate - the rotational angular velocity of Earth
|
|---|
| 113 | in radians/sec from IERS Conventions (2003). */
|
|---|
| 114 | static const double PAL__SR = 7.2921150e-5;
|
|---|
| 115 |
|
|---|
| 116 | /* Gaussian gravitational constant (exact) */
|
|---|
| 117 | static const double PAL__GCON = 0.01720209895;
|
|---|
| 118 |
|
|---|
| 119 | /* DINT(A) - truncate to nearest whole number towards zero (double) */
|
|---|
| 120 | #define DINT(A) ((A)<0.0?ceil(A):floor(A))
|
|---|
| 121 |
|
|---|
| 122 | /* DNINT(A) - round to nearest whole number (double) */
|
|---|
| 123 | #define DNINT(A) ((A)<0.0?ceil((A)-0.5):floor((A)+0.5))
|
|---|
| 124 |
|
|---|
| 125 | /* DMAX(A,B) - return maximum value - evaluates arguments multiple times */
|
|---|
| 126 | #define DMAX(A,B) ((A) > (B) ? (A) : (B) )
|
|---|
| 127 |
|
|---|
| 128 | /* DMIN(A,B) - return minimum value - evaluates arguments multiple times */
|
|---|
| 129 | #define DMIN(A,B) ((A) < (B) ? (A) : (B) )
|
|---|
| 130 |
|
|---|
| 131 | /* We actually prefer to use C99 copysign() but we define this here as a backup
|
|---|
| 132 | but it will not detect -0.0 so is not useful for palDfltin. */
|
|---|
| 133 | /* DSIGN(A,B) - magnitude of A with sign of B (double) */
|
|---|
| 134 | #define DSIGN(A,B) ((B)<0.0?-fabs(A):fabs(A))
|
|---|
| 135 |
|
|---|
| 136 | #endif
|
|---|
