\sin \left(x + \varepsilon\right) - \sin x
\sin x \cdot \frac{\sqrt[3]{{\left({\left(\cos \varepsilon\right)}^{3} - 1\right)}^{3}}}{\cos \varepsilon \cdot \left(\cos \varepsilon + 1\right) + 1} + \cos x \cdot \sin \varepsilondouble f(double x, double eps) {
double r127037 = x;
double r127038 = eps;
double r127039 = r127037 + r127038;
double r127040 = sin(r127039);
double r127041 = sin(r127037);
double r127042 = r127040 - r127041;
return r127042;
}
double f(double x, double eps) {
double r127043 = x;
double r127044 = sin(r127043);
double r127045 = eps;
double r127046 = cos(r127045);
double r127047 = 3.0;
double r127048 = pow(r127046, r127047);
double r127049 = 1.0;
double r127050 = r127048 - r127049;
double r127051 = pow(r127050, r127047);
double r127052 = cbrt(r127051);
double r127053 = r127046 + r127049;
double r127054 = r127046 * r127053;
double r127055 = r127054 + r127049;
double r127056 = r127052 / r127055;
double r127057 = r127044 * r127056;
double r127058 = cos(r127043);
double r127059 = sin(r127045);
double r127060 = r127058 * r127059;
double r127061 = r127057 + r127060;
return r127061;
}




Bits error versus x




Bits error versus eps
Results
| Original | 36.9 |
|---|---|
| Target | 15.3 |
| Herbie | 0.4 |
Initial program 36.9
rmApplied sin-sum21.4
Applied associate--l+21.4
Taylor expanded around inf 21.4
Simplified0.4
rmApplied flip3--0.4
Simplified0.4
Simplified0.4
rmApplied add-cbrt-cube0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2020065
(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)))