double code(double x, double y, double z, double t, double a) {
return ((double) (((double) (((double) (x * y)) - ((double) (z * t)))) / a));
}
double code(double x, double y, double z, double t, double a) {
double VAR;
if ((((double) (((double) (((double) (x * y)) - ((double) (z * t)))) / a)) <= -inf.0)) {
VAR = ((double) (((double) (x * ((double) (y / a)))) - ((double) (((double) sqrt(((double) (z * ((double) (t / a)))))) * ((double) sqrt(((double) (z * ((double) (t / a))))))))));
} else {
double VAR_1;
if ((((double) (((double) (((double) (x * y)) - ((double) (z * t)))) / a)) <= 7.935708403600279e+291)) {
VAR_1 = ((double) (((double) (((double) (x * y)) - ((double) (z * t)))) / a));
} else {
VAR_1 = ((double) (((double) (((double) cbrt(((double) (y / a)))) * ((double) (x * ((double) (((double) cbrt(((double) (y / a)))) * ((double) cbrt(((double) (y / a)))))))))) - ((double) (z * ((double) (t / a))))));
}
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
Results
| Original | 7.8 |
|---|---|
| Target | 6.3 |
| Herbie | 1.6 |
if (/ (- (* x y) (* z t)) a) < -inf.0Initial program 64.0
rmApplied div-sub64.0
Simplified31.6
Simplified0.3
rmApplied add-sqr-sqrt9.4
if -inf.0 < (/ (- (* x y) (* z t)) a) < 7.9357084036002786e291Initial program 0.9
if 7.9357084036002786e291 < (/ (- (* x y) (* z t)) a) Initial program 53.2
rmApplied div-sub53.2
Simplified30.3
Simplified4.2
rmApplied add-cube-cbrt4.7
Applied associate-*r*4.7
Final simplification1.6
herbie shell --seed 2020179
(FPCore (x y z t a)
:name "Data.Colour.Matrix:inverse from colour-2.3.3, B"
:precision binary64
:herbie-target
(if (< z -2.468684968699548e+170) (- (* (/ y a) x) (* (/ t a) z)) (if (< z 6.309831121978371e-71) (/ (- (* x y) (* z t)) a) (- (* (/ y a) x) (* (/ t a) z))))
(/ (- (* x y) (* z t)) a))