Average Error: 12.3 → 10.0
Time: 13.4s
Precision: binary64
\[\]
\[\]
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
	return ((double) (((double) (((double) (x * ((double) (((double) (y * z)) - ((double) (t * a)))))) - ((double) (b * ((double) (((double) (c * z)) - ((double) (i * a)))))))) + ((double) (j * ((double) (((double) (c * t)) - ((double) (i * y))))))));
}
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
	double VAR;
	if (((j <= -1.676145690950233e-32) || !(j <= 1.9025665083417535e-125))) {
		VAR = ((double) (((double) (((double) (((double) (((double) (z * x)) * y)) - ((double) (((double) (x * a)) * t)))) - ((double) (((double) cbrt(((double) (((double) (z * c)) - ((double) (a * i)))))) * ((double) (b * ((double) (((double) cbrt(((double) (((double) (z * c)) - ((double) (a * i)))))) * ((double) cbrt(((double) (((double) (z * c)) - ((double) (a * i)))))))))))))) + ((double) (j * ((double) (((double) (t * c)) - ((double) (y * i))))))));
	} else {
		VAR = ((double) (((double) (((double) (x * ((double) (((double) (z * y)) - ((double) (a * t)))))) + ((double) (b * ((double) (((double) (a * i)) - ((double) (z * c)))))))) + ((double) (((double) (t * ((double) (j * c)))) - ((double) (y * ((double) (j * i))))))));
	}
	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

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Target

Original12.3
Target16.2
Herbie10.0
\[\]

Derivation

  1. Split input into 2 regimes
  2. if j < -1.6761456909502331e-32 or 1.90256650834175349e-125 < j

    1. Initial program 8.9

      \[\]
    2. Using strategy rm
    3. Applied add-cube-cbrt9.1

      \[\leadsto \]
    4. Applied associate-*r*9.2

      \[\leadsto \]
    5. Simplified9.2

      \[\leadsto \]
    6. Using strategy rm
    7. Applied sub-neg9.2

      \[\leadsto \]
    8. Applied distribute-lft-in9.2

      \[\leadsto \]
    9. Simplified9.8

      \[\leadsto \]
    10. Simplified9.4

      \[\leadsto \]
    11. Using strategy rm
    12. Applied associate-*r*9.0

      \[\leadsto \]
    13. Using strategy rm
    14. Applied associate-*r*9.4

      \[\leadsto \]

    if -1.6761456909502331e-32 < j < 1.90256650834175349e-125

    1. Initial program 16.3

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

      \[\leadsto \]
    4. Applied distribute-lft-in16.3

      \[\leadsto \]
    5. Simplified13.6

      \[\leadsto \]
    6. Simplified10.7

      \[\leadsto \]
  3. Recombined 2 regimes into one program.
  4. Final simplification10.0

    \[\leadsto \]

Reproduce

herbie shell --seed 2020191 
(FPCore (x y z t a b c i j)
  :name "Linear.Matrix:det33 from linear-1.19.1.3"
  :precision binary64

  :herbie-target
  (if (< t -8.120978919195912e-33) (- (* x (- (* z y) (* a t))) (- (* b (- (* z c) (* a i))) (* (- (* c t) (* y i)) j))) (if (< t -4.712553818218485e-169) (+ (- (* x (- (* y z) (* t a))) (* b (- (* c z) (* i a)))) (/ (* j (- (pow (* c t) 2.0) (pow (* i y) 2.0))) (+ (* c t) (* i y)))) (if (< t -7.633533346031584e-308) (- (* x (- (* z y) (* a t))) (- (* b (- (* z c) (* a i))) (* (- (* c t) (* y i)) j))) (if (< t 1.0535888557455487e-139) (+ (- (* x (- (* y z) (* t a))) (* b (- (* c z) (* i a)))) (/ (* j (- (pow (* c t) 2.0) (pow (* i y) 2.0))) (+ (* c t) (* i y)))) (- (* x (- (* z y) (* a t))) (- (* b (- (* z c) (* a i))) (* (- (* c t) (* y i)) j)))))))

  (+ (- (* x (- (* y z) (* t a))) (* b (- (* c z) (* i a)))) (* j (- (* c t) (* i y)))))