[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: eigs arpack arpack++ and Debian
From: |
Christophe Prud'homme |
Subject: |
Re: eigs arpack arpack++ and Debian |
Date: |
Tue, 21 Mar 2006 15:47:30 +0100 |
User-agent: |
KMail/1.9.1 |
David,
| >I will have a look
| >I never wanted to use arpack directly but rather arpack++ which simplifies
| >things a lot. Moreover you have different strategies depending on the
| > nature of the problem which is easier for a human being to handle with
| > arpack++
|
| I tried to avoid arpack++ to try and avoid the need for an additional
| dependency. I suppose that arpack++ would be ok, if it made things
| simpler...
actually arpack++ uses mostly C++ templates
only a 2/3 files are actually compiled that all
| I handled all of the cases in eigs.cc that are handled by the same
| function in matlab, and that is what I meant by a "full
| implementation"... I'm missing cholB flags treatment, as I didn't have a
| sparse cholesky at the time I wrote eigs.cc. I do now and so that should
| be easy to add. There is also another small case that wasn't implemented
| for reasons that escape me at the moment.
yeah ufsparse support is of great help here otherwise you can forget about
doing efficient/reliable eigenvalue solves specially in the case of
generalized eigenvalue problems
| >| If you have experience with arpack, then perhaps you might like to take
| >| a look at this code and see if you can see what I did wrong...
| >
| >I have some experience with arpack/arpack++ which has been positive so far
| >As I said in my previous email I have a nice testsuite in arpack++ for a
| > wide range of problems that could be helpful to debug eigs() in octave.
|
| This is great news... The test I have are octave specific, and should
| all run (but don't). If your tests allow comparison with similar results
| in octave, that would be a great help.
I downloaded arpack++ from debian, I do
configure;make;make check
attached you have the logfile of "make check" in the examples/product/sym
where we require only access to matrix/vector product and the matrix is
symmetric. Only a few examples/tests among many others (band/(non)sym,
dense/(non)sym...)
there are also some matrices construction in examples/matrices/{(non)sym |
complex} that could be implemented in octave.
In any case I will see what I can do to help. First I need to understand your
code (eigs.cc)
Best regards
C.
--
Christophe Prud'homme
EPFL SB IACS CMCS
MA B2 534 (Bâtiment MA)
Station 8
CH-1015 Lausanne
Tel: +41 (0)21 693 25 47
Fax: +41 (0)21 693 43 03
log
Description: Text Data