diff --git a/src/jacobi.cpp b/src/jacobi.cpp index d597793..f40efc8 100644 --- a/src/jacobi.cpp +++ b/src/jacobi.cpp @@ -72,7 +72,7 @@ void jacobi_eigensolver(const arma::mat& A, // Use a do while so that max_offdiag gets a value before doing a comparison do { - max_offdiag = max_offdiag_symmetric(A_m, k, l); + max_offdiag = max_offdiag_symmetric(A_m, k, l); jacobi_rotate(A_m, R, k, l); } while (max_offdiag >= eps && ++iterations < maxiter);