[Help-gsl] Re: linear algebra doubt

[Rodney Sparapani]

Adhithya Ravichandran wrote:
> Hi
> I am using the linear algebra header for one of my simulations. I want to invert a matrix and i tried using the gsl LU invert function for real numbers. 
>
> Whatever be the order of the square matrix , only the last two rows of the output has valid 
> numbers. the rest of it is either "nan" or "inf" which is wrong.
>
> Am I using the function in a wrong manner? I have been using a permutation 
> matrix which is the identity matrix. do i need to change it?
>
> Adhithya

I had some trouble with LU also.  I'm not sure if there is a GSL bug or 
I'm doing something wrong.  I have had success with the SVD method.  So, 
I'd suggest giving that a try.

Rodney