Fastest code for Gamma function

Article View: comp.lang.fortran

Article #6560

#6560

From: Leonard Howell
Date: Sat, 06 May 2000 00:00

10 lines
1352 bytes

Hi all --

I need to compute the gamma probability density function  g(x;a,b)=C(a,b)
x**(a-1) exp(-x/b) on the order of 10**10 calls (on many different
occasions) and thus would like to use the FASTEST algorithm available.  It
must also be double precision with 14 or 15 decimal places of accuracy.
Advice/coded FORTRAN solutions greatly appreciated.  Other language
solutions also welcome, if the speed and accuracy is achieved.

Thanks, Leonard

Technical Details

Message-ID: <fGNQ4.51301$WF.2761743@bgtnsc04-news.ops.worldnet.att.net>

Path: arc1-ca.pugleaf.net!unknown!not-for-mail

Show Full Headers

Complete Article Headers

X-Google-Language: ENGLISH,ASCII-7-bit
X-Google-Thread: 10d632,8bad7910b3fe0dc6,start
X-Google-Attributes: gid10d632,public
X-Google-Thread: 104500,8bad7910b3fe0dc6,start
X-Google-Attributes: gid104500,public
X-Google-Thread: 1094ba,8bad7910b3fe0dc6,start
X-Google-Attributes: gid1094ba,public
From: "Leonard Howell" <leowood@worldnet.att.net>
Subject: Fastest code for Gamma function
Date: 2000/05/06
Message-ID: <fGNQ4.51301$WF.2761743@bgtnsc04-news.ops.worldnet.att.net>
X-Deja-AN: 619864555
X-Priority: 3
X-MimeOLE: Produced By Microsoft MimeOLE V5.00.2314.1300
X-Complaints-To: abuse@worldnet.att.net
X-Trace: bgtnsc04-news.ops.worldnet.att.net 957589259 63.10.52.223 (Sat, 06 May 2000 05:00:59 GMT)
Organization: AT&T Worldnet
X-MSMail-Priority: Normal
NNTP-Posting-Date: Sat, 06 May 2000 05:00:59 GMT
Newsgroups: comp.lang.fortran,computer42.statistik,sci.stat.math