From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 29742 invoked by alias); 14 Jun 2006 20:53:28 -0000 Received: (qmail 29733 invoked by uid 22791); 14 Jun 2006 20:53:27 -0000 X-Spam-Check-By: sourceware.org Received: from network-theory.com (HELO mail.network-theory.co.uk) (66.199.228.187) by sourceware.org (qpsmtpd/0.31) with ESMTP; Wed, 14 Jun 2006 20:53:25 +0000 From: "Brian Gough" MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit Message-ID: <17552.30530.701370.888894@hp2.network-theory.co.uk> Date: Thu, 15 Jun 2006 01:17:00 -0000 To: Patrick Alken Cc: gsl-discuss@sourceware.org Subject: Re: unsymmetric eigenvalue update In-Reply-To: <20060613155620.GA21626@hippogriff.physics.drexel.edu> References: <20060613155620.GA21626@hippogriff.physics.drexel.edu> X-Return-Hmac: ab32c3ff0f1ca6eaeeaf2885ebf3da5d Mailing-List: contact gsl-discuss-help@sourceware.org; run by ezmlm Precedence: bulk List-Subscribe: List-Archive: List-Post: List-Help: , Sender: gsl-discuss-owner@sourceware.org X-SW-Source: 2006-q2/txt/msg00102.txt.bz2 Patrick Alken writes: > I have run my code through Brian Gough's test program > on 12x12 and 200x200 random matrices, as well as some integer > matrix cases and find no significant disagreement with LAPACK. > ..... > I propose adding the double-shift method to gsl (assuming it > passes any further tests needed) which would be a perfectly > fine eigenvalue solver (EISPACK used it for a number of years) > until I (or someone else) can find the time to implement method 3. Excellent work! Useful summary. I will add the code at the next opportunity. -- Brian Gough