[Help-gsl] gsl_linalg_LU_invert

[Warren Weckesser]

Matthias,

Your variable m is a "matrix view" of the array data[].  When
you call gsl_linalg_LU_decomp, this data is overwritten.  To make
your code work, you can first save a copy of the matrix. I have
attached a modified version of your program that will print the
2x2 identity matrix, as expected.

Warren

gsl_inverse_test2.c
Description: Text Data