\tan \left(x + \varepsilon\right) - \tan x
\left(\left(\frac{\sin \varepsilon}{\left(1 - \frac{{\left(\sin x\right)}^{3} \cdot {\left(\sin \varepsilon\right)}^{3}}{{\left(\cos x\right)}^{3} \cdot {\left(\cos \varepsilon\right)}^{3}}\right) \cdot \cos \varepsilon} + \mathsf{fma}\left(\frac{{\left(\sin x\right)}^{2}}{{\left(\cos x\right)}^{2}}, \frac{{\left(\sin \varepsilon\right)}^{3}}{\left(1 - \frac{{\left(\sin x\right)}^{3} \cdot {\left(\sin \varepsilon\right)}^{3}}{{\left(\cos x\right)}^{3} \cdot {\left(\cos \varepsilon\right)}^{3}}\right) \cdot {\left(\cos \varepsilon\right)}^{3}}, \frac{{\left(\sin x\right)}^{2} \cdot \sin \varepsilon}{\cos \varepsilon \cdot \left(\left(1 - \frac{{\left(\sin x\right)}^{3} \cdot {\left(\sin \varepsilon\right)}^{3}}{{\left(\cos x\right)}^{3} \cdot {\left(\cos \varepsilon\right)}^{3}}\right) \cdot {\left(\cos x\right)}^{2}\right)}\right)\right) + \frac{{\left(\sin \varepsilon\right)}^{2}}{\left(1 - \frac{{\left(\sin x\right)}^{3} \cdot {\left(\sin \varepsilon\right)}^{3}}{{\left(\cos x\right)}^{3} \cdot {\left(\cos \varepsilon\right)}^{3}}\right) \cdot {\left(\cos \varepsilon\right)}^{2}} \cdot \left(\frac{\sin x}{\cos x} + \frac{{\left(\sin x\right)}^{3}}{{\left(\cos x\right)}^{3}}\right)\right) + \log \left(e^{\frac{\sin x}{\left(1 - \frac{{\left(\sin x\right)}^{3} \cdot {\left(\sin \varepsilon\right)}^{3}}{{\left(\cos x\right)}^{3} \cdot {\left(\cos \varepsilon\right)}^{3}}\right) \cdot \cos x} - \frac{\sin x}{\cos x}}\right)double code(double x, double eps) {
return ((double) (((double) tan(((double) (x + eps)))) - ((double) tan(x))));
}
double code(double x, double eps) {
return ((double) (((double) (((double) (((double) (((double) sin(eps)) / ((double) (((double) (1.0 - ((double) (((double) (((double) pow(((double) sin(x)), 3.0)) * ((double) pow(((double) sin(eps)), 3.0)))) / ((double) (((double) pow(((double) cos(x)), 3.0)) * ((double) pow(((double) cos(eps)), 3.0)))))))) * ((double) cos(eps)))))) + ((double) fma(((double) (((double) pow(((double) sin(x)), 2.0)) / ((double) pow(((double) cos(x)), 2.0)))), ((double) (((double) pow(((double) sin(eps)), 3.0)) / ((double) (((double) (1.0 - ((double) (((double) (((double) pow(((double) sin(x)), 3.0)) * ((double) pow(((double) sin(eps)), 3.0)))) / ((double) (((double) pow(((double) cos(x)), 3.0)) * ((double) pow(((double) cos(eps)), 3.0)))))))) * ((double) pow(((double) cos(eps)), 3.0)))))), ((double) (((double) (((double) pow(((double) sin(x)), 2.0)) * ((double) sin(eps)))) / ((double) (((double) cos(eps)) * ((double) (((double) (1.0 - ((double) (((double) (((double) pow(((double) sin(x)), 3.0)) * ((double) pow(((double) sin(eps)), 3.0)))) / ((double) (((double) pow(((double) cos(x)), 3.0)) * ((double) pow(((double) cos(eps)), 3.0)))))))) * ((double) pow(((double) cos(x)), 2.0)))))))))))) + ((double) (((double) (((double) pow(((double) sin(eps)), 2.0)) / ((double) (((double) (1.0 - ((double) (((double) (((double) pow(((double) sin(x)), 3.0)) * ((double) pow(((double) sin(eps)), 3.0)))) / ((double) (((double) pow(((double) cos(x)), 3.0)) * ((double) pow(((double) cos(eps)), 3.0)))))))) * ((double) pow(((double) cos(eps)), 2.0)))))) * ((double) (((double) (((double) sin(x)) / ((double) cos(x)))) + ((double) (((double) pow(((double) sin(x)), 3.0)) / ((double) pow(((double) cos(x)), 3.0)))))))))) + ((double) log(((double) exp(((double) (((double) (((double) sin(x)) / ((double) (((double) (1.0 - ((double) (((double) (((double) pow(((double) sin(x)), 3.0)) * ((double) pow(((double) sin(eps)), 3.0)))) / ((double) (((double) pow(((double) cos(x)), 3.0)) * ((double) pow(((double) cos(eps)), 3.0)))))))) * ((double) cos(x)))))) - ((double) (((double) sin(x)) / ((double) cos(x))))))))))));
}




Bits error versus x




Bits error versus eps
Results
| Original | 36.8 |
|---|---|
| Target | 14.9 |
| Herbie | 0.5 |
Initial program 36.8
rmApplied tan-sum21.8
Simplified21.8
rmApplied flip3--21.8
Applied associate-/r/21.8
Applied fma-neg21.8
Taylor expanded around inf 21.9
Simplified0.5
rmApplied add-log-exp0.5
Final simplification0.5
herbie shell --seed 2020114 +o rules:numerics
(FPCore (x eps)
:name "2tan (problem 3.3.2)"
:precision binary64
:herbie-target
(/ (sin eps) (* (cos x) (cos (+ x eps))))
(- (tan (+ x eps)) (tan x)))