[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