From owner-freebsd-numerics@freebsd.org Wed Sep 5 19:02:25 2018 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2610:1c1:1:606c::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id 7AD52FFA31A for ; Wed, 5 Sep 2018 19:02:25 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (troutmask.apl.washington.edu [128.95.76.21]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (Client CN "troutmask", Issuer "troutmask" (not verified)) by mx1.freebsd.org (Postfix) with ESMTPS id 11FF88EBA2 for ; Wed, 5 Sep 2018 19:02:25 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (localhost [127.0.0.1]) by troutmask.apl.washington.edu (8.15.2/8.15.2) with ESMTPS id w85J2N2b027979 (version=TLSv1.2 cipher=DHE-RSA-AES256-GCM-SHA384 bits=256 verify=NO); Wed, 5 Sep 2018 12:02:23 -0700 (PDT) (envelope-from sgk@troutmask.apl.washington.edu) Received: (from sgk@localhost) by troutmask.apl.washington.edu (8.15.2/8.15.2/Submit) id w85J2Njc027978; Wed, 5 Sep 2018 12:02:23 -0700 (PDT) (envelope-from sgk) Date: Wed, 5 Sep 2018 12:02:23 -0700 From: Steve Kargl To: Bruce Evans Cc: freebsd-numerics@freebsd.org Subject: Re: j0 (and y0) in the range 2 <= x < (p/2)*log(2) Message-ID: <20180905190223.GA27865@troutmask.apl.washington.edu> Reply-To: sgk@troutmask.apl.washington.edu References: <20180903235724.GA95333@troutmask.apl.washington.edu> <20180905201540.D1142@besplex.bde.org> <20180905152104.GA26453@troutmask.apl.washington.edu> <20180906034525.A2959@besplex.bde.org> MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Disposition: inline In-Reply-To: <20180906034525.A2959@besplex.bde.org> User-Agent: Mutt/1.10.1 (2018-07-13) X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.27 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Wed, 05 Sep 2018 19:02:25 -0000 On Thu, Sep 06, 2018 at 04:09:05AM +1000, Bruce Evans wrote: > On Wed, 5 Sep 2018, Steve Kargl wrote: > > > On Wed, Sep 05, 2018 at 10:06:29PM +1000, Bruce Evans wrote: > >> On Mon, 3 Sep 2018, Steve Kargl wrote: > >> > >>> Anyone know where the approximations for j0 (and y0) come from? > >> > >> I think they are ordinary minimax rational approximations for related > >> functions. As you noticed, the asymptotic expansion doesn't work below > >> about x = 8 (it is off by about 10% for j0(2). But we want to use the > >> single formula given by the asymptotic expansion for all the subintervals: > > > > I've scoured the literature and web for methods of computing > > Bessel functions. These functions are important to my real > > work. I have not found any paper, webpage, documentation, etc. > > that describes what "the related functions" are. > > They are just the functions in the asymptotic expansion with errors corrected > as I discussed. And as I noted, there is no documentation stating the approximations pzero(x) and qzero(x) aren't approximations for the asymptotic series P0(x) and Q0(x). If you are correct, then pzero(x) and qzero(x) are approximations to fudge*P0(x) and fudge*Q0(x). What fudge is and how it is determined is not documented. -- Steve