x + \frac{y \cdot \left(z - t\right)}{z - a}\begin{array}{l}
\mathbf{if}\;y \leq -3.099128616291226 \cdot 10^{-118}:\\
\;\;\;\;x + y \cdot \left(\frac{z}{z - a} - \frac{t}{z - a}\right)\\
\mathbf{elif}\;y \leq 7.043274856239975 \cdot 10^{-34}:\\
\;\;\;\;x + \frac{y \cdot \left(z - t\right)}{z - a}\\
\mathbf{else}:\\
\;\;\;\;x + y \cdot \frac{1}{\frac{z - a}{z - t}}\\
\end{array}(FPCore (x y z t a) :precision binary64 (+ x (/ (* y (- z t)) (- z a))))
(FPCore (x y z t a)
:precision binary64
(if (<= y -3.099128616291226e-118)
(+ x (* y (- (/ z (- z a)) (/ t (- z a)))))
(if (<= y 7.043274856239975e-34)
(+ x (/ (* y (- z t)) (- z a)))
(+ x (* y (/ 1.0 (/ (- z a) (- z t))))))))double code(double x, double y, double z, double t, double a) {
return ((double) (x + (((double) (y * ((double) (z - t)))) / ((double) (z - a)))));
}
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((y <= -3.099128616291226e-118)) {
tmp = ((double) (x + ((double) (y * ((double) ((z / ((double) (z - a))) - (t / ((double) (z - a)))))))));
} else {
double tmp_1;
if ((y <= 7.043274856239975e-34)) {
tmp_1 = ((double) (x + (((double) (y * ((double) (z - t)))) / ((double) (z - a)))));
} else {
tmp_1 = ((double) (x + ((double) (y * (1.0 / (((double) (z - a)) / ((double) (z - t))))))));
}
tmp = tmp_1;
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 10.2 |
|---|---|
| Target | 1.2 |
| Herbie | 0.5 |
if y < -3.09912861629122609e-118Initial program Error: 14.9 bits
SimplifiedError: 0.6 bits
rmApplied div-subError: 0.6 bits
if -3.09912861629122609e-118 < y < 7.04327485623997529e-34Initial program Error: 0.5 bits
if 7.04327485623997529e-34 < y Initial program Error: 20.1 bits
SimplifiedError: 0.4 bits
rmApplied clear-numError: 0.5 bits
Final simplificationError: 0.5 bits
herbie shell --seed 2020203
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisTicks from plot-0.2.3.4, A"
:precision binary64
:herbie-target
(+ x (/ y (/ (- z a) (- z t))))
(+ x (/ (* y (- z t)) (- z a))))