[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