[SciPy-User] generic_flapack.pyf and geqp3

[Skipper Seabold]

On Mon, Jul 11, 2011 at 8:44 PM, Collin Stocks <collinstocks at gmail.com> wrote:
> Hi, all.
>
> I am planning to try to add functionality to scipy.linalg.qr(),
> specifically to allow qr decomposition with pivoting. However, I have
> almost no knowledge of how to wrap the function in scipy/linalg/
> generic_flapack.pyf.
>
> Could somebody please point me in the correct direction?

Far from an expert, but I've used the 'smart way' in f2py to wrap some
LAPACK stuff. Basically, you run f2py on the fortran source, then edit
the .pyf file as you need to (use the f2py docs, the routine you're
wrapping's documentation, and the pyf examples in scipy for guidance),
and then run f2py again to build it (with -llapack or whatever if you
need to link against other libraries).

http://cens.ioc.ee/projects/f2py2e/usersguide/index.html#the-smart-way

I'd be interested to see an example of how to accomplish something
similar with fwrap.

Skipper