30 real& x, real& y, real& azi, real& rk)
const {
32 _earth.GenInverse(lat0, lon0, lat, lon,
35 t, azi0, azi, m, M, t, t);
47 real& lat, real& lon, real& azi, real& rk)
const {
51 s = _a * atan(rho/_a);
52 bool little = rho <= _a;
60 int count = numit_, trip = 0;
61 real lat1, lon1, azi1, M;
65 line.
Position(s, lat1, lon1, azi1, m, M, t);
70 real ds = little ? (m - rho * M) * M : (rho * m - M) * m;
73 if (!(fabs(ds) >= eps_ * _a))
77 lat = lat1; lon = lon1; azi = azi1; rk = M;
Math::real Position(real s12, real &lat2, real &lon2, real &azi2, real &m12, real &M12, real &M21, real &S12) const
void Forward(real lat0, real lon0, real lat, real lon, real &x, real &y, real &azi, real &rk) const
void Reverse(real lat0, real lon0, real x, real y, real &lat, real &lon, real &azi, real &rk) const