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;
}




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
Results
| Original | 12.3 |
|---|---|
| Target | 16.2 |
| Herbie | 10.0 |
if j < -1.6761456909502331e-32 or 1.90256650834175349e-125 < j Initial program 8.9
rmApplied add-cube-cbrt9.1
Applied associate-*r*9.2
Simplified9.2
rmApplied sub-neg9.2
Applied distribute-lft-in9.2
Simplified9.8
Simplified9.4
rmApplied associate-*r*9.0
rmApplied associate-*r*9.4
if -1.6761456909502331e-32 < j < 1.90256650834175349e-125Initial program 16.3
rmApplied sub-neg16.3
Applied distribute-lft-in16.3
Simplified13.6
Simplified10.7
Final simplification10.0
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)))))