Average Error: 12.1 → 9.5
Time: 13.0s
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 ((x <= -2.8015300348686148e-43)) {
		VAR = ((double) (((double) (((double) (x * ((double) (((double) (y * z)) - ((double) (t * a)))))) + ((double) (b * ((double) (((double) (a * i)) - ((double) (z * c)))))))) + ((double) (((double) (((double) cbrt(((double) (j * ((double) (((double) (t * c)) - ((double) (y * i)))))))) * ((double) cbrt(((double) (j * ((double) (((double) (t * c)) - ((double) (y * i)))))))))) * ((double) cbrt(((double) (((double) cbrt(((double) (((double) (t * c)) - ((double) (y * i)))))) * ((double) (j * ((double) (((double) cbrt(((double) (((double) (t * c)) - ((double) (y * i)))))) * ((double) cbrt(((double) (((double) (t * c)) - ((double) (y * i))))))))))))))))));
	} else {
		double VAR_1;
		if ((x <= 9.369818488516253e+54)) {
			VAR_1 = ((double) (((double) (j * ((double) (((double) (t * c)) - ((double) (y * i)))))) + ((double) (((double) (((double) (y * ((double) (x * z)))) - ((double) (t * ((double) (x * a)))))) + ((double) (((double) (i * ((double) (a * b)))) - ((double) (c * ((double) (z * b))))))))));
		} else {
			VAR_1 = ((double) (((double) (((double) (x * ((double) (((double) (y * z)) - ((double) (t * a)))))) + ((double) (b * ((double) (((double) (a * i)) - ((double) (z * c)))))))) + ((double) (((double) (((double) cbrt(((double) (j * ((double) (((double) (t * c)) - ((double) (y * i)))))))) * ((double) cbrt(((double) (j * ((double) (((double) (t * c)) - ((double) (y * i)))))))))) * ((double) cbrt(((double) (((double) (((double) cbrt(j)) * ((double) cbrt(j)))) * ((double) (((double) (((double) (t * c)) - ((double) (y * i)))) * ((double) cbrt(j))))))))))));
		}
		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
Try it out
Your Program's Arguments
Results
Enter valid numbers for all inputs
Target
Original12.1
Target16.2
Herbie9.5
\[\]
Derivation
  1. Split input into 3 regimes
  2. if  x  < -2.8015300348686148e-43
    1. Initial program 7.6
      \[\]
    2. Using strategy rm
    3. Applied add-cube-cbrt7.8
      \[\leadsto \]
    4. Using strategy rm
    5. Applied add-cube-cbrt7.8
      \[\leadsto \]
    6. Applied associate-*r*7.8
      \[\leadsto \]
    if -2.8015300348686148e-43 <  x  < 9.36981848851625302e54
    1. Initial program 15.0
      \[\]
    2. Using strategy rm
    3. Applied sub-neg15.0
      \[\leadsto \]
    4. Applied distribute-lft-in15.0
      \[\leadsto \]
    5. Simplified15.2
      \[\leadsto \]
    6. Simplified15.0
      \[\leadsto \]
    7. Using strategy rm
    8. Applied sub-neg15.0
      \[\leadsto \]
    9. Applied distribute-lft-in15.0
      \[\leadsto \]
    10. Simplified13.2
      \[\leadsto \]
    11. Simplified10.8
      \[\leadsto \]
    if 9.36981848851625302e54 <  x 
    1. Initial program 6.7
      \[\]
    2. Using strategy rm
    3. Applied add-cube-cbrt6.9
      \[\leadsto \]
    4. Using strategy rm
    5. Applied add-cube-cbrt6.9
      \[\leadsto \]
    6. Applied associate-*l*6.9
      \[\leadsto \]
    7. Simplified6.9
      \[\leadsto \]
  3. Recombined 3 regimes into one program.
  4. Final simplification9.5
    \[\leadsto \]
Reproduce
herbie shell --seed 2020190 
(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)))))