\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 r254537 = x;
double r254538 = y;
double r254539 = r254537 * r254538;
double r254540 = z;
double r254541 = t;
double r254542 = r254540 * r254541;
double r254543 = 16.0;
double r254544 = r254542 / r254543;
double r254545 = r254539 + r254544;
double r254546 = a;
double r254547 = b;
double r254548 = r254546 * r254547;
double r254549 = 4.0;
double r254550 = r254548 / r254549;
double r254551 = r254545 - r254550;
double r254552 = c;
double r254553 = r254551 + r254552;
return r254553;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r254554 = x;
double r254555 = y;
double r254556 = r254554 * r254555;
double r254557 = z;
double r254558 = t;
double r254559 = r254557 * r254558;
double r254560 = 16.0;
double r254561 = r254559 / r254560;
double r254562 = r254556 + r254561;
double r254563 = a;
double r254564 = b;
double r254565 = r254563 * r254564;
double r254566 = 4.0;
double r254567 = r254565 / r254566;
double r254568 = r254562 - r254567;
double r254569 = c;
double r254570 = r254568 + r254569;
return r254570;
}



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