\sin \left(x + \varepsilon\right) - \sin x
\sin \varepsilon \cdot \left(\cos x - \frac{\sin x}{\cos \varepsilon + 1} \cdot \sin \varepsilon\right)double f(double x, double eps) {
double r116364 = x;
double r116365 = eps;
double r116366 = r116364 + r116365;
double r116367 = sin(r116366);
double r116368 = sin(r116364);
double r116369 = r116367 - r116368;
return r116369;
}
double f(double x, double eps) {
double r116370 = eps;
double r116371 = sin(r116370);
double r116372 = x;
double r116373 = cos(r116372);
double r116374 = sin(r116372);
double r116375 = cos(r116370);
double r116376 = 1.0;
double r116377 = r116375 + r116376;
double r116378 = r116374 / r116377;
double r116379 = r116378 * r116371;
double r116380 = r116373 - r116379;
double r116381 = r116371 * r116380;
return r116381;
}




Bits error versus x




Bits error versus eps
Results
| Original | 37.5 |
|---|---|
| Target | 15.3 |
| Herbie | 0.4 |
Initial program 37.5
rmApplied sin-sum22.1
rmApplied *-un-lft-identity22.1
Applied *-un-lft-identity22.1
Applied distribute-lft-out--22.1
Simplified0.4
rmApplied flip--0.5
Applied associate-*r/0.5
Simplified0.5
Final simplification0.4
herbie shell --seed 2019294
(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)))