This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: multivariate gaussian distribution (Code)
Hi Brian,
I remember reading that the matrix decomposition through
eigensystems was more 'stable', while Cholesky was
probably faster. See the implementation in R:
http://rweb.stat.umn.edu/R/library/MASS/html/mvrnorm.html
I used the LU for determinant and matrix inverse only.
I'm not an expert in this field though...
Emmanuel
P.S.: in my opinion the multivariate_gaussian_pdf_LU is rather
useless, because for repeated calls, you want the inverse
to be computed once and for all.
On Mon, Dec 29, 2003 at 11:41:10AM +0000, Brian Gough wrote:
> Emmanuel Benazera writes:
> > Hi,
> >
> > here is a tarball. Hope it helps.
> >
>
> Is there any advantage of LU over Cholesky decomposition?
>
> --
> Brian Gough
>
> Network Theory Ltd
> 15 Royal Park
> Bristol BS8 3AL
> United Kingdom
>
> Tel: +44 (0)117 3179309
> Fax: +44 (0)117 9048108
> Web: http://www.network-theory.co.uk/