Commit 15eec065 authored by Philipp Arras's avatar Philipp Arras

Change naming

parent a6f727db
...@@ -29,28 +29,29 @@ from .linear_operator import LinearOperator ...@@ -29,28 +29,29 @@ from .linear_operator import LinearOperator
class MaskOperator(LinearOperator): class MaskOperator(LinearOperator):
"""Implementation of a mask response """Implementation of a mask response
Takes a field, applies a mask and returns the values of the field in a Takes a field, applies flags and returns the values of the field in a
UnstructuredDomain. It can be used as response operator. :class:`UnstructuredDomain`.
Parameters Parameters
---------- ----------
mask : Field flags : Field
Is converted to boolean. Where True, the input field is flagged.
""" """
def __init__(self, mask): def __init__(self, flags):
if not isinstance(mask, Field): if not isinstance(flags, Field):
raise TypeError raise TypeError
self._domain = DomainTuple.make(mask.domain) self._domain = DomainTuple.make(flags.domain)
self._mask = np.logical_not(mask.to_global_data()) self._flags = np.logical_not(flags.to_global_data())
self._target = DomainTuple.make(UnstructuredDomain(self._mask.sum())) self._target = DomainTuple.make(UnstructuredDomain(self._flags.sum()))
self._capability = self.TIMES | self.ADJOINT_TIMES self._capability = self.TIMES | self.ADJOINT_TIMES
def apply(self, x, mode): def apply(self, x, mode):
self._check_input(x, mode) self._check_input(x, mode)
x = x.to_global_data() x = x.to_global_data()
if mode == self.TIMES: if mode == self.TIMES:
res = x[self._mask] res = x[self._flags]
return Field.from_global_data(, res) return Field.from_global_data(, res)
res = np.empty(self.domain.shape, x.dtype) res = np.empty(self.domain.shape, x.dtype)
res[self._mask] = x res[self._flags] = x
res[~self._mask] = 0 res[~self._flags] = 0
return Field.from_global_data(self.domain, res) return Field.from_global_data(self.domain, res)
Markdown is supported
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