Re: R-alpha: R-0.15alpha gamma function bug

Rob Creecy (rcreecy@census.gov)
Tue, 14 Jan 1997 18:25:44 -0500


Message-Id: <32DC15F8.2B83@census.gov>
Date: Tue, 14 Jan 1997 18:25:44 -0500
From: Rob Creecy <rcreecy@census.gov>
To: r-testers@stat.math.ethz.ch
Subject: Re: R-alpha: R-0.15alpha gamma function bug

I have since discovered the bug fix I suggested is not
a fix after all, and I am not now sure why I have this bug
at all. I'm still trying to figure it out.

Rob Creecy

Gregory R. Warnes wrote:
> 
> No Problem here:
> 
> Solaris-gorn% uname -a
> SunOS gorn 5.5.1 Generic sun4u sparc SUNW,Ultra-1
> 
> Solaris-gorn% R
> 
> R : Copyright 1996, Robert Gentleman and Ross Ihaka
> Version 0.15  (Released: December 25, 1996)
> 
> R is free software and comes with ABSOLUTELY NO WARRANTY.
> You are welcome to redistribute it under certain conditions.
> Type "license()" for details.
> 
> > gamma(5)
> [1] 24
> 
> On Tue, 14 Jan 1997, Rob Creecy wrote:
> 
> > Has anyone else seen this?
> >
> > Version R-0.15alpha, compiled using cc on Solaris 2.3
> >
> > When first calling
> > >gamma(5)
> > [1] 0
> > >gamma(5)
> > [1] 24
> >
> > Apparently the variable
> >
> >  int signgam = 0; line 22 in lgamma.c
> > should be
> >  int signgam = 1;
> >
> > This seems to fix the problem.
> >
> > Rob Creecy
> > Census Bureau

> -------------------------------------------------------------------------------
>     Gregory R. Warnes          | It is high time that the ideal of success
> warnes@biostat.washington.edu  |  be replaced by the ideal of service.
>                                |                       Albert Einstein
> -------------------------------------------------------------------------------
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
r-testers mailing list -- For info or help, send "info" or "help",
To [un]subscribe, send "[un]subscribe"
(in the "body", not the subject !)  To: r-testers-request@stat.math.ethz.ch
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-