x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z}\begin{array}{l}
\mathbf{if}\;z \le -9.79529762078173593280535198605757013158 \cdot 10^{279}:\\
\;\;\;\;t + \left(\frac{x}{z} - \frac{t}{z}\right) \cdot y\\
\mathbf{elif}\;z \le 1.598379849628712025701715812034079117588 \cdot 10^{185}:\\
\;\;\;\;x + \frac{\sqrt[3]{y - z} \cdot \sqrt[3]{y - z}}{\sqrt[3]{a - z}} \cdot \left(\frac{t - x}{\sqrt[3]{a - z}} \cdot \frac{\sqrt[3]{y - z}}{\sqrt[3]{a - z}}\right)\\
\mathbf{else}:\\
\;\;\;\;t + \left(\frac{x}{z} - \frac{t}{z}\right) \cdot y\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r33524512 = x;
double r33524513 = y;
double r33524514 = z;
double r33524515 = r33524513 - r33524514;
double r33524516 = t;
double r33524517 = r33524516 - r33524512;
double r33524518 = r33524515 * r33524517;
double r33524519 = a;
double r33524520 = r33524519 - r33524514;
double r33524521 = r33524518 / r33524520;
double r33524522 = r33524512 + r33524521;
return r33524522;
}
double f(double x, double y, double z, double t, double a) {
double r33524523 = z;
double r33524524 = -9.795297620781736e+279;
bool r33524525 = r33524523 <= r33524524;
double r33524526 = t;
double r33524527 = x;
double r33524528 = r33524527 / r33524523;
double r33524529 = r33524526 / r33524523;
double r33524530 = r33524528 - r33524529;
double r33524531 = y;
double r33524532 = r33524530 * r33524531;
double r33524533 = r33524526 + r33524532;
double r33524534 = 1.598379849628712e+185;
bool r33524535 = r33524523 <= r33524534;
double r33524536 = r33524531 - r33524523;
double r33524537 = cbrt(r33524536);
double r33524538 = r33524537 * r33524537;
double r33524539 = a;
double r33524540 = r33524539 - r33524523;
double r33524541 = cbrt(r33524540);
double r33524542 = r33524538 / r33524541;
double r33524543 = r33524526 - r33524527;
double r33524544 = r33524543 / r33524541;
double r33524545 = r33524537 / r33524541;
double r33524546 = r33524544 * r33524545;
double r33524547 = r33524542 * r33524546;
double r33524548 = r33524527 + r33524547;
double r33524549 = r33524535 ? r33524548 : r33524533;
double r33524550 = r33524525 ? r33524533 : r33524549;
return r33524550;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 24.8 |
|---|---|
| Target | 12.0 |
| Herbie | 10.7 |
if z < -9.795297620781736e+279 or 1.598379849628712e+185 < z Initial program 51.1
rmApplied add-cube-cbrt51.4
Applied times-frac26.1
rmApplied add-cube-cbrt25.6
Applied times-frac25.6
Applied associate-*l*25.6
Taylor expanded around inf 23.6
Simplified13.6
if -9.795297620781736e+279 < z < 1.598379849628712e+185Initial program 20.4
rmApplied add-cube-cbrt20.8
Applied times-frac10.6
rmApplied add-cube-cbrt10.6
Applied times-frac10.6
Applied associate-*l*10.3
Final simplification10.7
herbie shell --seed 2019200
(FPCore (x y z t a)
:name "Graphics.Rendering.Chart.Axis.Types:invLinMap from Chart-1.5.3"
:herbie-target
(if (< z -1.2536131056095036e+188) (- t (* (/ y z) (- t x))) (if (< z 4.446702369113811e+64) (+ x (/ (- y z) (/ (- a z) (- t x)))) (- t (* (/ y z) (- t x)))))
(+ x (/ (* (- y z) (- t x)) (- a z))))