\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 r166425 = x;
double r166426 = y;
double r166427 = r166425 * r166426;
double r166428 = z;
double r166429 = t;
double r166430 = r166428 * r166429;
double r166431 = 16.0;
double r166432 = r166430 / r166431;
double r166433 = r166427 + r166432;
double r166434 = a;
double r166435 = b;
double r166436 = r166434 * r166435;
double r166437 = 4.0;
double r166438 = r166436 / r166437;
double r166439 = r166433 - r166438;
double r166440 = c;
double r166441 = r166439 + r166440;
return r166441;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r166442 = x;
double r166443 = y;
double r166444 = r166442 * r166443;
double r166445 = z;
double r166446 = t;
double r166447 = r166445 * r166446;
double r166448 = 16.0;
double r166449 = r166447 / r166448;
double r166450 = r166444 + r166449;
double r166451 = a;
double r166452 = b;
double r166453 = r166451 * r166452;
double r166454 = 4.0;
double r166455 = r166453 / r166454;
double r166456 = r166450 - r166455;
double r166457 = c;
double r166458 = r166456 + r166457;
return r166458;
}



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