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) (t * i)))))))) + ((double) (j * ((double) (((double) (c * a)) - ((double) (y * i))))))));
}
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double VAR;
if ((i <= -1.7241027508523938e+92)) {
VAR = ((double) (((double) (((double) (x * ((double) (((double) (y * z)) - ((double) (t * a)))))) + ((double) (b * ((double) (((double) (i * t)) - ((double) (z * c)))))))) + ((double) (((double) (c * ((double) (a * j)))) - ((double) (i * ((double) (y * j))))))));
} else {
double VAR_1;
if ((i <= -113366.18846745034)) {
VAR_1 = ((double) (((double) (((double) (((double) (y * ((double) (x * z)))) - ((double) (a * ((double) (x * t)))))) + ((double) (((double) (((double) cbrt(b)) * ((double) cbrt(b)))) * ((double) (((double) cbrt(b)) * ((double) (((double) (i * t)) - ((double) (z * c)))))))))) + ((double) (j * ((double) (((double) (a * c)) - ((double) (i * y))))))));
} else {
double VAR_2;
if ((i <= 6.209812177855934e-197)) {
VAR_2 = ((double) (((double) (j * ((double) (((double) (a * c)) - ((double) (i * y)))))) + ((double) (((double) (x * ((double) (((double) (y * z)) - ((double) (t * a)))))) + ((double) (((double) cbrt(((double) cbrt(b)))) * ((double) (((double) (((double) cbrt(b)) * ((double) (((double) cbrt(b)) * ((double) cbrt(((double) (((double) cbrt(b)) * ((double) cbrt(b)))))))))) * ((double) (((double) (i * t)) - ((double) (z * c))))))))))));
} else {
double VAR_3;
if ((i <= 1.3065923009124668e-46)) {
VAR_3 = ((double) (((double) (((double) (((double) (y * ((double) (x * z)))) - ((double) (a * ((double) (x * t)))))) + ((double) (((double) (((double) cbrt(b)) * ((double) cbrt(b)))) * ((double) (((double) cbrt(b)) * ((double) (((double) (i * t)) - ((double) (z * c)))))))))) + ((double) (j * ((double) (((double) (a * c)) - ((double) (i * y))))))));
} else {
VAR_3 = ((double) (((double) (((double) (x * ((double) (((double) (y * z)) - ((double) (t * a)))))) + ((double) (b * ((double) (((double) (i * t)) - ((double) (z * c)))))))) + ((double) (((double) (c * ((double) (a * j)))) - ((double) (i * ((double) (y * j))))))));
}
VAR_2 = VAR_3;
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
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.0 |
|---|---|
| Target | 20.0 |
| Herbie | 10.7 |
if i < -1.7241027508523938e92 or 1.30659230091246675e-46 < i Initial program 16.9
rmApplied sub-neg16.9
Applied distribute-lft-in16.9
Simplified16.9
Simplified12.7
if -1.7241027508523938e92 < i < -113366.188467450338 or 6.20981217785593431e-197 < i < 1.30659230091246675e-46Initial program 10.0
rmApplied add-cube-cbrt10.3
Applied associate-*l*10.3
Simplified10.3
rmApplied sub-neg10.3
Applied distribute-lft-in10.3
Simplified10.4
Simplified10.3
if -113366.188467450338 < i < 6.20981217785593431e-197Initial program 8.8
rmApplied add-cube-cbrt9.1
Applied associate-*l*9.1
Simplified9.1
rmApplied associate-*r*9.1
Simplified9.1
rmApplied add-cube-cbrt9.1
Applied cbrt-prod9.1
Applied associate-*r*9.1
Simplified9.1
Final simplification10.7
herbie shell --seed 2020190
(FPCore (x y z t a b c i j)
:name "Data.Colour.Matrix:determinant from colour-2.3.3, A"
:precision binary64
:herbie-target
(if (< x -1.469694296777705e-64) (+ (- (* x (- (* y z) (* t a))) (/ (* b (- (pow (* c z) 2.0) (pow (* t i) 2.0))) (+ (* c z) (* t i)))) (* j (- (* c a) (* y i)))) (if (< x 3.2113527362226803e-147) (- (* (- (* b i) (* x a)) t) (- (* z (* c b)) (* j (- (* c a) (* y i))))) (+ (- (* x (- (* y z) (* t a))) (/ (* b (- (pow (* c z) 2.0) (pow (* t i) 2.0))) (+ (* c z) (* t i)))) (* j (- (* c a) (* y i))))))
(+ (- (* x (- (* y z) (* t a))) (* b (- (* c z) (* t i)))) (* j (- (* c a) (* y i)))))