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

SciPy scipy-tickets at scipy.net
Fri Nov 17 15:06:01 CST 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 NeilenMarais):

 Replying to [comment:6 aric]:
 > I did some further looking into this.  I do believe I have a working
 interface now.

 > I made some tweaks to the pyf file and also discovered that one of the
 LAPACK routines
 > was the culprit in producing the strange results for <s,d>neupd.  If you
 compile in your
 > own version of <s,d>lahqr.f with suitable compiler options everything
 seems to work.
 >
 > e.g. -ffloat-store or
 > my default with setup.py is

 > /usr/bin/g77 -g -Wall -fno-second-underscore -fPIC -O3 -funroll-loops
 -march=pentium4 -mmmx -msse2 -msse -fomit-frame-pointer -malign-double

 Interestingly, I no longer get segfaults even with the default compiler
 options. I wonder if it wasn't just some F2PY issue that's been sorted out
 in the meantime (I try to keep up to date with svn scipy)

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


More information about the Scipy-tickets mailing list