\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 r271223 = x;
double r271224 = y;
double r271225 = r271223 * r271224;
double r271226 = z;
double r271227 = t;
double r271228 = r271226 * r271227;
double r271229 = 16.0;
double r271230 = r271228 / r271229;
double r271231 = r271225 + r271230;
double r271232 = a;
double r271233 = b;
double r271234 = r271232 * r271233;
double r271235 = 4.0;
double r271236 = r271234 / r271235;
double r271237 = r271231 - r271236;
double r271238 = c;
double r271239 = r271237 + r271238;
return r271239;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r271240 = x;
double r271241 = y;
double r271242 = r271240 * r271241;
double r271243 = z;
double r271244 = t;
double r271245 = r271243 * r271244;
double r271246 = 16.0;
double r271247 = r271245 / r271246;
double r271248 = r271242 + r271247;
double r271249 = a;
double r271250 = b;
double r271251 = r271249 * r271250;
double r271252 = 4.0;
double r271253 = r271251 / r271252;
double r271254 = r271248 - r271253;
double r271255 = c;
double r271256 = r271254 + r271255;
return r271256;
}



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