\sin \left(x + \varepsilon\right) - \sin x
\begin{array}{l}
\mathbf{if}\;\varepsilon \le -1.755630184432992 \cdot 10^{-08}:\\
\;\;\;\;\left(\sin x \cdot \cos \varepsilon + \cos x \cdot \sin \varepsilon\right) - \sin x\\
\mathbf{elif}\;\varepsilon \le 7.800526234462867 \cdot 10^{-25}:\\
\;\;\;\;2 \cdot \left(\sin \left(\frac{\varepsilon}{2}\right) \cdot \cos \left(\frac{x + \left(x + \varepsilon\right)}{2}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\cos x \cdot \sin \varepsilon - \sin x\right) + \sin x \cdot \cos \varepsilon\\
\end{array}double f(double x, double eps) {
double r5270634 = x;
double r5270635 = eps;
double r5270636 = r5270634 + r5270635;
double r5270637 = sin(r5270636);
double r5270638 = sin(r5270634);
double r5270639 = r5270637 - r5270638;
return r5270639;
}
double f(double x, double eps) {
double r5270640 = eps;
double r5270641 = -1.755630184432992e-08;
bool r5270642 = r5270640 <= r5270641;
double r5270643 = x;
double r5270644 = sin(r5270643);
double r5270645 = cos(r5270640);
double r5270646 = r5270644 * r5270645;
double r5270647 = cos(r5270643);
double r5270648 = sin(r5270640);
double r5270649 = r5270647 * r5270648;
double r5270650 = r5270646 + r5270649;
double r5270651 = r5270650 - r5270644;
double r5270652 = 7.800526234462867e-25;
bool r5270653 = r5270640 <= r5270652;
double r5270654 = 2.0;
double r5270655 = r5270640 / r5270654;
double r5270656 = sin(r5270655);
double r5270657 = r5270643 + r5270640;
double r5270658 = r5270643 + r5270657;
double r5270659 = r5270658 / r5270654;
double r5270660 = cos(r5270659);
double r5270661 = r5270656 * r5270660;
double r5270662 = r5270654 * r5270661;
double r5270663 = r5270649 - r5270644;
double r5270664 = r5270663 + r5270646;
double r5270665 = r5270653 ? r5270662 : r5270664;
double r5270666 = r5270642 ? r5270651 : r5270665;
return r5270666;
}




Bits error versus x




Bits error versus eps
Results
| Original | 37.1 |
|---|---|
| Target | 15.3 |
| Herbie | 0.7 |
if eps < -1.755630184432992e-08Initial program 29.9
rmApplied sin-sum0.6
if -1.755630184432992e-08 < eps < 7.800526234462867e-25Initial program 45.1
rmApplied diff-sin45.1
Simplified0.3
if 7.800526234462867e-25 < eps Initial program 29.9
rmApplied sin-sum1.6
Applied associate--l+1.7
Final simplification0.7
herbie shell --seed 2019135
(FPCore (x eps)
:name "2sin (example 3.3)"
:herbie-target
(* 2 (* (cos (+ x (/ eps 2))) (sin (/ eps 2))))
(- (sin (+ x eps)) (sin x)))