\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 r261070 = x;
double r261071 = y;
double r261072 = r261070 * r261071;
double r261073 = z;
double r261074 = t;
double r261075 = r261073 * r261074;
double r261076 = 16.0;
double r261077 = r261075 / r261076;
double r261078 = r261072 + r261077;
double r261079 = a;
double r261080 = b;
double r261081 = r261079 * r261080;
double r261082 = 4.0;
double r261083 = r261081 / r261082;
double r261084 = r261078 - r261083;
double r261085 = c;
double r261086 = r261084 + r261085;
return r261086;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r261087 = x;
double r261088 = y;
double r261089 = r261087 * r261088;
double r261090 = z;
double r261091 = t;
double r261092 = r261090 * r261091;
double r261093 = 16.0;
double r261094 = r261092 / r261093;
double r261095 = r261089 + r261094;
double r261096 = a;
double r261097 = b;
double r261098 = r261096 * r261097;
double r261099 = 4.0;
double r261100 = r261098 / r261099;
double r261101 = r261095 - r261100;
double r261102 = c;
double r261103 = r261101 + r261102;
return r261103;
}



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 2020034
(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))