\frac{x \cdot y - z \cdot t}{a}\frac{x \cdot y - t \cdot z}{a}double f(double x, double y, double z, double t, double a) {
double r52754483 = x;
double r52754484 = y;
double r52754485 = r52754483 * r52754484;
double r52754486 = z;
double r52754487 = t;
double r52754488 = r52754486 * r52754487;
double r52754489 = r52754485 - r52754488;
double r52754490 = a;
double r52754491 = r52754489 / r52754490;
return r52754491;
}
double f(double x, double y, double z, double t, double a) {
double r52754492 = x;
double r52754493 = y;
double r52754494 = r52754492 * r52754493;
double r52754495 = t;
double r52754496 = z;
double r52754497 = r52754495 * r52754496;
double r52754498 = r52754494 - r52754497;
double r52754499 = a;
double r52754500 = r52754498 / r52754499;
return r52754500;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 7.9 |
|---|---|
| Target | 5.8 |
| Herbie | 7.9 |
Initial program 7.9
rmApplied *-un-lft-identity7.9
Applied associate-/r*7.9
Simplified7.9
Final simplification7.9
herbie shell --seed 2019174 +o rules:numerics
(FPCore (x y z t a)
:name "Data.Colour.Matrix:inverse from colour-2.3.3, B"
: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))