test_lm_space.py 2.93 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
14
15
16
17
#
# Copyright(C) 2013-2017 Max-Planck-Society
#
# 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
Jait Dixit's avatar
Jait Dixit committed
20
21
22
import unittest
import numpy as np

23
from numpy.testing import assert_, assert_equal, assert_raises,\
24
        assert_almost_equal
Martin Reinecke's avatar
Martin Reinecke committed
25
from nifty2go import LMSpace
Jait Dixit's avatar
Jait Dixit committed
26
27
from test.common import expand

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
35
36
            'lmax': 5,
            'mmax': 5,
            'shape': (36,),
            'harmonic': True,
            'dim': 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,
42
            'dim': 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
73
    @expand([['lmax', int],
            ['mmax', int],
            ['dim', int]])
74
    def test_property_ret_type(self, attribute, expected_type):
75
        l = 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']):
84
                LMSpace(lmax, mmax)
Jait Dixit's avatar
Jait Dixit committed
85
        else:
86
            l = 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
91
    def test_dvol(self):
        assert_almost_equal(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 = LMSpace(lmax)
96
        assert_almost_equal(l.get_k_length_array().val, expected)