Average Error: 3.7 → 0.7
Time: 4.0s
Precision: 64
\[\left(x \cdot 2 - \left(\left(y \cdot 9\right) \cdot z\right) \cdot t\right) + \left(a \cdot 27\right) \cdot b\]
\[\begin{array}{l} \mathbf{if}\;\left(y \cdot 9\right) \cdot z \le -1.9974616173648538 \cdot 10^{116} \lor \neg \left(\left(y \cdot 9\right) \cdot z \le 7.2920232258789333 \cdot 10^{208}\right):\\ \;\;\;\;\left(x \cdot 2 - \left(y \cdot 9\right) \cdot \left(z \cdot t\right)\right) + \left(a \cdot 27\right) \cdot b\\ \mathbf{else}:\\ \;\;\;\;x \cdot 2 + \left(27 \cdot \left(a \cdot b\right) - 9 \cdot \left(t \cdot \left(z \cdot y\right)\right)\right)\\ \end{array}\]
\left(x \cdot 2 - \left(\left(y \cdot 9\right) \cdot z\right) \cdot t\right) + \left(a \cdot 27\right) \cdot b
\begin{array}{l}
\mathbf{if}\;\left(y \cdot 9\right) \cdot z \le -1.9974616173648538 \cdot 10^{116} \lor \neg \left(\left(y \cdot 9\right) \cdot z \le 7.2920232258789333 \cdot 10^{208}\right):\\
\;\;\;\;\left(x \cdot 2 - \left(y \cdot 9\right) \cdot \left(z \cdot t\right)\right) + \left(a \cdot 27\right) \cdot b\\

\mathbf{else}:\\
\;\;\;\;x \cdot 2 + \left(27 \cdot \left(a \cdot b\right) - 9 \cdot \left(t \cdot \left(z \cdot y\right)\right)\right)\\

\end{array}
double f(double x, double y, double z, double t, double a, double b) {
        double r676146 = x;
        double r676147 = 2.0;
        double r676148 = r676146 * r676147;
        double r676149 = y;
        double r676150 = 9.0;
        double r676151 = r676149 * r676150;
        double r676152 = z;
        double r676153 = r676151 * r676152;
        double r676154 = t;
        double r676155 = r676153 * r676154;
        double r676156 = r676148 - r676155;
        double r676157 = a;
        double r676158 = 27.0;
        double r676159 = r676157 * r676158;
        double r676160 = b;
        double r676161 = r676159 * r676160;
        double r676162 = r676156 + r676161;
        return r676162;
}

double f(double x, double y, double z, double t, double a, double b) {
        double r676163 = y;
        double r676164 = 9.0;
        double r676165 = r676163 * r676164;
        double r676166 = z;
        double r676167 = r676165 * r676166;
        double r676168 = -1.9974616173648538e+116;
        bool r676169 = r676167 <= r676168;
        double r676170 = 7.292023225878933e+208;
        bool r676171 = r676167 <= r676170;
        double r676172 = !r676171;
        bool r676173 = r676169 || r676172;
        double r676174 = x;
        double r676175 = 2.0;
        double r676176 = r676174 * r676175;
        double r676177 = t;
        double r676178 = r676166 * r676177;
        double r676179 = r676165 * r676178;
        double r676180 = r676176 - r676179;
        double r676181 = a;
        double r676182 = 27.0;
        double r676183 = r676181 * r676182;
        double r676184 = b;
        double r676185 = r676183 * r676184;
        double r676186 = r676180 + r676185;
        double r676187 = r676181 * r676184;
        double r676188 = r676182 * r676187;
        double r676189 = r676166 * r676163;
        double r676190 = r676177 * r676189;
        double r676191 = r676164 * r676190;
        double r676192 = r676188 - r676191;
        double r676193 = r676176 + r676192;
        double r676194 = r676173 ? r676186 : r676193;
        return r676194;
}

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

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Target

Original3.7
Target2.6
Herbie0.7
\[\begin{array}{l} \mathbf{if}\;y \lt 7.590524218811189 \cdot 10^{-161}:\\ \;\;\;\;\left(x \cdot 2 - \left(\left(y \cdot 9\right) \cdot z\right) \cdot t\right) + a \cdot \left(27 \cdot b\right)\\ \mathbf{else}:\\ \;\;\;\;\left(x \cdot 2 - 9 \cdot \left(y \cdot \left(t \cdot z\right)\right)\right) + \left(a \cdot 27\right) \cdot b\\ \end{array}\]

Derivation

  1. Split input into 2 regimes
  2. if (* (* y 9.0) z) < -1.9974616173648538e+116 or 7.292023225878933e+208 < (* (* y 9.0) z)

    1. Initial program 20.1

      \[\left(x \cdot 2 - \left(\left(y \cdot 9\right) \cdot z\right) \cdot t\right) + \left(a \cdot 27\right) \cdot b\]
    2. Using strategy rm
    3. Applied associate-*l*2.3

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

    if -1.9974616173648538e+116 < (* (* y 9.0) z) < 7.292023225878933e+208

    1. Initial program 0.5

      \[\left(x \cdot 2 - \left(\left(y \cdot 9\right) \cdot z\right) \cdot t\right) + \left(a \cdot 27\right) \cdot b\]
    2. Using strategy rm
    3. Applied associate-*l*0.4

      \[\leadsto \left(x \cdot 2 - \left(\left(y \cdot 9\right) \cdot z\right) \cdot t\right) + \color{blue}{a \cdot \left(27 \cdot b\right)}\]
    4. Using strategy rm
    5. Applied sub-neg0.4

      \[\leadsto \color{blue}{\left(x \cdot 2 + \left(-\left(\left(y \cdot 9\right) \cdot z\right) \cdot t\right)\right)} + a \cdot \left(27 \cdot b\right)\]
    6. Applied associate-+l+0.4

      \[\leadsto \color{blue}{x \cdot 2 + \left(\left(-\left(\left(y \cdot 9\right) \cdot z\right) \cdot t\right) + a \cdot \left(27 \cdot b\right)\right)}\]
    7. Simplified0.4

      \[\leadsto x \cdot 2 + \color{blue}{\left(27 \cdot \left(a \cdot b\right) - \left(\left(y \cdot 9\right) \cdot z\right) \cdot t\right)}\]
    8. Taylor expanded around inf 0.4

      \[\leadsto x \cdot 2 + \color{blue}{\left(27 \cdot \left(a \cdot b\right) - 9 \cdot \left(t \cdot \left(z \cdot y\right)\right)\right)}\]
  3. Recombined 2 regimes into one program.
  4. Final simplification0.7

    \[\leadsto \begin{array}{l} \mathbf{if}\;\left(y \cdot 9\right) \cdot z \le -1.9974616173648538 \cdot 10^{116} \lor \neg \left(\left(y \cdot 9\right) \cdot z \le 7.2920232258789333 \cdot 10^{208}\right):\\ \;\;\;\;\left(x \cdot 2 - \left(y \cdot 9\right) \cdot \left(z \cdot t\right)\right) + \left(a \cdot 27\right) \cdot b\\ \mathbf{else}:\\ \;\;\;\;x \cdot 2 + \left(27 \cdot \left(a \cdot b\right) - 9 \cdot \left(t \cdot \left(z \cdot y\right)\right)\right)\\ \end{array}\]

Reproduce

herbie shell --seed 2020083 
(FPCore (x y z t a b)
  :name "Diagrams.Solve.Polynomial:cubForm  from diagrams-solve-0.1, A"
  :precision binary64

  :herbie-target
  (if (< y 7.590524218811189e-161) (+ (- (* x 2) (* (* (* y 9) z) t)) (* a (* 27 b))) (+ (- (* x 2) (* 9 (* y (* t z)))) (* (* a 27) b)))

  (+ (- (* x 2) (* (* (* y 9) z) t)) (* (* a 27) b)))