10 #if !defined(GEOGRAPHICLIB_GEODESICLINE_HPP)
11 #define GEOGRAPHICLIB_GEODESICLINE_HPP 1
75 static const int nC1_ = Geodesic::nC1_;
76 static const int nC1p_ = Geodesic::nC1p_;
77 static const int nC2_ = Geodesic::nC2_;
78 static const int nC3_ = Geodesic::nC3_;
79 static const int nC4_ = Geodesic::nC4_;
82 real _lat1, _lon1, _azi1;
83 real _a, _f, _b, _c2, _f1, _salp0, _calp0, _k2,
84 _salp1, _calp1, _ssig1, _csig1, _dn1, _stau1, _ctau1, _somg1, _comg1,
85 _A1m1, _A2m1, _A3c, _B11, _B21, _B31, _A4, _B41;
88 real _C1a[nC1_ + 1], _C1pa[nC1p_ + 1], _C2a[nC2_ + 1], _C3a[nC3_],
94 real azi1, real salp1, real calp1,
98 real azi1, real salp1, real calp1,
99 unsigned caps,
bool arcmode, real s13_a13);
102 CAP_NONE = Geodesic::CAP_NONE,
103 CAP_C1 = Geodesic::CAP_C1,
104 CAP_C1p = Geodesic::CAP_C1p,
105 CAP_C2 = Geodesic::CAP_C2,
106 CAP_C3 = Geodesic::CAP_C3,
107 CAP_C4 = Geodesic::CAP_C4,
108 CAP_ALL = Geodesic::CAP_ALL,
109 CAP_MASK = Geodesic::CAP_MASK,
110 OUT_ALL = Geodesic::OUT_ALL,
111 OUT_MASK = Geodesic::OUT_MASK,
227 unsigned caps = ALL);
282 real& lat2, real& lon2, real& azi2,
283 real& m12, real& M12, real& M21,
286 return GenPosition(
false, s12,
287 LATITUDE | LONGITUDE | AZIMUTH |
288 REDUCEDLENGTH | GEODESICSCALE | AREA,
289 lat2, lon2, azi2, t, m12, M12, M21, S12);
297 return GenPosition(
false, s12,
298 LATITUDE | LONGITUDE,
299 lat2, lon2, t, t, t, t, t, t);
308 return GenPosition(
false, s12,
309 LATITUDE | LONGITUDE | AZIMUTH,
310 lat2, lon2, azi2, t, t, t, t, t);
317 real& azi2, real& m12)
const {
319 return GenPosition(
false, s12,
320 LATITUDE | LONGITUDE |
321 AZIMUTH | REDUCEDLENGTH,
322 lat2, lon2, azi2, t, m12, t, t, t);
329 real& azi2, real& M12, real& M21)
332 return GenPosition(
false, s12,
333 LATITUDE | LONGITUDE |
334 AZIMUTH | GEODESICSCALE,
335 lat2, lon2, azi2, t, t, M12, M21, t);
342 real& lat2, real& lon2, real& azi2,
343 real& m12, real& M12, real& M21)
346 return GenPosition(
false, s12,
347 LATITUDE | LONGITUDE | AZIMUTH |
348 REDUCEDLENGTH | GEODESICSCALE,
349 lat2, lon2, azi2, t, m12, M12, M21, t);
395 real& s12, real& m12, real& M12, real& M21,
397 GenPosition(
true, a12,
398 LATITUDE | LONGITUDE | AZIMUTH | DISTANCE |
399 REDUCEDLENGTH | GEODESICSCALE | AREA,
400 lat2, lon2, azi2, s12, m12, M12, M21, S12);
409 GenPosition(
true, a12,
410 LATITUDE | LONGITUDE,
411 lat2, lon2, t, t, t, t, t, t);
418 real& lat2, real& lon2, real& azi2)
421 GenPosition(
true, a12,
422 LATITUDE | LONGITUDE | AZIMUTH,
423 lat2, lon2, azi2, t, t, t, t, t);
432 GenPosition(
true, a12,
433 LATITUDE | LONGITUDE | AZIMUTH | DISTANCE,
434 lat2, lon2, azi2, s12, t, t, t, t);
441 real& s12, real& m12)
const {
443 GenPosition(
true, a12,
444 LATITUDE | LONGITUDE | AZIMUTH |
445 DISTANCE | REDUCEDLENGTH,
446 lat2, lon2, azi2, s12, m12, t, t, t);
453 real& s12, real& M12, real& M21)
456 GenPosition(
true, a12,
457 LATITUDE | LONGITUDE | AZIMUTH |
458 DISTANCE | GEODESICSCALE,
459 lat2, lon2, azi2, s12, t, M12, M21, t);
466 real& s12, real& m12, real& M12, real& M21)
469 GenPosition(
true, a12,
470 LATITUDE | LONGITUDE | AZIMUTH |
471 DISTANCE | REDUCEDLENGTH | GEODESICSCALE,
472 lat2, lon2, azi2, s12, m12, M12, M21, t);
537 Math::real GenPosition(
bool arcmode,
real s12_a12,
unsigned outmask,
556 void SetDistance(
real s13);
567 void SetArc(
real a13);
579 void GenSetDistance(
bool arcmode,
real s13_a13);
589 bool Init()
const {
return _caps != 0U; }
616 {
if (Init()) { sazi1 = _salp1; cazi1 = _calp1; } }
634 {
if (Init()) { sazi0 = _salp0; cazi0 = _calp0; } }
674 return (_caps & testcaps) == testcaps;
685 {
return Init() ? (arcmode ? _a13 : _s13) :
Math::NaN(); }
702 #endif // GEOGRAPHICLIB_GEODESICLINE_HPP