[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] About gsl_linalg_cholesky_svx(A,x)
From: |
Chuhu Yang |
Subject: |
[Bug-gsl] About gsl_linalg_cholesky_svx(A,x) |
Date: |
Sat, 16 Sep 2006 20:36:26 -0500 |
User-agent: |
Microsoft-Entourage/11.2.1.051004 |
Dear Mir or Madam,
Recently, I was using gsl. I found one problem with the function
gsl_linalg_cholesky_svx(A,x).
I want to obtain a solution to a linear function AX=b. I used
gsl_linalg_cholesky_decomp(A) first. I obtained a correct Cholesky factor L.
A after calling gsl_linalg_cholesky_decomp is a combination of L and L'.
Then I used gsl_linalg_cholesky_svx(A,X) with A from the first step. I
couldn't obtain correct solution for X. I guess there is a bug in
gsl_linalg_cholesky_svx(A,X).
However, I can get a correct solution by using
gsl_linalg_HH_svx(A,X) directly.
I hope this might help to make the function gsl_linalg_cholesky_svx better.
Thanks,
Chuhu
- [Bug-gsl] About gsl_linalg_cholesky_svx(A,x),
Chuhu Yang <=