\sin \left(x + \varepsilon\right) - \sin x
\sqrt[3]{\left(\left(\sin x \cdot \sin x\right) \cdot \sin x\right) \cdot \left(\left(\left(\cos \varepsilon - 1\right) \cdot \left(\cos \varepsilon - 1\right)\right) \cdot \left(\cos \varepsilon - 1\right)\right)} + \cos x \cdot \sin \varepsilondouble f(double x, double eps) {
double r236421 = x;
double r236422 = eps;
double r236423 = r236421 + r236422;
double r236424 = sin(r236423);
double r236425 = sin(r236421);
double r236426 = r236424 - r236425;
return r236426;
}
double f(double x, double eps) {
double r236427 = x;
double r236428 = sin(r236427);
double r236429 = r236428 * r236428;
double r236430 = r236429 * r236428;
double r236431 = eps;
double r236432 = cos(r236431);
double r236433 = 1.0;
double r236434 = r236432 - r236433;
double r236435 = r236434 * r236434;
double r236436 = r236435 * r236434;
double r236437 = r236430 * r236436;
double r236438 = cbrt(r236437);
double r236439 = cos(r236427);
double r236440 = sin(r236431);
double r236441 = r236439 * r236440;
double r236442 = r236438 + r236441;
return r236442;
}




Bits error versus x




Bits error versus eps
Results
| Original | 36.7 |
|---|---|
| Target | 14.8 |
| Herbie | 0.4 |
Initial program 36.7
rmApplied sin-sum21.8
rmApplied *-un-lft-identity21.8
Applied *-un-lft-identity21.8
Applied distribute-lft-out--21.8
Simplified0.4
rmApplied add-cbrt-cube0.4
Applied add-cbrt-cube0.4
Applied cbrt-unprod0.4
Simplified0.4
rmApplied add-cbrt-cube0.4
Applied add-cbrt-cube0.5
Applied cbrt-unprod0.4
Applied rem-cube-cbrt0.4
Final simplification0.4
herbie shell --seed 2019356
(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)))