test_lm_space.py 2.94 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program.  If not, see <http://www.gnu.org/licenses/>.
Theo Steininger's avatar
Theo Steininger committed
13
#
Martin Reinecke's avatar
Martin Reinecke committed
14
# Copyright(C) 2013-2018 Max-Planck-Society
Theo Steininger's avatar
Theo Steininger committed
15
16
17
#
# NIFTy is being developed at the Max-Planck-Institut fuer Astrophysik
# and financially supported by the Studienstiftung des deutschen Volkes.
18

Martin Reinecke's avatar
Martin Reinecke committed
19
from __future__ import division
20

Jait Dixit's avatar
Jait Dixit committed
21
22
23
import unittest
from test.common import expand

24
25
26
27
import nifty5 as ift
import numpy as np
from numpy.testing import assert_, assert_allclose, assert_equal, assert_raises

Martin Reinecke's avatar
Martin Reinecke committed
28
# [lmax, expected]
29
CONSTRUCTOR_CONFIGS = [
30
        [5, None, {
Jait Dixit's avatar
Jait Dixit committed
31
32
33
34
            'lmax': 5,
            'mmax': 5,
            'shape': (36,),
            'harmonic': True,
Martin Reinecke's avatar
Martin Reinecke committed
35
            'size': 36,
Jait Dixit's avatar
Jait Dixit committed
36
            }],
37
        [7, 4, {
Jait Dixit's avatar
Jait Dixit committed
38
            'lmax': 7,
39
40
            'mmax': 4,
            'shape': (52,),
Jait Dixit's avatar
Jait Dixit committed
41
            'harmonic': True,
Martin Reinecke's avatar
Martin Reinecke committed
42
            'size': 52,
Jait Dixit's avatar
Jait Dixit committed
43
            }],
44
        [-1, 28, {
Jait Dixit's avatar
Jait Dixit committed
45
46
47
48
            'error': ValueError
            }]
    ]

Theo Steininger's avatar
Theo Steininger committed
49

50
def _k_length_array_helper(index_arr, lmax):
Martin Reinecke's avatar
Martin Reinecke committed
51
52
53
54
    if index_arr <= lmax:
        index_half = index_arr
    else:
        if (index_arr - lmax) % 2 == 0:
Martin Reinecke's avatar
Martin Reinecke committed
55
            index_half = (index_arr + lmax)//2
Martin Reinecke's avatar
Martin Reinecke committed
56
        else:
Martin Reinecke's avatar
Martin Reinecke committed
57
            index_half = (index_arr + lmax + 1)//2
Martin Reinecke's avatar
Martin Reinecke committed
58

59
60
    m = np.ceil(((2*lmax + 1) - np.sqrt((2*lmax + 1)**2 -
                 8*(index_half - lmax)))/2).astype(int)
Martin Reinecke's avatar
Martin Reinecke committed
61

62
    return index_half - m*(2*lmax + 1 - m)//2
Martin Reinecke's avatar
Martin Reinecke committed
63

Theo Steininger's avatar
Theo Steininger committed
64

65
66
def get_k_length_array_configs():
    da_0 = [_k_length_array_helper(idx, 5) for idx in np.arange(36)]
Martin Reinecke's avatar
Martin Reinecke committed
67
    return [[5, da_0]]
68
69


70
class LMSpaceInterfaceTests(unittest.TestCase):
Jait Dixit's avatar
Jait Dixit committed
71
72
    @expand([['lmax', int],
            ['mmax', int],
Martin Reinecke's avatar
Martin Reinecke committed
73
            ['size', int]])
74
    def test_property_ret_type(self, attribute, expected_type):
Martin Reinecke's avatar
Martin Reinecke committed
75
        l = ift.LMSpace(7, 5)
76
        assert_(isinstance(getattr(l, attribute), expected_type))
Jait Dixit's avatar
Jait Dixit committed
77
78
79


class LMSpaceFunctionalityTests(unittest.TestCase):
80
    @expand(CONSTRUCTOR_CONFIGS)
81
    def test_constructor(self, lmax, mmax, expected):
82
83
        if 'error' in expected:
            with assert_raises(expected['error']):
Martin Reinecke's avatar
Martin Reinecke committed
84
                ift.LMSpace(lmax, mmax)
Jait Dixit's avatar
Jait Dixit committed
85
        else:
Martin Reinecke's avatar
Martin Reinecke committed
86
            l = ift.LMSpace(lmax, mmax)
Martin Reinecke's avatar
Martin Reinecke committed
87
            for key, value in expected.items():
88
                assert_equal(getattr(l, key), value)
Jait Dixit's avatar
Jait Dixit committed
89

90
    def test_dvol(self):
Martin Reinecke's avatar
Martin Reinecke committed
91
        assert_allclose(ift.LMSpace(5).dvol, 1.)
Jait Dixit's avatar
Jait Dixit committed
92

93
94
    @expand(get_k_length_array_configs())
    def test_k_length_array(self, lmax, expected):
Martin Reinecke's avatar
Martin Reinecke committed
95
        l = ift.LMSpace(lmax)
Martin Reinecke's avatar
Martin Reinecke committed
96
        assert_allclose(l.get_k_length_array().to_global_data(), expected)