\tan \left(x + \varepsilon\right) - \tan x
\left(\frac{\frac{\sin x}{\cos x}}{1 - \frac{\sqrt[3]{\sin x} \cdot \sqrt[3]{\sin x}}{\cos x \cdot \cos \varepsilon} \cdot \frac{\sqrt[3]{\sin x}}{\frac{1}{\sin \varepsilon}}} - \frac{\sin x}{\cos x}\right) + \frac{\frac{\sin \varepsilon}{\cos \varepsilon}}{1 - \frac{\sin x}{\cos x} \cdot \frac{\sin \varepsilon}{\cos \varepsilon}}double f(double x, double eps) {
double r3622319 = x;
double r3622320 = eps;
double r3622321 = r3622319 + r3622320;
double r3622322 = tan(r3622321);
double r3622323 = tan(r3622319);
double r3622324 = r3622322 - r3622323;
return r3622324;
}
double f(double x, double eps) {
double r3622325 = x;
double r3622326 = sin(r3622325);
double r3622327 = cos(r3622325);
double r3622328 = r3622326 / r3622327;
double r3622329 = 1.0;
double r3622330 = cbrt(r3622326);
double r3622331 = r3622330 * r3622330;
double r3622332 = eps;
double r3622333 = cos(r3622332);
double r3622334 = r3622327 * r3622333;
double r3622335 = r3622331 / r3622334;
double r3622336 = sin(r3622332);
double r3622337 = r3622329 / r3622336;
double r3622338 = r3622330 / r3622337;
double r3622339 = r3622335 * r3622338;
double r3622340 = r3622329 - r3622339;
double r3622341 = r3622328 / r3622340;
double r3622342 = r3622341 - r3622328;
double r3622343 = r3622336 / r3622333;
double r3622344 = r3622328 * r3622343;
double r3622345 = r3622329 - r3622344;
double r3622346 = r3622343 / r3622345;
double r3622347 = r3622342 + r3622346;
return r3622347;
}




Bits error versus x




Bits error versus eps
Results
| Original | 37.0 |
|---|---|
| Target | 15.4 |
| Herbie | 12.6 |
Initial program 37.0
rmApplied tan-sum21.6
Taylor expanded around -inf 21.7
Simplified12.6
Taylor expanded around inf 12.6
Simplified12.6
rmApplied div-inv12.6
Applied add-cube-cbrt12.6
Applied times-frac12.6
Final simplification12.6
herbie shell --seed 2019133 +o rules:numerics
(FPCore (x eps)
:name "2tan (problem 3.3.2)"
:herbie-target
(/ (sin eps) (* (cos x) (cos (+ x eps))))
(- (tan (+ x eps)) (tan x)))