[SciPy-dev] Wrapping band matrices

[Arnd Baecker]

Hi Pearu,

On Thu, 1 Jun 2006, Pearu Peterson wrote:

>
> Hi Arnd,
> Could you set up a scipy ticket for your request and assign it to me?

ok, I will.

> At the moment I don't have time for testing your wrappers

That's fine - actually, before you really start looking into them,
it would be better that we do some of the things in the TODO file before
(<s,d> stuff and add some more unit tests)
- I will let you know when we managed to do that.

At this point it would be great if could give me a pointer why
I don't have the `overwrite_ab` option anymore  - does one have to
change the wrapper somehow?

Many thanks,

Arnd


> but I'm
> interested integrating them to scipy/lib/lapack (the plan is to get rid
> of lapack wrappers from linalg module as other scipy packages might use
> lapack wrappers as well).
> Thanks,
> Pearu
>
> On Wed, 31 May 2006, Arnd Baecker wrote:
>
> > Hi,
> >
> > I would like to contribute f2py generated
> > wrappers for the LAPACK routines
> >  dgbrfs,  dgbtrf, dgbtrs, dsbev, dsbevd, dsbevx
> > which deal with symmetric band matrices.
> >
> > Originally they were developed with Numeric
> > and the old f2py.
> >
> > At this point I have two questions:
> >
> > a) before the function signature looked like:
> >
> >    dsbev - Function signature:
> >      w,z,info = dsbev(ab,[compute_v,lower,ldab,overwrite_ab])
> >
> >   but now it is
> >
> >    dsbev - Function signature:
> >      w,z,info = dsbev(ab,[compute_v,lower,ldab])
> >
> >   I.e.: `overwrite_ab` is missing.
> >
> >   What has to be done to get this again?
> >
> > b) What would be the best approach for the integration
> >   into scipy?
> >
> >   The present status is at
> >      http://www.physik.tu-dresden.de/~baecker/python/band.zip
> >
> >   The routines work "stand-alone"
> >     python setup.py build
> >     python tests/test_band.py
> >   and  (after more clean-up/documentation)_
> >   might_ be a good example/framework
> >   for others to add further LAPACK routines.
> >   In particular I think it is easier to test new stuff separately
> >   than to directly integrate it into
> >   scipy/Lib/linalg/generic_flapack.pyf
> >
> >   Would it make sense to add this to the scipy sandbox at some point?
> >   (I would prefer to leave the final integration
> >    into scipy/Lib/linalg/* to a real expert...)
> >
> >   In any case, before this, the next steps are:
> >     - do the <tchar=s,d>sbev stuff for all routines
> >       in generic_band.pyf
> >     - address FIXMEs in generic_band.pyf
> >     - add unit tests for all routines
> >       (only one example is given at the moment)
> >
> > Many thanks,
> >
> > Arnd
> >
> > _______________________________________________
> > Scipy-dev mailing list
> > Scipy-dev at scipy.net
> > http://www.scipy.net/mailman/listinfo/scipy-dev
> >
>
> _______________________________________________
> Scipy-dev mailing list
> Scipy-dev at scipy.net
> http://www.scipy.net/mailman/listinfo/scipy-dev
>
>