Average Error: 12.3 → 9.6
Time: 12.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 ((j <= -6.809710970042638e+27)) {
		VAR = ((double) (((double) (((double) (((double) (y * ((double) (x * z)))) - ((double) (t * ((double) (x * a)))))) + ((double) (((double) (((double) cbrt(b)) * ((double) cbrt(b)))) * ((double) (((double) cbrt(b)) * ((double) (((double) (a * i)) - ((double) (z * c)))))))))) + ((double) (j * ((double) (((double) (t * c)) - ((double) (y * i))))))));
	} else {
		double VAR_1;
		if ((j <= 7.11380026928101e-16)) {
			VAR_1 = ((double) (((double) (((double) (x * ((double) (((double) (y * z)) - ((double) (t * a)))))) + ((double) (b * ((double) (((double) (a * i)) - ((double) (z * c)))))))) + ((double) (((double) (c * ((double) (j * t)))) - ((double) (i * ((double) (j * y))))))));
		} else {
			VAR_1 = ((double) (((double) (j * ((double) (((double) (t * c)) - ((double) (y * i)))))) + ((double) (((double) (((double) (z * ((double) (y * x)))) - ((double) (a * ((double) (x * t)))))) + ((double) (b * ((double) (((double) (a * i)) - ((double) (z * c))))))))));
		}
		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.3
Target16.3
Herbie9.6
\[\]
Derivation
  1. Split input into 3 regimes
  2. if  j  < -6.8097109700426381e27
    1. Initial program 7.3
      \[\]
    2. Using strategy rm
    3. Applied add-cube-cbrt7.5
      \[\leadsto \]
    4. Applied associate-*l*7.5
      \[\leadsto \]
    5. Simplified7.5
      \[\leadsto \]
    6. Using strategy rm
    7. Applied sub-neg7.5
      \[\leadsto \]
    8. Applied distribute-lft-in7.5
      \[\leadsto \]
    9. Simplified8.2
      \[\leadsto \]
    10. Simplified8.6
      \[\leadsto \]
    if -6.8097109700426381e27 <  j  < 7.1138002692810102e-16
    1. Initial program 15.2
      \[\]
    2. Using strategy rm
    3. Applied sub-neg15.2
      \[\leadsto \]
    4. Applied distribute-lft-in15.2
      \[\leadsto \]
    5. Simplified12.7
      \[\leadsto \]
    6. Simplified10.1
      \[\leadsto \]
    if 7.1138002692810102e-16 <  j 
    1. Initial program 8.1
      \[\]
    2. Using strategy rm
    3. Applied sub-neg8.1
      \[\leadsto \]
    4. Applied distribute-lft-in8.1
      \[\leadsto \]
    5. Simplified7.8
      \[\leadsto \]
    6. Simplified8.8
      \[\leadsto \]
  3. Recombined 3 regimes into one program.
  4. Final simplification9.6
    \[\leadsto \]
Reproduce
herbie shell --seed 2020192 
(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)))))