\sin \left(x + \varepsilon\right) - \sin x
\left(\sin x \cdot \cos \varepsilon - \sin x\right) + \sin \varepsilon \cdot \cos x
double f(double x, double eps) {
double r13957823 = x;
double r13957824 = eps;
double r13957825 = r13957823 + r13957824;
double r13957826 = sin(r13957825);
double r13957827 = sin(r13957823);
double r13957828 = r13957826 - r13957827;
return r13957828;
}
double f(double x, double eps) {
double r13957829 = x;
double r13957830 = sin(r13957829);
double r13957831 = eps;
double r13957832 = cos(r13957831);
double r13957833 = r13957830 * r13957832;
double r13957834 = r13957833 - r13957830;
double r13957835 = sin(r13957831);
double r13957836 = cos(r13957829);
double r13957837 = r13957835 * r13957836;
double r13957838 = r13957834 + r13957837;
return r13957838;
}




Bits error versus x




Bits error versus eps
Results
| Original | 37.0 |
|---|---|
| Target | 15.2 |
| Herbie | 0.4 |
Initial program 37.0
rmApplied sin-sum21.7
Taylor expanded around -inf 21.7
Simplified0.4
rmApplied fma-udef0.4
Final simplification0.4
herbie shell --seed 2019107 +o rules:numerics
(FPCore (x eps)
:name "2sin (example 3.3)"
:herbie-target
(* 2 (* (cos (+ x (/ eps 2))) (sin (/ eps 2))))
(- (sin (+ x eps)) (sin x)))