\sin \left(x + \varepsilon\right) - \sin x
\left(\sin x \cdot \cos \varepsilon + \left(-\sin x\right)\right) + \cos x \cdot \sin \varepsilon
double f(double x, double eps) {
double r88355 = x;
double r88356 = eps;
double r88357 = r88355 + r88356;
double r88358 = sin(r88357);
double r88359 = sin(r88355);
double r88360 = r88358 - r88359;
return r88360;
}
double f(double x, double eps) {
double r88361 = x;
double r88362 = sin(r88361);
double r88363 = eps;
double r88364 = cos(r88363);
double r88365 = r88362 * r88364;
double r88366 = -r88362;
double r88367 = r88365 + r88366;
double r88368 = cos(r88361);
double r88369 = sin(r88363);
double r88370 = r88368 * r88369;
double r88371 = r88367 + r88370;
return r88371;
}




Bits error versus x




Bits error versus eps
Results
| Original | 37.5 |
|---|---|
| Target | 14.9 |
| Herbie | 0.5 |
Initial program 37.5
rmApplied sin-sum22.5
rmApplied *-un-lft-identity22.5
Applied *-un-lft-identity22.5
Applied distribute-lft-out--22.5
Simplified0.4
rmApplied fma-udef0.4
rmApplied sub-neg0.4
Applied distribute-lft-in0.5
Simplified0.5
Final simplification0.5
herbie shell --seed 2019353 +o rules:numerics
(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)))