|
From: | David Bateman |
Subject: | Re: Sparse-merge |
Date: | Wed, 09 Mar 2005 12:28:50 +0100 |
User-agent: | Mozilla Thunderbird 0.8 (X11/20040923) |
John W. Eaton wrote:
On 8-Mar-2005, Richard Hindmarsh <address@hidden> wrote:| Thanks. Will 2.9.x (or does the main trunk) have an interface to ARPACK | (e.g. through "eigs"), or is anyone working on this?Not yet. SEe the list of projects for sparse matrices in the file http://www.octave.org/cgi-bin/viewcvs.cgi/octave/PROJECTS?rev=1.60&view=auto David Bateman, who wrote much of the sprase code, may have some additional ideas about how to implement eigs. If you are interested in helping, please let us know. Thanks, jwe
ARPACK is definitely on the todo list. It is a complete lack of time at the moment that will prevent me from doing something here.
The priority for me at the moment is to complete the direct solver code, this means fix up the matrix type caching code, using TAUCS to implement the sparse Cholesky solver and identify the QR solver for the minimum norm solution. TAUCS has the advantage that it has drop tolerances and the cholinc function can also directly be implemented. I believe that Andy Adler wanted to look at this, so I won't touch this till he has had a chance at it.
If you want to do the ARPACK interfacing, I'd welcome the help and would be willing to vet the code and give pointers.
Regards David -- David Bateman address@hiddenMotorola Labs - Paris +33 1 69 35 48 04 (Ph) Parc Les Algorithmes, Commune de St Aubin +33 1 69 35 77 01 (Fax) 91193 Gif-Sur-Yvette FRANCE
The information contained in this communication has been classified as: [x] General Business Information [ ] Motorola Internal Use Only [ ] Motorola Confidential Proprietary
------------------------------------------------------------- Octave is freely available under the terms of the GNU GPL. Octave's home on the web: http://www.octave.org How to fund new projects: http://www.octave.org/funding.html Subscription information: http://www.octave.org/archive.html -------------------------------------------------------------
[Prev in Thread] | Current Thread | [Next in Thread] |