David Ronis wrote:
> Just updated from cvs/trunk on a i686-linux box, running clisp-2.47. I
> had the following testsuite failure:
>
>
> Running tests in rtest_gamma:
> ********************** Problem 134 ***************
> Input:
> closeto(double_factorial(10.0), 3840.0, 1.0E-12)
>
>
> Result:
> 1.364242052659392E-12
>
> This differed from the expected result:
> true
Can you compute the following values for me?
:lisp (expt 2d0 (/ #c(10d0 0d0) 2))
:lisp (gamma-lanczos (+ 1 (/ #c(10d0 0d0) 2)))
:lisp (expt (/ 2 (float pi)) (/ (- 1 (cos (* (float pi) #c(10d0 0d0)))) 4))
This is what I get with clisp:
#C(32.0 0.0)
#C(119.99999999999996 0.0)
#C(1.0 0.0)
So double_factorial should be the product of these, which is
#C(3839.9999999999986 0.0)
Ray