\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 r166214 = x;
double r166215 = eps;
double r166216 = r166214 + r166215;
double r166217 = sin(r166216);
double r166218 = sin(r166214);
double r166219 = r166217 - r166218;
return r166219;
}
double f(double x, double eps) {
double r166220 = x;
double r166221 = sin(r166220);
double r166222 = eps;
double r166223 = cos(r166222);
double r166224 = 1.0;
double r166225 = r166223 - r166224;
double r166226 = r166221 * r166225;
double r166227 = 3.0;
double r166228 = pow(r166226, r166227);
double r166229 = cbrt(r166228);
double r166230 = cos(r166220);
double r166231 = sin(r166222);
double r166232 = r166230 * r166231;
double r166233 = r166229 + r166232;
return r166233;
}




Bits error versus x




Bits error versus eps
Results
| Original | 36.5 |
|---|---|
| Target | 15.2 |
| Herbie | 0.4 |
Initial program 36.5
rmApplied sin-sum21.3
rmApplied *-un-lft-identity21.3
Applied *-un-lft-identity21.3
Applied distribute-lft-out--21.3
Simplified0.4
rmApplied sub-neg0.4
Applied distribute-lft-in0.4
Simplified0.4
rmApplied add-cbrt-cube0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2020034
(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)))