Fixed bug in two_sided_eigs.

nekStab · Nov 17, 2023 · 87651b9 · 87651b9
1 parent 3418b19
commit 87651b9
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/IterativeSolvers.f90 b/src/IterativeSolvers.f90
@@ -413,7 +413,7 @@ subroutine two_sided_eigs(A, V, W, rvecs, lvecs, eigvals, residuals, info, nev,
     end do lanczos
 
     call sort_index(residuals, indices, reverse=.false.)
-    eigvals = eigvals(indices) ; rvecs = rvecs(:, indices) ; lvecs = lvecs(indices, :)
+    eigvals = eigvals(indices) ; rvecs = rvecs(:, indices) ; lvecs = lvecs(:, indices)
     return
   end subroutine two_sided_eigs