[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

More information about the SciPy-user mailing list