\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 r267707 = x;
double r267708 = y;
double r267709 = r267707 * r267708;
double r267710 = z;
double r267711 = t;
double r267712 = r267710 * r267711;
double r267713 = 16.0;
double r267714 = r267712 / r267713;
double r267715 = r267709 + r267714;
double r267716 = a;
double r267717 = b;
double r267718 = r267716 * r267717;
double r267719 = 4.0;
double r267720 = r267718 / r267719;
double r267721 = r267715 - r267720;
double r267722 = c;
double r267723 = r267721 + r267722;
return r267723;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r267724 = x;
double r267725 = y;
double r267726 = r267724 * r267725;
double r267727 = z;
double r267728 = t;
double r267729 = r267727 * r267728;
double r267730 = 16.0;
double r267731 = r267729 / r267730;
double r267732 = r267726 + r267731;
double r267733 = a;
double r267734 = b;
double r267735 = r267733 * r267734;
double r267736 = 4.0;
double r267737 = r267735 / r267736;
double r267738 = r267732 - r267737;
double r267739 = c;
double r267740 = r267738 + r267739;
return r267740;
}



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.2
Final simplification0.2
herbie shell --seed 2020049
(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))