\frac{x \cdot y - z \cdot t}{a}\frac{x \cdot y - z \cdot t}{a}double f(double x, double y, double z, double t, double a) {
double r163822 = x;
double r163823 = y;
double r163824 = r163822 * r163823;
double r163825 = z;
double r163826 = t;
double r163827 = r163825 * r163826;
double r163828 = r163824 - r163827;
double r163829 = a;
double r163830 = r163828 / r163829;
return r163830;
}
double f(double x, double y, double z, double t, double a) {
double r163831 = x;
double r163832 = y;
double r163833 = r163831 * r163832;
double r163834 = z;
double r163835 = t;
double r163836 = r163834 * r163835;
double r163837 = r163833 - r163836;
double r163838 = a;
double r163839 = r163837 / r163838;
return r163839;
}




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.8 |
| Herbie | 7.5 |
Initial program 7.5
Final simplification7.5
herbie shell --seed 2020045 +o rules:numerics
(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))