\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 r318985 = x;
double r318986 = y;
double r318987 = r318985 * r318986;
double r318988 = z;
double r318989 = t;
double r318990 = r318988 * r318989;
double r318991 = 16.0;
double r318992 = r318990 / r318991;
double r318993 = r318987 + r318992;
double r318994 = a;
double r318995 = b;
double r318996 = r318994 * r318995;
double r318997 = 4.0;
double r318998 = r318996 / r318997;
double r318999 = r318993 - r318998;
double r319000 = c;
double r319001 = r318999 + r319000;
return r319001;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r319002 = x;
double r319003 = y;
double r319004 = r319002 * r319003;
double r319005 = z;
double r319006 = t;
double r319007 = r319005 * r319006;
double r319008 = 16.0;
double r319009 = r319007 / r319008;
double r319010 = r319004 + r319009;
double r319011 = a;
double r319012 = b;
double r319013 = r319011 * r319012;
double r319014 = 4.0;
double r319015 = r319013 / r319014;
double r319016 = r319010 - r319015;
double r319017 = c;
double r319018 = r319016 + r319017;
return r319018;
}



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