\frac{x \cdot y - z \cdot t}{a}\frac{x \cdot y + \left(-t \cdot z\right)}{a}double f(double x, double y, double z, double t, double a) {
double r500347 = x;
double r500348 = y;
double r500349 = r500347 * r500348;
double r500350 = z;
double r500351 = t;
double r500352 = r500350 * r500351;
double r500353 = r500349 - r500352;
double r500354 = a;
double r500355 = r500353 / r500354;
return r500355;
}
double f(double x, double y, double z, double t, double a) {
double r500356 = x;
double r500357 = y;
double r500358 = r500356 * r500357;
double r500359 = t;
double r500360 = z;
double r500361 = r500359 * r500360;
double r500362 = -r500361;
double r500363 = r500358 + r500362;
double r500364 = a;
double r500365 = r500363 / r500364;
return r500365;
}




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
Simplified7.8
Final simplification7.8
herbie shell --seed 2019303 +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.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))