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. _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
