From e89c69e109b11131d49950e59af192b556a2c425 Mon Sep 17 00:00:00 2001
From: Jens Renders <jens.renders@uantwerpen.be>
Date: Wed, 24 Mar 2021 22:44:51 +0100
Subject: Implemented _adjoint of OpTomo by equating it to _transpose
---
python/astra/optomo.py | 6 ++++++
1 file changed, 6 insertions(+)
diff --git a/python/astra/optomo.py b/python/astra/optomo.py
index ff67791..5ecff8f 100644
--- a/python/astra/optomo.py
+++ b/python/astra/optomo.py
@@ -95,6 +95,9 @@ class OpTomo(scipy.sparse.linalg.LinearOperator):
def _transpose(self):
return self.transposeOpTomo
+ # real operator
+ _adjoint = _transpose
+
def __checkArray(self, arr, shp):
if len(arr.shape)==1:
arr = arr.reshape(shp)
@@ -249,6 +252,9 @@ class OpTomoTranspose(scipy.sparse.linalg.LinearOperator):
def _transpose(self):
return self.parent
+ # real operator
+ _adjoint = _transpose
+
def __mul__(self,s):
# Catch the case of a backprojection of 2D/3D data
if isinstance(s, np.ndarray) and s.shape==self.parent.sshape:
--
cgit v1.2.3