[Help-gsl] linear system A x = b

[address@hidden]

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?


Thanks


--Laura