This is the mail archive of the gsl-discuss@sources.redhat.com mailing list for the GSL project.

Index Nav:	[Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav:	[Date Prev] [Date Next]	[Thread Prev] [Thread Next]
Other format:	[Raw text]

Re: Bug in gsl_linalg_LU_decomp

From: Brian Gough <bjg at network-theory dot co dot uk>
To: Andrew Jason Penner <ajpenner at Physics dot UManitoba dot CA>
Cc: gsl-discuss at sources dot redhat dot com
Date: Thu, 8 Aug 2002 21:53:12 +0100 (BST)
Subject: Re: Bug in gsl_linalg_LU_decomp
References: <Pine.OSF.4.21.0208071352020.21271-100000@aurora.physics.umanitoba.ca>

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.

References:
- Bug in gsl_linalg_LU_decomp
  - From: Andrew Jason Penner

Index Nav:	[Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav:	[Date Prev] [Date Next]	[Thread Prev] [Thread Next]