test_sugar.py 2.78 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
# 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/>.
#
# Copyright(C) 2013-2019 Max-Planck-Society
#
# NIFTy is being developed at the Max-Planck-Institut fuer Astrophysik.

import numpy as np
19
import pytest
Philipp Arras's avatar
Philipp Arras committed
20
from numpy.testing import assert_equal
21

Martin Reinecke's avatar
Martin Reinecke committed
22
import nifty7 as ift
23

24
from .common import setup_function, teardown_function
25

26
27
pmp = pytest.mark.parametrize

28
29
30
31

def test_get_signal_variance():
    space = ift.RGSpace(3)
    hspace = space.get_default_codomain()
Philipp Arras's avatar
Philipp Arras committed
32
33
    sv = ift.get_signal_variance(lambda x: np.ones_like(x), hspace)
    assert_equal(sv, 3.)
34
35
36

    space = ift.RGSpace(3, distances=1.)
    hspace = space.get_default_codomain()
Martin Reinecke's avatar
Martin Reinecke committed
37

38
39
    def spec2(k):
        t = np.zeros_like(k)
Martin Reinecke's avatar
Martin Reinecke committed
40
        t[k == 0] = 1.
41
        return t
42
    assert_equal(ift.get_signal_variance(spec2, hspace), 1/9.)
Philipp Arras's avatar
Philipp Arras committed
43
44
45
46
47


def test_exec_time():
    dom = ift.RGSpace(12, harmonic=True)
    op = ift.HarmonicTransformOperator(dom)
Martin Reinecke's avatar
Martin Reinecke committed
48
    op1 = op.ptw("exp")
Philipp Arras's avatar
Philipp Arras committed
49
    lh = ift.GaussianEnergy(domain=op.target, sampling_dtype=np.float64) @ op1
Philipp Arras's avatar
Philipp Arras committed
50
51
    ic = ift.GradientNormController(iteration_limit=2)
    ham = ift.StandardHamiltonian(lh, ic_samp=ic)
52
    kl = ift.MetricGaussianKL.make(ift.full(ham.domain, 0.), ham, 1, False)
Philipp Arras's avatar
Philipp Arras committed
53
54
55
56
57
58
    ops = [op, op1, lh, ham, kl]
    for oo in ops:
        for wm in [True, False]:
            ift.exec_time(oo, wm)


59
60
61
@pmp('mf', [False, True])
@pmp('cplx', [False, True])
def test_calc_pos(mf, cplx):
Philipp Arras's avatar
Philipp Arras committed
62
    dom = ift.RGSpace(12, harmonic=True)
Martin Reinecke's avatar
Martin Reinecke committed
63
    op = ift.HarmonicTransformOperator(dom).ptw("exp")
64
65
66
67
68
    if mf:
        op = op.ducktape_left('foo')
        dom = ift.makeDomain({'': dom})
    if cplx:
        op = op + 1j*op
69
    fld = op(0.1 * ift.from_random(op.domain, 'normal'))
Philipp Arras's avatar
Philipp Arras committed
70
71
    pos = ift.calculate_position(op, fld)
    ift.extra.assert_allclose(op(pos), fld, 1e-1, 1e-1)
72
73
74
75
76
77
78
79


def test_isinstance_helpers():
    dom = ift.RGSpace(12, harmonic=True)
    op = ift.ScalingOperator(dom, 12.)
    fld = ift.full(dom, 0.)
    lin = ift.Linearization.make_var(fld)
    assert not ift.is_fieldlike(op)
Philipp Arras's avatar
Philipp Arras committed
80
    assert ift.is_fieldlike(lin)
81
82
83
84
85
86
87
    assert ift.is_fieldlike(fld)
    assert not ift.is_linearization(op)
    assert ift.is_linearization(lin)
    assert not ift.is_linearization(fld)
    assert ift.is_operator(op)
    assert not ift.is_operator(lin)
    assert not ift.is_operator(fld)