Average Error: 27.2 → 27.7
Time: 50.0s
Precision: binary64
\[\]
\[\]
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k, double y0, double y1, double y2, double y3, double y4, double y5) {
	return ((double) (((double) (((double) (((double) (((double) (((double) (((double) (((double) (x * y)) - ((double) (z * t)))) * ((double) (((double) (a * b)) - ((double) (c * i)))))) - ((double) (((double) (((double) (x * j)) - ((double) (z * k)))) * ((double) (((double) (y0 * b)) - ((double) (y1 * i)))))))) + ((double) (((double) (((double) (x * y2)) - ((double) (z * y3)))) * ((double) (((double) (y0 * c)) - ((double) (y1 * a)))))))) + ((double) (((double) (((double) (t * j)) - ((double) (y * k)))) * ((double) (((double) (y4 * b)) - ((double) (y5 * i)))))))) - ((double) (((double) (((double) (t * y2)) - ((double) (y * y3)))) * ((double) (((double) (y4 * c)) - ((double) (y5 * a)))))))) + ((double) (((double) (((double) (k * y2)) - ((double) (j * y3)))) * ((double) (((double) (y4 * y1)) - ((double) (y5 * y0))))))));
}
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k, double y0, double y1, double y2, double y3, double y4, double y5) {
	double VAR;
	if ((y5 <= -4.808809147439703e+18)) {
		VAR = ((double) (((double) (((double) (((double) (((double) (((double) (((double) (b * ((double) (a * ((double) (((double) (x * y)) - ((double) (z * t)))))))) + ((double) (c * ((double) (i * ((double) (((double) (z * t)) - ((double) (x * y)))))))))) + ((double) (((double) (((double) (x * j)) - ((double) (z * k)))) * ((double) (((double) (i * y1)) - ((double) (b * y0)))))))) + ((double) (((double) (((double) (x * y2)) - ((double) (z * y3)))) * ((double) (((double) (c * y0)) - ((double) (a * y1)))))))) + ((double) (((double) (((double) (t * j)) - ((double) (y * k)))) * ((double) (((double) (b * y4)) - ((double) (y5 * i)))))))) + ((double) (((double) (y5 * ((double) (a * ((double) (((double) (t * y2)) - ((double) (y * y3)))))))) + ((double) (y4 * ((double) (c * ((double) (((double) (y * y3)) - ((double) (t * y2)))))))))))) + ((double) (((double) (((double) (k * y2)) - ((double) (j * y3)))) * ((double) (((double) (y1 * y4)) - ((double) (y5 * y0))))))));
	} else {
		double VAR_1;
		if ((y5 <= -1.9307705767884762e-78)) {
			VAR_1 = ((double) (((double) (((double) (((double) (k * y2)) - ((double) (j * y3)))) * ((double) (((double) (y1 * y4)) - ((double) (y5 * y0)))))) + ((double) (((double) (((double) (((double) (((double) (t * j)) - ((double) (y * k)))) * ((double) (((double) (b * y4)) - ((double) (y5 * i)))))) + ((double) (((double) (((double) (((double) (x * y2)) - ((double) (z * y3)))) * ((double) (((double) (c * y0)) - ((double) (a * y1)))))) + ((double) (((double) (((double) (((double) (x * y)) - ((double) (z * t)))) * ((double) (((double) (b * a)) - ((double) (c * i)))))) + ((double) (((double) (((double) (j * ((double) (i * ((double) (x * y1)))))) + ((double) (z * ((double) (y0 * ((double) (b * k)))))))) - ((double) (k * ((double) (z * ((double) (i * y1)))))))))))))) + ((double) (((double) (((double) (t * y2)) - ((double) (y * y3)))) * ((double) (((double) (y5 * a)) - ((double) (c * y4))))))))));
		} else {
			double VAR_2;
			if ((y5 <= -3.4297215494461836e-116)) {
				VAR_2 = ((double) (((double) (((double) (((double) (k * y2)) - ((double) (j * y3)))) * ((double) (((double) (y1 * y4)) - ((double) (y5 * y0)))))) + ((double) (((double) (((double) (((double) (((double) (t * j)) - ((double) (y * k)))) * ((double) (((double) (b * y4)) - ((double) (y5 * i)))))) + ((double) (((double) (((double) (b * ((double) (a * ((double) (((double) (x * y)) - ((double) (z * t)))))))) + ((double) (c * ((double) (i * ((double) (((double) (z * t)) - ((double) (x * y)))))))))) + ((double) (((double) (((double) (x * y2)) - ((double) (z * y3)))) * ((double) (((double) (c * y0)) - ((double) (a * y1)))))))))) + ((double) (((double) (((double) (t * y2)) - ((double) (y * y3)))) * ((double) (((double) (y5 * a)) - ((double) (c * y4))))))))));
			} else {
				double VAR_3;
				if ((y5 <= -1.013671810830034e-187)) {
					VAR_3 = ((double) (((double) (((double) (((double) (k * y2)) - ((double) (j * y3)))) * ((double) (((double) (y1 * y4)) - ((double) (y5 * y0)))))) + ((double) (((double) (((double) (((double) (((double) (((double) (x * y2)) - ((double) (z * y3)))) * ((double) (((double) (c * y0)) - ((double) (a * y1)))))) + ((double) (((double) (((double) (((double) (x * y)) - ((double) (z * t)))) * ((double) (((double) (b * a)) - ((double) (c * i)))))) + ((double) (((double) (((double) (x * j)) - ((double) (z * k)))) * ((double) (((double) (i * y1)) - ((double) (b * y0)))))))))) + ((double) (((double) (k * ((double) (i * ((double) (y5 * y)))))) - ((double) (((double) (i * ((double) (t * ((double) (y5 * j)))))) + ((double) (k * ((double) (b * ((double) (y * y4)))))))))))) + ((double) (((double) (((double) (t * y2)) - ((double) (y * y3)))) * ((double) (((double) (y5 * a)) - ((double) (c * y4))))))))));
				} else {
					double VAR_4;
					if ((y5 <= -2.7574949932241726e-232)) {
						VAR_4 = ((double) (((double) (((double) (((double) (((double) (((double) (((double) (b * ((double) (a * ((double) (((double) (x * y)) - ((double) (z * t)))))))) + ((double) (c * ((double) (i * ((double) (((double) (z * t)) - ((double) (x * y)))))))))) + ((double) (((double) (((double) (x * j)) - ((double) (z * k)))) * ((double) (((double) (i * y1)) - ((double) (b * y0)))))))) + ((double) (((double) (((double) (x * y2)) - ((double) (z * y3)))) * ((double) (((double) (c * y0)) - ((double) (a * y1)))))))) + ((double) (((double) (((double) (t * j)) - ((double) (y * k)))) * ((double) (((double) (b * y4)) - ((double) (y5 * i)))))))) + ((double) (((double) cbrt(((double) (((double) (c * y4)) - ((double) (y5 * a)))))) * ((double) (((double) (((double) cbrt(((double) (((double) (c * y4)) - ((double) (y5 * a)))))) * ((double) cbrt(((double) (((double) (c * y4)) - ((double) (y5 * a)))))))) * ((double) (((double) (y * y3)) - ((double) (t * y2)))))))))) + ((double) (((double) (y1 * ((double) (y4 * ((double) (((double) (k * y2)) - ((double) (j * y3)))))))) + ((double) (y5 * ((double) (y0 * ((double) (((double) (j * y3)) - ((double) (k * y2))))))))))));
					} else {
						double VAR_5;
						if ((y5 <= 6.874962237785278e-283)) {
							VAR_5 = ((double) (((double) (((double) (((double) (k * y2)) - ((double) (j * y3)))) * ((double) (((double) (y1 * y4)) - ((double) (y5 * y0)))))) + ((double) (((double) (((double) (((double) (((double) (((double) (x * y2)) - ((double) (z * y3)))) * ((double) (((double) (c * y0)) - ((double) (a * y1)))))) + ((double) (((double) (((double) (((double) (x * y)) - ((double) (z * t)))) * ((double) (((double) (b * a)) - ((double) (c * i)))))) + ((double) (((double) (((double) (x * j)) - ((double) (z * k)))) * ((double) (((double) (i * y1)) - ((double) (b * y0)))))))))) + ((double) (((double) (k * ((double) (i * ((double) (y5 * y)))))) - ((double) (((double) (i * ((double) (t * ((double) (y5 * j)))))) + ((double) (k * ((double) (b * ((double) (y * y4)))))))))))) + ((double) (((double) (((double) (t * y2)) - ((double) (y * y3)))) * ((double) (((double) (y5 * a)) - ((double) (c * y4))))))))));
						} else {
							double VAR_6;
							if ((y5 <= 1.194052064162617e-181)) {
								VAR_6 = ((double) (((double) (((double) (((double) (k * y2)) - ((double) (j * y3)))) * ((double) (((double) (y1 * y4)) - ((double) (y5 * y0)))))) + ((double) (((double) (((double) (((double) (((double) (t * j)) - ((double) (y * k)))) * ((double) (((double) (b * y4)) - ((double) (y5 * i)))))) + ((double) (((double) (((double) (((double) (((double) (x * y)) - ((double) (z * t)))) * ((double) (((double) (b * a)) - ((double) (c * i)))))) + ((double) (((double) (((double) (x * j)) - ((double) (z * k)))) * ((double) (((double) (i * y1)) - ((double) (b * y0)))))))) + ((double) (((double) (y0 * ((double) (c * ((double) (((double) (x * y2)) - ((double) (z * y3)))))))) + ((double) (y1 * ((double) (a * ((double) (((double) (z * y3)) - ((double) (x * y2)))))))))))))) + ((double) (((double) (((double) (t * y2)) - ((double) (y * y3)))) * ((double) (((double) (y5 * a)) - ((double) (c * y4))))))))));
							} else {
								double VAR_7;
								if ((y5 <= 4.810834636764084e-104)) {
									VAR_7 = ((double) (((double) (((double) (((double) (k * y2)) - ((double) (j * y3)))) * ((double) (((double) (y1 * y4)) - ((double) (y5 * y0)))))) + ((double) (((double) (((double) (((double) (((double) (b * ((double) (a * ((double) (((double) (x * y)) - ((double) (z * t)))))))) + ((double) (c * ((double) (i * ((double) (((double) (z * t)) - ((double) (x * y)))))))))) + ((double) (((double) (((double) (x * j)) - ((double) (z * k)))) * ((double) (((double) (i * y1)) - ((double) (b * y0)))))))) + ((double) (((double) (((double) (x * y2)) - ((double) (z * y3)))) * ((double) (((double) (c * y0)) - ((double) (a * y1)))))))) + ((double) (((double) cbrt(((double) (((double) (c * y4)) - ((double) (y5 * a)))))) * ((double) (((double) (((double) cbrt(((double) (((double) (c * y4)) - ((double) (y5 * a)))))) * ((double) cbrt(((double) (((double) (c * y4)) - ((double) (y5 * a)))))))) * ((double) (((double) (y * y3)) - ((double) (t * y2))))))))))));
								} else {
									VAR_7 = ((double) (((double) (((double) (((double) (((double) (((double) (((double) (b * ((double) (a * ((double) (((double) (x * y)) - ((double) (z * t)))))))) + ((double) (c * ((double) (i * ((double) (((double) (z * t)) - ((double) (x * y)))))))))) + ((double) (((double) (((double) (x * j)) - ((double) (z * k)))) * ((double) (((double) (i * y1)) - ((double) (b * y0)))))))) + ((double) (((double) (((double) (x * y2)) - ((double) (z * y3)))) * ((double) (((double) (c * y0)) - ((double) (a * y1)))))))) + ((double) (((double) (((double) (t * j)) - ((double) (y * k)))) * ((double) (((double) (b * y4)) - ((double) (y5 * i)))))))) + ((double) (((double) cbrt(((double) (((double) (c * y4)) - ((double) (y5 * a)))))) * ((double) (((double) (((double) cbrt(((double) (((double) (c * y4)) - ((double) (y5 * a)))))) * ((double) cbrt(((double) (((double) (c * y4)) - ((double) (y5 * a)))))))) * ((double) (((double) (y * y3)) - ((double) (t * y2)))))))))) + ((double) (((double) (y1 * ((double) (y4 * ((double) (((double) (k * y2)) - ((double) (j * y3)))))))) + ((double) (y5 * ((double) (y0 * ((double) (((double) (j * y3)) - ((double) (k * y2))))))))))));
								}
								VAR_6 = VAR_7;
							}
							VAR_5 = VAR_6;
						}
						VAR_4 = VAR_5;
					}
					VAR_3 = VAR_4;
				}
				VAR_2 = VAR_3;
			}
			VAR_1 = VAR_2;
		}
		VAR = VAR_1;
	}
	return VAR;
}

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

Bits error versus y0

Bits error versus y1

Bits error versus y2

Bits error versus y3

Bits error versus y4

Bits error versus y5

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Target

Original27.2
Target30.9
Herbie27.7
\[\]

Derivation

  1. Split input into 7 regimes
  2. if y5 < -4808809147439703040

    1. Initial program 29.8

      \[\]
    2. Using strategy rm
    3. Applied sub-neg29.8

      \[\leadsto \]
    4. Applied distribute-lft-in29.8

      \[\leadsto \]
    5. Simplified29.6

      \[\leadsto \]
    6. Simplified29.4

      \[\leadsto \]
    7. Using strategy rm
    8. Applied sub-neg29.4

      \[\leadsto \]
    9. Applied distribute-lft-in29.4

      \[\leadsto \]
    10. Simplified29.3

      \[\leadsto \]
    11. Simplified27.2

      \[\leadsto \]

    if -4808809147439703040 < y5 < -1.93077057678847617e-78

    1. Initial program 24.9

      \[\]
    2. Taylor expanded around inf 27.4

      \[\leadsto \]
    3. Simplified28.1

      \[\leadsto \]

    if -1.93077057678847617e-78 < y5 < -3.42972154944618364e-116

    1. Initial program 29.9

      \[\]
    2. Using strategy rm
    3. Applied sub-neg29.9

      \[\leadsto \]
    4. Applied distribute-lft-in29.9

      \[\leadsto \]
    5. Simplified30.7

      \[\leadsto \]
    6. Simplified31.0

      \[\leadsto \]
    7. Taylor expanded around 0 34.6

      \[\leadsto \]

    if -3.42972154944618364e-116 < y5 < -1.013671810830034e-187 or -2.7574949932241726e-232 < y5 < 6.87496223778527841e-283

    1. Initial program 26.2

      \[\]
    2. Taylor expanded around inf 27.8

      \[\leadsto \]
    3. Simplified27.4

      \[\leadsto \]

    if -1.013671810830034e-187 < y5 < -2.7574949932241726e-232 or 4.8108346367640839e-104 < y5

    1. Initial program 27.2

      \[\]
    2. Using strategy rm
    3. Applied sub-neg27.2

      \[\leadsto \]
    4. Applied distribute-lft-in27.2

      \[\leadsto \]
    5. Simplified27.3

      \[\leadsto \]
    6. Simplified27.4

      \[\leadsto \]
    7. Using strategy rm
    8. Applied add-cube-cbrt27.5

      \[\leadsto \]
    9. Applied associate-*r*27.5

      \[\leadsto \]
    10. Using strategy rm
    11. Applied sub-neg27.5

      \[\leadsto \]
    12. Applied distribute-lft-in27.5

      \[\leadsto \]
    13. Simplified27.8

      \[\leadsto \]
    14. Simplified26.9

      \[\leadsto \]

    if 6.87496223778527841e-283 < y5 < 1.194052064162617e-181

    1. Initial program 26.9

      \[\]
    2. Using strategy rm
    3. Applied sub-neg26.9

      \[\leadsto \]
    4. Applied distribute-lft-in26.9

      \[\leadsto \]
    5. Simplified27.1

      \[\leadsto \]
    6. Simplified26.6

      \[\leadsto \]

    if 1.194052064162617e-181 < y5 < 4.8108346367640839e-104

    1. Initial program 26.0

      \[\]
    2. Using strategy rm
    3. Applied sub-neg26.0

      \[\leadsto \]
    4. Applied distribute-lft-in26.0

      \[\leadsto \]
    5. Simplified26.9

      \[\leadsto \]
    6. Simplified27.0

      \[\leadsto \]
    7. Using strategy rm
    8. Applied add-cube-cbrt27.1

      \[\leadsto \]
    9. Applied associate-*r*27.1

      \[\leadsto \]
    10. Taylor expanded around 0 30.1

      \[\leadsto \]
  3. Recombined 7 regimes into one program.
  4. Final simplification27.7

    \[\leadsto \]

Reproduce

herbie shell --seed 2020192 
(FPCore (x y z t a b c i j k y0 y1 y2 y3 y4 y5)
  :name "Linear.Matrix:det44 from linear-1.19.1.3"
  :precision binary64

  :herbie-target
  (if (< y4 -7.206256231996481e+60) (- (- (* (- (* b a) (* i c)) (- (* y x) (* t z))) (- (* (- (* j x) (* k z)) (- (* y0 b) (* i y1))) (* (- (* j t) (* k y)) (- (* y4 b) (* y5 i))))) (- (/ (- (* y2 t) (* y3 y)) (/ 1.0 (- (* y4 c) (* y5 a)))) (* (- (* y2 k) (* y3 j)) (- (* y4 y1) (* y5 y0))))) (if (< y4 -3.364603505246317e-66) (+ (- (- (- (* (* t c) (* i z)) (* (* a t) (* b z))) (* (* y c) (* i x))) (* (- (* b y0) (* i y1)) (- (* j x) (* k z)))) (- (* (- (* y0 c) (* a y1)) (- (* x y2) (* z y3))) (- (* (- (* t y2) (* y y3)) (- (* y4 c) (* a y5))) (* (- (* y1 y4) (* y5 y0)) (- (* k y2) (* j y3)))))) (if (< y4 -1.2000065055686116e-105) (+ (+ (- (* (- (* j t) (* k y)) (- (* y4 b) (* y5 i))) (* (* y3 y) (- (* y5 a) (* y4 c)))) (+ (* (* y5 a) (* t y2)) (* (- (* k y2) (* j y3)) (- (* y4 y1) (* y5 y0))))) (- (* (- (* x y2) (* z y3)) (- (* c y0) (* a y1))) (- (* (- (* b y0) (* i y1)) (- (* j x) (* k z))) (* (- (* y x) (* z t)) (- (* b a) (* i c)))))) (if (< y4 6.718963124057495e-279) (+ (- (- (- (* (* k y) (* y5 i)) (* (* y b) (* y4 k))) (* (* y5 t) (* i j))) (- (* (- (* y2 t) (* y3 y)) (- (* y4 c) (* y5 a))) (* (- (* y2 k) (* y3 j)) (- (* y4 y1) (* y5 y0))))) (- (* (- (* b a) (* i c)) (- (* y x) (* t z))) (- (* (- (* j x) (* k z)) (- (* y0 b) (* i y1))) (* (- (* y2 x) (* y3 z)) (- (* c y0) (* y1 a)))))) (if (< y4 4.77962681403792e-222) (+ (+ (- (* (- (* j t) (* k y)) (- (* y4 b) (* y5 i))) (* (* y3 y) (- (* y5 a) (* y4 c)))) (+ (* (* y5 a) (* t y2)) (* (- (* k y2) (* j y3)) (- (* y4 y1) (* y5 y0))))) (- (* (- (* x y2) (* z y3)) (- (* c y0) (* a y1))) (- (* (- (* b y0) (* i y1)) (- (* j x) (* k z))) (* (- (* y x) (* z t)) (- (* b a) (* i c)))))) (if (< y4 2.2852241541266835e-175) (+ (- (- (- (* (* k y) (* y5 i)) (* (* y b) (* y4 k))) (* (* y5 t) (* i j))) (- (* (- (* y2 t) (* y3 y)) (- (* y4 c) (* y5 a))) (* (- (* y2 k) (* y3 j)) (- (* y4 y1) (* y5 y0))))) (- (* (- (* b a) (* i c)) (- (* y x) (* t z))) (- (* (- (* j x) (* k z)) (- (* y0 b) (* i y1))) (* (- (* y2 x) (* y3 z)) (- (* c y0) (* y1 a)))))) (+ (- (+ (+ (- (* (- (* x y) (* z t)) (- (* a b) (* c i))) (- (* k (* i (* z y1))) (+ (* j (* i (* x y1))) (* y0 (* k (* z b)))))) (- (* z (* y3 (* a y1))) (+ (* y2 (* x (* a y1))) (* y0 (* z (* c y3)))))) (* (- (* t j) (* y k)) (- (* y4 b) (* y5 i)))) (* (- (* t y2) (* y y3)) (- (* y4 c) (* y5 a)))) (* (- (* k y2) (* j y3)) (- (* y4 y1) (* y5 y0))))))))))

  (+ (- (+ (+ (- (* (- (* x y) (* z t)) (- (* a b) (* c i))) (* (- (* x j) (* z k)) (- (* y0 b) (* y1 i)))) (* (- (* x y2) (* z y3)) (- (* y0 c) (* y1 a)))) (* (- (* t j) (* y k)) (- (* y4 b) (* y5 i)))) (* (- (* t y2) (* y y3)) (- (* y4 c) (* y5 a)))) (* (- (* k y2) (* j y3)) (- (* y4 y1) (* y5 y0)))))