\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 r186352 = x;
double r186353 = y;
double r186354 = r186352 * r186353;
double r186355 = z;
double r186356 = t;
double r186357 = r186355 * r186356;
double r186358 = 16.0;
double r186359 = r186357 / r186358;
double r186360 = r186354 + r186359;
double r186361 = a;
double r186362 = b;
double r186363 = r186361 * r186362;
double r186364 = 4.0;
double r186365 = r186363 / r186364;
double r186366 = r186360 - r186365;
double r186367 = c;
double r186368 = r186366 + r186367;
return r186368;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r186369 = x;
double r186370 = y;
double r186371 = r186369 * r186370;
double r186372 = z;
double r186373 = t;
double r186374 = r186372 * r186373;
double r186375 = 16.0;
double r186376 = r186374 / r186375;
double r186377 = r186371 + r186376;
double r186378 = a;
double r186379 = b;
double r186380 = r186378 * r186379;
double r186381 = 4.0;
double r186382 = r186380 / r186381;
double r186383 = r186377 - r186382;
double r186384 = c;
double r186385 = r186383 + r186384;
return r186385;
}



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