\left(x + y\right) - \frac{\left(z - t\right) \cdot y}{a - t}\begin{array}{l}
\mathbf{if}\;\left(x + y\right) - \frac{\left(z - t\right) \cdot y}{a - t} \le -5.14917433408142322 \cdot 10^{-200}:\\
\;\;\;\;\left(x + y\right) - \frac{\frac{\sqrt[3]{z - t} \cdot \sqrt[3]{z - t}}{\sqrt[3]{a - t}}}{\sqrt[3]{\sqrt[3]{a - t}}} \cdot \left(\frac{\frac{\sqrt[3]{z - t}}{\sqrt[3]{a - t}}}{\sqrt[3]{\sqrt[3]{a - t}}} \cdot \frac{y}{\sqrt[3]{\sqrt[3]{a - t}}}\right)\\
\mathbf{elif}\;\left(x + y\right) - \frac{\left(z - t\right) \cdot y}{a - t} \le 0.0:\\
\;\;\;\;\frac{z \cdot y}{t} + x\\
\mathbf{else}:\\
\;\;\;\;\left(x + y\right) - \frac{\sqrt[3]{\frac{z - t}{\sqrt[3]{a - t} \cdot \sqrt[3]{a - t}}} \cdot \sqrt[3]{\frac{z - t}{\sqrt[3]{a - t} \cdot \sqrt[3]{a - t}}}}{\sqrt[3]{\sqrt[3]{a - t}}} \cdot \left(\frac{\sqrt[3]{\frac{z - t}{\sqrt[3]{a - t} \cdot \sqrt[3]{a - t}}}}{\sqrt[3]{\sqrt[3]{a - t}}} \cdot \frac{y}{\sqrt[3]{\sqrt[3]{a - t}}}\right)\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r639942 = x;
double r639943 = y;
double r639944 = r639942 + r639943;
double r639945 = z;
double r639946 = t;
double r639947 = r639945 - r639946;
double r639948 = r639947 * r639943;
double r639949 = a;
double r639950 = r639949 - r639946;
double r639951 = r639948 / r639950;
double r639952 = r639944 - r639951;
return r639952;
}
double f(double x, double y, double z, double t, double a) {
double r639953 = x;
double r639954 = y;
double r639955 = r639953 + r639954;
double r639956 = z;
double r639957 = t;
double r639958 = r639956 - r639957;
double r639959 = r639958 * r639954;
double r639960 = a;
double r639961 = r639960 - r639957;
double r639962 = r639959 / r639961;
double r639963 = r639955 - r639962;
double r639964 = -5.149174334081423e-200;
bool r639965 = r639963 <= r639964;
double r639966 = cbrt(r639958);
double r639967 = r639966 * r639966;
double r639968 = cbrt(r639961);
double r639969 = r639967 / r639968;
double r639970 = cbrt(r639968);
double r639971 = r639969 / r639970;
double r639972 = r639966 / r639968;
double r639973 = r639972 / r639970;
double r639974 = r639954 / r639970;
double r639975 = r639973 * r639974;
double r639976 = r639971 * r639975;
double r639977 = r639955 - r639976;
double r639978 = 0.0;
bool r639979 = r639963 <= r639978;
double r639980 = r639956 * r639954;
double r639981 = r639980 / r639957;
double r639982 = r639981 + r639953;
double r639983 = r639968 * r639968;
double r639984 = r639958 / r639983;
double r639985 = cbrt(r639984);
double r639986 = r639985 * r639985;
double r639987 = r639986 / r639970;
double r639988 = r639985 / r639970;
double r639989 = r639988 * r639974;
double r639990 = r639987 * r639989;
double r639991 = r639955 - r639990;
double r639992 = r639979 ? r639982 : r639991;
double r639993 = r639965 ? r639977 : r639992;
return r639993;
}




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.5 |
| Herbie | 8.5 |
if (- (+ x y) (/ (* (- z t) y) (- a t))) < -5.149174334081423e-200Initial program 13.0
rmApplied add-cube-cbrt13.2
Applied times-frac7.4
rmApplied add-cube-cbrt7.4
Applied cbrt-prod7.4
Applied *-un-lft-identity7.4
Applied times-frac7.5
Applied associate-*r*7.5
Simplified7.5
rmApplied cbrt-prod7.6
Applied add-cube-cbrt7.6
Applied times-frac7.6
Applied times-frac7.6
Applied associate-*l*6.8
if -5.149174334081423e-200 < (- (+ x y) (/ (* (- z t) y) (- a t))) < 0.0Initial program 54.9
Taylor expanded around inf 20.5
if 0.0 < (- (+ x y) (/ (* (- z t) y) (- a t))) Initial program 12.9
rmApplied add-cube-cbrt13.1
Applied times-frac7.9
rmApplied add-cube-cbrt7.9
Applied cbrt-prod7.9
Applied *-un-lft-identity7.9
Applied times-frac7.9
Applied associate-*r*7.9
Simplified7.9
rmApplied cbrt-prod8.0
Applied add-cube-cbrt8.0
Applied times-frac8.0
Applied associate-*l*7.8
Final simplification8.5
herbie shell --seed 2020025
(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 (- 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 (- a t))) y))))
(- (+ x y) (/ (* (- z t) y) (- a t))))