linear system

linear system

[lauradiara]

Hi all,

can someone help me?

I don't know how to solve linear system (Ax=b) which have det A = 0.

In this case the linear system can have nothing or infinite solutions.

How can I solve this kinkd of system with gsl?


Thanks

--Laura

linear system

[lauradiara]

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