\sin \left(x + \varepsilon\right) - \sin x
\left(\sin x \cdot \cos \varepsilon + \left(-\sin x\right)\right) + \cos x \cdot \sin \varepsilon
double f(double x, double eps) {
double r118232 = x;
double r118233 = eps;
double r118234 = r118232 + r118233;
double r118235 = sin(r118234);
double r118236 = sin(r118232);
double r118237 = r118235 - r118236;
return r118237;
}
double f(double x, double eps) {
double r118238 = x;
double r118239 = sin(r118238);
double r118240 = eps;
double r118241 = cos(r118240);
double r118242 = r118239 * r118241;
double r118243 = -r118239;
double r118244 = r118242 + r118243;
double r118245 = cos(r118238);
double r118246 = sin(r118240);
double r118247 = r118245 * r118246;
double r118248 = r118244 + r118247;
return r118248;
}




Bits error versus x




Bits error versus eps
Results
| Original | 37.2 |
|---|---|
| Target | 15.1 |
| Herbie | 0.4 |
Initial program 37.2
rmApplied sin-sum22.0
Applied associate--l+22.0
Taylor expanded around inf 22.0
Simplified0.4
rmApplied sub-neg0.4
Applied distribute-lft-in0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2020018
(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)))