double code(double x) {
return ((double) (((double) (1.0 / 2.0)) * ((double) log(((double) (((double) (1.0 + x)) / ((double) (1.0 - x))))))));
}
double code(double x) {
return ((double) (((double) (1.0 / 2.0)) * ((double) (((double) (0.6666666666666666 * ((double) pow(((double) (x / 1.0)), 3.0)))) + ((double) (((double) (2.0 * x)) + ((double) log(((double) pow(((double) exp(0.4)), ((double) (((double) pow(x, 5.0)) / ((double) pow(1.0, 5.0))))))))))))));
}



Bits error versus x
Results
Initial program 58.5
rmApplied log-div58.5
Taylor expanded around 0 0.3
Simplified0.3
rmApplied add-log-exp0.3
Simplified0.3
Final simplification0.3
herbie shell --seed 2020191
(FPCore (x)
:name "Hyperbolic arc-(co)tangent"
:precision binary64
(* (/ 1.0 2.0) (log (/ (+ 1.0 x) (- 1.0 x)))))