\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 r213595 = x;
double r213596 = y;
double r213597 = r213595 * r213596;
double r213598 = z;
double r213599 = t;
double r213600 = r213598 * r213599;
double r213601 = 16.0;
double r213602 = r213600 / r213601;
double r213603 = r213597 + r213602;
double r213604 = a;
double r213605 = b;
double r213606 = r213604 * r213605;
double r213607 = 4.0;
double r213608 = r213606 / r213607;
double r213609 = r213603 - r213608;
double r213610 = c;
double r213611 = r213609 + r213610;
return r213611;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r213612 = x;
double r213613 = y;
double r213614 = r213612 * r213613;
double r213615 = z;
double r213616 = t;
double r213617 = r213615 * r213616;
double r213618 = 16.0;
double r213619 = r213617 / r213618;
double r213620 = r213614 + r213619;
double r213621 = a;
double r213622 = b;
double r213623 = r213621 * r213622;
double r213624 = 4.0;
double r213625 = r213623 / r213624;
double r213626 = r213620 - r213625;
double r213627 = c;
double r213628 = r213626 + r213627;
return r213628;
}



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