[Scipy-tickets] [SciPy] #231: A wrapper for ARPACK

SciPy scipy-tickets at scipy.net
Mon Oct 9 22:52:58 CDT 2006


#231: A wrapper for ARPACK
--------------------------+-------------------------------------------------
 Reporter:  nils          |        Owner:  somebody
     Type:  enhancement   |       Status:  new     
 Priority:  normal        |    Milestone:          
Component:  scipy.sparse  |      Version:          
 Severity:  normal        |   Resolution:          
 Keywords:                |  
--------------------------+-------------------------------------------------
Comment (by aric):

 I have attached a first cut at an ARPACK wrapper.
 It uses Travis' suggestion and follows a style similar
 to interative.py - all you need to provide is a numpy
 array or matrix or other object with a matvec method
 (like the sparse matrices).

 There are at least two glaring defects:

 1) Only "mode1" type eigenproblems are implemented (Ax=wx)
    Generalized eigenproblems (Ax=wMX) and various inverse and shifted
    modes would be a great addition.

 2) There is some issue with the single and double precision non-symmetric
    methods.  Either they give the wrong answer or I don't understand the
 output.
    The interface is high style 20th century FORTRAN and I might have made
 an error
    in workspace sizes or the like.  The symmetric and complex solvers seem
 to
    work (some tests are provided).

-- 
Ticket URL: <http://projects.scipy.org/scipy/scipy/ticket/231#comment:1>
SciPy <http://www.scipy.org/>
SciPy is open-source software for mathematics, science, and engineering.


More information about the Scipy-tickets mailing list