coursepointer 0.3.4

Converts waypoints into Garmin FIT course points
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126

/**
 * \file ProjTest.cpp
 * \brief Command line utility for testing transverse Mercator projections
 *
 * Copyright (c) Charles Karney (2008-2022) <karney@alum.mit.edu> and licensed
 * under the MIT/X11 License.  For more information, see
 * https://geographiclib.sourceforge.io/
 **********************************************************************/

#include "GeographicLib/EllipticFunction.hpp"
#include "GeographicLib/TransverseMercator.hpp"
#include "GeographicLib/TransverseMercatorExact.hpp"
#include "GeographicLib/PolarStereographic.hpp"
#include <iostream>
#include <cassert>

using namespace GeographicLib;
using namespace std;

int main() {
  {
    Math::real x, y, gamma, k;
    x = y = gamma = k = 0;
    PolarStereographic::UPS().Forward(true, Math::NaN<Math::real>(), 30.0,
                                    x, y, gamma, k);
    cout << x << " " << y << " " << k << "\n";
    x = y = gamma = k = 0;
    PolarStereographic::UPS().Forward(false, -80.0, Math::NaN<Math::real>(),
                                    x, y, gamma, k);
    cout << x << " " << y << " " << gamma << "\n";
  }
  {
    Math::real lat, lon, gamma, k;
    lat = lon = gamma = k = 0;
    PolarStereographic::UPS().Reverse(true, Math::NaN<Math::real>(), 0.0,
                                    lat, lon, gamma, k);
    cout << lat << " " << lon << " " << gamma << " " << k << "\n";
    lat = lon = gamma = k = 0;
    PolarStereographic::UPS().Reverse(false,0.0, Math::NaN<Math::real>(),
                                    lat, lon, gamma, k);
    cout << lat << " " << lon << " " << gamma << " " << k << "\n";
  }
  {
    Math::real x, y, gamma, k;
    x = y = gamma = k = 0;
    TransverseMercator::UTM().Forward(0.0, Math::NaN<Math::real>(), 0.0,
                                    x, y, gamma, k);
    cout << x << " " << y << " " << gamma << " " << k << "\n";
    x = y = gamma = k = 0;
    TransverseMercator::UTM().Forward(0.0, 0.0, Math::NaN<Math::real>(),
                                    x, y, gamma, k);
    cout << x << " " << y << " " << gamma << " " << k << "\n";
    x = y = gamma = k = 0;
    TransverseMercator::UTM().Forward(Math::NaN<Math::real>(), 0.0, 0.0,
                                    x, y, gamma, k);
    cout << x << " " << y << " " << gamma << " " << k << "\n";
  }
  {
    Math::real lat, lon, gamma, k;
    lat = lon = gamma = k = 0;
    TransverseMercator::UTM().Reverse(0.0, Math::NaN<Math::real>(), 0.0,
                                    lat, lon, gamma, k);
    cout << lat << " " << lon << " " << gamma << " " << k << "\n";
    lat = lon = gamma = k = 0;
    TransverseMercator::UTM().Reverse(0.0, 0.0, Math::NaN<Math::real>(),
                                    lat, lon, gamma, k);
    cout << lat << " " << lon << " " << gamma << " " << k << "\n";
    lat = lon = gamma = k = 0;
    TransverseMercator::UTM().Reverse(Math::NaN<Math::real>(), 0.0, 0.0,
                                    lat, lon, gamma, k);
    cout << lon << "\n";
  }
  {
    Math::real x, y, gamma, k;
    x = y = gamma = k = 0;
    TransverseMercatorExact::UTM().Forward(0.0, Math::NaN<Math::real>(), 0.0,
                                         x, y, gamma, k);
    cout << x << " " << y << " " << gamma << " " << k << "\n";
    x = y = gamma = k = 0;
    TransverseMercatorExact::UTM().Forward(0.0, 0.0, Math::NaN<Math::real>(),
                                         x, y, gamma, k);
    cout << x << " " << y << " " << gamma << " " << k << "\n";
    x = y = gamma = k = 0;
    TransverseMercatorExact::UTM().Forward(Math::NaN<Math::real>(), 0.0, 0.0,
                                         x, y, gamma, k);
    cout << x << " " << y << " " << gamma << " " << k << "\n";
  }
  {
    Math::real lat, lon, gamma, k;
    lat = lon = gamma = k = 0;
    TransverseMercatorExact::UTM().Reverse(0.0, Math::NaN<Math::real>(), 0.0,
                                         lat, lon, gamma, k);
    cout << lat << " " << lon << " " << gamma << " " << k << "\n";
    lat = lon = gamma = k = 0;
    TransverseMercatorExact::UTM().Reverse(0.0, 0.0, Math::NaN<Math::real>(),
                                         lat, lon, gamma, k);
    cout << lat << " " << lon << " " << gamma << " " << k << "\n";
    lat = lon = gamma = k = 0;
    TransverseMercatorExact::UTM().Reverse(Math::NaN<Math::real>(), 0.0, 0.0,
                                         lat, lon, gamma, k);
    cout << lon << "\n";
  }
  {
    EllipticFunction e(Math::NaN<Math::real>());
    Math::real sn, cn, dn;
    cout << " k2 " << e.k2()
         << " kp2 " << e.kp2()
         << " K " << e.K()
         << " E " << e.E()
         << " KE " << e.KE() << "\n";
    sn = cn = dn = 0;
    e.sncndn(Math::real(0.1), sn, cn, dn);
    cout << " sncndn " << sn << " " << cn << " " << dn
         << " E(phi) " << e.E(Math::real(0.1))
         << " E(sncndn) " << e.E(sn, cn, dn) << "\n";
  }
  {
    EllipticFunction e(Math::real(0.1));
    Math::real sn, cn, dn;
    sn = cn = dn = 0;
    e.sncndn(Math::NaN<Math::real>(), sn, cn, dn);
    cout << " sncndn " << sn << " " << cn << " " << dn
         << " E(phi) " << e.E(Math::NaN<Math::real>())
         << " E(sncndn) " << e.E(sn, cn, dn) << "\n";
  }
}