From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 4937 invoked by alias); 12 Dec 2007 16:28:59 -0000 Received: (qmail 4917 invoked by uid 22791); 12 Dec 2007 16:28:57 -0000 X-Spam-Check-By: sourceware.org Received: from network-theory.com (HELO mail.network-theory.co.uk) (66.199.228.187) by sourceware.org (qpsmtpd/0.31) with ESMTP; Wed, 12 Dec 2007 16:28:53 +0000 Date: Wed, 12 Dec 2007 16:28:00 -0000 Message-ID: <87fxy73mil.wl%bjg@network-theory.co.uk> From: Brian Gough To: Richard Mathar Cc: gsl-discuss@sourceware.org Subject: Re: Hypergeometric 2F1() specfunc: cases x=1 In-Reply-To: <200712112114.lBBLERNe005280@amer.strw.leidenuniv.nl> References: <200712112114.lBBLERNe005280@amer.strw.leidenuniv.nl> User-Agent: Wanderlust/2.14.0 (Africa) Emacs/22.1 Mule/5.0 (SAKAKI) MIME-Version: 1.0 (generated by SEMI 1.14.6 - "Maruoka") Content-Type: text/plain; charset=US-ASCII X-Message-Mac: 87016ad6893dd434d138828a9a8115ff Mailing-List: contact gsl-discuss-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: gsl-discuss-owner@sourceware.org X-SW-Source: 2007-q4/txt/msg00040.txt.bz2 At Tue, 11 Dec 2007 22:14:27 +0100, Richard Mathar wrote: > It would be nice if the implementation of the Gaussian hypergeometric > functions would catch the special case of unit argument and compute it > via the standard formula with the Gamma functions (Abramowitz & Stegun 15.1.20). Sounds reasonable --- can you write some test cases with high precision values for specfunc/test_hyperg.c to make sure it's working. -- Brian Gough