This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
linear system
- From: "lauradiara\ at libero\ dot it" <lauradiara at libero dot it>
- To: "gsl-discuss" <gsl-discuss at sources dot redhat dot com>
- Date: Thu, 27 Nov 2003 10:14:03 +0100
- Subject: linear system
Hi all,
I've implemented a function which solves linear system A x = b with the user input.
Reading the docs I have noticed that in gsl_linalg.h there are two very similar
methods which differ for the type of input parameters: gsl_linalg_LU_solve and
gsl_linalg_complex_LU_solve.
What are the differences? (I know that the second one works on complex numbers.)
Could I use only gsl_linalg_complex_LU_solve? If the user insert double numbers,
I insert the img part of complex number=0. Is it true?
What are the real usages of gsl_linalg_complex_LU_solve?
Other question: how can I solve linear system A x = b with indinite solutions?
For example: we suppose to have:
x+2y+3z=3
2x+3y+8z=4
3x+2y+17z=1
Thi system have infinite solutions:
x=-1-7t
y=2+2t
z=t t in R
Can i use gsl_linalg_LU_solve or gsl_linalg_complex_LU_solve to solve this kind
of linear system?
Cheers and thanks for your support
--Laura