\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 r240535 = x;
double r240536 = y;
double r240537 = r240535 * r240536;
double r240538 = z;
double r240539 = t;
double r240540 = r240538 * r240539;
double r240541 = 16.0;
double r240542 = r240540 / r240541;
double r240543 = r240537 + r240542;
double r240544 = a;
double r240545 = b;
double r240546 = r240544 * r240545;
double r240547 = 4.0;
double r240548 = r240546 / r240547;
double r240549 = r240543 - r240548;
double r240550 = c;
double r240551 = r240549 + r240550;
return r240551;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r240552 = x;
double r240553 = y;
double r240554 = r240552 * r240553;
double r240555 = z;
double r240556 = t;
double r240557 = r240555 * r240556;
double r240558 = 16.0;
double r240559 = r240557 / r240558;
double r240560 = r240554 + r240559;
double r240561 = a;
double r240562 = b;
double r240563 = r240561 * r240562;
double r240564 = 4.0;
double r240565 = r240563 / r240564;
double r240566 = r240560 - r240565;
double r240567 = c;
double r240568 = r240566 + r240567;
return r240568;
}



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.2
Final simplification0.2
herbie shell --seed 2020001
(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))