Commit 680a22ee authored by Philipp Arras's avatar Philipp Arras
Browse files

Remove ValueInserter again since it is not needed

parent 15617c75
......@@ -47,7 +47,6 @@ from .operators.outer_product_operator import OuterProduct
from .operators.simple_linear_operators import (
VdotOperator, ConjugationOperator, Realizer,
FieldAdapter, GeometryRemover, NullOperator)
from .operators.value_inserter import ValueInserter
from .operators.energy_operators import (
EnergyOperator, GaussianEnergy, PoissonianEnergy, InverseGammaLikelihood,
BernoulliEnergy, Hamiltonian, SampledKullbachLeiblerDivergence)
# 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
# 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 <>.
# Copyright(C) 2013-2018 Max-Planck-Society
# NIFTy is being developed at the Max-Planck-Institut fuer Astrophysik
# and financially supported by the Studienstiftung des deutschen Volkes.
from __future__ import absolute_import, division, print_function
import numpy as np
from ..compat import *
from ..domain_tuple import DomainTuple
from import UnstructuredDomain
from ..field import Field
from .linear_operator import LinearOperator
class ValueInserter(LinearOperator):
"""Operator which inserts one value into a field.
target : Domain, tuple of Domain or DomainTuple
index : tuple
The index of the target into which the value of the domain shall be
def __init__(self, target, index):
from ..sugar import makeDomain
self._domain = makeDomain(UnstructuredDomain(1))
self._target = DomainTuple.make(target)
if not isinstance(index, tuple):
raise TypeError
self._index = index
self._capability = self.TIMES | self.ADJOINT_TIMES
# Check if index is in bounds
def apply(self, x, mode):
self._check_input(x, mode)
x = x.to_global_data()
if mode == self.TIMES:
res = np.zeros(, dtype=x.dtype)
res[self._index] = x
res = np.full((1,),x[self._index])
return Field.from_global_data(self._tgt(mode), res)
......@@ -44,14 +44,6 @@ class Consistency_Tests(unittest.TestCase):
op = ift.LOSResponse(sp, starts, ends, sigma_low, sigma_ups)
ift.extra.consistency_check(op, dtype, dtype)
def testValueInserter(self):
op = ift.ValueInserter(ift.RGSpace([23, 44]), (2, 43))
lambda: ift.ValueInserter(ift.RGSpace(3), (7,)))
lambda: ift.ValueInserter(ift.RGSpace(3), 2))
@expand(product(_h_spaces + _p_spaces + _pow_spaces,
[np.float64, np.complex128]))
def testOperatorCombinations(self, sp, dtype):
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment