bigfloats of special functions (e.g. Bessel).. wasRe: integrate bessel_j errcatch?
Subject: bigfloats of special functions (e.g. Bessel).. wasRe: integrate bessel_j errcatch?
From: Richard Fateman
Date: Thu, 08 Dec 2011 09:30:11 -0800
On 11/9/2011 12:44 PM, Barton Willis wrote:
> In addition to being slow, only the 2F1 hypergeometric function analytically continues to outside the convergence disk.
>
> --Barton
>
I have been looking around at the literature on "unrestricted
algorithms" for elementary and special functions,
esp. work by Richard Brent.
These are algorithms appropriate for bigfloats in that instead of
compute sin(x), or bessel_j(0,x),
for a given x, take another parameter n, where n is the number of bits
of precision required.
Given this amount of (fairly readable) stuff including code, it seems
some student project
could be handed this task of upgrading numerical evaluation to work for
bigfloats.
I note that bessel_j(0,3.12b0) remains unevaluated, right now.
Any takers?
RJF