\sin \left(x + \varepsilon\right) - \sin x
\sqrt[3]{{\left(\sin x \cdot \left(\cos \varepsilon - 1\right)\right)}^{3}} + \cos x \cdot \sin \varepsilondouble f(double x, double eps) {
double r131701 = x;
double r131702 = eps;
double r131703 = r131701 + r131702;
double r131704 = sin(r131703);
double r131705 = sin(r131701);
double r131706 = r131704 - r131705;
return r131706;
}
double f(double x, double eps) {
double r131707 = x;
double r131708 = sin(r131707);
double r131709 = eps;
double r131710 = cos(r131709);
double r131711 = 1.0;
double r131712 = r131710 - r131711;
double r131713 = r131708 * r131712;
double r131714 = 3.0;
double r131715 = pow(r131713, r131714);
double r131716 = cbrt(r131715);
double r131717 = cos(r131707);
double r131718 = sin(r131709);
double r131719 = r131717 * r131718;
double r131720 = r131716 + r131719;
return r131720;
}




Bits error versus x




Bits error versus eps
Results
| Original | 37.4 |
|---|---|
| Target | 15.1 |
| Herbie | 0.5 |
Initial program 37.4
rmApplied sin-sum22.2
Taylor expanded around inf 22.2
Simplified0.4
rmApplied add-cbrt-cube0.5
Applied add-cbrt-cube0.5
Applied cbrt-unprod0.5
Simplified0.5
Final simplification0.5
herbie shell --seed 2020047
(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)))