[Numpy-discussion] Suggestions for GSoC Projects

Jennifer stone jenny.stone125 at gmail.com
Sun Feb 23 11:05:07 EST 2014


>
> Typically, you need to tell mpmath to use an appropriate precision for
> the evaluation:
>
> >>> mpmath.mp.dps <http://bl-1.com/click/load/AzZebQZnUGVfPQNgU2Q-b0231>= 300
> >>> float(mpmath.hyp2f1('10','5','-300.5','0.5'))
> -3.8520270815239185e+32
>
> Oh k! I didn't pay that much heed to it and set mp.dps<http://bl-1.com/click/load/UGVbaAdmUmdXNVU2U2c-b0231>to 100. Wonder why the
difference was so drastic. Anyways thanks a lot.

>
>     ***
>
> What can help in hyp2f1 for large values of a,b,c is use of recurrence
> relations. These are typically stable in one direction only. [1] This
> seems to be still a partially open research question...
>
> Our current hyp2f1 implementation does use recurrences (hyp2f1ra), but
> perhaps they are not invoked for this case. The problem here can be the
> accurate determination of the convergence region for each parameter value.
>
> [1] http://www.ams.org/journals/mcom/2007-76-259/S0025-5718-07-01918-7/<http://bl-1.com/click/load/BzIAM1c2UmcCYFU2U2Y-b0231>
>
> --
> Pauli Virtanen
>
>
> _______________________________________________
> NumPy-Discussion mailing list
> NumPy-Discussion at scipy.org
> http://mail.scipy.org/mailman/listinfo/numpy-discussion<http://bl-1.com/click/load/XWgMP1Y3BTAEZlQ3BDY-b0231>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/numpy-discussion/attachments/20140223/14ed5e40/attachment.html>


More information about the NumPy-Discussion mailing list