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):\\
\;\;\;\;t + y \cdot \left(\frac{x}{z} - \frac{t}{z}\right)\\
\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 r509747 = x;
double r509748 = y;
double r509749 = z;
double r509750 = r509748 - r509749;
double r509751 = t;
double r509752 = r509751 - r509747;
double r509753 = r509750 * r509752;
double r509754 = a;
double r509755 = r509754 - r509749;
double r509756 = r509753 / r509755;
double r509757 = r509747 + r509756;
return r509757;
}
double f(double x, double y, double z, double t, double a) {
double r509758 = z;
double r509759 = -5.229039851220961e+234;
bool r509760 = r509758 <= r509759;
double r509761 = 1.150967563391517e+227;
bool r509762 = r509758 <= r509761;
double r509763 = !r509762;
bool r509764 = r509760 || r509763;
double r509765 = t;
double r509766 = y;
double r509767 = x;
double r509768 = r509767 / r509758;
double r509769 = r509765 / r509758;
double r509770 = r509768 - r509769;
double r509771 = r509766 * r509770;
double r509772 = r509765 + r509771;
double r509773 = r509766 - r509758;
double r509774 = a;
double r509775 = r509774 - r509758;
double r509776 = cbrt(r509775);
double r509777 = r509776 * r509776;
double r509778 = r509773 / r509777;
double r509779 = r509765 - r509767;
double r509780 = r509779 / r509776;
double r509781 = r509778 * r509780;
double r509782 = r509767 + r509781;
double r509783 = r509764 ? r509772 : r509782;
return r509783;
}




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))))