X-Recipient: archive-cygwin@delorie.com
X-Spam-Check-By: sourceware.org
Message-ID: <ba40711f0710120545m7ea72cc5oa0604a4a5ea25072@mail.gmail.com>
Date: Fri, 12 Oct 2007 08:45:45 -0400
From: "Lev Bishop" <lev.bishop@gmail.com>
To: cygwin@cygwin.com
Subject: Re: A problem with tgamma function
In-Reply-To: <Pine.OSF.4.58.0710120952130.4055@ax0rm1.roma1.infn.it>
MIME-Version: 1.0
Content-Type: text/plain; charset=ISO-8859-1
Content-Transfer-Encoding: 7bit
Content-Disposition: inline
References: <Pine.OSF.4.58.0710120952130.4055@ax0rm1.roma1.infn.it>
Mailing-List: contact cygwin-help@cygwin.com; run by ezmlm
Precedence: bulk
List-Id: <cygwin.cygwin.com>
List-Unsubscribe: <mailto:cygwin-unsubscribe-archive-cygwin=delorie.com@cygwin.com>
List-Subscribe: <mailto:cygwin-subscribe@cygwin.com>
List-Archive: <http://sourceware.org/ml/cygwin/>
List-Post: <mailto:cygwin@cygwin.com>
List-Help: <mailto:cygwin-help@cygwin.com>, <http://sourceware.org/ml/#faqs>
Sender: cygwin-owner@cygwin.com
Mail-Followup-To: cygwin@cygwin.com
Delivered-To: mailing list cygwin@cygwin.com

On 10/12/07, Angelo Graziosi wrote:
>
> Lev Bishop wrote:
>
> > I'm sure glibc and newlib would both appreciate a good algorithm for
> > tgamma(), if you felt like submitting one...
>
> It seems that a good algorithm has yet been coded in the GAMMA
> implementation of CERNLIB
> (http://wwwasdoc.web.cern.ch/wwwasdoc/shortwrupsdir/c302/top.html,
> http://cernlib.web.cern.ch/cernlib/download/2006_source/src/mathlib/gen/c/gamma64.F),
> because with that implementation the test I posted passes...

One immediately obvious problem with that implementation is that it
doesn't handle negative parameter.

Lev

--
Unsubscribe info:      http://cygwin.com/ml/#unsubscribe-simple
Problem reports:       http://cygwin.com/problems.html
Documentation:         http://cygwin.com/docs.html
FAQ:                   http://cygwin.com/faq/