[Numpy-discussion] some work on arpack
Nils Wagner
nwagner at iam.uni-stuttgart.de
Wed Aug 16 12:50:05 EDT 2006
On Wed, 16 Aug 2006 08:08:16 -0700
"David Grant" <davidgrant at gmail.com> wrote:
> On 8/16/06, Nils Wagner <nwagner at iam.uni-stuttgart.de>
>wrote:
>>
>> Keith Goodman wrote:
>> > On 8/15/06, David Grant <davidgrant at gmail.com> wrote:
>> >
>> >
>> >> My idea is (if I have time) to write an eigs-like
>>function in python
>> >> that will only perform a subset of what Matlab's eigs
>>does for. It
>> >> will, for example, compute a certain number of
>>eigenvalues and
>> >> eigenvectors for a real, sparse, symmetric matrix
>>(the case I'm
>> >> interested in)
>> >>
>> >
>> >
>> AFAIK, pysparse (in the sandbox) includes a module that
>>implements a
>> Jacobi-Davidson
>> eigenvalue solver for the symmetric, generalised matrix
>>eigenvalue
>> problem (JDSYM).
>> Did someone test pysparse ?
>>
>> I did try pysparse a few years ago (I think right before
>>sparse stuff came
> into scipy). I think there is probably an old post
>asking the list about
> sparse stuff and I think Travis had just written it and
>told me about it...
> can't remember. Can JDSYM just return the k lowest
>eigenvalues/eigenvectors?
>
> --
> David Grant
> http://www.davidgrant.ca
Yes. See
http://people.web.psi.ch/geus/pyfemax/pysparse_examples.html
for details.
Nils
More information about the NumPy-Discussion
mailing list