[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Help-gsl] problem in using LU decomp
|
From: |
Brian Gough |
|
Subject: |
Re: [Help-gsl] problem in using LU decomp |
|
Date: |
Fri, 22 May 2009 16:03:06 +0100 |
|
User-agent: |
Wanderlust/2.14.0 (Africa) Emacs/22.1 Mule/5.0 (SAKAKI) |
At Wed, 20 May 2009 22:38:01 +0200,
Serle Shuman wrote:
> Can anyone see what I am doing wrong. I seem to get and infinite v[0] value
> when it should be -0.0000015531.
Hi,
The matrix is singular so there's no unique solution here. The blas
triangular solver routines return nan/inf in this case.
> Also can one use LU decomposition solver on a non-square matrix or does one
> need to pad it up with zeros.
Requirements on the matrix should be in the LU documentation in the
manual, if not please file a bug report to address@hidden, thanks.
--
Brian Gough
(GSL Maintainer)
Support freedom by joining the FSF
http://www.fsf.org/associate/support_freedom/join_fsf?referrer=37