\frac{x \cdot y - z \cdot t}{a}\frac{x \cdot y}{a} - \frac{t \cdot z}{a}double f(double x, double y, double z, double t, double a) {
double r849830 = x;
double r849831 = y;
double r849832 = r849830 * r849831;
double r849833 = z;
double r849834 = t;
double r849835 = r849833 * r849834;
double r849836 = r849832 - r849835;
double r849837 = a;
double r849838 = r849836 / r849837;
return r849838;
}
double f(double x, double y, double z, double t, double a) {
double r849839 = x;
double r849840 = y;
double r849841 = r849839 * r849840;
double r849842 = a;
double r849843 = r849841 / r849842;
double r849844 = t;
double r849845 = z;
double r849846 = r849844 * r849845;
double r849847 = r849846 / r849842;
double r849848 = r849843 - r849847;
return r849848;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 8.0 |
|---|---|
| Target | 6.2 |
| Herbie | 8.0 |
Initial program 8.0
rmApplied div-inv8.1
rmApplied add-cube-cbrt8.9
Taylor expanded around 0 8.0
Final simplification8.0
herbie shell --seed 2020039
(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))