[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Bug-gsl] bug in file eigen/jacobi.c

From: Brian Gough
Subject: Re: [Bug-gsl] bug in file eigen/jacobi.c
Date: Tue, 13 Jan 2009 21:33:24 +0000
User-agent: Wanderlust/2.14.0 (Africa) Emacs/22.2 Mule/5.0 (SAKAKI)

At Wed, 24 Dec 2008 14:59:29 +0100,
stefan wrote:
> I have found a bug in the function:
> int gsl_eigen_jacobi (gsl_matrix * a, gsl_vector * eval, gsl_matrix * 
> evec, unsigned int max_rot, unsigned int *nrot)
> The Problem is that the routine does not stop if a good solution is 
> found! That means if I set the parameter max_rot for example to 1000 and 
> i call the function then nrot is set to 1000, so nrot is always max_rot. 
> I figured out the the function:


This function only works for symmetric matrices.  The example program
has an nonsymmetric matrix so it will not converge.  For non-symmetric
matrices use gsl_eigen_nonsymm instead.

best regards

Brian Gough
(GSL Maintainer)

Support freedom by joining the FSF!

reply via email to

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