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

Average Error: 5.5 → 1.6

Time: 12.4s

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 \le -2.03888388659653177 \cdot 10^{295}:\\ \;\;\;\;\left(\left(\left(\left(x \cdot 18\right) \cdot \left(y \cdot \left(z \cdot t\right)\right) - \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\\ \mathbf{elif}\;\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 2.661920863197471 \cdot 10^{300}:\\ \;\;\;\;\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) - j \cdot \left(27 \cdot k\right)\\ \mathbf{else}:\\ \;\;\;\;\left(\left(\left(\left(x \cdot 18\right) \cdot \left(\left(y \cdot z\right) \cdot t\right) - \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 r153948 = x;
        double r153949 = 18.0;
        double r153950 = r153948 * r153949;
        double r153951 = y;
        double r153952 = r153950 * r153951;
        double r153953 = z;
        double r153954 = r153952 * r153953;
        double r153955 = t;
        double r153956 = r153954 * r153955;
        double r153957 = a;
        double r153958 = 4.0;
        double r153959 = r153957 * r153958;
        double r153960 = r153959 * r153955;
        double r153961 = r153956 - r153960;
        double r153962 = b;
        double r153963 = c;
        double r153964 = r153962 * r153963;
        double r153965 = r153961 + r153964;
        double r153966 = r153948 * r153958;
        double r153967 = i;
        double r153968 = r153966 * r153967;
        double r153969 = r153965 - r153968;
        double r153970 = j;
        double r153971 = 27.0;
        double r153972 = r153970 * r153971;
        double r153973 = k;
        double r153974 = r153972 * r153973;
        double r153975 = r153969 - r153974;
        return r153975;
}

↓

double f(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
        double r153976 = x;
        double r153977 = 18.0;
        double r153978 = r153976 * r153977;
        double r153979 = y;
        double r153980 = r153978 * r153979;
        double r153981 = z;
        double r153982 = r153980 * r153981;
        double r153983 = t;
        double r153984 = r153982 * r153983;
        double r153985 = a;
        double r153986 = 4.0;
        double r153987 = r153985 * r153986;
        double r153988 = r153987 * r153983;
        double r153989 = r153984 - r153988;
        double r153990 = b;
        double r153991 = c;
        double r153992 = r153990 * r153991;
        double r153993 = r153989 + r153992;
        double r153994 = r153976 * r153986;
        double r153995 = i;
        double r153996 = r153994 * r153995;
        double r153997 = r153993 - r153996;
        double r153998 = -2.0388838865965318e+295;
        bool r153999 = r153997 <= r153998;
        double r154000 = r153981 * r153983;
        double r154001 = r153979 * r154000;
        double r154002 = r153978 * r154001;
        double r154003 = r154002 - r153988;
        double r154004 = r154003 + r153992;
        double r154005 = r154004 - r153996;
        double r154006 = j;
        double r154007 = 27.0;
        double r154008 = r154006 * r154007;
        double r154009 = k;
        double r154010 = r154008 * r154009;
        double r154011 = r154005 - r154010;
        double r154012 = 2.661920863197471e+300;
        bool r154013 = r153997 <= r154012;
        double r154014 = r154007 * r154009;
        double r154015 = r154006 * r154014;
        double r154016 = r153997 - r154015;
        double r154017 = r153979 * r153981;
        double r154018 = r154017 * r153983;
        double r154019 = r153978 * r154018;
        double r154020 = r154019 - r153988;
        double r154021 = r154020 + r153992;
        double r154022 = r154021 - r153996;
        double r154023 = r154022 - r154010;
        double r154024 = r154013 ? r154016 : r154023;
        double r154025 = r153999 ? r154011 : r154024;
        return r154025;
}

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

Derivation

1. Split input into 3 regimes

2. if (- (+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c)) (* (* x 4.0) i)) < -2.0388838865965318e+295

   1. Initial program 44.3

      \[\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* 28.1

      \[\leadsto \left(\left(\left(\color{blue}{\left(\left(x \cdot 18\right) \cdot \left(y \cdot z\right)\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\]
   4. Using strategy rm

   5. Applied associate-*l* 16.2

      \[\leadsto \left(\left(\left(\color{blue}{\left(x \cdot 18\right) \cdot \left(\left(y \cdot z\right) \cdot t\right)} - \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. Using strategy rm

   7. Applied associate-*l* 7.9

      \[\leadsto \left(\left(\left(\left(x \cdot 18\right) \cdot \color{blue}{\left(y \cdot \left(z \cdot t\right)\right)} - \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 -2.0388838865965318e+295 < (- (+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c)) (* (* x 4.0) i)) < 2.661920863197471e+300

   1. Initial program 0.3

      \[\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(\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) - \color{blue}{j \cdot \left(27 \cdot k\right)}\]

   if 2.661920863197471e+300 < (- (+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c)) (* (* x 4.0) i))

   1. Initial program 51.2

      \[\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* 33.3

      \[\leadsto \left(\left(\left(\color{blue}{\left(\left(x \cdot 18\right) \cdot \left(y \cdot z\right)\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\]
   4. Using strategy rm

   5. Applied associate-*l* 15.0

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

4. Final simplification 1.6

   \[\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 \le -2.03888388659653177 \cdot 10^{295}:\\ \;\;\;\;\left(\left(\left(\left(x \cdot 18\right) \cdot \left(y \cdot \left(z \cdot t\right)\right) - \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\\ \mathbf{elif}\;\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 2.661920863197471 \cdot 10^{300}:\\ \;\;\;\;\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) - j \cdot \left(27 \cdot k\right)\\ \mathbf{else}:\\ \;\;\;\;\left(\left(\left(\left(x \cdot 18\right) \cdot \left(\left(y \cdot z\right) \cdot t\right) - \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 2020056 +o rules:numerics
(FPCore (x y z t a b c i j k)
  :name "Diagrams.Solve.Polynomial:cubForm  from diagrams-solve-0.1"
  :precision binary64
  (- (- (+ (- (* (* (* (* x 18) y) z) t) (* (* a 4) t)) (* b c)) (* (* x 4) i)) (* (* j 27) k)))