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



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