[Bug-gsl] Output of the gsl_linalg_complex_LU_solve

From:

Ho, Don T.

Subject:

[Bug-gsl] Output of the gsl_linalg_complex_LU_solve

Date:

Thu, 22 Mar 2007 15:56:16 -0400

The result of solving for A*x = B using the gsl_linalg_complex_LU_sovle
is wrong.
But if I do the inverse of A using gsl_linalg_complex_LU_invert() and
then take B*(A)^-1 then the output x is correct.
I didn't have the source code to play with.