1 // Example of using the GeographicLib::SphericalEngine class
2
3 #include <iostream>
4 #include <exception>
5 #include <vector>
6 #include <GeographicLib/SphericalEngine.hpp>
7
8 using namespace std;
9 using namespace GeographicLib;
10
main()11 int main() {
12 // See also example-SphericHarmonic.cpp
13 try {
14 int N = 3; // The maxium degree
15 double ca[] = {10, 9, 8, 7, 6, 5, 4, 3, 2, 1}; // cosine coefficients
16 vector<double> C(ca, ca + (N + 1) * (N + 2) / 2);
17 double sa[] = {6, 5, 4, 3, 2, 1}; // sine coefficients
18 vector<double> S(sa, sa + N * (N + 1) / 2);
19 SphericalEngine::coeff c[1];
20 c[0] = SphericalEngine::coeff(C, S, N);
21 double f[] = {1};
22 double x = 2, y = 3, z = 1, a = 1;
23 double v, vx, vy, vz;
24 v = SphericalEngine::Value<true, SphericalEngine::FULL, 1>
25 (c, f, x, y, z, a, vx, vy, vz);
26 cout << v << " " << vx << " " << vy << " " << vz << "\n";
27 }
28 catch (const exception& e) {
29 cerr << "Caught exception: " << e.what() << "\n";
30 return 1;
31 }
32 }
33