\frac{1}{x} - \frac{1}{\tan x}\left(0.0222222222222222231 \cdot {x}^{3} + 0.00211640211640211654 \cdot {x}^{5}\right) + 0.333333333333333315 \cdot xdouble f(double x) {
double r128324 = 1.0;
double r128325 = x;
double r128326 = r128324 / r128325;
double r128327 = tan(r128325);
double r128328 = r128324 / r128327;
double r128329 = r128326 - r128328;
return r128329;
}
double f(double x) {
double r128330 = 0.022222222222222223;
double r128331 = x;
double r128332 = 3.0;
double r128333 = pow(r128331, r128332);
double r128334 = r128330 * r128333;
double r128335 = 0.0021164021164021165;
double r128336 = 5.0;
double r128337 = pow(r128331, r128336);
double r128338 = r128335 * r128337;
double r128339 = r128334 + r128338;
double r128340 = 0.3333333333333333;
double r128341 = r128340 * r128331;
double r128342 = r128339 + r128341;
return r128342;
}




Bits error versus x
Results
| Original | 59.9 |
|---|---|
| Target | 0.1 |
| Herbie | 0.3 |
Initial program 59.9
Taylor expanded around 0 0.3
rmApplied associate-+r+0.3
Final simplification0.3
herbie shell --seed 2020020
(FPCore (x)
:name "invcot (example 3.9)"
:precision binary64
:pre (and (< -0.026 x) (< x 0.026))
:herbie-target
(if (< (fabs x) 0.026) (* (/ x 3) (+ 1 (/ (* x x) 15))) (- (/ 1 x) (/ 1 (tan x))))
(- (/ 1 x) (/ 1 (tan x))))