\sin \left(x + \varepsilon\right) - \sin x
\mathsf{fma}\left(\sin \varepsilon, \cos x, \mathsf{fma}\left(\sin x, \cos \varepsilon, -\sin x\right)\right)double f(double x, double eps) {
double r4548898 = x;
double r4548899 = eps;
double r4548900 = r4548898 + r4548899;
double r4548901 = sin(r4548900);
double r4548902 = sin(r4548898);
double r4548903 = r4548901 - r4548902;
return r4548903;
}
double f(double x, double eps) {
double r4548904 = eps;
double r4548905 = sin(r4548904);
double r4548906 = x;
double r4548907 = cos(r4548906);
double r4548908 = sin(r4548906);
double r4548909 = cos(r4548904);
double r4548910 = -r4548908;
double r4548911 = fma(r4548908, r4548909, r4548910);
double r4548912 = fma(r4548905, r4548907, r4548911);
return r4548912;
}




Bits error versus x




Bits error versus eps
| Original | 36.8 |
|---|---|
| Target | 15.0 |
| Herbie | 0.3 |
Initial program 36.8
rmApplied sin-sum21.8
rmApplied add-cube-cbrt21.9
Applied associate-*l*21.9
Taylor expanded around inf 21.8
Simplified0.4
rmApplied fma-neg0.3
Final simplification0.3
herbie shell --seed 2019172 +o rules:numerics
(FPCore (x eps)
:name "2sin (example 3.3)"
:herbie-target
(* 2.0 (* (cos (+ x (/ eps 2.0))) (sin (/ eps 2.0))))
(- (sin (+ x eps)) (sin x)))