some work on arpack
Building an arpack extension turned out to be surprisingly simple. For example for dsaupd: f2py c dsaupd.f m dsaupd L/usr/lib/blas/atlas:/usr/lib/lapack/atlas llapack lblas larpack It took me a long time to get the command down to something that simple. Took me a while even to figure out I could just use the arpack library on my computer rather than relinking all of arpack! I was able to import the dsaupd.so python module just fine and I was also able to call it just fine. I'll have to tweak the pyf file in order to get some proper output. But this gives me confidence that arpack is easy to hook into which is what others have said in the past, but without any experience with f2py I had no idea myself. f2py is awesome, for anyone who doesn't know. Matlab has interfaces for the arpack functions like dsaupd, dseupd, dnaupd, znaupd, zneupd (the mex file documentation claims those are the only ones, but they have more). Matlab has a C interface to these functions in arpackc.mex* and the script eigs.m does the grunt work, providing a very highlevel interface as well as doing some linear algebra (the same type of stuff that is done in arpack's examples directory I gather) and various other things. My idea is (if I have time) to write an eigslike 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)... I hope that this subsetofmatlab'seigs function will not be too hard to write. Then more functionality can be added on to eigs.py later... Does this make sense? Has anyone else started work on arpack integration at all?  David Grant http://www.davidgrant.ca
On 8/15/06, David Grant <davidgrant@gmail.com> wrote:
My idea is (if I have time) to write an eigslike 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)
Will it also work for a real, dense, symmetric matrix? That's the case I'm interested in. But even if it doesn't, your work is great news for numpy.
Keith Goodman wrote:
On 8/15/06, David Grant <davidgrant@gmail.com> wrote:
My idea is (if I have time) to write an eigslike 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 JacobiDavidson eigenvalue solver for the symmetric, generalised matrix eigenvalue problem (JDSYM). Did someone test pysparse ? Nils
Will it also work for a real, dense, symmetric matrix? That's the case I'm interested in. But even if it doesn't, your work is great news for numpy.
 Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with preintegrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.asus.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642 _______________________________________________ Numpydiscussion mailing list Numpydiscussion@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/numpydiscussion
Nils Wagner wrote:
Keith Goodman wrote:
On 8/15/06, David Grant <davidgrant@gmail.com> wrote:
My idea is (if I have time) to write an eigslike 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 JacobiDavidson eigenvalue solver for the symmetric, generalised matrix eigenvalue problem (JDSYM). Did someone test pysparse ?
Nils
Will it also work for a real, dense, symmetric matrix? That's the case I'm interested in. But even if it doesn't, your work is great news for numpy.
 Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with preintegrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.asus.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642 _______________________________________________ Numpydiscussion mailing list Numpydiscussion@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/numpydiscussion
 Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with preintegrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.asus.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642 _______________________________________________ Numpydiscussion mailing list Numpydiscussion@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/numpydiscussion
Ok it's not ready... gcc: Lib/sandbox/pysparse/src/spmatrixmodule.c In file included from Lib/sandbox/pysparse/src/spmatrixmodule.c:17: Lib/sandbox/pysparse/src/ll_mat.c: In function ‘LLMat_matvec_transp’: Lib/sandbox/pysparse/src/ll_mat.c:760: error: ‘CONTIGUOUS’ undeclared (first use in this function) Lib/sandbox/pysparse/src/ll_mat.c:760: error: (Each undeclared identifier is reported only once Lib/sandbox/pysparse/src/ll_mat.c:760: error: for each function it appears in.) Lib/sandbox/pysparse/src/ll_mat.c: In function ‘LLMat_matvec’: Lib/sandbox/pysparse/src/ll_mat.c:797: error: ‘CONTIGUOUS’ undeclared (first use in this function) In file included from Lib/sandbox/pysparse/src/spmatrixmodule.c:18: Lib/sandbox/pysparse/src/csr_mat.c: In function ‘CSRMat_matvec_transp’: Lib/sandbox/pysparse/src/csr_mat.c:119: error: ‘CONTIGUOUS’ undeclared (first use in this function) Lib/sandbox/pysparse/src/csr_mat.c: In function ‘CSRMat_matvec’: Lib/sandbox/pysparse/src/csr_mat.c:146: error: ‘CONTIGUOUS’ undeclared (first use in this function) In file included from Lib/sandbox/pysparse/src/spmatrixmodule.c:19: Lib/sandbox/pysparse/src/sss_mat.c: In function ‘SSSMat_matvec’: Lib/sandbox/pysparse/src/sss_mat.c:83: error: ‘CONTIGUOUS’ undeclared (first use in this function) In file included from Lib/sandbox/pysparse/src/spmatrixmodule.c:17: Lib/sandbox/pysparse/src/ll_mat.c: In function ‘LLMat_matvec_transp’: Lib/sandbox/pysparse/src/ll_mat.c:760: error: ‘CONTIGUOUS’ undeclared (first use in this function) Lib/sandbox/pysparse/src/ll_mat.c:760: error: (Each undeclared identifier is reported only once Lib/sandbox/pysparse/src/ll_mat.c:760: error: for each function it appears in.) Lib/sandbox/pysparse/src/ll_mat.c: In function ‘LLMat_matvec’: Lib/sandbox/pysparse/src/ll_mat.c:797: error: ‘CONTIGUOUS’ undeclared (first use in this function) In file included from Lib/sandbox/pysparse/src/spmatrixmodule.c:18: Lib/sandbox/pysparse/src/csr_mat.c: In function ‘CSRMat_matvec_transp’: Lib/sandbox/pysparse/src/csr_mat.c:119: error: ‘CONTIGUOUS’ undeclared (first use in this function) Lib/sandbox/pysparse/src/csr_mat.c: In function ‘CSRMat_matvec’: Lib/sandbox/pysparse/src/csr_mat.c:146: error: ‘CONTIGUOUS’ undeclared (first use in this function) In file included from Lib/sandbox/pysparse/src/spmatrixmodule.c:19: Lib/sandbox/pysparse/src/sss_mat.c: In function ‘SSSMat_matvec’: Lib/sandbox/pysparse/src/sss_mat.c:83: error: ‘CONTIGUOUS’ undeclared (first use in this function) error: Command "gcc pthread fnostrictaliasing DNDEBUG O2 fmessagelength=0 Wall D_FORTIFY_SOURCE=2 g fPIC ILib/sandbox/pysparse/include/ I/usr/lib64/python2.4/sitepackages/numpy/core/include I/usr/include/python2.4 c Lib/sandbox/pysparse/src/spmatrixmodule.c o build/temp.linuxx86_642.4/Lib/sandbox/pysparse/src/spmatrixmodule.o" failed with exit status 1 Nils
On 8/16/06, Nils Wagner <nwagner@iam.unistuttgart.de> wrote:
Keith Goodman wrote:
On 8/15/06, David Grant <davidgrant@gmail.com> wrote:
My idea is (if I have time) to write an eigslike 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 JacobiDavidson 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
On Wed, 16 Aug 2006 08:08:16 0700 "David Grant" <davidgrant@gmail.com> wrote:
On 8/16/06, Nils Wagner <nwagner@iam.unistuttgart.de> wrote:
Keith Goodman wrote:
On 8/15/06, David Grant <davidgrant@gmail.com> wrote:
My idea is (if I have time) to write an eigslike 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 JacobiDavidson 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
On 8/16/06, Keith Goodman <kwgoodman@gmail.com> wrote:
On 8/15/06, David Grant <davidgrant@gmail.com> wrote:
My idea is (if I have time) to write an eigslike 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)
Will it also work for a real, dense, symmetric matrix? That's the case I'm interested in. But even if it doesn't, your work is great news for numpy.
Real, dense, symmetric, well doesn't scipy already have something for this? I'm honestly not sure on the arpack side of things, I thought arpack was only useful (over other tools) for sparse matrices, I could be wrong.  David Grant http://www.davidgrant.ca
Hello all
Original Message From: numpydiscussionbounces@lists.sourceforge.net [mailto:numpy discussionbounces@lists.sourceforge.net] On Behalf Of David Grant Sent: 16 August 2006 17:11 To: Discussion of Numerical Python Subject: Re: [Numpydiscussion] some work on arpack
On 8/16/06, Keith Goodman <kwgoodman@gmail.com> wrote:
On 8/15/06, David Grant <davidgrant@gmail.com> wrote:
My idea is (if I have time) to write an eigslike 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)
Will it also work for a real, dense, symmetric matrix? That's the case I'm interested in. But even if it doesn't, your work is great news for numpy.
Real, dense, symmetric, well doesn't scipy already have something for this? I'm honestly not sure on the arpack side of things, I thought arpack was only useful (over other tools) for sparse matrices, I could be wrong.
Maybe SciPy can also do this, but what makes ARPACK useful is that it can get you a few eigenvalues and eigenvectors of a massive matrix without having to have the whole thing in memory. Instead, you provide ARPACK with a function that does A*x on your matrix. ARPACK passes a few x's to your function and a few eigenvalues and eigenvectors fall out. I recently used MATLAB's eigs to do exactly this. I had a dense matrix A with dimensions m x n, where m >> n. I wanted the eigenvalues of A'A (which has dimensions m x m, which is too large to keep in memory). But I could keep A and A' in memory I could quickly calculate A'A*x, which is what ARPACK needs. Cheers, Albert
On 8/16/06, Albert Strasheim <fullung@gmail.com> wrote:
Hello all
Original Message From: numpydiscussionbounces@lists.sourceforge.net [mailto:numpy discussionbounces@lists.sourceforge.net] On Behalf Of David Grant Sent: 16 August 2006 17:11 To: Discussion of Numerical Python Subject: Re: [Numpydiscussion] some work on arpack
On 8/16/06, Keith Goodman <kwgoodman@gmail.com> wrote:
On 8/15/06, David Grant <davidgrant@gmail.com> wrote:
> My idea is (if I have time) to write an eigslike 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)
Will it also work for a real, dense, symmetric matrix? That's the case I'm interested in. But even if it doesn't, your work is great news for numpy.
Real, dense, symmetric, well doesn't scipy already have something for this? I'm honestly not sure on the arpack side of things, I thought arpack was only useful (over other tools) for sparse matrices, I could be wrong.
Maybe SciPy can also do this, but what makes ARPACK useful is that it can get you a few eigenvalues and eigenvectors of a massive matrix without having to have the whole thing in memory. Instead, you provide ARPACK with a function that does A*x on your matrix. ARPACK passes a few x's to your function and a few eigenvalues and eigenvectors fall out.
Cool, thanks for the info.  David Grant http://www.davidgrant.ca
On 8/15/06, David Grant <davidgrant@gmail.com> wrote:
My idea is (if I have time) to write an eigslike 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)... I hope that this subsetofmatlab'seigs function will not be too hard to write. Then more functionality can be added on to eigs.py later... Does this make sense?
Did you, or anybody else on the list, have any luck making a numpy version of eigs?  Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with preintegrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.asus.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
Here's something I've been wondering: is it ok to port Matlab functions over to python? If so, then it's maybe an afternoon's work to get eigs working, given eigs.m and python Arpack wrappers. I'm 99% sure it would *not* be kosher to copypaste their code directly, but it's not a copy paste job because of the differences between Python and .m. So it seems a grayer area to me. bb On 10/22/06, Keith Goodman <kwgoodman@gmail.com> wrote:
On 8/15/06, David Grant <davidgrant@gmail.com> wrote:
My idea is (if I have time) to write an eigslike 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)... I hope that this subsetofmatlab'seigs function will not be too hard to write. Then more functionality can be added on to eigs.py later... Does this make sense?
Did you, or anybody else on the list, have any luck making a numpy version of eigs?
 Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with preintegrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.asus.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642 _______________________________________________ Numpydiscussion mailing list Numpydiscussion@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/numpydiscussion
 Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with preintegrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.asus.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
On 10/21/06, Bill Baxter <wbaxter@gmail.com> wrote:
Here's something I've been wondering: is it ok to port Matlab functions over to python? If so, then it's maybe an afternoon's work to get eigs working, given eigs.m and python Arpack wrappers.
I'm 99% sure it would *not* be kosher to copypaste their code directly, but it's not a copy paste job because of the differences between Python and .m. So it seems a grayer area to me.
bb
I haven't the faintest idea. What did the EULA say? What about a copyright? This sounds like a job for a lawyer, but I would guess it would be a bad idea unless the code is from a third party with a free license. Chuck  Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with preintegrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.asus.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642 _______________________________________________ Numpydiscussion mailing list Numpydiscussion@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/numpydiscussion
On 10/22/06, Charles R Harris <charlesr.harris@gmail.com> wrote:
On 10/21/06, Bill Baxter <wbaxter@gmail.com> wrote:
Here's something I've been wondering: is it ok to port Matlab functions over to python? If so, then it's maybe an afternoon's work to get eigs working, given eigs.m and python Arpack wrappers. [...]
I haven't the faintest idea. What did the EULA say? What about a copyright? This sounds like a job for a lawyer, but I would guess it would be a bad idea unless the code is from a third party with a free license.
Gee, there's an idea. :)  from License.txt  "Except as expressly provided by this Agreement, including the attached Addendum, Licensee may not adapt, translate, or convert "Mfiles", "MDLfiles" or "Pcode" contained in the Programs in order to create software, a principal purpose of which is to perform the same or similar functions as Programs licensed by MathWorks or which is intended to replace any component of the Programs. The Licensee may not incorporate or use "Mfiles", "Pcode", source code, or any other part of the Programs in or as part of another computer program without the consent of MathWorks. A Licensed User may modify pieces of MathWorks' code for the Licensed User's own use. A Licensed User may share such modified code with others provided each recipient is also a Licensed User for the original form of the code. "  So. If you're a Matlab licensee it would be ok to use it in your own stuff, but putting it in Numpy is definitely out of the question. bb  Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with preintegrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.asus.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
On Sun, Oct 22, 2006 at 09:27:39AM +0900, Bill Baxter wrote:
On 10/22/06, Charles R Harris <charlesr.harris@gmail.com> wrote:
On 10/21/06, Bill Baxter <wbaxter@gmail.com> wrote:
Here's something I've been wondering: is it ok to port Matlab functions over to python? If so, then it's maybe an afternoon's work to get eigs working, given eigs.m and python Arpack wrappers. [...]
I haven't the faintest idea. What did the EULA say? What about a copyright? This sounds like a job for a lawyer, but I would guess it would be a bad idea unless the code is from a third party with a free license.
Gee, there's an idea. :)
 from License.txt  [...] 
So. If you're a Matlab licensee it would be ok to use it in your own stuff, but putting it in Numpy is definitely out of the question.
A lot of the code in the Octave Forge is licensed under the BSD license, so you can use that. Unfortunately, the main Octave code is released under the GPL. Cheers Stéfan  Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with preintegrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.asus.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
On Sat, Oct 21, 2006 at 02:05:42PM 0700, Keith Goodman wrote:
Did you, or anybody else on the list, have any luck making a numpy version of eigs?
I made a start at an ARPACK wrapper, see http://projects.scipy.org/scipy/scipy/ticket/231 and the short thread at scipydev http://thread.gmane.org/gmane.comp.python.scientific.devel/5166/focus=5175 In addition to the wrapper there is a Python interface (and some tests). I don't know if the interface is like "eigs"  I don't use Matlab. It will give you a few eigenvalues and eigenvectors for the standard eigenproblem (Ax=lx) for any type of A (symmetric/nonsymmetric, real/complex, single/double, sparse/nonsparse). The generalized and shifted modes are not implemented. I need to find some time (or some help) to get it finished. Regards, Aric  Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with preintegrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.asus.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
Hello all
Original Message From: numpydiscussionbounces@lists.sourceforge.net [mailto:numpy discussionbounces@lists.sourceforge.net] On Behalf Of Aric Hagberg Sent: Sunday, October 22, 2006 4:09 PM To: Discussion of Numerical Python Subject: Re: [Numpydiscussion] some work on arpack
On Sat, Oct 21, 2006 at 02:05:42PM 0700, Keith Goodman wrote:
Did you, or anybody else on the list, have any luck making a numpy version of eigs?
I made a start at an ARPACK wrapper, see http://projects.scipy.org/scipy/scipy/ticket/231 and the short thread at scipydev http://thread.gmane.org/gmane.comp.python.scientific.devel/5166/focus=5175
In addition to the wrapper there is a Python interface (and some tests). I don't know if the interface is like "eigs"  I don't use Matlab.
If you're interested in making an interface that looks like MATLAB's, you can take a peek at:
It will give you a few eigenvalues and eigenvectors for the standard eigenproblem (Ax=lx) for any type of A (symmetric/nonsymmetric, real/complex, single/double, sparse/nonsparse).
The generalized and shifted modes are not implemented. I need to find some time (or some help) to get it finished.
Regards, Aric
 Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with preintegrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.asus.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642 _______________________________________________ Numpydiscussion mailing list Numpydiscussion@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/numpydiscussion
 Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with preintegrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.asus.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
Argh.
<snip>
In addition to the wrapper there is a Python interface (and some tests). I don't know if the interface is like "eigs"  I don't use Matlab.
http://www.mathworks.de/access/helpdesk/help/techdoc/ref/index.html?/access/ helpdesk/help/techdoc/ref/eigs.html should give you some idea of how MATLAB's eigs works. Cheers, Albert  Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with preintegrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.asus.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
On 10/22/06, Aric Hagberg <hagberg@lanl.gov> wrote:
On Sat, Oct 21, 2006 at 02:05:42PM 0700, Keith Goodman wrote:
Did you, or anybody else on the list, have any luck making a numpy version of eigs?
I made a start at an ARPACK wrapper, see http://projects.scipy.org/scipy/scipy/ticket/231 and the short thread at scipydev http://thread.gmane.org/gmane.comp.python.scientific.devel/5166/focus=5175
That looks very promising!  Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with preintegrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.asus.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
Don't know if this is any use, but to me (not knowing Fortran nearly so well as C++) this looks pretty useful: http://www.caam.rice.edu/software/ARPACK/arpack++.html http://www.ime.unicamp.br/~chico/arpack++/ It provides a nice highlevel interface on top of ARPACK. I could see it being useful on a number of levels: 1) actually use its wrappers directly instead of calling the f2py'ed fortran code  at least as a stop gap measure to cover the holes ("shifted modes" etc)  this could potentially even be faster since the reversecommunication interface requires some iteration, and this way the loops would run in compiled C++ instead of python. 2) use it's code just as a reference to help write more wrappers to the fortran code 3) if nothing else, it has pretty decent documentation. For instance it includes a nice table of the amount of storage you need to reserve for various different calling modes. (It seems ARPACK's main documentation is not freely available, so I'm not sure what's in it, but if you have that then ARPACK++'s docs may not be much help). As far as I can tell, it must be under the same licensing terms as ARPACK itself. It doesn't actually specify anything as far as I could see. But it seems to be linked prominently from the ARPACK website and appears as an offshoot project. bb On 10/22/06, Aric Hagberg <hagberg@lanl.gov> wrote:
On Sat, Oct 21, 2006 at 02:05:42PM 0700, Keith Goodman wrote:
Did you, or anybody else on the list, have any luck making a numpy version of eigs?
I made a start at an ARPACK wrapper, see http://projects.scipy.org/scipy/scipy/ticket/231 and the short thread at scipydev http://thread.gmane.org/gmane.comp.python.scientific.devel/5166/focus=5175
In addition to the wrapper there is a Python interface (and some tests). I don't know if the interface is like "eigs"  I don't use Matlab.
It will give you a few eigenvalues and eigenvectors for the standard eigenproblem (Ax=lx) for any type of A (symmetric/nonsymmetric, real/complex, single/double, sparse/nonsparse).
The generalized and shifted modes are not implemented. I need to find some time (or some help) to get it finished.
Regards, Aric
participants (8)

Albert Strasheim

Aric Hagberg

Bill Baxter

Charles R Harris

David Grant

Keith Goodman

Nils Wagner

Stefan van der Walt