SF [2159499] Full bigfloat precision for Gamma after the second call
Subject: SF [2159499] Full bigfloat precision for Gamma after the second call
From: Raymond Toy
Date: Sat, 18 Oct 2008 10:27:12 -0400
van Nek wrote:
>
> (defun fprt18231_3 nil
> (let ((a 1823176476672000))
> (setq a (ash a (ash fpprec 1)))
> (fpquotient (intofp (isqrt a)) (intofp (ash 1 fpprec))) ))
Can you try this version? I think it's correct, and gets rid of the
fpquotient. My tests with cmucl indicate that this is still quite a bit
slower than the original.
(defun fprt18231_3a nil
(let ((a 1823176476672000))
(setq a (ash a (ash fpprec 1)))
(destructuring-bind (mantissa exp)
(intofp (isqrt a))
(list mantissa (- exp fpprec)))))
Ray