X-Recipient: archive-cygwin AT delorie DOT com X-Spam-Check-By: sourceware.org Message-ID: Date: Fri, 12 Oct 2007 08:45:45 -0400 From: "Lev Bishop" To: cygwin AT cygwin DOT com Subject: Re: A problem with tgamma function In-Reply-To: MIME-Version: 1.0 Content-Type: text/plain; charset=ISO-8859-1 Content-Transfer-Encoding: 7bit Content-Disposition: inline References: Mailing-List: contact cygwin-help AT cygwin DOT com; run by ezmlm Precedence: bulk List-Id: List-Unsubscribe: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: cygwin-owner AT cygwin DOT com Mail-Followup-To: cygwin AT cygwin DOT com Delivered-To: mailing list cygwin AT cygwin DOT 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/