This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: product of a matrix by the inverse of a symmetric one
- From: Brian Gough <bjg at network-theory dot co dot uk>
- To: fjmolina at lbl dot gov
- Cc: <gsl-discuss at sources dot redhat dot com>
- Date: Wed, 16 Jul 2003 16:29:17 +0100
- Subject: Re: product of a matrix by the inverse of a symmetric one
- References: <16147.585.145727.634096@dhcp-63-193.cse.ucsc.edu>
Francisco J Molina writes:
>
> I am trying to multiply A by (B^(-1)) where B^(-1) is the inverse of B and B
> is symmetric.
> I think that the best way to do this is to use
>
> gsl_linalg_LU_decomp
>
> to get an LU decomposition of B and then use gsl_blas_dtrsm several times.
> Am I right?
Yes, that's what I would do.
--
Brian