Average Error: 20.9 → 8.9
Time: 5.2s
Precision: 64
\[\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}\]
\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;
}

Error

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

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Target

Original20.9
Target14.5
Herbie8.9
\[\begin{array}{l} \mathbf{if}\;\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} \lt -1.10015674080410512 \cdot 10^{-171}:\\ \;\;\;\;\frac{\left(\left(x \cdot 9\right) \cdot y - \left(z \cdot 4\right) \cdot \left(t \cdot a\right)\right) + b}{z \cdot c}\\ \mathbf{elif}\;\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} \lt -0.0:\\ \;\;\;\;\frac{\frac{\left(\left(x \cdot 9\right) \cdot y - \left(\left(z \cdot 4\right) \cdot t\right) \cdot a\right) + b}{z}}{c}\\ \mathbf{elif}\;\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} \lt 1.17088779117474882 \cdot 10^{-53}:\\ \;\;\;\;\frac{\left(\left(x \cdot 9\right) \cdot y - \left(z \cdot 4\right) \cdot \left(t \cdot a\right)\right) + b}{z \cdot c}\\ \mathbf{elif}\;\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} \lt 2.8768236795461372 \cdot 10^{130}:\\ \;\;\;\;\left(\left(9 \cdot \frac{y}{c}\right) \cdot \frac{x}{z} + \frac{b}{c \cdot z}\right) - 4 \cdot \frac{a \cdot t}{c}\\ \mathbf{elif}\;\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} \lt 1.3838515042456319 \cdot 10^{158}:\\ \;\;\;\;\frac{\left(\left(x \cdot 9\right) \cdot y - \left(z \cdot 4\right) \cdot \left(t \cdot a\right)\right) + b}{z \cdot c}\\ \mathbf{else}:\\ \;\;\;\;\left(9 \cdot \left(\frac{y}{c \cdot z} \cdot x\right) + \frac{b}{c \cdot z}\right) - 4 \cdot \frac{a \cdot t}{c}\\ \end{array}\]

Derivation

  1. Split input into 5 regimes
  2. 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)

    1. Initial program 34.7

      \[\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}\]
    2. Taylor expanded around 0 27.0

      \[\leadsto \color{blue}{\left(\frac{b}{z \cdot c} + 9 \cdot \frac{x \cdot y}{z \cdot c}\right) - 4 \cdot \frac{a \cdot t}{c}}\]
    3. Using strategy rm
    4. Applied associate-/l*15.4

      \[\leadsto \left(\frac{b}{z \cdot c} + 9 \cdot \color{blue}{\frac{x}{\frac{z \cdot c}{y}}}\right) - 4 \cdot \frac{a \cdot t}{c}\]
    5. Using strategy rm
    6. Applied *-un-lft-identity15.4

      \[\leadsto \left(\frac{b}{z \cdot c} + 9 \cdot \frac{x}{\frac{z \cdot c}{y}}\right) - 4 \cdot \frac{a \cdot t}{\color{blue}{1 \cdot c}}\]
    7. Applied times-frac13.5

      \[\leadsto \left(\frac{b}{z \cdot c} + 9 \cdot \frac{x}{\frac{z \cdot c}{y}}\right) - 4 \cdot \color{blue}{\left(\frac{a}{1} \cdot \frac{t}{c}\right)}\]
    8. Simplified13.5

      \[\leadsto \left(\frac{b}{z \cdot c} + 9 \cdot \frac{x}{\frac{z \cdot c}{y}}\right) - 4 \cdot \left(\color{blue}{a} \cdot \frac{t}{c}\right)\]
    9. Using strategy rm
    10. Applied *-un-lft-identity13.5

      \[\leadsto \left(\frac{b}{z \cdot c} + 9 \cdot \frac{x}{\frac{z \cdot c}{\color{blue}{1 \cdot y}}}\right) - 4 \cdot \left(a \cdot \frac{t}{c}\right)\]
    11. Applied times-frac9.4

      \[\leadsto \left(\frac{b}{z \cdot c} + 9 \cdot \frac{x}{\color{blue}{\frac{z}{1} \cdot \frac{c}{y}}}\right) - 4 \cdot \left(a \cdot \frac{t}{c}\right)\]
    12. Simplified9.4

      \[\leadsto \left(\frac{b}{z \cdot c} + 9 \cdot \frac{x}{\color{blue}{z} \cdot \frac{c}{y}}\right) - 4 \cdot \left(a \cdot \frac{t}{c}\right)\]

    if -1.0598385219903945e+153 < (* (* x 9.0) y) < -6.963166086682921e+62 or 1.0650715768046093e-240 < (* (* x 9.0) y) < 7.544481860240967e+188

    1. Initial program 17.3

      \[\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}\]
    2. Taylor expanded around 0 7.2

      \[\leadsto \color{blue}{\left(\frac{b}{z \cdot c} + 9 \cdot \frac{x \cdot y}{z \cdot c}\right) - 4 \cdot \frac{a \cdot t}{c}}\]
    3. Using strategy rm
    4. Applied associate-*r/7.2

      \[\leadsto \left(\frac{b}{z \cdot c} + \color{blue}{\frac{9 \cdot \left(x \cdot y\right)}{z \cdot c}}\right) - 4 \cdot \frac{a \cdot t}{c}\]

    if -6.963166086682921e+62 < (* (* x 9.0) y) < -4.236715216656786e-261

    1. Initial program 16.0

      \[\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}\]
    2. Taylor expanded around 0 7.0

      \[\leadsto \color{blue}{\left(\frac{b}{z \cdot c} + 9 \cdot \frac{x \cdot y}{z \cdot c}\right) - 4 \cdot \frac{a \cdot t}{c}}\]
    3. Using strategy rm
    4. Applied associate-/l*9.9

      \[\leadsto \left(\frac{b}{z \cdot c} + 9 \cdot \color{blue}{\frac{x}{\frac{z \cdot c}{y}}}\right) - 4 \cdot \frac{a \cdot t}{c}\]
    5. Using strategy rm
    6. Applied *-un-lft-identity9.9

      \[\leadsto \left(\frac{b}{z \cdot c} + 9 \cdot \frac{x}{\frac{z \cdot c}{y}}\right) - 4 \cdot \frac{a \cdot t}{\color{blue}{1 \cdot c}}\]
    7. Applied times-frac9.1

      \[\leadsto \left(\frac{b}{z \cdot c} + 9 \cdot \frac{x}{\frac{z \cdot c}{y}}\right) - 4 \cdot \color{blue}{\left(\frac{a}{1} \cdot \frac{t}{c}\right)}\]
    8. Simplified9.1

      \[\leadsto \left(\frac{b}{z \cdot c} + 9 \cdot \frac{x}{\frac{z \cdot c}{y}}\right) - 4 \cdot \left(\color{blue}{a} \cdot \frac{t}{c}\right)\]
    9. Using strategy rm
    10. Applied associate-/r*11.3

      \[\leadsto \left(\color{blue}{\frac{\frac{b}{z}}{c}} + 9 \cdot \frac{x}{\frac{z \cdot c}{y}}\right) - 4 \cdot \left(a \cdot \frac{t}{c}\right)\]

    if -4.236715216656786e-261 < (* (* x 9.0) y) < 2.1541262158678e-321

    1. Initial program 17.8

      \[\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}\]
    2. Taylor expanded around 0 8.1

      \[\leadsto \color{blue}{\left(\frac{b}{z \cdot c} + 9 \cdot \frac{x \cdot y}{z \cdot c}\right) - 4 \cdot \frac{a \cdot t}{c}}\]
    3. Using strategy rm
    4. Applied times-frac8.1

      \[\leadsto \left(\frac{b}{z \cdot c} + 9 \cdot \color{blue}{\left(\frac{x}{z} \cdot \frac{y}{c}\right)}\right) - 4 \cdot \frac{a \cdot t}{c}\]
    5. Applied associate-*r*8.1

      \[\leadsto \left(\frac{b}{z \cdot c} + \color{blue}{\left(9 \cdot \frac{x}{z}\right) \cdot \frac{y}{c}}\right) - 4 \cdot \frac{a \cdot t}{c}\]

    if 2.1541262158678e-321 < (* (* x 9.0) y) < 3.4351592935597035e-305

    1. Initial program 16.6

      \[\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}\]
    2. Using strategy rm
    3. Applied *-un-lft-identity16.6

      \[\leadsto \frac{\color{blue}{1 \cdot \left(\left(\left(x \cdot 9\right) \cdot y - \left(\left(z \cdot 4\right) \cdot t\right) \cdot a\right) + b\right)}}{z \cdot c}\]
    4. Applied times-frac13.9

      \[\leadsto \color{blue}{\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}}\]
  3. Recombined 5 regimes into one program.
  4. Final simplification8.9

    \[\leadsto \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}\]

Reproduce

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)))