From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 4441 invoked by alias); 4 Dec 2003 08:46:24 -0000 Mailing-List: contact gsl-discuss-help@sources.redhat.com; run by ezmlm Precedence: bulk List-Subscribe: List-Archive: List-Post: List-Help: , Sender: gsl-discuss-owner@sources.redhat.com Received: (qmail 4426 invoked from network); 4 Dec 2003 08:46:23 -0000 Received: from unknown (HELO pasmtp.tele.dk) (193.162.159.95) by sources.redhat.com with SMTP; 4 Dec 2003 08:46:23 -0000 Received: from herodot (unknown [62.242.210.149]) by pasmtp.tele.dk (Postfix) with ESMTP id 7246F1EC3EE for ; Thu, 4 Dec 2003 09:46:22 +0100 (CET) Received: from geisshirt.dk (unknown [62.242.210.146]) by herodot (Postfix) with ESMTP id 0772F27AD2 for ; Thu, 4 Dec 2003 09:46:21 +0100 (CET) Message-ID: <3FCEF411.4020501@geisshirt.dk> Date: Thu, 04 Dec 2003 08:46:00 -0000 From: Kenneth Geisshirt User-Agent: Mozilla/5.0 (X11; U; Linux i686; en-US; rv:1.5) Gecko/20031107 Debian/1.5-2 StumbleUpon/1.87 X-Accept-Language: da, en MIME-Version: 1.0 To: gsl-discuss@sources.redhat.com Subject: Re: Nonsymmetric Eigenproblems References: <3FCE0A5B.777CE630@mecha.uni-stuttgart.de> In-Reply-To: <3FCE0A5B.777CE630@mecha.uni-stuttgart.de> Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: 8bit X-SW-Source: 2003-q4/txt/msg00116.txt.bz2 Nils Wagner wrote: > I am wondering if there is any effort concerning the handling of > standard and generalized eigenvalue problems with > nonsymmetric matrices. > > Afaik, there are only routines for Hermitian (complex) > symmetric (real case) matrices in GSL-1.4. > > Any comments or suggestions ? Back in 1994/5 I ported some routines from Pascal (implementation due to my thesis supervisor who ported it from Algol). I haven't used it for years. The description in the header file says: Eigen is a library for computing eigenvalues and eigenvectors of general matrices. There is only one routine exported, namely Eigen. The meaning of the arguments to Eigen is: 1. The dimension of the general matrix (n). 2. A general matrix (A). 3. The maximal number of iterations. 4. The precision. 5. A vector with the eigenvalues. 6. A matrix with the eigenvectors. The source code is around 1000 lines of ISO C. I would be happy to contribute it to GSL. Kneth -- Kenneth Geisshirt, M.Sc., Ph.D. http://kenneth.geisshirt.dk Grøndals Parkvej 2A, 3. sal kenneth@geisshirt.dk DK-2720 Vanløse +45 38 87 78 38