[Help-gsl] Question on gsl_blas_strsm

[Johannes Wohlmuth]

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);


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?

I show you an example: 
A=
198.417         0.654522
0.654522        202.07
0.775645        -0.0124067
0.345305        1.88726
B=
198.517         0.654522
0.654522        202.17
C= (calculated with lower part)
0.999496        1.63098e-06
0.00329705      0.999495
0.00390719      -7.4017e-05
0.00173942      0.00932939
C= (calculated with upper part)
0.999486        0.00323749
1.63092e-06     0.999505
0.00390739      -6.13675e-05
0.00170864      0.00933502
C=(calculate with octave)
   9.9950e-01    1.6309e-06
    1.6309e-06    9.9951e-01
    3.9074e-03   -7.4018e-05
    1.7087e-03    9.3295e-03


vielen Dank für die Hilfe
Johannes
_______________________________________________________________________
Viren-Scan für Ihren PC! Jetzt für jeden. Sofort, online und kostenlos.
Gleich testen! http://www.pc-sicherheit.web.de/freescan/?mc=022222