From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 13738 invoked by alias); 14 Aug 2006 14:24:19 -0000 Received: (qmail 13730 invoked by uid 22791); 14 Aug 2006 14:24:18 -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; Mon, 14 Aug 2006 15:24:13 +0100 Date: Mon, 14 Aug 2006 14:24:00 -0000 Message-ID: <874pwftmyf.wl%bjg@network-theory.co.uk> From: Brian Gough 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> User-Agent: Wanderlust/2.14.0 (Africa) Emacs/21.3 Mule/5.0 (SAKAKI) Reply-To: gsl-discuss@sourceware.org MIME-Version: 1.0 (generated by SEMI 1.14.6 - "Maruoka") Content-Type: text/plain; charset=US-ASCII X-Message-Mac: 34ed8117e2361a1780803f301889275b 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-q3/txt/msg00014.txt.bz2 At Tue, 13 Jun 2006 09:56:20 -0600, > 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. > > Attached is the latest (final? :)) patch. The latest patch is in sources.redhat.com CVS with some extra cleanups. It works well. I modified the balance() routine to return the diagonal scaling elements D (for computing the eigenvectors later) and moved the memory allocation up out of hessenberg() into the workspace. -- Brian Gough