\frac{x \cdot y - z \cdot t}{a}\begin{array}{l}
\mathbf{if}\;x \cdot y - z \cdot t \le -1.4289436252700815 \cdot 10^{209}:\\
\;\;\;\;x \cdot \frac{y}{a} - \frac{t}{\frac{a}{z}}\\
\mathbf{elif}\;x \cdot y - z \cdot t \le 3.05233143825895052 \cdot 10^{218}:\\
\;\;\;\;\frac{1}{\frac{a}{x \cdot y - z \cdot t}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\frac{a}{y}} - \frac{t}{\frac{a}{z}}\\
\end{array}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) (x * y)) - ((double) (z * t)))) <= -1.4289436252700815e+209)) {
VAR = ((double) (((double) (x * ((double) (y / a)))) - ((double) (t / ((double) (a / z))))));
} else {
double VAR_1;
if ((((double) (((double) (x * y)) - ((double) (z * t)))) <= 3.0523314382589505e+218)) {
VAR_1 = ((double) (1.0 / ((double) (a / ((double) (((double) (x * y)) - ((double) (z * t))))))));
} else {
VAR_1 = ((double) (((double) (x / ((double) (a / y)))) - ((double) (t / ((double) (a / z))))));
}
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.5 |
|---|---|
| Target | 5.5 |
| Herbie | 1.1 |
if (- (* x y) (* z t)) < -1.4289436252700815e+209Initial program 29.1
rmApplied div-sub29.1
Simplified29.1
rmApplied associate-/l*14.1
rmApplied *-un-lft-identity14.1
Applied times-frac1.1
Simplified1.1
if -1.4289436252700815e+209 < (- (* x y) (* z t)) < 3.0523314382589505e+218Initial program 0.8
rmApplied clear-num1.2
if 3.0523314382589505e+218 < (- (* x y) (* z t)) Initial program 32.6
rmApplied div-sub32.6
Simplified32.6
rmApplied associate-/l*17.3
rmApplied associate-/l*1.1
Final simplification1.1
herbie shell --seed 2020130
(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))