Problems with the Gamma function

• Subject: Problems with the Gamma function
• From: Robert Dodier
• Date: Sat, 13 Sep 2008 12:55:53 -0600

On 9/13/08, Dieter Kaiser <drdieterkaiser at web.de> wrote:

>                  (ignore-errors
>                    (* (sqrt (float (* 2 pi)))
>                       (+ ss (aref c 0))
>                       (exp (+ (- zgh) (* 2 lnzp)))))))
>             (cond ((null result)
>                    ;; No result. Overflow.
>                    (merror "Overflow in `gamma-lanczos'."))
>                   ((or (float-nan-p (realpart result))
>                        (float-inf-p (realpart result)))
>                    ;; Result, but beyond extreme values. Overflow.
>                    (merror "Overflow in `gamma-lanczos'."))
>                   (t result)))))))

Seems OK to me. I am assuming that REALPART is just
fetching something i.e. there is no computation which could
cause an error even if there were no error before that point.

Thanks for working on this,

Robert Dodier

• Prev by Date: segmentation fault with Maxima/ECL
• Next by Date: segmentation fault with Maxima/ECL
• Previous by thread: Problems with the Gamma function
• Next by thread: Problems with the Gamma function
• Index(es):
  • Date
  • Thread