There is a maintenance of MPCDF Gitlab on Thursday, April 22st 2020, 9:00 am CEST - Expect some service interruptions during this time

test_adjoint.py 11.3 KB
Newer Older
Martin Reinecke's avatar
Martin Reinecke committed
1 2 3 4 5 6 7 8 9 10 11 12 13
# 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/>.
#
Philipp Arras's avatar
Philipp Arras committed
14
# Copyright(C) 2013-2020 Max-Planck-Society
Martin Reinecke's avatar
Martin Reinecke committed
15
#
16
# NIFTy is being developed at the Max-Planck-Institut fuer Astrophysik.
Martin Reinecke's avatar
Martin Reinecke committed
17

Philipp Arras's avatar
Philipp Arras committed
18 19
import numpy as np
import pytest
20

Martin Reinecke's avatar
5->6  
Martin Reinecke committed
21
import nifty6 as ift
22

23
from ..common import list2fixture, setup_function, teardown_function
Martin Reinecke's avatar
Martin Reinecke committed
24

Philipp Arras's avatar
Philipp Arras committed
25 26
_h_RG_spaces = [ift.RGSpace(7, distances=0.2, harmonic=True),
                ift.RGSpace((12, 46), distances=(.2, .3), harmonic=True)]
Philipp Arras's avatar
Philipp Arras committed
27
_h_spaces = _h_RG_spaces + [ift.LMSpace(17)]
Philipp Arras's avatar
Philipp Arras committed
28 29
_p_RG_spaces = [ift.RGSpace(19, distances=0.7),
                ift.RGSpace((1, 2, 3, 6), distances=(0.2, 0.25, 0.34, .8))]
30
_p_spaces = _p_RG_spaces + [ift.HPSpace(17), ift.GLSpace(8, 13)]
Philipp Arras's avatar
Philipp Arras committed
31 32
_pow_spaces = [ift.PowerSpace(ift.RGSpace((17, 38), (0.99, 1340), harmonic=True)),
               ift.PowerSpace(ift.LMSpace(18), ift.PowerSpace.useful_binbounds(ift.LMSpace(18), False))]
33

Philipp Arras's avatar
Philipp Arras committed
34 35 36
pmp = pytest.mark.parametrize
dtype = list2fixture([np.float64, np.complex128])

Martin Reinecke's avatar
Martin Reinecke committed
37

Philipp Arras's avatar
Philipp Arras committed
38 39
@pmp('sp', _p_RG_spaces)
def testLOSResponse(sp, dtype):
Martin Reinecke's avatar
Martin Reinecke committed
40 41 42 43
    starts = ift.random.current_rng().standard_normal((len(sp.shape), 10))
    ends = ift.random.current_rng().standard_normal((len(sp.shape), 10))
    sigma_low = 1e-4*ift.random.current_rng().standard_normal(10)
    sigma_ups = 1e-5*ift.random.current_rng().standard_normal(10)
Philipp Arras's avatar
Philipp Arras committed
44 45 46 47 48 49
    op = ift.LOSResponse(sp, starts, ends, sigma_low, sigma_ups)
    ift.extra.consistency_check(op, dtype, dtype)


@pmp('sp', _h_spaces + _p_spaces + _pow_spaces)
def testOperatorCombinations(sp, dtype):
50 51
    a = ift.DiagonalOperator(ift.Field.from_random(sp, "normal", dtype=dtype))
    b = ift.DiagonalOperator(ift.Field.from_random(sp, "normal", dtype=dtype))
Philipp Arras's avatar
Philipp Arras committed
52 53 54 55 56 57
    op = ift.SandwichOperator.make(a, b)
    ift.extra.consistency_check(op, dtype, dtype)
    op = a(b)
    ift.extra.consistency_check(op, dtype, dtype)
    op = a + b
    ift.extra.consistency_check(op, dtype, dtype)
58 59
    op = a - b
    ift.extra.consistency_check(op, dtype, dtype)
Philipp Arras's avatar
Philipp Arras committed
60 61 62 63


def testLinearInterpolator():
    sp = ift.RGSpace((10, 8), distances=(0.1, 3.5))
Martin Reinecke's avatar
Martin Reinecke committed
64
    pos = ift.random.current_rng().random((2, 23))
Philipp Arras's avatar
Philipp Arras committed
65 66 67 68 69 70
    pos[0, :] *= 0.9
    pos[1, :] *= 7*3.5
    op = ift.LinearInterpolator(sp, pos)
    ift.extra.consistency_check(op)


71 72 73 74
@pmp('sp', _h_spaces + _p_spaces + _pow_spaces)
def testRealizer(sp):
    op = ift.Realizer(sp)
    ift.extra.consistency_check(op, np.complex128, np.float64,
Martin Reinecke's avatar
Martin Reinecke committed
75
                                only_r_linear=True)
76 77


Philipp Arras's avatar
Philipp Arras committed
78 79 80 81 82 83 84
@pmp('sp', _h_spaces + _p_spaces + _pow_spaces)
def testImaginizer(sp):
    op = ift.Imaginizer(sp)
    ift.extra.consistency_check(op, np.complex128, np.float64,
                                only_r_linear=True)


85 86 87 88
@pmp('sp', _h_spaces + _p_spaces + _pow_spaces)
def testConjugationOperator(sp):
    op = ift.ConjugationOperator(sp)
    ift.extra.consistency_check(op, np.complex128, np.complex128,
Martin Reinecke's avatar
Martin Reinecke committed
89
                                only_r_linear=True)
90 91


Philipp Arras's avatar
Philipp Arras committed
92 93
@pmp('sp', _h_spaces + _p_spaces + _pow_spaces)
def testOperatorAdaptor(sp, dtype):
94
    op = ift.DiagonalOperator(ift.Field.from_random(sp, "normal", dtype=dtype))
Philipp Arras's avatar
Philipp Arras committed
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
    ift.extra.consistency_check(op.adjoint, dtype, dtype)
    ift.extra.consistency_check(op.inverse, dtype, dtype)
    ift.extra.consistency_check(op.inverse.adjoint, dtype, dtype)
    ift.extra.consistency_check(op.adjoint.inverse, dtype, dtype)


@pmp('sp1', _h_spaces + _p_spaces + _pow_spaces)
@pmp('sp2', _h_spaces + _p_spaces + _pow_spaces)
def testNullOperator(sp1, sp2, dtype):
    op = ift.NullOperator(sp1, sp2)
    ift.extra.consistency_check(op, dtype, dtype)
    mdom1 = ift.MultiDomain.make({'a': sp1})
    mdom2 = ift.MultiDomain.make({'b': sp2})
    op = ift.NullOperator(mdom1, mdom2)
    ift.extra.consistency_check(op, dtype, dtype)
    op = ift.NullOperator(sp1, mdom2)
    ift.extra.consistency_check(op, dtype, dtype)
    op = ift.NullOperator(mdom1, sp2)
    ift.extra.consistency_check(op, dtype, dtype)


@pmp('sp', _p_RG_spaces)
def testHarmonicSmoothingOperator(sp, dtype):
    op = ift.HarmonicSmoothingOperator(sp, 0.1)
    ift.extra.consistency_check(op, dtype, dtype)


@pmp('sp', _h_spaces + _p_spaces + _pow_spaces)
def testDOFDistributor(sp, dtype):
    # TODO: Test for DomainTuple
    if sp.size < 4:
        return
    dofdex = np.arange(sp.size).reshape(sp.shape) % 3
Martin Reinecke's avatar
Martin Reinecke committed
128
    dofdex = ift.Field.from_raw(sp, dofdex)
Philipp Arras's avatar
Philipp Arras committed
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 162 163 164 165 166 167 168 169 170
    op = ift.DOFDistributor(dofdex)
    ift.extra.consistency_check(op, dtype, dtype)


@pmp('sp', _h_spaces)
def testPPO(sp, dtype):
    op = ift.PowerDistributor(target=sp)
    ift.extra.consistency_check(op, dtype, dtype)
    ps = ift.PowerSpace(
        sp, ift.PowerSpace.useful_binbounds(sp, logarithmic=False, nbin=3))
    op = ift.PowerDistributor(target=sp, power_space=ps)
    ift.extra.consistency_check(op, dtype, dtype)
    ps = ift.PowerSpace(
        sp, ift.PowerSpace.useful_binbounds(sp, logarithmic=True, nbin=3))
    op = ift.PowerDistributor(target=sp, power_space=ps)
    ift.extra.consistency_check(op, dtype, dtype)


@pmp('sp', _h_RG_spaces + _p_RG_spaces)
def testFFT(sp, dtype):
    op = ift.FFTOperator(sp)
    ift.extra.consistency_check(op, dtype, dtype)
    op = ift.FFTOperator(sp.get_default_codomain())
    ift.extra.consistency_check(op, dtype, dtype)


@pmp('sp', _h_RG_spaces + _p_RG_spaces)
def testHartley(sp, dtype):
    op = ift.HartleyOperator(sp)
    ift.extra.consistency_check(op, dtype, dtype)
    op = ift.HartleyOperator(sp.get_default_codomain())
    ift.extra.consistency_check(op, dtype, dtype)


@pmp('sp', _h_spaces)
def testHarmonic(sp, dtype):
    op = ift.HarmonicTransformOperator(sp)
    ift.extra.consistency_check(op, dtype, dtype)


@pmp('sp', _p_spaces)
def testMask(sp, dtype):
171
    f = ift.from_random(sp).val
Philipp Arras's avatar
Philipp Arras committed
172 173
    mask = np.zeros_like(f)
    mask[f > 0] = 1
Martin Reinecke's avatar
Martin Reinecke committed
174
    mask = ift.Field.from_raw(sp, mask)
Philipp Arras's avatar
Philipp Arras committed
175 176 177 178 179 180
    op = ift.MaskOperator(mask)
    ift.extra.consistency_check(op, dtype, dtype)


@pmp('sp', _h_spaces + _p_spaces)
def testDiagonal(sp, dtype):
181
    op = ift.DiagonalOperator(ift.Field.from_random(sp, dtype=dtype))
Philipp Arras's avatar
Philipp Arras committed
182 183 184 185 186 187 188 189 190 191 192 193
    ift.extra.consistency_check(op, dtype, dtype)


@pmp('sp', _h_spaces + _p_spaces + _pow_spaces)
def testGeometryRemover(sp, dtype):
    op = ift.GeometryRemover(sp)
    ift.extra.consistency_check(op, dtype, dtype)


@pmp('spaces', [0, 1, 2, 3, (0, 1), (0, 2), (0, 1, 2), (0, 2, 3), (1, 3)])
@pmp('wgt', [0, 1, 2, -1])
def testContractionOperator(spaces, wgt, dtype):
Philipp Arras's avatar
Philipp Arras committed
194
    dom = (ift.RGSpace(1), ift.RGSpace(2), ift.GLSpace(3), ift.HPSpace(2))
Philipp Arras's avatar
Philipp Arras committed
195 196 197 198 199
    op = ift.ContractionOperator(dom, spaces, wgt)
    ift.extra.consistency_check(op, dtype, dtype)


def testDomainTupleFieldInserter():
200 201
    target = ift.DomainTuple.make((ift.UnstructuredDomain([3, 2]),
                                   ift.UnstructuredDomain(7),
Philipp Arras's avatar
Philipp Arras committed
202
                                   ift.RGSpace([4, 22])))
203
    op = ift.DomainTupleFieldInserter(target, 1, (5,))
Philipp Arras's avatar
Philipp Arras committed
204 205 206 207 208 209 210
    ift.extra.consistency_check(op)


@pmp('space', [0, 2])
@pmp('factor', [1, 2, 2.7])
@pmp('central', [False, True])
def testZeroPadder(space, factor, dtype, central):
Philipp Arras's avatar
Philipp Arras committed
211 212
    dom = (ift.RGSpace(4), ift.UnstructuredDomain(5), ift.RGSpace(3, 4),
           ift.HPSpace(2))
Philipp Arras's avatar
Philipp Arras committed
213
    newshape = [int(factor*ll) for ll in dom[space].shape]
Philipp Arras's avatar
Philipp Arras committed
214 215 216 217
    op = ift.FieldZeroPadder(dom, newshape, space, central)
    ift.extra.consistency_check(op, dtype, dtype)


Philipp Arras's avatar
Philipp Arras committed
218 219 220 221
@pmp('args', [[ift.RGSpace((13, 52, 40)), (4, 6, 25), None],
              [ift.RGSpace((128, 128)), (45, 48), 0],
              [ift.RGSpace(13), (7,), None],
              [(ift.HPSpace(3), ift.RGSpace((12, 24), distances=0.3)), (12, 12), 1]])
Philipp Arras's avatar
Philipp Arras committed
222 223 224 225
def testRegridding(args):
    op = ift.RegriddingOperator(*args)
    ift.extra.consistency_check(op)

Martin Reinecke's avatar
Martin Reinecke committed
226

Philipp Arras's avatar
Philipp Arras committed
227 228 229 230
@pmp('fdomain', [(ift.RGSpace((3, 5, 4)), ift.RGSpace((16,), distances=(7.,))),
                 ift.HPSpace(12)])
@pmp('domain', [(ift.RGSpace(2), ift.GLSpace(10)),
                ift.RGSpace((10, 12), distances=(0.1, 1.))])
Philipp Arras's avatar
Philipp Arras committed
231
def testOuter(fdomain, domain):
232 233
    f = ift.from_random(fdomain, 'normal')
    op = ift.OuterProduct(domain, f)
Philipp Arras's avatar
Philipp Arras committed
234
    ift.extra.consistency_check(op)
Philipp Arras's avatar
Philipp Arras committed
235 236 237


@pmp('sp', _h_spaces + _p_spaces + _pow_spaces)
Philipp Arras's avatar
Philipp Arras committed
238 239
@pmp('seed', [12, 3])
def testValueInserter(sp, seed):
Martin Reinecke's avatar
Martin Reinecke committed
240 241 242 243 244 245 246 247 248
    with ift.random.Context(seed):
        ind = []
        for ss in sp.shape:
            if ss == 1:
                ind.append(0)
            else:
                ind.append(int(ift.random.current_rng().integers(0, ss-1)))
        op = ift.ValueInserter(sp, ind)
        ift.extra.consistency_check(op)
249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266


@pmp('sp', _pow_spaces)
def testSlopeRemover(sp):
    op = ift.library.correlated_fields._SlopeRemover(sp)
    ift.extra.consistency_check(op)


@pmp('sp', _pow_spaces)
def testTwoLogIntegrations(sp):
    op = ift.library.correlated_fields._TwoLogIntegrations(sp)
    ift.extra.consistency_check(op)


@pmp('sp', _h_spaces + _p_spaces + _pow_spaces)
def testSpecialSum(sp):
    op = ift.library.correlated_fields._SpecialSum(sp)
    ift.extra.consistency_check(op)
267

Philipp Arras's avatar
Philipp Arras committed
268

Philipp Arras's avatar
Reorder  
Philipp Arras committed
269 270 271 272 273 274 275 276
@pmp('dofdex', [(0,), (1,), (0, 1), (1, 0)])
def testCorFldDistr(dofdex):
    tgt = ift.UnstructuredDomain(len(dofdex))
    dom = ift.UnstructuredDomain(2)
    op = ift.library.correlated_fields._Distributor(dofdex, dom, tgt)
    ift.extra.consistency_check(op)


277
def metatestMatrixProductOperator(sp, mat_shape, seed, **kwargs):
Martin Reinecke's avatar
Martin Reinecke committed
278 279 280 281 282 283 284
    with ift.random.Context(seed):
        mat = ift.random.current_rng().standard_normal(mat_shape)
        op = ift.MatrixProductOperator(sp, mat, **kwargs)
        ift.extra.consistency_check(op)
        mat = mat + 1j*ift.random.current_rng().standard_normal(mat_shape)
        op = ift.MatrixProductOperator(sp, mat, **kwargs)
        ift.extra.consistency_check(op)
Philipp Arras's avatar
Philipp Arras committed
285

Philipp Arras's avatar
Philipp Arras committed
286

287 288 289 290 291 292
@pmp('sp', [ift.RGSpace(10)])
@pmp('spaces', [None, (0,)])
@pmp('seed', [12, 3])
def testMatrixProductOperator_1d(sp, spaces, seed):
    mat_shape = sp.shape * 2
    metatestMatrixProductOperator(sp, mat_shape, seed, spaces=spaces)
Philipp Arras's avatar
Philipp Arras committed
293

Philipp Arras's avatar
Philipp Arras committed
294

295 296
@pmp('sp', [ift.DomainTuple.make((ift.RGSpace((2)), ift.RGSpace((10))))])
@pmp('spaces', [(0,), (1,), (0, 1)])
297
@pmp('seed', [12, 3])
298 299 300 301 302 303 304
def testMatrixProductOperator_2d_spaces(sp, spaces, seed):
    appl_shape = []
    for sp_idx in spaces:
        appl_shape += sp[sp_idx].shape
    appl_shape = tuple(appl_shape)
    mat_shape = appl_shape * 2
    metatestMatrixProductOperator(sp, mat_shape, seed, spaces=spaces)
305

Philipp Arras's avatar
Philipp Arras committed
306

307 308 309 310 311 312
@pmp('sp', [ift.RGSpace((2, 10))])
@pmp('seed', [12, 3])
def testMatrixProductOperator_2d_flatten(sp, seed):
    appl_shape = (ift.utilities.my_product(sp.shape),)
    mat_shape = appl_shape * 2
    metatestMatrixProductOperator(sp, mat_shape, seed, flatten=True)
313

Philipp Arras's avatar
Philipp Arras committed
314

Philipp Arras's avatar
Philipp Arras committed
315 316
@pmp('seed', [12, 3])
def testPartialExtractor(seed):
Martin Reinecke's avatar
Martin Reinecke committed
317 318 319 320 321 322 323 324
    with ift.random.Context(seed):
        tgt = {'a': ift.RGSpace(1), 'b': ift.RGSpace(2)}
        dom = tgt.copy()
        dom['c'] = ift.RGSpace(3)
        dom = ift.MultiDomain.make(dom)
        tgt = ift.MultiDomain.make(tgt)
        op = ift.PartialExtractor(dom, tgt)
        ift.extra.consistency_check(op)
Philipp Arras's avatar
Philipp Arras committed
325

Philipp Arras's avatar
Philipp Arras committed
326

Philipp Arras's avatar
Philipp Arras committed
327 328 329 330 331
@pmp('seed', [12, 3])
def testSlowFieldAdapter(seed):
    dom = {'a': ift.RGSpace(1), 'b': ift.RGSpace(2)}
    op = ift.operators.simple_linear_operators._SlowFieldAdapter(dom, 'a')
    ift.extra.consistency_check(op)