\frac{x \cdot y - \left(z \cdot 9\right) \cdot t}{a \cdot 2}\begin{array}{l}
\mathbf{if}\;x \cdot y - \left(z \cdot 9\right) \cdot t \le -1.378384659417442190831979112913996094847 \cdot 10^{220}:\\
\;\;\;\;\frac{x}{\frac{a}{y}} \cdot 0.5 - \frac{z}{a} \cdot \left(t \cdot 4.5\right)\\
\mathbf{elif}\;x \cdot y - \left(z \cdot 9\right) \cdot t \le 1.623858966667658017012234116157301497207 \cdot 10^{298}:\\
\;\;\;\;\frac{1}{2 \cdot a} \cdot \left(x \cdot y - \left(z \cdot 9\right) \cdot t\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\frac{a}{y}} \cdot 0.5 - \left(\frac{z}{a} \cdot t\right) \cdot 4.5\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r34839664 = x;
double r34839665 = y;
double r34839666 = r34839664 * r34839665;
double r34839667 = z;
double r34839668 = 9.0;
double r34839669 = r34839667 * r34839668;
double r34839670 = t;
double r34839671 = r34839669 * r34839670;
double r34839672 = r34839666 - r34839671;
double r34839673 = a;
double r34839674 = 2.0;
double r34839675 = r34839673 * r34839674;
double r34839676 = r34839672 / r34839675;
return r34839676;
}
double f(double x, double y, double z, double t, double a) {
double r34839677 = x;
double r34839678 = y;
double r34839679 = r34839677 * r34839678;
double r34839680 = z;
double r34839681 = 9.0;
double r34839682 = r34839680 * r34839681;
double r34839683 = t;
double r34839684 = r34839682 * r34839683;
double r34839685 = r34839679 - r34839684;
double r34839686 = -1.3783846594174422e+220;
bool r34839687 = r34839685 <= r34839686;
double r34839688 = a;
double r34839689 = r34839688 / r34839678;
double r34839690 = r34839677 / r34839689;
double r34839691 = 0.5;
double r34839692 = r34839690 * r34839691;
double r34839693 = r34839680 / r34839688;
double r34839694 = 4.5;
double r34839695 = r34839683 * r34839694;
double r34839696 = r34839693 * r34839695;
double r34839697 = r34839692 - r34839696;
double r34839698 = 1.623858966667658e+298;
bool r34839699 = r34839685 <= r34839698;
double r34839700 = 1.0;
double r34839701 = 2.0;
double r34839702 = r34839701 * r34839688;
double r34839703 = r34839700 / r34839702;
double r34839704 = r34839703 * r34839685;
double r34839705 = r34839693 * r34839683;
double r34839706 = r34839705 * r34839694;
double r34839707 = r34839692 - r34839706;
double r34839708 = r34839699 ? r34839704 : r34839707;
double r34839709 = r34839687 ? r34839697 : r34839708;
return r34839709;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 7.4 |
|---|---|
| Target | 5.5 |
| Herbie | 0.9 |
if (- (* x y) (* (* z 9.0) t)) < -1.3783846594174422e+220Initial program 30.4
Taylor expanded around 0 30.2
rmApplied associate-/l*15.8
rmApplied *-un-lft-identity15.8
Applied times-frac0.9
Applied associate-*r*0.9
Simplified0.9
if -1.3783846594174422e+220 < (- (* x y) (* (* z 9.0) t)) < 1.623858966667658e+298Initial program 0.8
rmApplied div-inv0.9
if 1.623858966667658e+298 < (- (* x y) (* (* z 9.0) t)) Initial program 58.9
Taylor expanded around 0 58.4
rmApplied associate-/l*32.7
rmApplied *-un-lft-identity32.7
Applied times-frac0.3
Simplified0.3
Final simplification0.9
herbie shell --seed 2019170
(FPCore (x y z t a)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, I"
: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.0 t))) (* a 2.0)) (- (* (/ y a) (* x 0.5)) (* (/ t a) (* z 4.5)))))
(/ (- (* x y) (* (* z 9.0) t)) (* a 2.0)))