Re: [SciPy-dev] Current status of arpack.speigs.eigvals
by Neilen Marais other posts by this author
Nov 20 2006 8:36AM messages near this date
[SciPy-dev] Current status of arpack.speigs.eigvals
|
[SciPy-dev] ARPACK wrapper
You'll find it works better for big matrices, since that's what the ARPACK
"tunables" were chosen for. Anyway, no, it's far from ready for primetime.
I've used it successfully to solve my problems though. Patches appreciated
(though wait for my update later today).
Regards
Neilen
On Mon, 20 Nov 2006 14:07:03 +0100, Nils Wagner wrote:
> Hi all,
>
> I am wondering if arpack.speigs.eigvals is ready to compute eigenvalues.
>
> a = random.rand(10,10)
> matvec = lambda x:dot(a,x)
> nev = 4
> ncv = 2*nev
> w1,v1 = arpack.speigs.eigvals(matvec,a.shape[0],nev=nev,ncv=ncv)
>
> Traceback (most recent call last):
> File "test_eigs.py", line 34, in ?
> w1,v1 = arpack.speigs.eigvals(matvec,a.shape[0],nev=nev,ncv=ncv)
> File
> "/usr/lib64/python2.4/site-packages/scipy/sandbox/arpack/speigs.py",
> line 71, in eigvals
> if info != 0: raise "Hell" # Indicates some error during the
> Arnouldi iterations
> Hell
>
>
> Am I missing something ?
>
>
> Nils
--
you know its kind of tragic
we live in the new world
but we've lost the magic
-- Battery 9 (www.battery9.co.za)
_______________________________________________
Scipy-dev mailing list
Scipy-dev@[...].org
http://projects.scipy.org/mailman/listinfo/scipy-dev
Thread:
Nils Wagner
Neilen Marais
|