\frac{1}{2} \cdot \log \left(\frac{1 + x}{1 - x}\right)0.5 \cdot \left(\log 1 + \left(x + x\right)\right)
(FPCore (x) :precision binary64 (* (/ 1.0 2.0) (log (/ (+ 1.0 x) (- 1.0 x)))))
(FPCore (x) :precision binary64 (* 0.5 (+ (log 1.0) (+ x x))))
double code(double x) {
return (1.0 / 2.0) * log((1.0 + x) / (1.0 - x));
}
double code(double x) {
return 0.5 * (log(1.0) + (x + x));
}



Bits error versus x
Results
Initial program 58.6
Simplified58.6
Taylor expanded around 0 0.6
Simplified0.6
Final simplification0.6
herbie shell --seed 2020277
(FPCore (x)
:name "Hyperbolic arc-(co)tangent"
:precision binary64
(* (/ 1.0 2.0) (log (/ (+ 1.0 x) (- 1.0 x)))))