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
|
/*
* (C) Copyright 1996- ECMWF.
*
* This software is licensed under the terms of the Apache Licence Version 2.0
* which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
*
* In applying this licence, ECMWF does not waive the privileges and immunities
* granted to it by virtue of its status as an intergovernmental organisation nor
* does it submit to any jurisdiction.
*/
#include <memory>
#include "eckit/geo/Projection.h"
#include "eckit/geo/figure/UnitSphere.h"
#include "eckit/geo/projection/LonLatToXYZ.h" // to test Reverse
#include "eckit/geo/projection/Reverse.h"
#include "eckit/spec/Custom.h"
#include "eckit/testing/Test.h"
namespace eckit::geo::test {
CASE("projection: none") {
Point p = PointLonLat{1, 1};
std::unique_ptr<Projection> projection(ProjectionFactoryType::instance().get("none").create(spec::Custom{}));
EXPECT(points_equal(p, projection->inv(p)));
EXPECT(points_equal(p, projection->fwd(p)));
}
CASE("projection: reverse") {
projection::LonLatToXYZ ab(new figure::UnitSphere);
projection::Reverse<projection::LonLatToXYZ> ba(new figure::UnitSphere);
PointLonLat p = NORTH_POLE;
PointXYZ q{0., 0., 1.};
ASSERT(points_equal(q, ab.fwd(p)));
ASSERT(points_equal(p, ab.inv(q)));
// ensure fwd(PointXYZ) -> PointLonLat, inv(PointLonLat) -> PointXYZ
EXPECT(points_equal(p, ba.fwd(q)));
EXPECT(points_equal(q, ba.inv(p)));
ASSERT(ab.spec().get_string("type") == "ll-to-xyz");
EXPECT(ba.spec().get_string("type") == "reverse-ll-to-xyz");
}
} // namespace eckit::geo::test
int main(int argc, char** argv) {
return eckit::testing::run_tests(argc, argv);
}
|
