\left(x + y\right) - \frac{\left(z - t\right) \cdot y}{a - t}\begin{array}{l}
\mathbf{if}\;t \le -4.1831138401846353 \cdot 10^{161} \lor \neg \left(t \le 4.07362133015771308 \cdot 10^{220}\right):\\
\;\;\;\;\frac{z \cdot y}{t} + x\\
\mathbf{else}:\\
\;\;\;\;\left(x + y\right) - \left(\frac{z - t}{\sqrt[3]{a - t} \cdot \sqrt[3]{a - t}} \cdot \frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\sqrt[3]{\sqrt[3]{a - t} \cdot \sqrt[3]{a - t}}}\right) \cdot \frac{\sqrt[3]{y}}{\sqrt[3]{\sqrt[3]{a - t}}}\\
\end{array}double code(double x, double y, double z, double t, double a) {
return ((double) (((double) (x + y)) - ((double) (((double) (((double) (z - t)) * y)) / ((double) (a - t))))));
}
double code(double x, double y, double z, double t, double a) {
double VAR;
if (((t <= -4.1831138401846353e+161) || !(t <= 4.073621330157713e+220))) {
VAR = ((double) (((double) (((double) (z * y)) / t)) + x));
} else {
VAR = ((double) (((double) (x + y)) - ((double) (((double) (((double) (((double) (z - t)) / ((double) (((double) cbrt(((double) (a - t)))) * ((double) cbrt(((double) (a - t)))))))) * ((double) (((double) (((double) cbrt(y)) * ((double) cbrt(y)))) / ((double) cbrt(((double) (((double) cbrt(((double) (a - t)))) * ((double) cbrt(((double) (a - t)))))))))))) * ((double) (((double) cbrt(y)) / ((double) cbrt(((double) cbrt(((double) (a - t))))))))))));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 16.7 |
|---|---|
| Target | 8.3 |
| Herbie | 9.9 |
if t < -4.1831138401846353e161 or 4.07362133015771308e220 < t Initial program 34.5
Taylor expanded around inf 16.7
if -4.1831138401846353e161 < t < 4.07362133015771308e220Initial program 12.0
rmApplied add-cube-cbrt12.2
Applied times-frac8.4
rmApplied add-cube-cbrt8.5
Applied cbrt-prod8.5
Applied add-cube-cbrt8.5
Applied times-frac8.5
Applied associate-*r*8.2
Final simplification9.9
herbie shell --seed 2020152
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisTick from plot-0.2.3.4, B"
:precision binary64
:herbie-target
(if (< (- (+ x y) (/ (* (- z t) y) (- a t))) -1.3664970889390727e-07) (- (+ y x) (* (* (- z t) (/ 1.0 (- a t))) y)) (if (< (- (+ x y) (/ (* (- z t) y) (- a t))) 1.4754293444577233e-239) (/ (- (* y (- a z)) (* x t)) (- a t)) (- (+ y x) (* (* (- z t) (/ 1.0 (- a t))) y))))
(- (+ x y) (/ (* (- z t) y) (- a t))))