\frac{\log \left(1 - x\right)}{\log \left(1 + x\right)}\frac{1}{\frac{1 \cdot x + \log 1}{\log 1 - \left(1 \cdot x + \frac{1}{2} \cdot \frac{{x}^{2}}{{1}^{2}}\right)} - \frac{\frac{1}{2} \cdot \frac{{x}^{2}}{{1}^{2}}}{\log 1 - \left(1 \cdot x + \frac{1}{2} \cdot \frac{{x}^{2}}{{1}^{2}}\right)}}double f(double x) {
double r102453 = 1.0;
double r102454 = x;
double r102455 = r102453 - r102454;
double r102456 = log(r102455);
double r102457 = r102453 + r102454;
double r102458 = log(r102457);
double r102459 = r102456 / r102458;
return r102459;
}
double f(double x) {
double r102460 = 1.0;
double r102461 = 1.0;
double r102462 = x;
double r102463 = r102461 * r102462;
double r102464 = log(r102461);
double r102465 = r102463 + r102464;
double r102466 = 0.5;
double r102467 = 2.0;
double r102468 = pow(r102462, r102467);
double r102469 = pow(r102461, r102467);
double r102470 = r102468 / r102469;
double r102471 = r102466 * r102470;
double r102472 = r102463 + r102471;
double r102473 = r102464 - r102472;
double r102474 = r102465 / r102473;
double r102475 = r102471 / r102473;
double r102476 = r102474 - r102475;
double r102477 = r102460 / r102476;
return r102477;
}




Bits error versus x
Results
| Original | 61.4 |
|---|---|
| Target | 0.3 |
| Herbie | 0.4 |
Initial program 61.4
Taylor expanded around 0 60.5
Taylor expanded around 0 0.4
rmApplied clear-num0.4
rmApplied div-sub0.4
Final simplification0.4
herbie shell --seed 2020035
(FPCore (x)
:name "qlog (example 3.10)"
:precision binary64
:pre (and (< -1 x) (< x 1))
:herbie-target
(- (+ (+ (+ 1 x) (/ (* x x) 2)) (* 0.4166666666666667 (pow x 3))))
(/ (log (- 1 x)) (log (+ 1 x))))