x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z}\begin{array}{l}
\mathbf{if}\;x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z} \le -2.32490640665928128674486260054726540816 \cdot 10^{-294} \lor \neg \left(x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z} \le 0.0\right):\\
\;\;\;\;x + \left(\frac{y - z}{\sqrt[3]{a - z} \cdot \sqrt[3]{a - z}} \cdot \frac{\sqrt[3]{t - x} \cdot \sqrt[3]{t - x}}{\sqrt[3]{1}}\right) \cdot \frac{\sqrt[3]{t - x}}{\sqrt[3]{a - z}}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{x \cdot y}{z} + t\right) - \frac{t \cdot y}{z}\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r435361 = x;
double r435362 = y;
double r435363 = z;
double r435364 = r435362 - r435363;
double r435365 = t;
double r435366 = r435365 - r435361;
double r435367 = r435364 * r435366;
double r435368 = a;
double r435369 = r435368 - r435363;
double r435370 = r435367 / r435369;
double r435371 = r435361 + r435370;
return r435371;
}
double f(double x, double y, double z, double t, double a) {
double r435372 = x;
double r435373 = y;
double r435374 = z;
double r435375 = r435373 - r435374;
double r435376 = t;
double r435377 = r435376 - r435372;
double r435378 = r435375 * r435377;
double r435379 = a;
double r435380 = r435379 - r435374;
double r435381 = r435378 / r435380;
double r435382 = r435372 + r435381;
double r435383 = -2.3249064066592813e-294;
bool r435384 = r435382 <= r435383;
double r435385 = 0.0;
bool r435386 = r435382 <= r435385;
double r435387 = !r435386;
bool r435388 = r435384 || r435387;
double r435389 = cbrt(r435380);
double r435390 = r435389 * r435389;
double r435391 = r435375 / r435390;
double r435392 = cbrt(r435377);
double r435393 = r435392 * r435392;
double r435394 = 1.0;
double r435395 = cbrt(r435394);
double r435396 = r435393 / r435395;
double r435397 = r435391 * r435396;
double r435398 = r435392 / r435389;
double r435399 = r435397 * r435398;
double r435400 = r435372 + r435399;
double r435401 = r435372 * r435373;
double r435402 = r435401 / r435374;
double r435403 = r435402 + r435376;
double r435404 = r435376 * r435373;
double r435405 = r435404 / r435374;
double r435406 = r435403 - r435405;
double r435407 = r435388 ? r435400 : r435406;
return r435407;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 24.5 |
|---|---|
| Target | 11.8 |
| Herbie | 8.5 |
if (+ x (/ (* (- y z) (- t x)) (- a z))) < -2.3249064066592813e-294 or 0.0 < (+ x (/ (* (- y z) (- t x)) (- a z))) Initial program 21.3
rmApplied add-cube-cbrt21.8
Applied times-frac8.3
rmApplied *-un-lft-identity8.3
Applied cbrt-prod8.3
Applied add-cube-cbrt8.4
Applied times-frac8.4
Applied associate-*r*7.8
if -2.3249064066592813e-294 < (+ x (/ (* (- y z) (- t x)) (- a z))) < 0.0Initial program 60.7
Taylor expanded around inf 16.7
Final simplification8.5
herbie shell --seed 2019291
(FPCore (x y z t a)
:name "Graphics.Rendering.Chart.Axis.Types:invLinMap from Chart-1.5.3"
:precision binary64
:herbie-target
(if (< z -1.25361310560950359e188) (- t (* (/ y z) (- t x))) (if (< z 4.44670236911381103e64) (+ x (/ (- y z) (/ (- a z) (- t x)))) (- t (* (/ y z) (- t x)))))
(+ x (/ (* (- y z) (- t x)) (- a z))))