\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 -5.81880234558190669 \cdot 10^{134} \lor \neg \left(x \cdot y - \left(z \cdot 9\right) \cdot t \le 1.95624520253703639 \cdot 10^{220}\right):\\
\;\;\;\;\left(x \cdot 0.5\right) \cdot \frac{y}{a} - \left(t \cdot 4.5\right) \cdot \frac{z}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{a} \cdot \frac{x \cdot y - \left(z \cdot 9\right) \cdot t}{2}\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r684227 = x;
double r684228 = y;
double r684229 = r684227 * r684228;
double r684230 = z;
double r684231 = 9.0;
double r684232 = r684230 * r684231;
double r684233 = t;
double r684234 = r684232 * r684233;
double r684235 = r684229 - r684234;
double r684236 = a;
double r684237 = 2.0;
double r684238 = r684236 * r684237;
double r684239 = r684235 / r684238;
return r684239;
}
double f(double x, double y, double z, double t, double a) {
double r684240 = x;
double r684241 = y;
double r684242 = r684240 * r684241;
double r684243 = z;
double r684244 = 9.0;
double r684245 = r684243 * r684244;
double r684246 = t;
double r684247 = r684245 * r684246;
double r684248 = r684242 - r684247;
double r684249 = -5.818802345581907e+134;
bool r684250 = r684248 <= r684249;
double r684251 = 1.9562452025370364e+220;
bool r684252 = r684248 <= r684251;
double r684253 = !r684252;
bool r684254 = r684250 || r684253;
double r684255 = 0.5;
double r684256 = r684240 * r684255;
double r684257 = a;
double r684258 = r684241 / r684257;
double r684259 = r684256 * r684258;
double r684260 = 4.5;
double r684261 = r684246 * r684260;
double r684262 = r684243 / r684257;
double r684263 = r684261 * r684262;
double r684264 = r684259 - r684263;
double r684265 = 1.0;
double r684266 = r684265 / r684257;
double r684267 = 2.0;
double r684268 = r684248 / r684267;
double r684269 = r684266 * r684268;
double r684270 = r684254 ? r684264 : r684269;
return r684270;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 8.0 |
|---|---|
| Target | 5.6 |
| Herbie | 1.7 |
if (- (* x y) (* (* z 9.0) t)) < -5.818802345581907e+134 or 1.9562452025370364e+220 < (- (* x y) (* (* z 9.0) t)) Initial program 24.6
Taylor expanded around 0 24.2
rmApplied *-un-lft-identity24.2
Applied times-frac14.1
Applied associate-*r*14.2
Simplified14.2
rmApplied *-un-lft-identity14.2
Applied times-frac2.9
Applied associate-*r*2.8
Simplified2.8
if -5.818802345581907e+134 < (- (* x y) (* (* z 9.0) t)) < 1.9562452025370364e+220Initial program 1.1
rmApplied *-un-lft-identity1.1
Applied times-frac1.2
Final simplification1.7
herbie shell --seed 2020100
(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)))