\left(x + y\right) - \frac{\left(z - t\right) \cdot y}{a - t}\begin{array}{l}
\mathbf{if}\;a \le -2.78725408054180311 \cdot 10^{-108}:\\
\;\;\;\;\left(x + y\right) - \left(\frac{\frac{z - t}{\sqrt[3]{a - t} \cdot \sqrt[3]{a - t}}}{\sqrt[3]{\sqrt[3]{a - t} \cdot \sqrt[3]{a - t}}} \cdot \frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\sqrt[3]{\sqrt[3]{\sqrt[3]{a - t}}} \cdot \sqrt[3]{\sqrt[3]{\sqrt[3]{a - t}}}}\right) \cdot \frac{\sqrt[3]{y}}{\sqrt[3]{\sqrt[3]{\sqrt[3]{a - t}}}}\\
\mathbf{elif}\;a \le 1.1307683724034909 \cdot 10^{-188}:\\
\;\;\;\;\frac{z \cdot y}{t} + x\\
\mathbf{else}:\\
\;\;\;\;\left(x + y\right) - \frac{\frac{\frac{z - t}{\sqrt[3]{a - t} \cdot \sqrt[3]{a - t}}}{\sqrt[3]{\sqrt[3]{a - t} \cdot \sqrt[3]{a - t}}}}{\sqrt[3]{\sqrt[3]{\sqrt[3]{a - t}} \cdot \sqrt[3]{\sqrt[3]{a - t}}}} \cdot \frac{y}{\sqrt[3]{\sqrt[3]{\sqrt[3]{a - t}}}}\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r691576 = x;
double r691577 = y;
double r691578 = r691576 + r691577;
double r691579 = z;
double r691580 = t;
double r691581 = r691579 - r691580;
double r691582 = r691581 * r691577;
double r691583 = a;
double r691584 = r691583 - r691580;
double r691585 = r691582 / r691584;
double r691586 = r691578 - r691585;
return r691586;
}
double f(double x, double y, double z, double t, double a) {
double r691587 = a;
double r691588 = -2.787254080541803e-108;
bool r691589 = r691587 <= r691588;
double r691590 = x;
double r691591 = y;
double r691592 = r691590 + r691591;
double r691593 = z;
double r691594 = t;
double r691595 = r691593 - r691594;
double r691596 = r691587 - r691594;
double r691597 = cbrt(r691596);
double r691598 = r691597 * r691597;
double r691599 = r691595 / r691598;
double r691600 = cbrt(r691598);
double r691601 = r691599 / r691600;
double r691602 = cbrt(r691591);
double r691603 = r691602 * r691602;
double r691604 = cbrt(r691597);
double r691605 = cbrt(r691604);
double r691606 = r691605 * r691605;
double r691607 = r691603 / r691606;
double r691608 = r691601 * r691607;
double r691609 = r691602 / r691605;
double r691610 = r691608 * r691609;
double r691611 = r691592 - r691610;
double r691612 = 1.1307683724034909e-188;
bool r691613 = r691587 <= r691612;
double r691614 = r691593 * r691591;
double r691615 = r691614 / r691594;
double r691616 = r691615 + r691590;
double r691617 = r691604 * r691604;
double r691618 = cbrt(r691617);
double r691619 = r691601 / r691618;
double r691620 = r691591 / r691605;
double r691621 = r691619 * r691620;
double r691622 = r691592 - r691621;
double r691623 = r691613 ? r691616 : r691622;
double r691624 = r691589 ? r691611 : r691623;
return r691624;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 16.5 |
|---|---|
| Target | 8.5 |
| Herbie | 9.2 |
if a < -2.787254080541803e-108Initial program 15.5
rmApplied add-cube-cbrt15.6
Applied times-frac8.8
rmApplied add-cube-cbrt8.9
Applied cbrt-prod8.9
Applied *-un-lft-identity8.9
Applied times-frac8.9
Applied associate-*r*8.8
Simplified8.8
rmApplied add-cube-cbrt8.9
Applied add-cube-cbrt8.9
Applied times-frac8.9
Applied associate-*r*8.8
if -2.787254080541803e-108 < a < 1.1307683724034909e-188Initial program 20.0
Taylor expanded around inf 9.7
if 1.1307683724034909e-188 < a Initial program 15.3
rmApplied add-cube-cbrt15.4
Applied times-frac9.2
rmApplied add-cube-cbrt9.2
Applied cbrt-prod9.2
Applied *-un-lft-identity9.2
Applied times-frac9.2
Applied associate-*r*9.2
Simplified9.2
rmApplied add-cube-cbrt9.2
Applied cbrt-prod9.2
Applied *-un-lft-identity9.2
Applied times-frac9.2
Applied associate-*r*9.3
Simplified9.3
Final simplification9.2
herbie shell --seed 2020021
(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))))