[SciPy-dev] ARPACK wrapper

Nils Wagner nwagner at iam.uni-stuttgart.de
Mon Nov 20 02:45:47 CST 2006


Neilen Marais wrote:
> Hi 
>
> I've commited the beginnings of an  ARPACK
> wrapper partially adressing
> http://projects.scipy.org/scipy/scipy/ticket/231 in r2323
>
> I based it on the arpack-0.10.tar.bz2 source that aric attached to #231. I
> added an ARPACK driver that demonstrates how to use ARPACK's shift-invert
> mode to solve a generalised eigensystem. Unfortunately I haven't generated
> any compact test data for that yet (I was using matrices created by my own
> FEM app) but will soon.
>
> As things stand now a lot of work still remains to be done. I can think of
> (off the top of my head)
>
> 1) Merging Aric and my drivers
> 2) Handling more cases with an easy-to-use interface
> 3) Using the symmetric/complex/etc. modes
> 4) Many more, no doubt ;)
>
> I'm away for the rest of the weekend, but I look forward to comments on
> Monday :)
>
> Regards
> Neilen
>
>   
Hi Neilen,

Thank you very much for your hard work on the ARPACK wrapper.
Just now I have installed the sandbox package and tried to solve a
random eigenvalue
problem.

In order to compare the results with the workhorse eig I have used a
very small
order n. The number of desired eigenpairs is equal to k=4 in my example,
but the
shape of the array of eigenvectors is (n,k+1) and for the eigenvalues it
is (k+1,).
The eigenvectors returned by arpack.eigen are zero.

Nils


-------------- next part --------------
A non-text attachment was scrubbed...
Name: test_eigs.py
Type: text/x-python
Size: 133 bytes
Desc: not available
Url : http://projects.scipy.org/pipermail/scipy-dev/attachments/20061120/c55b2b12/attachment.py 


More information about the Scipy-dev mailing list