\sin \left(x + \varepsilon\right) - \sin x
\frac{\sin x \cdot \left(-\sin \varepsilon \cdot \sin \varepsilon\right)}{\cos \varepsilon + 1} + \cos x \cdot \sin \varepsilondouble f(double x, double eps) {
double r177736 = x;
double r177737 = eps;
double r177738 = r177736 + r177737;
double r177739 = sin(r177738);
double r177740 = sin(r177736);
double r177741 = r177739 - r177740;
return r177741;
}
double f(double x, double eps) {
double r177742 = x;
double r177743 = sin(r177742);
double r177744 = eps;
double r177745 = sin(r177744);
double r177746 = r177745 * r177745;
double r177747 = -r177746;
double r177748 = r177743 * r177747;
double r177749 = cos(r177744);
double r177750 = 1.0;
double r177751 = r177749 + r177750;
double r177752 = r177748 / r177751;
double r177753 = cos(r177742);
double r177754 = r177753 * r177745;
double r177755 = r177752 + r177754;
return r177755;
}




Bits error versus x




Bits error versus eps
Results
| Original | 37.1 |
|---|---|
| Target | 14.9 |
| Herbie | 0.4 |
Initial program 37.1
rmApplied sin-sum22.0
Taylor expanded around inf 22.0
Simplified0.4
rmApplied flip--0.5
Applied associate-*r/0.5
Simplified0.5
rmApplied sub-1-cos0.4
Final simplification0.4
herbie shell --seed 2019347
(FPCore (x eps)
:name "2sin (example 3.3)"
:precision binary64
:herbie-target
(* 2 (* (cos (+ x (/ eps 2))) (sin (/ eps 2))))
(- (sin (+ x eps)) (sin x)))