jdqr: Pass the interface to orthonormalize().

BIMAU · Nov 12, 2023 · b3a0c93 · b3a0c93
1 parent c6643b0
commit b3a0c93
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/jadapy/jdqr.py b/jadapy/jdqr.py
@@ -117,7 +117,7 @@ def jdqr(A, num=5, target=Target.SmallestMagnitude, tol=1e-8, lock_tol=None, M=N
                     A, prec, Q[:, 0:k+nev], Y[:, 0:k+nev], H[0:k+nev, 0:k+nev],
                     sigma, r[:, 0:nev], solver_tolerance, solver_maxit, interface)
 
-            orthonormalize(V[:, 0:m], V[:, m:m+nev], M=M, MV=None if MV is None else MV[:, 0:m])
+            orthonormalize(V[:, 0:m], V[:, m:m+nev], M=M, MV=None if MV is None else MV[:, 0:m], interface=interface)
 
         AV[:, m:m+nev] = A @ V[:, m:m+nev]
         if M is not None: