[Numpy-discussion] Suggestions for GSoC Projects
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'))
> 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.  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.
>  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
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the NumPy-Discussion