\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 r156498 = x;
double r156499 = y;
double r156500 = r156498 * r156499;
double r156501 = z;
double r156502 = t;
double r156503 = r156501 * r156502;
double r156504 = 16.0;
double r156505 = r156503 / r156504;
double r156506 = r156500 + r156505;
double r156507 = a;
double r156508 = b;
double r156509 = r156507 * r156508;
double r156510 = 4.0;
double r156511 = r156509 / r156510;
double r156512 = r156506 - r156511;
double r156513 = c;
double r156514 = r156512 + r156513;
return r156514;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r156515 = x;
double r156516 = y;
double r156517 = r156515 * r156516;
double r156518 = z;
double r156519 = t;
double r156520 = r156518 * r156519;
double r156521 = 16.0;
double r156522 = r156520 / r156521;
double r156523 = r156517 + r156522;
double r156524 = a;
double r156525 = b;
double r156526 = r156524 * r156525;
double r156527 = 4.0;
double r156528 = r156526 / r156527;
double r156529 = r156523 - r156528;
double r156530 = c;
double r156531 = r156529 + r156530;
return r156531;
}



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