\left(x + y\right) - \frac{\left(z - t\right) \cdot y}{a - t}\begin{array}{l}
\mathbf{if}\;a \le -1.2521412398963187 \cdot 10^{-142}:\\
\;\;\;\;\sqrt[3]{y - \left(\frac{y}{\sqrt[3]{a - t}} \cdot \frac{\sqrt[3]{z - t}}{\sqrt[3]{a - t}}\right) \cdot \frac{\sqrt[3]{z - t} \cdot \sqrt[3]{z - t}}{\sqrt[3]{a - t}}} \cdot \left(\sqrt[3]{y - \left(\frac{y}{\sqrt[3]{a - t}} \cdot \frac{\sqrt[3]{z - t}}{\sqrt[3]{a - t}}\right) \cdot \frac{\sqrt[3]{z - t} \cdot \sqrt[3]{z - t}}{\sqrt[3]{a - t}}} \cdot \sqrt[3]{y - \left(\frac{y}{\sqrt[3]{a - t}} \cdot \frac{\sqrt[3]{z - t}}{\sqrt[3]{a - t}}\right) \cdot \frac{\sqrt[3]{z - t} \cdot \sqrt[3]{z - t}}{\sqrt[3]{a - t}}}\right) + x\\
\mathbf{elif}\;a \le 9.465845835877107 \cdot 10^{-212}:\\
\;\;\;\;x + \frac{y \cdot z}{t}\\
\mathbf{else}:\\
\;\;\;\;\left(y - \frac{\sqrt[3]{z - t} \cdot \sqrt[3]{z - t}}{\sqrt[3]{a - t}} \cdot \left(\sqrt[3]{\sqrt[3]{z - t} \cdot \sqrt[3]{z - t}} \cdot \left(\frac{y}{\sqrt[3]{a - t}} \cdot \frac{\sqrt[3]{\sqrt[3]{z - t}}}{\sqrt[3]{a - t}}\right)\right)\right) + x\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r29791556 = x;
double r29791557 = y;
double r29791558 = r29791556 + r29791557;
double r29791559 = z;
double r29791560 = t;
double r29791561 = r29791559 - r29791560;
double r29791562 = r29791561 * r29791557;
double r29791563 = a;
double r29791564 = r29791563 - r29791560;
double r29791565 = r29791562 / r29791564;
double r29791566 = r29791558 - r29791565;
return r29791566;
}
double f(double x, double y, double z, double t, double a) {
double r29791567 = a;
double r29791568 = -1.2521412398963187e-142;
bool r29791569 = r29791567 <= r29791568;
double r29791570 = y;
double r29791571 = t;
double r29791572 = r29791567 - r29791571;
double r29791573 = cbrt(r29791572);
double r29791574 = r29791570 / r29791573;
double r29791575 = z;
double r29791576 = r29791575 - r29791571;
double r29791577 = cbrt(r29791576);
double r29791578 = r29791577 / r29791573;
double r29791579 = r29791574 * r29791578;
double r29791580 = r29791577 * r29791577;
double r29791581 = r29791580 / r29791573;
double r29791582 = r29791579 * r29791581;
double r29791583 = r29791570 - r29791582;
double r29791584 = cbrt(r29791583);
double r29791585 = r29791584 * r29791584;
double r29791586 = r29791584 * r29791585;
double r29791587 = x;
double r29791588 = r29791586 + r29791587;
double r29791589 = 9.465845835877107e-212;
bool r29791590 = r29791567 <= r29791589;
double r29791591 = r29791570 * r29791575;
double r29791592 = r29791591 / r29791571;
double r29791593 = r29791587 + r29791592;
double r29791594 = cbrt(r29791580);
double r29791595 = cbrt(r29791577);
double r29791596 = r29791595 / r29791573;
double r29791597 = r29791574 * r29791596;
double r29791598 = r29791594 * r29791597;
double r29791599 = r29791581 * r29791598;
double r29791600 = r29791570 - r29791599;
double r29791601 = r29791600 + r29791587;
double r29791602 = r29791590 ? r29791593 : r29791601;
double r29791603 = r29791569 ? r29791588 : r29791602;
return r29791603;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 16.3 |
|---|---|
| Target | 8.5 |
| Herbie | 8.3 |
if a < -1.2521412398963187e-142Initial program 14.9
rmApplied add-cube-cbrt15.0
Applied times-frac8.8
rmApplied add-cube-cbrt8.9
Applied times-frac8.9
Applied associate-*l*8.8
rmApplied associate--l+6.5
rmApplied add-cube-cbrt6.8
if -1.2521412398963187e-142 < a < 9.465845835877107e-212Initial program 20.3
rmApplied add-cube-cbrt20.6
Applied times-frac18.4
rmApplied add-cube-cbrt18.5
Applied times-frac18.4
Applied associate-*l*17.2
rmApplied associate--l+10.6
Taylor expanded around inf 9.1
if 9.465845835877107e-212 < a Initial program 15.7
rmApplied add-cube-cbrt15.8
Applied times-frac10.4
rmApplied add-cube-cbrt10.5
Applied times-frac10.5
Applied associate-*l*10.1
rmApplied associate--l+7.3
rmApplied *-un-lft-identity7.3
Applied add-cube-cbrt8.1
Applied cbrt-prod9.1
Applied times-frac8.9
Applied associate-*l*9.1
Final simplification8.3
herbie shell --seed 2019165
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisTick from plot-0.2.3.4, B"
: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))))