22 #if defined (PI) 
23 #undef PI
201 const double CalendarAstronomer::PI = 3.14159265358979323846;  member in CalendarAstronomer
203 #define CalendarAstronomer_PI2  (CalendarAstronomer::PI*2.0)
204 #define RAD_HOUR  ( 12 / CalendarAstronomer::PI )     // radians -> hours
205 #define DEG_RAD ( CalendarAstronomer::PI / 180 )      // degrees -> radians
206 #define RAD_DEG  ( 180 / CalendarAstronomer::PI )     // radians -> degrees
222     return normalize(angle, CalendarAstronomer::PI * 2.0);  in norm2PI()
229 …return normalize(angle + CalendarAstronomer::PI, CalendarAstronomer::PI * 2.0) - CalendarAstronome…  in normPI()
503 …  double H = getLocalSidereal()*CalendarAstronomer::PI/12 - equatorial.ascension;     // Hour-angle  in eclipticToHorizon()
530 #define SUN_ETA_G    (279.403303 * CalendarAstronomer::PI/180) // Ecliptic longitude at epoch
531 #define SUN_OMEGA_G  (282.768422 * CalendarAstronomer::PI/180) // Ecliptic longitude of perigee
688     return  (CalendarAstronomer::PI/2);  in SUMMER_SOLSTICE()
710     return  ((CalendarAstronomer::PI*3)/2);  in WINTER_SOLSTICE()
1047 #define moonL0  (318.351648 * CalendarAstronomer::PI/180 )   // Mean long. at epoch
1048 #define moonP0 ( 36.340410 * CalendarAstronomer::PI/180 )   // Mean long. of perigee
1049 #define moonN0 ( 318.510107 * CalendarAstronomer::PI/180 )   // Mean long. of node
1050 #define moonI  (   5.145366 * CalendarAstronomer::PI/180 )   // Inclination of orbit
1055 #define moonT0 (   0.5181 * CalendarAstronomer::PI/180 )     // Angular size at distance A
1056 #define moonPi (   0.9507 * CalendarAstronomer::PI/180 )     // Parallax at distance A
1083         double meanLongitude = norm2PI(13.1763966*PI/180*day + moonL0);  in getMoonPosition()
1084         meanAnomalyMoon = norm2PI(meanLongitude - 0.1114041*PI/180 * day - moonP0);  in getMoonPosition()
1092         double evection = 1.2739*PI/180 * ::sin(2 * (meanLongitude - sunLongitude)  in getMoonPosition()
1094         double annual   = 0.1858*PI/180 * ::sin(meanAnomalySun);  in getMoonPosition()
1095         double a3       = 0.3700*PI/180 * ::sin(meanAnomalySun);  in getMoonPosition()
1106         double center = 6.2886*PI/180 * ::sin(meanAnomalyMoon);  in getMoonPosition()
1107         double a4 =     0.2140*PI/180 * ::sin(2 * meanAnomalyMoon);  in getMoonPosition()
1117 …    double variation = 0.6583*CalendarAstronomer::PI/180 * ::sin(2*(moonLongitude - sunLongitude));  in getMoonPosition()
1128         double nodeLongitude = norm2PI(moonN0 - 0.0529539*PI/180 * day);  in getMoonPosition()
1130         nodeLongitude -= 0.16*PI/180 * ::sin(meanAnomalySun);  in getMoonPosition()
1214     return   CalendarAstronomer::MoonAge(CalendarAstronomer::PI);  in FULL_MOON()