[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] gsl_linalg_LU_decomp and singular matrices
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] gsl_linalg_LU_decomp and singular matrices |
Date: |
Mon, 15 Jun 2009 10:53:59 +0100 |
User-agent: |
Wanderlust/2.14.0 (Africa) Emacs/22.2 Mule/5.0 (SAKAKI) |
At Fri, 12 Jun 2009 18:57:40 -0400,
Erwin Kalvelagen wrote:
> gsl_linalg_LU_decomp should return an error code when the input matrix is
> singular. See any numerical library for examples. I believe ignoring this
> makes this routine of little value.
Thanks for your email. The LU decomposition is valid for a singular
matrix so we don't return an error in that case.
It would make sense to return an error code from gsl_linalg_LU_svx
though, I will add one there.
--
Brian Gough
(GSL Maintainer)
Support freedom by joining the FSF
http://www.fsf.org/associate/support_freedom/join_fsf?referrer=37