Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, E

Average Error: 5.6 → 5.6

Time: 19.5s

Precision: 64

Math

\[\left(\left(\left(\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z\right) \cdot t - \left(a \cdot 4\right) \cdot t\right) + b \cdot c\right) - \left(x \cdot 4\right) \cdot i\right) - \left(j \cdot 27\right) \cdot k\]

↓

\[t \cdot \left(\left(x \cdot \left(18 \cdot y\right)\right) \cdot z - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \cdot i + j \cdot \left(27 \cdot k\right)\right)\right)\]

C

double f(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r620104 = x;
        double r620105 = 18.0;
        double r620106 = r620104 * r620105;
        double r620107 = y;
        double r620108 = r620106 * r620107;
        double r620109 = z;
        double r620110 = r620108 * r620109;
        double r620111 = t;
        double r620112 = r620110 * r620111;
        double r620113 = a;
        double r620114 = 4.0;
        double r620115 = r620113 * r620114;
        double r620116 = r620115 * r620111;
        double r620117 = r620112 - r620116;
        double r620118 = b;
        double r620119 = c;
        double r620120 = r620118 * r620119;
        double r620121 = r620117 + r620120;
        double r620122 = r620104 * r620114;
        double r620123 = i;
        double r620124 = r620122 * r620123;
        double r620125 = r620121 - r620124;
        double r620126 = j;
        double r620127 = 27.0;
        double r620128 = r620126 * r620127;
        double r620129 = k;
        double r620130 = r620128 * r620129;
        double r620131 = r620125 - r620130;
        return r620131;
}

↓

double f(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r620132 = t;
        double r620133 = x;
        double r620134 = 18.0;
        double r620135 = y;
        double r620136 = r620134 * r620135;
        double r620137 = r620133 * r620136;
        double r620138 = z;
        double r620139 = r620137 * r620138;
        double r620140 = a;
        double r620141 = 4.0;
        double r620142 = r620140 * r620141;
        double r620143 = r620139 - r620142;
        double r620144 = r620132 * r620143;
        double r620145 = b;
        double r620146 = c;
        double r620147 = r620145 * r620146;
        double r620148 = r620133 * r620141;
        double r620149 = i;
        double r620150 = r620148 * r620149;
        double r620151 = j;
        double r620152 = 27.0;
        double r620153 = k;
        double r620154 = r620152 * r620153;
        double r620155 = r620151 * r620154;
        double r620156 = r620150 + r620155;
        double r620157 = r620147 - r620156;
        double r620158 = r620144 + r620157;
        return r620158;
}

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

Bits error versus i

Bits error versus j

Bits error versus k

Target

Original	5.6
Target	1.7
Herbie	5.6

\[\begin{array}{l} \mathbf{if}\;t \lt -1.62108153975413982700795070153457058168 \cdot 10^{-69}:\\ \;\;\;\;\left(\left(18 \cdot t\right) \cdot \left(\left(x \cdot y\right) \cdot z\right) - \left(a \cdot t + i \cdot x\right) \cdot 4\right) - \left(\left(k \cdot j\right) \cdot 27 - c \cdot b\right)\\ \mathbf{elif}\;t \lt 165.6802794380522243500308832153677940369:\\ \;\;\;\;\left(\left(18 \cdot y\right) \cdot \left(x \cdot \left(z \cdot t\right)\right) - \left(a \cdot t + i \cdot x\right) \cdot 4\right) + \left(c \cdot b - 27 \cdot \left(k \cdot j\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\left(\left(18 \cdot t\right) \cdot \left(\left(x \cdot y\right) \cdot z\right) - \left(a \cdot t + i \cdot x\right) \cdot 4\right) - \left(\left(k \cdot j\right) \cdot 27 - c \cdot b\right)\\ \end{array}\]

Derivation

1. Split input into 2 regimes

2. if z < -4.053702415850972e+42 or 1.205035209333235e-88 < z

   1. Initial program 6.6

      \[\left(\left(\left(\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z\right) \cdot t - \left(a \cdot 4\right) \cdot t\right) + b \cdot c\right) - \left(x \cdot 4\right) \cdot i\right) - \left(j \cdot 27\right) \cdot k\]

   2. Simplified 6.6

      \[\leadsto \color{blue}{t \cdot \left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \cdot i + \left(j \cdot 27\right) \cdot k\right)\right)}\]
   3. Using strategy rm

   4. Applied associate-*l* 6.6

      \[\leadsto t \cdot \left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \cdot i + \color{blue}{j \cdot \left(27 \cdot k\right)}\right)\right)\]
   5. Using strategy rm

   6. Applied associate-*l* 6.6

      \[\leadsto t \cdot \left(\color{blue}{\left(x \cdot \left(18 \cdot y\right)\right)} \cdot z - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \cdot i + j \cdot \left(27 \cdot k\right)\right)\right)\]

   if -4.053702415850972e+42 < z < 1.205035209333235e-88

   1. Initial program 4.7

      \[\left(\left(\left(\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z\right) \cdot t - \left(a \cdot 4\right) \cdot t\right) + b \cdot c\right) - \left(x \cdot 4\right) \cdot i\right) - \left(j \cdot 27\right) \cdot k\]

   2. Simplified 4.7

      \[\leadsto \color{blue}{t \cdot \left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \cdot i + \left(j \cdot 27\right) \cdot k\right)\right)}\]
   3. Using strategy rm

   4. Applied associate-*l* 4.7

      \[\leadsto t \cdot \left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \cdot i + \color{blue}{j \cdot \left(27 \cdot k\right)}\right)\right)\]
   5. Using strategy rm

   6. Applied pow1 4.7

      \[\leadsto t \cdot \left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot \color{blue}{{z}^{1}} - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \cdot i + j \cdot \left(27 \cdot k\right)\right)\right)\]

   7. Applied pow1 4.7

      \[\leadsto t \cdot \left(\left(\left(x \cdot 18\right) \cdot \color{blue}{{y}^{1}}\right) \cdot {z}^{1} - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \cdot i + j \cdot \left(27 \cdot k\right)\right)\right)\]

   8. Applied pow1 4.7

      \[\leadsto t \cdot \left(\left(\left(x \cdot \color{blue}{{18}^{1}}\right) \cdot {y}^{1}\right) \cdot {z}^{1} - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \cdot i + j \cdot \left(27 \cdot k\right)\right)\right)\]

   9. Applied pow1 4.7

      \[\leadsto t \cdot \left(\left(\left(\color{blue}{{x}^{1}} \cdot {18}^{1}\right) \cdot {y}^{1}\right) \cdot {z}^{1} - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \cdot i + j \cdot \left(27 \cdot k\right)\right)\right)\]

   10. Applied pow-prod-down 4.7

       \[\leadsto t \cdot \left(\left(\color{blue}{{\left(x \cdot 18\right)}^{1}} \cdot {y}^{1}\right) \cdot {z}^{1} - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \cdot i + j \cdot \left(27 \cdot k\right)\right)\right)\]

   11. Applied pow-prod-down 4.7

       \[\leadsto t \cdot \left(\color{blue}{{\left(\left(x \cdot 18\right) \cdot y\right)}^{1}} \cdot {z}^{1} - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \cdot i + j \cdot \left(27 \cdot k\right)\right)\right)\]

   12. Applied pow-prod-down 4.7

       \[\leadsto t \cdot \left(\color{blue}{{\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z\right)}^{1}} - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \cdot i + j \cdot \left(27 \cdot k\right)\right)\right)\]

   13. Simplified 1.2

       \[\leadsto t \cdot \left({\color{blue}{\left(18 \cdot \left(x \cdot \left(z \cdot y\right)\right)\right)}}^{1} - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \cdot i + j \cdot \left(27 \cdot k\right)\right)\right)\]
3. Recombined 2 regimes into one program.

4. Final simplification 5.6

   \[\leadsto t \cdot \left(\left(x \cdot \left(18 \cdot y\right)\right) \cdot z - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \cdot i + j \cdot \left(27 \cdot k\right)\right)\right)\]

Reproduce

herbie shell --seed 2019297 
(FPCore (x y z t a b c i j k)
  :name "Diagrams.Solve.Polynomial:cubForm  from diagrams-solve-0.1, E"
  :precision binary64

  :herbie-target
  (if (< t -1.6210815397541398e-69) (- (- (* (* 18 t) (* (* x y) z)) (* (+ (* a t) (* i x)) 4)) (- (* (* k j) 27) (* c b))) (if (< t 165.680279438052224) (+ (- (* (* 18 y) (* x (* z t))) (* (+ (* a t) (* i x)) 4)) (- (* c b) (* 27 (* k j)))) (- (- (* (* 18 t) (* (* x y) z)) (* (+ (* a t) (* i x)) 4)) (- (* (* k j) 27) (* c b)))))

  (- (- (+ (- (* (* (* (* x 18) y) z) t) (* (* a 4) t)) (* b c)) (* (* x 4) i)) (* (* j 27) k)))