\left(x + y\right) - \frac{\left(z - t\right) \cdot y}{a - t}\begin{array}{l}
\mathbf{if}\;a \le -3.022383678732151115662024570927261203627 \cdot 10^{-117}:\\
\;\;\;\;\left(x + y\right) - \left(\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}}}\right) \cdot \left(\sqrt[3]{\frac{z - t}{\sqrt[3]{a - t} \cdot \sqrt[3]{a - t}}} \cdot \frac{y}{\sqrt[3]{a - t}}\right)\\
\mathbf{elif}\;a \le 1.132547627968652530752459014111872971988 \cdot 10^{-26}:\\
\;\;\;\;\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 \left(\sqrt[3]{y} \cdot \sqrt[3]{y}\right)\right) \cdot \frac{\sqrt[3]{y}}{\sqrt[3]{a - t}}\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r439228 = x;
double r439229 = y;
double r439230 = r439228 + r439229;
double r439231 = z;
double r439232 = t;
double r439233 = r439231 - r439232;
double r439234 = r439233 * r439229;
double r439235 = a;
double r439236 = r439235 - r439232;
double r439237 = r439234 / r439236;
double r439238 = r439230 - r439237;
return r439238;
}
double f(double x, double y, double z, double t, double a) {
double r439239 = a;
double r439240 = -3.022383678732151e-117;
bool r439241 = r439239 <= r439240;
double r439242 = x;
double r439243 = y;
double r439244 = r439242 + r439243;
double r439245 = z;
double r439246 = t;
double r439247 = r439245 - r439246;
double r439248 = r439239 - r439246;
double r439249 = cbrt(r439248);
double r439250 = r439249 * r439249;
double r439251 = r439247 / r439250;
double r439252 = cbrt(r439251);
double r439253 = r439252 * r439252;
double r439254 = r439243 / r439249;
double r439255 = r439252 * r439254;
double r439256 = r439253 * r439255;
double r439257 = r439244 - r439256;
double r439258 = 1.1325476279686525e-26;
bool r439259 = r439239 <= r439258;
double r439260 = r439245 * r439243;
double r439261 = r439260 / r439246;
double r439262 = r439261 + r439242;
double r439263 = cbrt(r439243);
double r439264 = r439263 * r439263;
double r439265 = r439251 * r439264;
double r439266 = r439263 / r439249;
double r439267 = r439265 * r439266;
double r439268 = r439244 - r439267;
double r439269 = r439259 ? r439262 : r439268;
double r439270 = r439241 ? r439257 : r439269;
return r439270;
}




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.8 |
| Herbie | 10.0 |
if a < -3.022383678732151e-117Initial program 14.6
rmApplied add-cube-cbrt14.7
Applied times-frac8.4
rmApplied add-cube-cbrt8.5
Applied associate-*l*8.5
if -3.022383678732151e-117 < a < 1.1325476279686525e-26Initial program 20.4
Taylor expanded around inf 13.7
if 1.1325476279686525e-26 < a Initial program 14.8
rmApplied add-cube-cbrt14.8
Applied times-frac7.1
rmApplied *-un-lft-identity7.1
Applied cbrt-prod7.1
Applied add-cube-cbrt7.2
Applied times-frac7.2
Applied associate-*r*7.1
Simplified7.1
Final simplification10.0
herbie shell --seed 2019306
(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-7) (- (+ y x) (* (* (- z t) (/ 1 (- a t))) y)) (if (< (- (+ x y) (/ (* (- z t) y) (- a t))) 1.47542934445772333e-239) (/ (- (* y (- a z)) (* x t)) (- a t)) (- (+ y x) (* (* (- z t) (/ 1 (- a t))) y))))
(- (+ x y) (/ (* (- z t) y) (- a t))))