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 -5.94158773471745396 \cdot 10^{-307} \lor \neg \left(x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z} \le 0.0\right):\\
\;\;\;\;x + \left(\frac{\sqrt[3]{y - z} \cdot \sqrt[3]{y - z}}{\sqrt[3]{a - z}} \cdot \frac{\sqrt[3]{y - z}}{\frac{\sqrt[3]{a - z}}{\frac{\sqrt[3]{t - x} \cdot {\left(\sqrt[3]{\sqrt[3]{t - x}}\right)}^{3}}{\sqrt[3]{\sqrt[3]{a - z} \cdot \sqrt[3]{a - z}}}}}\right) \cdot \frac{\sqrt[3]{t - x}}{\sqrt[3]{\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 r781204 = x;
double r781205 = y;
double r781206 = z;
double r781207 = r781205 - r781206;
double r781208 = t;
double r781209 = r781208 - r781204;
double r781210 = r781207 * r781209;
double r781211 = a;
double r781212 = r781211 - r781206;
double r781213 = r781210 / r781212;
double r781214 = r781204 + r781213;
return r781214;
}
double f(double x, double y, double z, double t, double a) {
double r781215 = x;
double r781216 = y;
double r781217 = z;
double r781218 = r781216 - r781217;
double r781219 = t;
double r781220 = r781219 - r781215;
double r781221 = r781218 * r781220;
double r781222 = a;
double r781223 = r781222 - r781217;
double r781224 = r781221 / r781223;
double r781225 = r781215 + r781224;
double r781226 = -5.941587734717454e-307;
bool r781227 = r781225 <= r781226;
double r781228 = 0.0;
bool r781229 = r781225 <= r781228;
double r781230 = !r781229;
bool r781231 = r781227 || r781230;
double r781232 = cbrt(r781218);
double r781233 = r781232 * r781232;
double r781234 = cbrt(r781223);
double r781235 = r781233 / r781234;
double r781236 = cbrt(r781220);
double r781237 = cbrt(r781236);
double r781238 = 3.0;
double r781239 = pow(r781237, r781238);
double r781240 = r781236 * r781239;
double r781241 = r781234 * r781234;
double r781242 = cbrt(r781241);
double r781243 = r781240 / r781242;
double r781244 = r781234 / r781243;
double r781245 = r781232 / r781244;
double r781246 = r781235 * r781245;
double r781247 = cbrt(r781234);
double r781248 = r781236 / r781247;
double r781249 = r781246 * r781248;
double r781250 = r781215 + r781249;
double r781251 = r781215 * r781216;
double r781252 = r781251 / r781217;
double r781253 = r781252 + r781219;
double r781254 = r781219 * r781216;
double r781255 = r781254 / r781217;
double r781256 = r781253 - r781255;
double r781257 = r781231 ? r781250 : r781256;
return r781257;
}




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 | 8.6 |
if (+ x (/ (* (- y z) (- t x)) (- a z))) < -5.941587734717454e-307 or 0.0 < (+ x (/ (* (- y z) (- t x)) (- a z))) Initial program 21.3
rmApplied add-cube-cbrt21.8
Applied times-frac8.2
rmApplied add-cube-cbrt8.3
Applied cbrt-prod8.3
Applied add-cube-cbrt8.5
Applied times-frac8.5
Applied associate-*r*8.0
rmApplied add-cube-cbrt8.1
rmApplied add-cube-cbrt8.1
Applied times-frac8.1
Applied associate-*l*7.7
Simplified7.8
if -5.941587734717454e-307 < (+ x (/ (* (- y z) (- t x)) (- a z))) < 0.0Initial program 61.2
Taylor expanded around inf 17.8
Final simplification8.6
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))))