On Wed, 28 Jun 2006 03:22:28 -0500 Robert Kern <[EMAIL PROTECTED]> wrote:
> [EMAIL PROTECTED] wrote: > > Hi, > > > > [TO]: NumPy uses Numeric's old wrapper to lapack algorithms. > > [TO]: > > [TO]: SciPy uses it's own f2py-generated wrapper (it doesn't rely on the > > [TO]: NumPy wrapper). > > [TO]: > > [TO]: The numpy.dual library exists so you can use the SciPy calls if > > the [TO]: person has SciPy installed or the NumPy ones otherwise. It > > exists [TO]: precisely for the purpose of seamlessly taking advantage of > > [TO]: algorithms/interfaces that exist in NumPy but are improved in > > SciPy. > > > > This strikes me as a little bit odd. Why not just provide the > > best-performing function to both SciPy and NumPy? Would NumPy be more > > difficult to install if the SciPy algorithm for inv() was incorporated? > > That's certainly the case for the FFT algorithms. Scipy wraps more (and > more complicated) FFT libraries that are faster than FFTPACK. > > Most of the linalg functionality should probably be wrapping the same > routines if an optimized LAPACK is available. However, changing the routine > used in numpy in the absence of an optimized LAPACK would require > reconstructing the f2c'ed lapack_lite library that we include with the > numpy source. That hasn't been touched in so long that I would hesitate to > do so. If you are willing to do the work and the testing to ensure that it > still works everywhere, we'd probably accept the change. Annoying to redo (as tracking down *good* LAPACK sources is a chore), but hardly as bad as it was. I added the scripts I used to generated lapack_lite.c et al to numpy/linalg/lapack_lite in svn. These are the same things that were used to generate those files in recent versions of Numeric (which numpy uses). You only need to specify the top-level routines; the scripts find the dependencies. I'd suggest using the source for LAPACK that Debian uses; the maintainer, Camm Maguire, has done a bunch of work adding patches to fix routines that have been floating around. For instance, eigenvalues works better than before (lot fewer segfaults). With this, the hard part is writing the wrapper routines. If someone wants to wrap extra routines, I can do the the lapack_lite generation for them. -- |>|\/|< /--------------------------------------------------------------------------\ |David M. Cooke http://arbutus.physics.mcmaster.ca/dmc/ |[EMAIL PROTECTED] Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with pre-integrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642 _______________________________________________ Numpy-discussion mailing list Numpy-discussion@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/numpy-discussion
