x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z}\begin{array}{l}
\mathbf{if}\;z \le -5.229039851220961 \cdot 10^{234} \lor \neg \left(z \le 1.15096756339151705 \cdot 10^{227}\right):\\
\;\;\;\;y \cdot \left(\frac{x}{z} - \frac{t}{z}\right) + t\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y - z}{\sqrt[3]{a - z} \cdot \sqrt[3]{a - z}} \cdot \frac{t - x}{\sqrt[3]{a - z}}\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r806573 = x;
double r806574 = y;
double r806575 = z;
double r806576 = r806574 - r806575;
double r806577 = t;
double r806578 = r806577 - r806573;
double r806579 = r806576 * r806578;
double r806580 = a;
double r806581 = r806580 - r806575;
double r806582 = r806579 / r806581;
double r806583 = r806573 + r806582;
return r806583;
}
double f(double x, double y, double z, double t, double a) {
double r806584 = z;
double r806585 = -5.229039851220961e+234;
bool r806586 = r806584 <= r806585;
double r806587 = 1.150967563391517e+227;
bool r806588 = r806584 <= r806587;
double r806589 = !r806588;
bool r806590 = r806586 || r806589;
double r806591 = y;
double r806592 = x;
double r806593 = r806592 / r806584;
double r806594 = t;
double r806595 = r806594 / r806584;
double r806596 = r806593 - r806595;
double r806597 = r806591 * r806596;
double r806598 = r806597 + r806594;
double r806599 = r806591 - r806584;
double r806600 = a;
double r806601 = r806600 - r806584;
double r806602 = cbrt(r806601);
double r806603 = r806602 * r806602;
double r806604 = r806599 / r806603;
double r806605 = r806594 - r806592;
double r806606 = r806605 / r806602;
double r806607 = r806604 * r806606;
double r806608 = r806592 + r806607;
double r806609 = r806590 ? r806598 : r806608;
return r806609;
}




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 | 12.0 |
| Herbie | 10.3 |
if z < -5.229039851220961e+234 or 1.150967563391517e+227 < z Initial program 52.2
rmApplied add-cube-cbrt52.4
Applied times-frac27.7
rmApplied add-cube-cbrt27.7
Applied cbrt-prod27.7
Applied associate-*r*27.7
rmApplied add-cube-cbrt27.9
Applied associate-*r*27.9
Taylor expanded around inf 23.2
Simplified13.1
if -5.229039851220961e+234 < z < 1.150967563391517e+227Initial program 19.8
rmApplied add-cube-cbrt20.3
Applied times-frac9.9
Final simplification10.3
herbie shell --seed 2020046
(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.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))))