Resolve incorrect handling of BlockArray reduction (see #408)

lanl · Apr 20, 2023 · 92b56f3 · 92b56f3
1 parent dd43264
commit 92b56f3
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/scico/linop/_util.py b/scico/linop/_util.py
@@ -154,8 +154,8 @@ def valid_adjoint(
 
     u = A(x)
     v = AT(y)
-    yTu = snp.dot(y.ravel().conj(), u.ravel())  # type: ignore
-    vTx = snp.dot(v.ravel().conj(), x.ravel())  # type: ignore
+    yTu = snp.sum(y.conj() * u)  # type: ignore
+    vTx = snp.sum(v.conj() * x)  # type: ignore
     err = snp.abs(yTu - vTx) / max(snp.abs(yTu), snp.abs(vTx))
     if eps is None:
         return err