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
|
/*
* (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 <numeric>
#include <vector>
#include "eckit/geo/Exceptions.h"
#include "eckit/geo/order/HEALPix.h"
#include "eckit/geo/order/Scan.h"
#include "eckit/spec/Custom.h"
#include "eckit/testing/Test.h"
namespace eckit::geo::test {
CASE("order=scan") {
for (std::string scan : {"i+j+", "i+j-"}) {
order::Scan order(scan);
EXPECT(scan == order.order());
}
}
CASE("healpix") {
using order::HEALPix;
auto is_power_of_two = [](size_t n) { return (n & (n - 1)) == 0; };
SECTION("order=ring/nested") {
size_t Nside = 4;
EXPECT(HEALPix::order_default() == HEALPix(spec::Custom{}).order());
EXPECT(HEALPix::RING == HEALPix(spec::Custom{{"order", HEALPix::RING}}).order());
EXPECT(HEALPix::NESTED == HEALPix(spec::Custom{{"order", HEALPix::NESTED}}).order());
for (const auto& type : {HEALPix::order_default(), HEALPix::RING, HEALPix::NESTED}) {
HEALPix order(type);
auto no_reorder = order.reorder(type, Nside);
std::vector<size_t> expected(no_reorder.size());
std::iota(expected.begin(), expected.end(), 0);
EXPECT(expected == no_reorder);
}
}
SECTION("exceptions") {
EXPECT_THROWS_AS(HEALPix(spec::Custom{{"order", "?"}}), exception::OrderError);
EXPECT_THROWS_AS(HEALPix("?"), exception::OrderError);
for (size_t Nside : {1, 2, 3}) {
const auto size = 12 * Nside * Nside;
EXPECT_NO_THROW(HEALPix{});
EXPECT_NO_THROW(HEALPix{HEALPix::order_default()});
}
}
}
} // namespace eckit::geo::test
int main(int argc, char** argv) {
return eckit::testing::run_tests(argc, argv);
}
|
