# define M_E 2.7182818284590452354
#endif
-// ONE_SECOND is pi/180/60/60, or about 100 feet at earths' equator
-#define ONE_SECOND 4.848136811E-6
+// SG_ONE_SECOND is pi/180/60/60, or about 100 feet at earths' equator
+#define SG_ONE_SECOND 4.848136811E-6
// Radius of Earth in kilometers at the equator. Another source had
// 6378.165 but this is probably close enough
-#define EARTH_RAD 6378.155
+#define SG_EARTH_RAD 6378.155
// Earth parameters for WGS 84, taken from LaRCsim/ls_constants.h
}
template <class T>
-inline const T SG_MIN(const T a, const T b) {
+inline const T SG_MIN2(const T a, const T b) {
return a < b ? a : b;
}
// return the minimum of three values
template <class T>
inline const T SG_MIN3( const T a, const T b, const T c) {
- return (a < b ? SG_MIN (a, c) : SG_MIN (b, c));
+ return (a < b ? SG_MIN2 (a, c) : SG_MIN2 (b, c));
}
template <class T>
-inline const T SG_MAX(const T a, const T b) {
+inline const T SG_MAX2(const T a, const T b) {
return a > b ? a : b;
}
// return the maximum of three values
template <class T>
inline const T SG_MAX3 (const T a, const T b, const T c) {
- return (a > b ? SG_MAX (a, c) : SG_MAX (b, c));
+ return (a > b ? SG_MAX2 (a, c) : SG_MAX2 (b, c));
}
//
lat_geoc = SGD_PI_2 - atan2( sqrt(cp.x()*cp.x() + cp.y()*cp.y()), cp.z() );
radius = sqrt( cp.x()*cp.x() + cp.y()*cp.y() + cp.z()*cp.z() );
- if( ( (SGD_PI_2 - lat_geoc) < ONE_SECOND ) // near North pole
- || ( (SGD_PI_2 + lat_geoc) < ONE_SECOND ) ) // near South pole
+ if( ( (SGD_PI_2 - lat_geoc) < SG_ONE_SECOND ) // near North pole
+ || ( (SGD_PI_2 + lat_geoc) < SG_ONE_SECOND ) ) // near South pole
{
result = radius - EQUATORIAL_RADIUS_M*E;
} else {
SG_USING_STD(cout);
#endif
-// ONE_SECOND is pi/180/60/60, or about 100 feet at earths' equator
-#define ONE_SECOND 4.848136811E-6
-
#define DOMAIN_ERR_DEBUG 1
double t_lat, x_alpha, mu_alpha, delt_mu, r_alpha, l_point, rho_alpha;
double sin_mu_a, denom,delt_lambda, lambda_sl, sin_lambda_sl;
- if( ( (SGD_PI_2 - lat_geoc) < ONE_SECOND ) // near North pole
- || ( (SGD_PI_2 + lat_geoc) < ONE_SECOND ) ) // near South pole
+ if( ( (SGD_PI_2 - lat_geoc) < SG_ONE_SECOND ) // near North pole
+ || ( (SGD_PI_2 + lat_geoc) < SG_ONE_SECOND ) ) // near South pole
{
*lat_geod = lat_geoc;
*sea_level_r = EQUATORIAL_RADIUS_M*E;
CMetarStation::for_each( dispMetarStation );
glEnd();
glColor3d(0.0, 0.0, 1.0);
- glutWireSphere(EARTH_RAD,32,16);
+ glutWireSphere(SG_EARTH_RAD,32,16);
glutSwapBuffers();
ii++;
glPopMatrix();