Re: eigs arpack arpack++ and Debian

[Christophe Prud'homme]

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