-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathring_test.go
161 lines (145 loc) · 4.72 KB
/
ring_test.go
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
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
package healpix
import (
"math"
"testing"
)
func TestFirstIndex(t *testing.T) {
testCases := []struct {
name string
order HealpixOrder
ringId int
fstIndex int
}{
{"Ring 0 @ order 0 = 0 index", 0, 0, 0},
{"Ring 1 @ order 0 = 4 index", 0, 1, 4},
{"Ring 2 @ order 0 = 8 index", 0, 2, 8},
{"Ring 0 @ order 1 = 0 index", 1, 0, 0},
{"Ring 1 @ order 1 = 4 index", 1, 1, 4},
{"Ring 2 @ order 1 = 12 index", 1, 2, 12},
{"Ring 3 @ order 1 = 20 index", 1, 3, 20},
{"Ring 4 @ order 1 = 28 index", 1, 4, 28},
{"Ring 5 @ order 1 = 36 index", 1, 5, 36},
{"Ring 6 @ order 1 = 44 index", 1, 6, 44},
{"Ring 0 @ order 2 = 0 index", 2, 0, 0},
{"Ring 1 @ order 2 = 4 index", 2, 1, 4},
{"Ring 2 @ order 2 = 8 index", 2, 2, 12},
{"Ring 3 @ order 2 = 24 index", 2, 3, 24},
{"Ring 4 @ order 2 = 40 index", 2, 4, 40},
{"Ring 10 @ order 2 = 136 index", 2, 10, 136},
{"Ring 11 @ order 2 = 152 index", 2, 11, 152},
{"Ring 12 @ order 2 = 168 index", 2, 12, 168},
{"Ring 13 @ order 2 = 180 index", 2, 13, 180},
{"Ring 14 @ order 2 = 188 index", 2, 14, 188},
}
for _, tc := range testCases {
t.Run(tc.name, func(t *testing.T) {
ring := NewRing(tc.order, tc.ringId)
if ring.FirstIndex() != tc.fstIndex {
t.Errorf("Ring %v @ order %v first index expected %v, got %v instead", tc.ringId, tc.order, tc.fstIndex, ring.FirstIndex())
}
})
}
}
func TestRingPixels(t *testing.T) {
testCases := []struct {
name string
order HealpixOrder
ringId int
offset bool
}{
{"Ring 0 @ order 0 = true", 0, 0, true},
{"Ring 1 @ order 0 = false", 0, 1, false},
{"Ring 2 @ order 0 = true", 0, 2, true},
{"Ring 0 @ order 1 = true", 1, 0, true},
{"Ring 1 @ order 1 = true", 1, 1, true},
{"Ring 2 @ order 1 = false", 1, 2, false},
{"Ring 3 @ order 1 = true", 1, 3, true},
{"Ring 4 @ order 1 = false", 1, 4, false},
{"Ring 5 @ order 1 = true", 1, 5, true},
{"Ring 6 @ order 1 = true", 1, 6, true},
{"Ring 0 @ order 2 = true", 2, 0, true},
{"Ring 1 @ order 2 = true", 2, 1, true},
{"Ring 2 @ order 2 = true", 2, 2, true},
{"Ring 3 @ order 2 = true", 2, 3, true},
{"Ring 4 @ order 2 = false", 2, 4, false},
{"Ring 10 @ order 2 = false", 2, 10, false},
{"Ring 11 @ order 2 = true", 2, 11, true},
{"Ring 12 @ order 2 = true", 2, 12, true},
{"Ring 13 @ order 2 = true", 2, 13, true},
{"Ring 14 @ order 2 = true", 2, 14, true},
}
for _, tc := range testCases {
t.Run(tc.name, func(t *testing.T) {
ring := NewRing(tc.order, tc.ringId)
if ring.IsOffset() != tc.offset {
t.Errorf("Ring %v @ order %v offset expected %v, got %v instead", tc.ringId, tc.order, tc.offset, ring.IsOffset())
}
})
}
}
func TestRingOffset(t *testing.T) {
testCases := []struct {
name string
order HealpixOrder
ringId int
pixels int
}{
{"Ring 0 @ order 0 = 4 pixels", 0, 0, 4},
{"Ring 1 @ order 0 = 4 pixels", 0, 1, 4},
{"Ring 2 @ order 0 = 4 pixels", 0, 2, 4},
{"Ring 1 @ order 1 = 8 pixels", 1, 1, 8},
{"Ring 2 @ order 1 = 8 pixels", 1, 2, 8},
{"Ring 3 @ order 1 = 8 pixels", 1, 3, 8},
{"Ring 4 @ order 1 = 8 pixels", 1, 4, 8},
{"Ring 5 @ order 1 = 8 pixels", 1, 5, 8},
{"Ring 6 @ order 1 = 4 pixels", 1, 6, 4},
{"Ring 0 @ order 2 = 0 pixels", 2, 0, 4},
{"Ring 1 @ order 2 = 4 pixels", 2, 1, 8},
{"Ring 2 @ order 2 = 8 pixels", 2, 2, 12},
{"Ring 3 @ order 2 = 24 pixels", 2, 3, 16},
{"Ring 4 @ order 2 = 40 pixels", 2, 4, 16},
{"Ring 10 @ order 2 = 136 pixels", 2, 10, 16},
{"Ring 11 @ order 2 = 152 pixels", 2, 11, 16},
{"Ring 12 @ order 2 = 168 pixels", 2, 12, 12},
{"Ring 13 @ order 2 = 180 pixels", 2, 13, 8},
{"Ring 14 @ order 2 = 188 pixels", 2, 14, 4},
}
for _, tc := range testCases {
t.Run(tc.name, func(t *testing.T) {
ring := NewRing(tc.order, tc.ringId)
if ring.Pixels() != tc.pixels {
t.Errorf("Ring %v @ order %v pixels expected %v, got %v instead", tc.ringId, tc.order, tc.pixels, ring.Pixels())
}
})
}
}
func withinTolerance(n1, n2, tolerance float64) bool {
if n1 == n2 {
return true
}
diff := math.Abs(n1 - n2)
if n2 == 0 {
return diff < tolerance
}
return (diff / math.Abs(n2)) < tolerance
}
func TestRingLatitude(t *testing.T) {
testCases := []struct {
name string
order HealpixOrder
ringId int
latitude float64
}{
{"Ring 0 @ order 0 = Pi/4", 0, 0, math.Pi/2 - math.Acos(float64(2)/float64(3))},
{"Ring 1 @ order 0 = 0", 0, 1, 0},
{"Ring 2 @ order 0 = -Pi/4", 0, 2, math.Pi/2 - math.Acos(float64(-2)/float64(3))},
}
for _, tc := range testCases {
t.Run(tc.name, func(t *testing.T) {
ring := NewRing(tc.order, tc.ringId)
if !withinTolerance(ring.Latitude(), tc.latitude, 1e-10) {
t.Errorf("Ring %v @ order %v latitude expected %v, got %v instead", tc.ringId, tc.order, tc.latitude, ring.Latitude())
}
})
}
}