[SciPy-dev] hyper1f1, ticket #659

[Moritz Helias]

Hello,

now you will also find a test for hyp1f1 attached to ticket #659:

1.) test_hyp1f1.nb
       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

2.) hyp1f1_test.py
       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  
relative error

3.) scipy_4393_hyp1f1_chgm.patch
      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  
4.41274416139e-05


I hope this will facilitate testing.

Greetings,

Moritz



On May 29, 2008, at 10:24 AM, Moritz Helias wrote:

> Hello,
>
> 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?
>
> Greetings,
>
> Moritz
>
> On May 28, 2008, at 10:16 PM, Stéfan van der Walt wrote:
>
>> 2008/5/28 Moritz Helias <helias@bccn.uni-freiburg.de>:
>>> special/specfun/specfun.f contains an implementation of M(a,b,x)
>>> which
>>> is called chgm.
>>> This works fine for me. It just has to be wrapped (which I did on my
>>> machine).
>>
>> Would you mind attaching the patch to a ticket for review?  Please
>> also include some tests to verify its behaviour.
>>
>> Thanks
>> Stéfan
>> _______________________________________________
>> Scipy-dev mailing list
>> Scipy-dev@scipy.org
>> http://projects.scipy.org/mailman/listinfo/scipy-dev
>
> _______________________________________________
> Scipy-dev mailing list
> Scipy-dev@scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-dev

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://projects.scipy.org/pipermail/scipy-dev/attachments/20080529/1e6442ba/attachment.html