\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 r249867 = x;
double r249868 = y;
double r249869 = r249867 * r249868;
double r249870 = z;
double r249871 = t;
double r249872 = r249870 * r249871;
double r249873 = 16.0;
double r249874 = r249872 / r249873;
double r249875 = r249869 + r249874;
double r249876 = a;
double r249877 = b;
double r249878 = r249876 * r249877;
double r249879 = 4.0;
double r249880 = r249878 / r249879;
double r249881 = r249875 - r249880;
double r249882 = c;
double r249883 = r249881 + r249882;
return r249883;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r249884 = x;
double r249885 = y;
double r249886 = r249884 * r249885;
double r249887 = z;
double r249888 = t;
double r249889 = r249887 * r249888;
double r249890 = 16.0;
double r249891 = r249889 / r249890;
double r249892 = r249886 + r249891;
double r249893 = a;
double r249894 = b;
double r249895 = r249893 * r249894;
double r249896 = 4.0;
double r249897 = r249895 / r249896;
double r249898 = r249892 - r249897;
double r249899 = c;
double r249900 = r249898 + r249899;
return r249900;
}



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