\frac{x \cdot y - \left(z \cdot 9\right) \cdot t}{a \cdot 2}\begin{array}{l}
\mathbf{if}\;a \le -1.228510223266635066615353468569416308744 \cdot 10^{-99}:\\
\;\;\;\;0.5 \cdot \left(x \cdot \frac{y}{a}\right) - 4.5 \cdot \frac{t}{\frac{a}{z}}\\
\mathbf{elif}\;a \le 64104884799497751389513187328:\\
\;\;\;\;\frac{0.5 \cdot \left(x \cdot y\right) - 4.5 \cdot \left(t \cdot z\right)}{a}\\
\mathbf{else}:\\
\;\;\;\;0.5 \cdot \left(\frac{x}{\sqrt[3]{a} \cdot \sqrt[3]{a}} \cdot \frac{y}{\sqrt[3]{a}}\right) - 4.5 \cdot \frac{t}{\frac{a}{z}}\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r487268 = x;
double r487269 = y;
double r487270 = r487268 * r487269;
double r487271 = z;
double r487272 = 9.0;
double r487273 = r487271 * r487272;
double r487274 = t;
double r487275 = r487273 * r487274;
double r487276 = r487270 - r487275;
double r487277 = a;
double r487278 = 2.0;
double r487279 = r487277 * r487278;
double r487280 = r487276 / r487279;
return r487280;
}
double f(double x, double y, double z, double t, double a) {
double r487281 = a;
double r487282 = -1.228510223266635e-99;
bool r487283 = r487281 <= r487282;
double r487284 = 0.5;
double r487285 = x;
double r487286 = y;
double r487287 = r487286 / r487281;
double r487288 = r487285 * r487287;
double r487289 = r487284 * r487288;
double r487290 = 4.5;
double r487291 = t;
double r487292 = z;
double r487293 = r487281 / r487292;
double r487294 = r487291 / r487293;
double r487295 = r487290 * r487294;
double r487296 = r487289 - r487295;
double r487297 = 6.410488479949775e+28;
bool r487298 = r487281 <= r487297;
double r487299 = r487285 * r487286;
double r487300 = r487284 * r487299;
double r487301 = r487291 * r487292;
double r487302 = r487290 * r487301;
double r487303 = r487300 - r487302;
double r487304 = r487303 / r487281;
double r487305 = cbrt(r487281);
double r487306 = r487305 * r487305;
double r487307 = r487285 / r487306;
double r487308 = r487286 / r487305;
double r487309 = r487307 * r487308;
double r487310 = r487284 * r487309;
double r487311 = r487310 - r487295;
double r487312 = r487298 ? r487304 : r487311;
double r487313 = r487283 ? r487296 : r487312;
return r487313;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 7.7 |
|---|---|
| Target | 5.2 |
| Herbie | 4.2 |
if a < -1.228510223266635e-99Initial program 9.0
Taylor expanded around 0 8.9
rmApplied associate-/l*7.4
rmApplied *-un-lft-identity7.4
Applied times-frac5.2
Simplified5.2
if -1.228510223266635e-99 < a < 6.410488479949775e+28Initial program 1.5
Taylor expanded around 0 1.6
rmApplied associate-*r/1.6
Applied associate-*r/1.5
Applied sub-div1.5
if 6.410488479949775e+28 < a Initial program 12.3
Taylor expanded around 0 12.2
rmApplied associate-/l*9.7
rmApplied add-cube-cbrt10.0
Applied times-frac5.6
Final simplification4.2
herbie shell --seed 2019326 +o rules:numerics
(FPCore (x y z t a)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, I"
:precision binary64
:herbie-target
(if (< a -2.090464557976709e+86) (- (* 0.5 (/ (* y x) a)) (* 4.5 (/ t (/ a z)))) (if (< a 2.144030707833976e+99) (/ (- (* x y) (* z (* 9 t))) (* a 2)) (- (* (/ y a) (* x 0.5)) (* (/ t a) (* z 4.5)))))
(/ (- (* x y) (* (* z 9) t)) (* a 2)))