This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: Bug in gsl_linalg_LU_decomp
Andrew Jason Penner writes:
>
> 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.
Can you send a small example program which shows the problem, thanks.