Changeset 14665


Ignore:
Timestamp:
Mar 7, 2011, 6:53:21 PM (9 years ago)
Author:
rme
Message:

In ATANH, in the complex argument case, use the formula from the spec
exactly. We were using the identity log(x/y) = log(x) - log(y), but
this appears to alter the branch cut.

Fixes ticket:829.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/source/lib/numbers.lisp

    r14451 r14665  
    745745      #+64-bit-target
    746746      (%single-float-atanh x))
    747     (/ (log (/ (+ 1 x) (- 1 x))) 2)))
    748 
     747    (/ (- (log (+ 1 x)) (log (- 1 x))) 2)))
    749748
    750749(defun ffloor (number &optional divisor)
Note: See TracChangeset for help on using the changeset viewer.