[SciPy-user] Logarithm of the determinant

[Nils Wagner]

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 ?

Nils