|
From: | Raymond Toy |
Subject: | [Gcl-devel] Re: acosh issues |
Date: | Wed, 16 Oct 2002 23:28:59 -0400 |
User-agent: | Mozilla/5.0 (X11; U; Linux i686; en-US; rv:1.2b) Gecko/20021008 |
Camm Maguire wrote:
Greetings, and thanks for these pointers! There are probably other functions that need attention as well. In any case, we should be encouraged that the original version basically had the scaling issuedown.Here is what I've just committed: (defun acosh (x) (* 2 (log (+ (sqrt (/ (1+ x) 2)) (sqrt (/ (1- x) 2))))))
Although this doesn't help in the case where x is real, but you might be interested in looking at the routines in CMUCL (see src/irrat.lisp) for these functions for complex args. They're translations of the routines in Kahan's article that Barton references.
Ray
[Prev in Thread] | Current Thread | [Next in Thread] |