\left(x + y\right) - \frac{\left(z - t\right) \cdot y}{a - t}\begin{array}{l}
\mathbf{if}\;a \le -1.607251433835726116796618339295063866918 \cdot 10^{-162}:\\
\;\;\;\;\left(x + y\right) - \frac{z - t}{\sqrt[3]{a - t} \cdot \left(\sqrt[3]{\sqrt[3]{a - t} \cdot \sqrt[3]{a - t}} \cdot \sqrt[3]{\sqrt[3]{a - t}}\right)} \cdot \frac{y}{\sqrt[3]{a - t}}\\
\mathbf{elif}\;a \le 4.889529153193783306974873529274732222518 \cdot 10^{-116}:\\
\;\;\;\;\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]{\frac{y}{\sqrt[3]{a - t}}} \cdot \sqrt[3]{\frac{y}{\sqrt[3]{a - t}}}\right)\right) \cdot \sqrt[3]{\frac{y}{\sqrt[3]{a - t}}}\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r411255 = x;
double r411256 = y;
double r411257 = r411255 + r411256;
double r411258 = z;
double r411259 = t;
double r411260 = r411258 - r411259;
double r411261 = r411260 * r411256;
double r411262 = a;
double r411263 = r411262 - r411259;
double r411264 = r411261 / r411263;
double r411265 = r411257 - r411264;
return r411265;
}
double f(double x, double y, double z, double t, double a) {
double r411266 = a;
double r411267 = -1.607251433835726e-162;
bool r411268 = r411266 <= r411267;
double r411269 = x;
double r411270 = y;
double r411271 = r411269 + r411270;
double r411272 = z;
double r411273 = t;
double r411274 = r411272 - r411273;
double r411275 = r411266 - r411273;
double r411276 = cbrt(r411275);
double r411277 = r411276 * r411276;
double r411278 = cbrt(r411277);
double r411279 = cbrt(r411276);
double r411280 = r411278 * r411279;
double r411281 = r411276 * r411280;
double r411282 = r411274 / r411281;
double r411283 = r411270 / r411276;
double r411284 = r411282 * r411283;
double r411285 = r411271 - r411284;
double r411286 = 4.8895291531937833e-116;
bool r411287 = r411266 <= r411286;
double r411288 = r411272 * r411270;
double r411289 = r411288 / r411273;
double r411290 = r411289 + r411269;
double r411291 = r411274 / r411277;
double r411292 = cbrt(r411283);
double r411293 = r411292 * r411292;
double r411294 = r411291 * r411293;
double r411295 = r411294 * r411292;
double r411296 = r411271 - r411295;
double r411297 = r411287 ? r411290 : r411296;
double r411298 = r411268 ? r411285 : r411297;
return r411298;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 16.0 |
|---|---|
| Target | 8.1 |
| Herbie | 9.1 |
if a < -1.607251433835726e-162Initial program 14.6
rmApplied add-cube-cbrt14.7
Applied times-frac8.8
rmApplied add-cube-cbrt8.8
Applied cbrt-prod8.9
if -1.607251433835726e-162 < a < 4.8895291531937833e-116Initial program 20.3
Taylor expanded around inf 10.5
if 4.8895291531937833e-116 < a Initial program 14.5
rmApplied add-cube-cbrt14.6
Applied times-frac8.4
rmApplied add-cube-cbrt8.5
Applied associate-*r*8.5
Final simplification9.1
herbie shell --seed 2019326
(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))))