\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 r450986 = x;
double r450987 = y;
double r450988 = r450986 * r450987;
double r450989 = z;
double r450990 = t;
double r450991 = r450989 * r450990;
double r450992 = 16.0;
double r450993 = r450991 / r450992;
double r450994 = r450988 + r450993;
double r450995 = a;
double r450996 = b;
double r450997 = r450995 * r450996;
double r450998 = 4.0;
double r450999 = r450997 / r450998;
double r451000 = r450994 - r450999;
double r451001 = c;
double r451002 = r451000 + r451001;
return r451002;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r451003 = x;
double r451004 = y;
double r451005 = r451003 * r451004;
double r451006 = z;
double r451007 = t;
double r451008 = r451006 * r451007;
double r451009 = 16.0;
double r451010 = r451008 / r451009;
double r451011 = r451005 + r451010;
double r451012 = a;
double r451013 = b;
double r451014 = r451012 * r451013;
double r451015 = 4.0;
double r451016 = r451014 / r451015;
double r451017 = r451011 - r451016;
double r451018 = c;
double r451019 = r451017 + r451018;
return r451019;
}



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