\frac{x \cdot y - z \cdot t}{a}\frac{x \cdot y + \left(-z \cdot t\right)}{a}double f(double x, double y, double z, double t, double a) {
double r521614 = x;
double r521615 = y;
double r521616 = r521614 * r521615;
double r521617 = z;
double r521618 = t;
double r521619 = r521617 * r521618;
double r521620 = r521616 - r521619;
double r521621 = a;
double r521622 = r521620 / r521621;
return r521622;
}
double f(double x, double y, double z, double t, double a) {
double r521623 = x;
double r521624 = y;
double r521625 = r521623 * r521624;
double r521626 = z;
double r521627 = t;
double r521628 = r521626 * r521627;
double r521629 = -r521628;
double r521630 = r521625 + r521629;
double r521631 = a;
double r521632 = r521630 / r521631;
return r521632;
}




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 | 5.8 |
| Herbie | 7.8 |
Initial program 7.8
rmApplied sub-neg7.8
Final simplification7.8
herbie shell --seed 2019303
(FPCore (x y z t a)
:name "Data.Colour.Matrix:inverse from colour-2.3.3, B"
:precision binary64
:herbie-target
(if (< z -2.46868496869954822e170) (- (* (/ y a) x) (* (/ t a) z)) (if (< z 6.30983112197837121e-71) (/ (- (* x y) (* z t)) a) (- (* (/ y a) x) (* (/ t a) z))))
(/ (- (* x y) (* z t)) a))