Fix givens rotation.

This commit is contained in:
Antonio Sánchez 2025-03-14 17:15:57 +00:00
parent 33b04fe518
commit 14f845a1a8

View File

@ -75,7 +75,7 @@ void r1updt(Matrix<Scalar, Dynamic, Dynamic> &s, const Matrix<Scalar, Dynamic, 1
/* givens rotation. */ /* givens rotation. */
w_givens[j] = givens; w_givens[j] = givens;
} else } else
v_givens[j] = IdentityRotation; w_givens[j] = IdentityRotation;
/* test for zero diagonal elements in the output s. */ /* test for zero diagonal elements in the output s. */
if (s(j, j) == 0.) { if (s(j, j) == 0.) {