\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 r163864 = x;
double r163865 = y;
double r163866 = r163864 * r163865;
double r163867 = z;
double r163868 = t;
double r163869 = r163867 * r163868;
double r163870 = 16.0;
double r163871 = r163869 / r163870;
double r163872 = r163866 + r163871;
double r163873 = a;
double r163874 = b;
double r163875 = r163873 * r163874;
double r163876 = 4.0;
double r163877 = r163875 / r163876;
double r163878 = r163872 - r163877;
double r163879 = c;
double r163880 = r163878 + r163879;
return r163880;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r163881 = x;
double r163882 = y;
double r163883 = r163881 * r163882;
double r163884 = z;
double r163885 = t;
double r163886 = r163884 * r163885;
double r163887 = 16.0;
double r163888 = r163886 / r163887;
double r163889 = r163883 + r163888;
double r163890 = a;
double r163891 = b;
double r163892 = r163890 * r163891;
double r163893 = 4.0;
double r163894 = r163892 / r163893;
double r163895 = r163889 - r163894;
double r163896 = c;
double r163897 = r163895 + r163896;
return r163897;
}



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