\sin \left(x + \varepsilon\right) - \sin x
\sin x \cdot \log \left(e^{\cos \varepsilon - 1}\right) + \cos x \cdot \sin \varepsilondouble f(double x, double eps) {
double r97833 = x;
double r97834 = eps;
double r97835 = r97833 + r97834;
double r97836 = sin(r97835);
double r97837 = sin(r97833);
double r97838 = r97836 - r97837;
return r97838;
}
double f(double x, double eps) {
double r97839 = x;
double r97840 = sin(r97839);
double r97841 = eps;
double r97842 = cos(r97841);
double r97843 = 1.0;
double r97844 = r97842 - r97843;
double r97845 = exp(r97844);
double r97846 = log(r97845);
double r97847 = r97840 * r97846;
double r97848 = cos(r97839);
double r97849 = sin(r97841);
double r97850 = r97848 * r97849;
double r97851 = r97847 + r97850;
return r97851;
}




Bits error versus x




Bits error versus eps
Results
| Original | 36.6 |
|---|---|
| Target | 15.3 |
| Herbie | 0.4 |
Initial program 36.6
rmApplied sin-sum21.2
rmApplied *-un-lft-identity21.2
Applied *-un-lft-identity21.2
Applied distribute-lft-out--21.2
Simplified0.4
rmApplied add-log-exp0.4
Applied add-log-exp0.4
Applied diff-log0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2020056
(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)))