Bug in gsl_linalg_LU_decomp

[Andrew Jason Penner]

	The bug is not major, I accidentially used a singular matrix for
the matrix to LU decompose, it created a LU decomposed matrix, and then I
used that matrix in gsl_linalg_LU_solver, and I did get a result.

	To the best of my knowledge, this should produce a warning, since
singular matricies should not decompose under LU techniques. If I was
more of a computer scientist I would attempt to fix this myself, however I
am not, so I leave the glory of the patch to someone more skilled.


Andrew Jason Penner
ajpenner@physics.umanitoba.ca

*********************************
Physics is the art of life. 
*********************************