[SciPy-dev] hyper1f1, ticket #659
Thu May 29 08:45:15 CDT 2008
now you will also find a test for hyp1f1 attached to ticket #659:
Mathematica script to create reference data at randomly drawn
points (a,b,x) in [-30,30] x [-30,30] x [-30,30] using Mathematica's
implementation of hyp1f1
uses reference data created from test_hyp1f1.nb and compares it
to the result of scipy's implementation.
calculates root mean square error relative error and maximum
is a patch (build based on revision 4393) to be applied to scipy
root directory to use chgm (from specfun.f) instead of hyperg (from
hyperg.c) for hyp1f1.
this patch applied, the maximum relative error on my machine is
I hope this will facilitate testing.
On May 29, 2008, at 10:24 AM, Moritz Helias wrote:
> I just attached a patch to ticket #659 which replaces the
> implementation of hyp1f1 by the fortran routine
> chgm (scipy/special/specfun/specfun.f).
> The file hyp1f1_ticket.py (found in the same ticket) may serve as a
> first test. I'll try to write a test, which compares the function to
> matlab's implementation.
> Maybe somebody has already done so for the old implementation?
> On May 28, 2008, at 10:16 PM, Stéfan van der Walt wrote:
>> 2008/5/28 Moritz Helias <email@example.com>:
>>> special/specfun/specfun.f contains an implementation of M(a,b,x)
>>> is called chgm.
>>> This works fine for me. It just has to be wrapped (which I did on my
>> Would you mind attaching the patch to a ticket for review? Please
>> also include some tests to verify its behaviour.
>> Scipy-dev mailing list
> Scipy-dev mailing list
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Scipy-dev