\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 r240284 = x;
double r240285 = y;
double r240286 = r240284 * r240285;
double r240287 = z;
double r240288 = t;
double r240289 = r240287 * r240288;
double r240290 = 16.0;
double r240291 = r240289 / r240290;
double r240292 = r240286 + r240291;
double r240293 = a;
double r240294 = b;
double r240295 = r240293 * r240294;
double r240296 = 4.0;
double r240297 = r240295 / r240296;
double r240298 = r240292 - r240297;
double r240299 = c;
double r240300 = r240298 + r240299;
return r240300;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r240301 = x;
double r240302 = y;
double r240303 = r240301 * r240302;
double r240304 = z;
double r240305 = t;
double r240306 = r240304 * r240305;
double r240307 = 16.0;
double r240308 = r240306 / r240307;
double r240309 = r240303 + r240308;
double r240310 = a;
double r240311 = b;
double r240312 = r240310 * r240311;
double r240313 = 4.0;
double r240314 = r240312 / r240313;
double r240315 = r240309 - r240314;
double r240316 = c;
double r240317 = r240315 + r240316;
return r240317;
}



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