This is the mail archive of the
gsl-discuss@sourceware.cygnus.com
mailing list for the GSL project.
Re: determinant.
- To: GSL discussion list <gsl-discuss at sourceware dot cygnus dot com>
- Subject: Re: determinant.
- From: Gerard Jungman <jungman at lanl dot gov>
- Date: Sun, 23 Jan 2000 01:12:34 -0700
- Organization: LANL T-8
- References: <20000123080636.26432.qmail@scam.xcf.berkeley.edu>
ali@scam.xcf.berkeley.edu wrote:
>
> Hi. I can't seem to find a general routine for computing determinants or
> inverses of a matrix in gsl. I'm also shocked to see that BLAS doesn't
> provide a determinant.
The determinant is a nearly useless object for numerical linear algebra.
libgsleigen provides one method for inversion, jacobi iteration:
int
gsl_eigen_invert_jacobi_impl(const gsl_matrix * matrix,
gsl_matrix * ainv,
unsigned int max_rot);
It may not have made it into the documentation yet.
--
G. Jungman