\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 r165450 = x;
double r165451 = y;
double r165452 = r165450 * r165451;
double r165453 = z;
double r165454 = t;
double r165455 = r165453 * r165454;
double r165456 = 16.0;
double r165457 = r165455 / r165456;
double r165458 = r165452 + r165457;
double r165459 = a;
double r165460 = b;
double r165461 = r165459 * r165460;
double r165462 = 4.0;
double r165463 = r165461 / r165462;
double r165464 = r165458 - r165463;
double r165465 = c;
double r165466 = r165464 + r165465;
return r165466;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r165467 = x;
double r165468 = y;
double r165469 = r165467 * r165468;
double r165470 = z;
double r165471 = t;
double r165472 = r165470 * r165471;
double r165473 = 16.0;
double r165474 = r165472 / r165473;
double r165475 = r165469 + r165474;
double r165476 = a;
double r165477 = b;
double r165478 = r165476 * r165477;
double r165479 = 4.0;
double r165480 = r165478 / r165479;
double r165481 = r165475 - r165480;
double r165482 = c;
double r165483 = r165481 + r165482;
return r165483;
}



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