This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
linear system problem
- From: <e70443 at estudiante dot uam dot es>
- To: gsl-discuss at sources dot redhat dot com
- Date: Thu, 29 Jan 2004 15:39:36 Europe/Madrid
- Subject: linear system problem
Dear sirs,
I am a student of Computer Engineering at Madrid
University, Spain.
First of all I want to congratulate you for the
very useful "Scientific Library GSL" you have.
The matter that concerns me is:
I need to solve the linear system A x = b, where A
is a rectangular M-by-N matrix, and M > N. Looking the
functions of the library, I decided to use de QR descomposition.
Before any calculation, I create the A matrix and b vector view, the 'x'
and 'tau' vectors of N size, and the 'residual' vector of M size.
First, I call the function gsl_linalg_QR_decomp() with the matrix A and
the vector 'tau'. Then I call the gsl_linalg_QR_lssolve() function with
the arguments A (with the QR descomposition created by the previous
function), tau, b, x, and residual, to find the solution to the linear
system.
I´m almost sure that these equations are correct but
the solution is not.
What should I expect ?
My question is; am I either using the right functions, or puting them in
the right order?
Shall I use another kind of descomposition?
Thank you very much for your cooperation
Sincerely,
Jaime Carsi Jaudenes.
---------------------------------------------
This message was sent using Endymion MailMan.
http://www.endymion.com/products/mailman/