\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 r131481 = x;
double r131482 = y;
double r131483 = r131481 * r131482;
double r131484 = z;
double r131485 = t;
double r131486 = r131484 * r131485;
double r131487 = 16.0;
double r131488 = r131486 / r131487;
double r131489 = r131483 + r131488;
double r131490 = a;
double r131491 = b;
double r131492 = r131490 * r131491;
double r131493 = 4.0;
double r131494 = r131492 / r131493;
double r131495 = r131489 - r131494;
double r131496 = c;
double r131497 = r131495 + r131496;
return r131497;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r131498 = x;
double r131499 = y;
double r131500 = r131498 * r131499;
double r131501 = z;
double r131502 = t;
double r131503 = r131501 * r131502;
double r131504 = 16.0;
double r131505 = r131503 / r131504;
double r131506 = r131500 + r131505;
double r131507 = a;
double r131508 = b;
double r131509 = r131507 * r131508;
double r131510 = 4.0;
double r131511 = r131509 / r131510;
double r131512 = r131506 - r131511;
double r131513 = c;
double r131514 = r131512 + r131513;
return r131514;
}



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 2019304 +o rules:numerics
(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))