\frac{1}{2} \cdot \log \left(\frac{1 + x}{1 - x}\right)\left(\mathsf{log1p}\left(x\right) - \mathsf{log1p}\left(\left(-x\right)\right)\right) \cdot \frac{1}{2}double f(double x) {
double r1343350 = 1.0;
double r1343351 = 2.0;
double r1343352 = r1343350 / r1343351;
double r1343353 = x;
double r1343354 = r1343350 + r1343353;
double r1343355 = r1343350 - r1343353;
double r1343356 = r1343354 / r1343355;
double r1343357 = log(r1343356);
double r1343358 = r1343352 * r1343357;
return r1343358;
}
double f(double x) {
double r1343359 = x;
double r1343360 = log1p(r1343359);
double r1343361 = -r1343359;
double r1343362 = log1p(r1343361);
double r1343363 = r1343360 - r1343362;
double r1343364 = 0.5;
double r1343365 = r1343363 * r1343364;
return r1343365;
}



Bits error versus x
Results
Initial program 58.5
Simplified58.5
rmApplied log-div58.5
Simplified50.5
rmApplied log1p-expm1-u50.5
Simplified0.0
Final simplification0.0
herbie shell --seed 2019129 +o rules:numerics
(FPCore (x)
:name "Hyperbolic arc-(co)tangent"
(* (/ 1 2) (log (/ (+ 1 x) (- 1 x)))))