[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.

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

reply via email to

[Prev in Thread] Current Thread [Next in Thread]