/
conformables.jl
40 lines (31 loc) · 1.54 KB
/
conformables.jl
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
################################################################################
conformables(map1::HealpixMap{T,RingOrder,AA1}, map2::HealpixMap{S,RingOrder,AA2}) where {T,S,AA1,AA2} =
map1.resolution.nside == map2.resolution.nside
conformables(
map1::HealpixMap{T,NestedOrder,AA1},
map2::HealpixMap{S,NestedOrder,AA2},
) where {T,S,AA1,AA2} = map1.resolution.nside == map2.resolution.nside
conformables(map1::HealpixMap{T,O1,AA1}, map2::HealpixMap{S,O2,AA2}) where {T,S,O1,O2,AA1,AA2} = false
################################################################################
conformables(
map1::PolarizedHealpixMap{T,RingOrder,AA1},
map2::PolarizedHealpixMap{S,RingOrder,AA2},
) where {T,S,AA1,AA2} = map1.i.resolution.nside == map2.i.resolution.nside
conformables(
map1::PolarizedHealpixMap{T,NestedOrder,AA1},
map2::PolarizedHealpixMap{S,NestedOrder,AA2},
) where {T,S,AA1,AA2} = map1.i.resolution.nside == map2.i.resolution.nside
conformables(
map1::PolarizedHealpixMap{T,O1,AA1},
map2::PolarizedHealpixMap{S,O2,AA2},
) where {T,S,O1,O2,AA1,AA2} = false
"""
conformables{T, S, O1, O2}(map1::HealpixMap{T, O1, AA1},
map2::HealpixMap{S, O2, AA2}) -> Bool
conformables{T, S, O1, O2}(map1::PolarizedHealpixMap{T, O1, AA1},
map2::PolarizedHealpixMap{S, O2, AA2}) -> Bool
Determine if two Healpix maps are "conformables", i.e., if their shape
and ordering are the same. The array types `AA1` and `AA2` are not considered
in testing conformability.
"""
conformables