\left(\left(x \cdot y + \frac{z \cdot t}{16}\right) - \frac{a \cdot b}{4}\right) + c\left(\left(x \cdot y + \frac{z \cdot t}{16}\right) - \frac{a \cdot b}{4}\right) + cdouble f(double x, double y, double z, double t, double a, double b, double c) {
double r216731 = x;
double r216732 = y;
double r216733 = r216731 * r216732;
double r216734 = z;
double r216735 = t;
double r216736 = r216734 * r216735;
double r216737 = 16.0;
double r216738 = r216736 / r216737;
double r216739 = r216733 + r216738;
double r216740 = a;
double r216741 = b;
double r216742 = r216740 * r216741;
double r216743 = 4.0;
double r216744 = r216742 / r216743;
double r216745 = r216739 - r216744;
double r216746 = c;
double r216747 = r216745 + r216746;
return r216747;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r216748 = x;
double r216749 = y;
double r216750 = r216748 * r216749;
double r216751 = z;
double r216752 = t;
double r216753 = r216751 * r216752;
double r216754 = 16.0;
double r216755 = r216753 / r216754;
double r216756 = r216750 + r216755;
double r216757 = a;
double r216758 = b;
double r216759 = r216757 * r216758;
double r216760 = 4.0;
double r216761 = r216759 / r216760;
double r216762 = r216756 - r216761;
double r216763 = c;
double r216764 = r216762 + r216763;
return r216764;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t



Bits error versus a



Bits error versus b



Bits error versus c
Results
Initial program 0.1
Final simplification0.1
herbie shell --seed 2020060
(FPCore (x y z t a b c)
:name "Diagrams.Solve.Polynomial:quartForm from diagrams-solve-0.1, C"
:precision binary64
(+ (- (+ (* x y) (/ (* z t) 16)) (/ (* a b) 4)) c))