Re: [Numpy-discussion] some work on arpack

[Nils W. <nw...@ia...>]

On Wed, 16 Aug 2006 08:08:16 -0700
  "David Grant" <dav...@gm...> wrote:
> On 8/16/06, Nils Wagner <nw...@ia...>=20
>wrote:
>>
>> Keith Goodman wrote:
>> > On 8/15/06, David Grant <dav...@gm...> wrote:
>> >
>> >
>> >> My idea is (if I have time) to write an eigs-like=20
>>function in python
>> >> that will only perform a subset of what Matlab's eigs=20
>>does for. It
>> >> will, for example, compute a certain number of=20
>>eigenvalues and
>> >> eigenvectors for a real, sparse, symmetric matrix=20
>>(the case I'm
>> >> interested in)
>> >>
>> >
>> >
>> AFAIK, pysparse (in the sandbox) includes a module that=20
>>implements a
>> Jacobi-Davidson
>> eigenvalue solver for the symmetric, generalised matrix=20
>>eigenvalue
>> problem (JDSYM).
>> Did someone test pysparse ?
>>
>> I did try pysparse a few years ago (I think right before=20
>>sparse stuff came
> into scipy). I think there is probably an old post=20
>asking the list about
> sparse stuff and I think Travis had just written it and=20
>told me about it...
> can't remember. Can JDSYM just return the k lowest=20
>eigenvalues/eigenvectors?
>=20
> --=20
> David Grant
> http://www.davidgrant.ca

Yes. See
  http://people.web.psi.ch/geus/pyfemax/pysparse_examples.html

for details.

Nils