\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 r325051 = x;
double r325052 = y;
double r325053 = r325051 * r325052;
double r325054 = z;
double r325055 = t;
double r325056 = r325054 * r325055;
double r325057 = 16.0;
double r325058 = r325056 / r325057;
double r325059 = r325053 + r325058;
double r325060 = a;
double r325061 = b;
double r325062 = r325060 * r325061;
double r325063 = 4.0;
double r325064 = r325062 / r325063;
double r325065 = r325059 - r325064;
double r325066 = c;
double r325067 = r325065 + r325066;
return r325067;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r325068 = x;
double r325069 = y;
double r325070 = r325068 * r325069;
double r325071 = z;
double r325072 = t;
double r325073 = r325071 * r325072;
double r325074 = 16.0;
double r325075 = r325073 / r325074;
double r325076 = r325070 + r325075;
double r325077 = a;
double r325078 = b;
double r325079 = r325077 * r325078;
double r325080 = 4.0;
double r325081 = r325079 / r325080;
double r325082 = r325076 - r325081;
double r325083 = c;
double r325084 = r325082 + r325083;
return r325084;
}



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