Commit 7b445aed authored by Pumpe, Daniel (dpumpe)'s avatar Pumpe, Daniel (dpumpe)
Browse files

Update response_operator.py according to Martins hints and insights

parent 2938e54d
Pipeline #12431 passed with stages
in 11 minutes and 3 seconds
......@@ -80,27 +80,22 @@ class ResponseOperator(LinearOperator):
shape_target = np.append(shape_target, self._domain[ii].shape)
self._target = self._parse_domain(FieldArray(shape_target))
self._sigma = sigma
self._exposure = exposure
self._kernel = len(self._domain)*[None]
for ii in xrange(len(self._kernel)):
self._kernel[ii] = SmoothingOperator(self._domain[ii],
sigma=self._sigma[ii])
self._composed_kernel = ComposedOperator(self._kernel)
self._exposure_op = len(self._domain)*[None]
if len(self._exposure_op) != len(self._kernel):
raise ValueError("Definition of kernel and exposure do not suit "
"each other")
else:
for ii in xrange(len(self._exposure_op)):
self._exposure_op[ii] = DiagonalOperator(
self._domain[ii],
diagonal=self._exposure[ii])
self._composed_exposure = ComposedOperator(self._exposure_op)
kernel_smoothing = len(self._domain)*[None]
kernel_exposure = len(self._domain)*[None]
if len(sigma)!= len(exposure):
raise ValueError("Length of smoothing kernel and length of"
"exposure do not match")
for ii in xrange(len(kernel_smoothing)):
kernel_smoothing[ii] = SmoothingOperator(self._domain[ii],
sigma=sigma[ii])
kernel_exposure[ii] = DiagonalOperator(self._domain[ii],
diagonal=exposure[ii])
self._composed_kernel = ComposedOperator(kernel_smoothing)
self._composed_exposure = ComposedOperator(kernel_exposure)
@property
def domain(self):
......
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