\frac{x \cdot y + z \cdot \left(t - a\right)}{y + z \cdot \left(b - y\right)}\left(x \cdot y + z \cdot \left(t - a\right)\right) \cdot \frac{1}{y + z \cdot \left(b - y\right)}double f(double x, double y, double z, double t, double a, double b) {
double r1306610 = x;
double r1306611 = y;
double r1306612 = r1306610 * r1306611;
double r1306613 = z;
double r1306614 = t;
double r1306615 = a;
double r1306616 = r1306614 - r1306615;
double r1306617 = r1306613 * r1306616;
double r1306618 = r1306612 + r1306617;
double r1306619 = b;
double r1306620 = r1306619 - r1306611;
double r1306621 = r1306613 * r1306620;
double r1306622 = r1306611 + r1306621;
double r1306623 = r1306618 / r1306622;
return r1306623;
}
double f(double x, double y, double z, double t, double a, double b) {
double r1306624 = x;
double r1306625 = y;
double r1306626 = r1306624 * r1306625;
double r1306627 = z;
double r1306628 = t;
double r1306629 = a;
double r1306630 = r1306628 - r1306629;
double r1306631 = r1306627 * r1306630;
double r1306632 = r1306626 + r1306631;
double r1306633 = 1.0;
double r1306634 = b;
double r1306635 = r1306634 - r1306625;
double r1306636 = r1306627 * r1306635;
double r1306637 = r1306625 + r1306636;
double r1306638 = r1306633 / r1306637;
double r1306639 = r1306632 * r1306638;
return r1306639;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a




Bits error versus b
Results
| Original | 23.5 |
|---|---|
| Target | 18.2 |
| Herbie | 23.6 |
Initial program 23.5
rmApplied clear-num23.6
rmApplied *-un-lft-identity23.6
Applied *-un-lft-identity23.6
Applied times-frac23.6
Applied add-cube-cbrt23.6
Applied times-frac23.6
Simplified23.6
Simplified23.5
rmApplied div-inv23.6
Final simplification23.6
herbie shell --seed 2020065
(FPCore (x y z t a b)
:name "Development.Shake.Progress:decay from shake-0.15.5"
:precision binary64
:herbie-target
(- (/ (+ (* z t) (* y x)) (+ y (* z (- b y)))) (/ a (+ (- b y) (/ y z))))
(/ (+ (* x y) (* z (- t a))) (+ y (* z (- b y)))))