x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z}\begin{array}{l}
\mathbf{if}\;z \le -1.488227998811165 \cdot 10^{232}:\\
\;\;\;\;\left(\frac{x \cdot y}{z} + t\right) - \frac{t \cdot y}{z}\\
\mathbf{elif}\;z \le 2.17941341566278614 \cdot 10^{139}:\\
\;\;\;\;x + \frac{y - z}{\frac{a - z}{t - x}}\\
\mathbf{elif}\;z \le 3.37471355853467873 \cdot 10^{210}:\\
\;\;\;\;x + \frac{\frac{y - z}{\sqrt[3]{a - z}}}{\sqrt[3]{a - z}} \cdot \frac{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 r807104 = x;
double r807105 = y;
double r807106 = z;
double r807107 = r807105 - r807106;
double r807108 = t;
double r807109 = r807108 - r807104;
double r807110 = r807107 * r807109;
double r807111 = a;
double r807112 = r807111 - r807106;
double r807113 = r807110 / r807112;
double r807114 = r807104 + r807113;
return r807114;
}
double f(double x, double y, double z, double t, double a) {
double r807115 = z;
double r807116 = -1.4882279988111647e+232;
bool r807117 = r807115 <= r807116;
double r807118 = x;
double r807119 = y;
double r807120 = r807118 * r807119;
double r807121 = r807120 / r807115;
double r807122 = t;
double r807123 = r807121 + r807122;
double r807124 = r807122 * r807119;
double r807125 = r807124 / r807115;
double r807126 = r807123 - r807125;
double r807127 = 2.1794134156627861e+139;
bool r807128 = r807115 <= r807127;
double r807129 = r807119 - r807115;
double r807130 = a;
double r807131 = r807130 - r807115;
double r807132 = r807122 - r807118;
double r807133 = r807131 / r807132;
double r807134 = r807129 / r807133;
double r807135 = r807118 + r807134;
double r807136 = 3.374713558534679e+210;
bool r807137 = r807115 <= r807136;
double r807138 = cbrt(r807131);
double r807139 = r807129 / r807138;
double r807140 = r807139 / r807138;
double r807141 = r807132 / r807138;
double r807142 = r807140 * r807141;
double r807143 = r807118 + r807142;
double r807144 = r807137 ? r807143 : r807126;
double r807145 = r807128 ? r807135 : r807144;
double r807146 = r807117 ? r807126 : r807145;
return r807146;
}




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.7 |
| Herbie | 13.1 |
if z < -1.4882279988111647e+232 or 3.374713558534679e+210 < z Initial program 51.7
Taylor expanded around inf 24.2
if -1.4882279988111647e+232 < z < 2.1794134156627861e+139Initial program 17.5
rmApplied associate-/l*10.4
if 2.1794134156627861e+139 < z < 3.374713558534679e+210Initial program 41.9
rmApplied add-cube-cbrt42.2
Applied times-frac18.2
rmApplied *-un-lft-identity18.2
Applied cbrt-prod18.2
Applied *-un-lft-identity18.2
Applied times-frac18.2
Applied associate-*r*18.2
Simplified18.2
Final simplification13.1
herbie shell --seed 2020047
(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))))