help-gsl
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Help-gsl] Question on gsl_blas_strsm


From: Brian Gough
Subject: Re: [Help-gsl] Question on gsl_blas_strsm
Date: Tue, 24 Oct 2006 09:51:17 +0100
User-agent: Mozilla/5.0 (X11; U; Linux i686; en-US; rv:1.8.0.6) Gecko/20060815 SeaMonkey/1.0.4

Johannes Wohlmuth wrote:
Hello,

I need some help with a Matrixproblem:
I want to solve C= A/B! B is a symmetric and positive definite
Matrix. I don't understand why there isn't any function like gsl_blas_ssysm! gsl_blas_strsm uses only the upper or the lower part of B? Is the following still the right function for my problem?
gsl_blas_strsm (CblasRight, CblasUpper, CblasTrans, CblasNonUnit,1.0, B, A);

Hello,

The blas functions are part of a standard, so we just implement what is in the standard. You can use the Cholesky decomposition to convert a symmetric positive definite system to a triangular system.

Another thing is: Why is there a difference when I use the upper or the lower part of the matrix B, when that
is symmetric?

If you still have this problem can you submit it as a bug report with an example program to address@hidden Thanks.

--
best regards,

Brian Gough
(GSL Maintainer)

Network Theory Ltd,
Publishing the GSL Manual - http://www.network-theory.co.uk/gsl/manual/




reply via email to

[Prev in Thread] Current Thread [Next in Thread]