[SciPy-Dev] Reimplementing hypergeometric functions

[Ted Pudlik]

Hello,

The current implementation of the confluent hypergeometric function in
SciPy has room for improvement---see gh-5349
<https://github.com/scipy/scipy/issues/5349> and the issues referenced
there.  But as has been discussed before
<https://github.com/scipy/scipy/issues/1561#issuecomment-130488352>, a recent
paper <http://arxiv.org/abs/1407.7786> reviews different algorithms for
this problem, and its authors even provide implementations
<http://datashare.is.ed.ac.uk/handle/10283/607> (albeit in MATLAB and with
an unclear license).

I'm planning to get in touch with the authors and see if they would be
willing to provide a copy of their MATLAB code with a BSD license.  If so,
I'll try to write a Cython translation and attempt to stitch the algorithms
together into a general-purpose routine to replace hyp1f1.

Does this sound reasonable?  Is anyone else working on reimplementing
hyp1f1 at the moment?

Best wishes,
Ted
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/scipy-dev/attachments/20160107/b3b70d9c/attachment.html>