diff --git a/src/afgl/util/lanczos.py b/src/afgl/util/lanczos.py
index 01c1f44..5dcc323 100644
--- a/src/afgl/util/lanczos.py
+++ b/src/afgl/util/lanczos.py
@@ -31,9 +31,9 @@ def lanczos(L, s, M):
         w = L @ V[:, j]
         alp[j] = np.dot(V[:, j], w)
 
-        w = w - V[:, j] * alp[j]
+        w = w - alp[j] * V[:, j]
         if j > 0:
-            w = w - V[:, j - 1] * beta[j - 1]
+            w = w - beta[j - 1] * V[:, j - 1]
 
         if j < M - 1:
             beta[j] = LA.norm(w)