\frac{1}{x} - \frac{1}{\tan x}{x}^{5} \cdot \frac{2}{945} + \frac{x}{\frac{\frac{1}{9} + \left(\left(\frac{1}{45} \cdot x\right) \cdot x - \frac{1}{3}\right) \cdot \left(\left(\frac{1}{45} \cdot x\right) \cdot x\right)}{\frac{1}{91125} \cdot \left(\left(x \cdot \left(x \cdot x\right)\right) \cdot \left(x \cdot \left(x \cdot x\right)\right)\right) + \frac{1}{27}}}double f(double x) {
double r2141989 = 1.0;
double r2141990 = x;
double r2141991 = r2141989 / r2141990;
double r2141992 = tan(r2141990);
double r2141993 = r2141989 / r2141992;
double r2141994 = r2141991 - r2141993;
return r2141994;
}
double f(double x) {
double r2141995 = x;
double r2141996 = 5.0;
double r2141997 = pow(r2141995, r2141996);
double r2141998 = 0.0021164021164021165;
double r2141999 = r2141997 * r2141998;
double r2142000 = 0.1111111111111111;
double r2142001 = 0.022222222222222223;
double r2142002 = r2142001 * r2141995;
double r2142003 = r2142002 * r2141995;
double r2142004 = 0.3333333333333333;
double r2142005 = r2142003 - r2142004;
double r2142006 = r2142005 * r2142003;
double r2142007 = r2142000 + r2142006;
double r2142008 = 1.0973936899862826e-05;
double r2142009 = r2141995 * r2141995;
double r2142010 = r2141995 * r2142009;
double r2142011 = r2142010 * r2142010;
double r2142012 = r2142008 * r2142011;
double r2142013 = 0.037037037037037035;
double r2142014 = r2142012 + r2142013;
double r2142015 = r2142007 / r2142014;
double r2142016 = r2141995 / r2142015;
double r2142017 = r2141999 + r2142016;
return r2142017;
}




Bits error versus x
Results
| Original | 59.8 |
|---|---|
| Target | 0.1 |
| Herbie | 0.0 |
Initial program 59.8
Taylor expanded around 0 0.3
Simplified0.3
rmApplied flip3-+1.2
Applied associate-*l/1.1
Simplified0.3
rmApplied associate-/l*0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2019133
(FPCore (x)
:name "invcot (example 3.9)"
: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))))