Problem with matrices bigger than 2x2 using gsl_linalg_LU_invert

[Brian Gough]

Philippe Huber writes:
 > First of all, thank you for having created such a great library. I have
 > a small problem with gsl_linalg_LU_invert. I am using a 4x4 matrix which
 > is not singular. I made a LU decomposition using gsl_linalg_LU_decomp.
 > Then, I send it to gsl_linalg_LU_invert and it crashes. I tried with a
 > 2x2 matrix and it works perfectly well.  I verified that the program also
 > crashes with a 3x3 matrix. Is it a bug or did I make something wrong.
 > Please apologize if I don't send my code, it is well too big.

Please send a small example which reproduces the bug, with details of 
your operating system, compiler, etc.
regards
Brian Gough