# [SciPy-user] Logarithm of the determinant

Robert Kern rkern at ucsd.edu
Thu Jul 14 18:22:24 CDT 2005

Nils Wagner wrote:
> Hi all,
>
> The GSL library has a function
> gsl_linalg_LU_lndet and gsl_linalg_complex_LU_lndet.
>
> These functions compute the logarithm of the absolute value of the
> determinant of a matrix A, \ln|det(A)|, from its LU decomposition, LU.
> This function may be useful if the direct computation of the determinant
> would overflow or underflow.
>
> In scipy we have only
> det(a, overwrite_a=0)
>     det(a, overwrite_a=0) -> d
>
>     Return determinant of a square matrix.
>
> Is it possible to add the logarithm
> of a determinant in scipy ?

As soon as someone writes a suitable implementation. We can't directly
use GSL's because of its GPL license.

--
Robert Kern
rkern at ucsd.edu

"In the fields of hell where the grass grows high
Are the graves of dreams allowed to die."
-- Richard Harter