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

Average Error: 5.1 → 1.0

Time: 23.6s

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\]

↓

\[\begin{array}{l} \mathbf{if}\;\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 = -\infty \lor \neg \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 \le 5.751099589137780109846545636083658425196 \cdot 10^{297}\right):\\ \;\;\;\;\left(\left(b \cdot c + \left(x \cdot \left(\left(18 \cdot y\right) \cdot \left(z \cdot t\right)\right) - \left(a \cdot 4\right) \cdot t\right)\right) - \left(x \cdot 4\right) \cdot i\right) - \left(j \cdot 27\right) \cdot k\\ \mathbf{else}:\\ \;\;\;\;\left(\left(\left(\left(\left(x \cdot \left(18 \cdot y\right)\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\\ \end{array}\]

C

double f(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r1113220 = x;
        double r1113221 = 18.0;
        double r1113222 = r1113220 * r1113221;
        double r1113223 = y;
        double r1113224 = r1113222 * r1113223;
        double r1113225 = z;
        double r1113226 = r1113224 * r1113225;
        double r1113227 = t;
        double r1113228 = r1113226 * r1113227;
        double r1113229 = a;
        double r1113230 = 4.0;
        double r1113231 = r1113229 * r1113230;
        double r1113232 = r1113231 * r1113227;
        double r1113233 = r1113228 - r1113232;
        double r1113234 = b;
        double r1113235 = c;
        double r1113236 = r1113234 * r1113235;
        double r1113237 = r1113233 + r1113236;
        double r1113238 = r1113220 * r1113230;
        double r1113239 = i;
        double r1113240 = r1113238 * r1113239;
        double r1113241 = r1113237 - r1113240;
        double r1113242 = j;
        double r1113243 = 27.0;
        double r1113244 = r1113242 * r1113243;
        double r1113245 = k;
        double r1113246 = r1113244 * r1113245;
        double r1113247 = r1113241 - r1113246;
        return r1113247;
}

↓

double f(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r1113248 = x;
        double r1113249 = 18.0;
        double r1113250 = r1113248 * r1113249;
        double r1113251 = y;
        double r1113252 = r1113250 * r1113251;
        double r1113253 = z;
        double r1113254 = r1113252 * r1113253;
        double r1113255 = t;
        double r1113256 = r1113254 * r1113255;
        double r1113257 = a;
        double r1113258 = 4.0;
        double r1113259 = r1113257 * r1113258;
        double r1113260 = r1113259 * r1113255;
        double r1113261 = r1113256 - r1113260;
        double r1113262 = b;
        double r1113263 = c;
        double r1113264 = r1113262 * r1113263;
        double r1113265 = r1113261 + r1113264;
        double r1113266 = r1113248 * r1113258;
        double r1113267 = i;
        double r1113268 = r1113266 * r1113267;
        double r1113269 = r1113265 - r1113268;
        double r1113270 = -inf.0;
        bool r1113271 = r1113269 <= r1113270;
        double r1113272 = 5.75109958913778e+297;
        bool r1113273 = r1113269 <= r1113272;
        double r1113274 = !r1113273;
        bool r1113275 = r1113271 || r1113274;
        double r1113276 = r1113249 * r1113251;
        double r1113277 = r1113253 * r1113255;
        double r1113278 = r1113276 * r1113277;
        double r1113279 = r1113248 * r1113278;
        double r1113280 = r1113279 - r1113260;
        double r1113281 = r1113264 + r1113280;
        double r1113282 = r1113281 - r1113268;
        double r1113283 = j;
        double r1113284 = 27.0;
        double r1113285 = r1113283 * r1113284;
        double r1113286 = k;
        double r1113287 = r1113285 * r1113286;
        double r1113288 = r1113282 - r1113287;
        double r1113289 = r1113248 * r1113276;
        double r1113290 = r1113289 * r1113253;
        double r1113291 = r1113290 * r1113255;
        double r1113292 = r1113291 - r1113260;
        double r1113293 = r1113292 + r1113264;
        double r1113294 = r1113293 - r1113268;
        double r1113295 = r1113294 - r1113287;
        double r1113296 = r1113275 ? r1113288 : r1113295;
        return r1113296;
}

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.1
Target	1.6
Herbie	1.0

\[\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 (- (+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c)) (* (* x 4.0) i)) < -inf.0 or 5.75109958913778e+297 < (- (+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c)) (* (* x 4.0) i))

   1. Initial program 53.9

      \[\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. Using strategy rm

   3. Applied pow1 53.9

      \[\leadsto \left(\left(\left(\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z\right) \cdot \color{blue}{{t}^{1}} - \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\]

   4. Applied pow1 53.9

      \[\leadsto \left(\left(\left(\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot \color{blue}{{z}^{1}}\right) \cdot {t}^{1} - \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\]

   5. Applied pow1 53.9

      \[\leadsto \left(\left(\left(\left(\left(\left(x \cdot 18\right) \cdot \color{blue}{{y}^{1}}\right) \cdot {z}^{1}\right) \cdot {t}^{1} - \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\]

   6. Applied pow1 53.9

      \[\leadsto \left(\left(\left(\left(\left(\left(x \cdot \color{blue}{{18}^{1}}\right) \cdot {y}^{1}\right) \cdot {z}^{1}\right) \cdot {t}^{1} - \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\]

   7. Applied pow1 53.9

      \[\leadsto \left(\left(\left(\left(\left(\left(\color{blue}{{x}^{1}} \cdot {18}^{1}\right) \cdot {y}^{1}\right) \cdot {z}^{1}\right) \cdot {t}^{1} - \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\]

   8. Applied pow-prod-down 53.9

      \[\leadsto \left(\left(\left(\left(\left(\color{blue}{{\left(x \cdot 18\right)}^{1}} \cdot {y}^{1}\right) \cdot {z}^{1}\right) \cdot {t}^{1} - \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\]

   9. Applied pow-prod-down 53.9

      \[\leadsto \left(\left(\left(\left(\color{blue}{{\left(\left(x \cdot 18\right) \cdot y\right)}^{1}} \cdot {z}^{1}\right) \cdot {t}^{1} - \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\]

   10. Applied pow-prod-down 53.9

       \[\leadsto \left(\left(\left(\color{blue}{{\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z\right)}^{1}} \cdot {t}^{1} - \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\]

   11. Applied pow-prod-down 53.9

       \[\leadsto \left(\left(\left(\color{blue}{{\left(\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z\right) \cdot t\right)}^{1}} - \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\]

   12. Simplified 7.1

       \[\leadsto \left(\left(\left({\color{blue}{\left(x \cdot \left(\left(18 \cdot y\right) \cdot \left(z \cdot t\right)\right)\right)}}^{1} - \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\]

   if -inf.0 < (- (+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c)) (* (* x 4.0) i)) < 5.75109958913778e+297

   1. Initial program 0.4

      \[\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. Using strategy rm

   3. Applied associate-*l* 0.4

      \[\leadsto \left(\left(\left(\left(\color{blue}{\left(x \cdot \left(18 \cdot y\right)\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\]
3. Recombined 2 regimes into one program.

4. Final simplification 1.0

   \[\leadsto \begin{array}{l} \mathbf{if}\;\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 = -\infty \lor \neg \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 \le 5.751099589137780109846545636083658425196 \cdot 10^{297}\right):\\ \;\;\;\;\left(\left(b \cdot c + \left(x \cdot \left(\left(18 \cdot y\right) \cdot \left(z \cdot t\right)\right) - \left(a \cdot 4\right) \cdot t\right)\right) - \left(x \cdot 4\right) \cdot i\right) - \left(j \cdot 27\right) \cdot k\\ \mathbf{else}:\\ \;\;\;\;\left(\left(\left(\left(\left(x \cdot \left(18 \cdot y\right)\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\\ \end{array}\]

Reproduce

herbie shell --seed 2019209 
(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)))