Displaying 1 result from an estimated 1 matches for "det_ge_real".
2022 Nov 09
1
det(diag(c(NaN, 1))) should be NaN, not 0
...0
[2,] NaN 0
> det(C)
[1] NaN
> (D <- t(C))
[,1] [,2]
[1,] NaN NaN
[2,] 0 0
> det(D)
[1] 0
Furthermore, the proposed change should _not_ decrease the performance
of determinant(A) for nonsingular 'A' ...
For those looking, the relevant C-level function is det_ge_real(),
defined in R-devel/src/modules/lapack/Lapack.c (at line 1260 in r83320).
Mikael
[1] https://github.com/Reference-LAPACK/lapack/blob/master/SRC/dgetrf.f