\frac{\left(\left(x \cdot 9\right) \cdot y - \left(\left(z \cdot 4\right) \cdot t\right) \cdot a\right) + b}{z \cdot c}\begin{array}{l}
\mathbf{if}\;\left(x \cdot 9\right) \cdot y \le -1.0598385219903945 \cdot 10^{153}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + 9 \cdot \frac{x}{z \cdot \frac{c}{y}}\right) - 4 \cdot \left(a \cdot \frac{t}{c}\right)\\
\mathbf{elif}\;\left(x \cdot 9\right) \cdot y \le -6.9631660866829206 \cdot 10^{62}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + \frac{9 \cdot \left(x \cdot y\right)}{z \cdot c}\right) - 4 \cdot \frac{a \cdot t}{c}\\
\mathbf{elif}\;\left(x \cdot 9\right) \cdot y \le -4.23671521665678558 \cdot 10^{-261}:\\
\;\;\;\;\left(\frac{\frac{b}{z}}{c} + 9 \cdot \frac{x}{\frac{z \cdot c}{y}}\right) - 4 \cdot \left(a \cdot \frac{t}{c}\right)\\
\mathbf{elif}\;\left(x \cdot 9\right) \cdot y \le 2.15413 \cdot 10^{-321}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + \left(9 \cdot \frac{x}{z}\right) \cdot \frac{y}{c}\right) - 4 \cdot \frac{a \cdot t}{c}\\
\mathbf{elif}\;\left(x \cdot 9\right) \cdot y \le 3.4351592935597035 \cdot 10^{-305}:\\
\;\;\;\;\frac{1}{z} \cdot \frac{\left(\left(x \cdot 9\right) \cdot y - \left(\left(z \cdot 4\right) \cdot t\right) \cdot a\right) + b}{c}\\
\mathbf{elif}\;\left(x \cdot 9\right) \cdot y \le 1.0650715768046093 \cdot 10^{-240}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + 9 \cdot \frac{x}{z \cdot \frac{c}{y}}\right) - 4 \cdot \left(a \cdot \frac{t}{c}\right)\\
\mathbf{elif}\;\left(x \cdot 9\right) \cdot y \le 7.54448186024096692 \cdot 10^{188}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + \frac{9 \cdot \left(x \cdot y\right)}{z \cdot c}\right) - 4 \cdot \frac{a \cdot t}{c}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + 9 \cdot \frac{x}{z \cdot \frac{c}{y}}\right) - 4 \cdot \left(a \cdot \frac{t}{c}\right)\\
\end{array}double f(double x, double y, double z, double t, double a, double b, double c) {
double r804500 = x;
double r804501 = 9.0;
double r804502 = r804500 * r804501;
double r804503 = y;
double r804504 = r804502 * r804503;
double r804505 = z;
double r804506 = 4.0;
double r804507 = r804505 * r804506;
double r804508 = t;
double r804509 = r804507 * r804508;
double r804510 = a;
double r804511 = r804509 * r804510;
double r804512 = r804504 - r804511;
double r804513 = b;
double r804514 = r804512 + r804513;
double r804515 = c;
double r804516 = r804505 * r804515;
double r804517 = r804514 / r804516;
return r804517;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r804518 = x;
double r804519 = 9.0;
double r804520 = r804518 * r804519;
double r804521 = y;
double r804522 = r804520 * r804521;
double r804523 = -1.0598385219903945e+153;
bool r804524 = r804522 <= r804523;
double r804525 = b;
double r804526 = z;
double r804527 = c;
double r804528 = r804526 * r804527;
double r804529 = r804525 / r804528;
double r804530 = r804527 / r804521;
double r804531 = r804526 * r804530;
double r804532 = r804518 / r804531;
double r804533 = r804519 * r804532;
double r804534 = r804529 + r804533;
double r804535 = 4.0;
double r804536 = a;
double r804537 = t;
double r804538 = r804537 / r804527;
double r804539 = r804536 * r804538;
double r804540 = r804535 * r804539;
double r804541 = r804534 - r804540;
double r804542 = -6.963166086682921e+62;
bool r804543 = r804522 <= r804542;
double r804544 = r804518 * r804521;
double r804545 = r804519 * r804544;
double r804546 = r804545 / r804528;
double r804547 = r804529 + r804546;
double r804548 = r804536 * r804537;
double r804549 = r804548 / r804527;
double r804550 = r804535 * r804549;
double r804551 = r804547 - r804550;
double r804552 = -4.236715216656786e-261;
bool r804553 = r804522 <= r804552;
double r804554 = r804525 / r804526;
double r804555 = r804554 / r804527;
double r804556 = r804528 / r804521;
double r804557 = r804518 / r804556;
double r804558 = r804519 * r804557;
double r804559 = r804555 + r804558;
double r804560 = r804559 - r804540;
double r804561 = 2.1541262158678e-321;
bool r804562 = r804522 <= r804561;
double r804563 = r804518 / r804526;
double r804564 = r804519 * r804563;
double r804565 = r804521 / r804527;
double r804566 = r804564 * r804565;
double r804567 = r804529 + r804566;
double r804568 = r804567 - r804550;
double r804569 = 3.4351592935597035e-305;
bool r804570 = r804522 <= r804569;
double r804571 = 1.0;
double r804572 = r804571 / r804526;
double r804573 = r804526 * r804535;
double r804574 = r804573 * r804537;
double r804575 = r804574 * r804536;
double r804576 = r804522 - r804575;
double r804577 = r804576 + r804525;
double r804578 = r804577 / r804527;
double r804579 = r804572 * r804578;
double r804580 = 1.0650715768046093e-240;
bool r804581 = r804522 <= r804580;
double r804582 = 7.544481860240967e+188;
bool r804583 = r804522 <= r804582;
double r804584 = r804583 ? r804551 : r804541;
double r804585 = r804581 ? r804541 : r804584;
double r804586 = r804570 ? r804579 : r804585;
double r804587 = r804562 ? r804568 : r804586;
double r804588 = r804553 ? r804560 : r804587;
double r804589 = r804543 ? r804551 : r804588;
double r804590 = r804524 ? r804541 : r804589;
return r804590;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a




Bits error versus b




Bits error versus c
Results
| Original | 20.9 |
|---|---|
| Target | 14.5 |
| Herbie | 8.9 |
if (* (* x 9.0) y) < -1.0598385219903945e+153 or 3.4351592935597035e-305 < (* (* x 9.0) y) < 1.0650715768046093e-240 or 7.544481860240967e+188 < (* (* x 9.0) y) Initial program 34.7
Taylor expanded around 0 27.0
rmApplied associate-/l*15.4
rmApplied *-un-lft-identity15.4
Applied times-frac13.5
Simplified13.5
rmApplied *-un-lft-identity13.5
Applied times-frac9.4
Simplified9.4
if -1.0598385219903945e+153 < (* (* x 9.0) y) < -6.963166086682921e+62 or 1.0650715768046093e-240 < (* (* x 9.0) y) < 7.544481860240967e+188Initial program 17.3
Taylor expanded around 0 7.2
rmApplied associate-*r/7.2
if -6.963166086682921e+62 < (* (* x 9.0) y) < -4.236715216656786e-261Initial program 16.0
Taylor expanded around 0 7.0
rmApplied associate-/l*9.9
rmApplied *-un-lft-identity9.9
Applied times-frac9.1
Simplified9.1
rmApplied associate-/r*11.3
if -4.236715216656786e-261 < (* (* x 9.0) y) < 2.1541262158678e-321Initial program 17.8
Taylor expanded around 0 8.1
rmApplied times-frac8.1
Applied associate-*r*8.1
if 2.1541262158678e-321 < (* (* x 9.0) y) < 3.4351592935597035e-305Initial program 16.6
rmApplied *-un-lft-identity16.6
Applied times-frac13.9
Final simplification8.9
herbie shell --seed 2020065
(FPCore (x y z t a b c)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, J"
:precision binary64
:herbie-target
(if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) -1.1001567408041051e-171) (/ (+ (- (* (* x 9) y) (* (* z 4) (* t a))) b) (* z c)) (if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) -0.0) (/ (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) z) c) (if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) 1.1708877911747488e-53) (/ (+ (- (* (* x 9) y) (* (* z 4) (* t a))) b) (* z c)) (if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) 2.876823679546137e+130) (- (+ (* (* 9 (/ y c)) (/ x z)) (/ b (* c z))) (* 4 (/ (* a t) c))) (if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) 1.3838515042456319e+158) (/ (+ (- (* (* x 9) y) (* (* z 4) (* t a))) b) (* z c)) (- (+ (* 9 (* (/ y (* c z)) x)) (/ b (* c z))) (* 4 (/ (* a t) c))))))))
(/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)))