\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 r202884 = x;
double r202885 = y;
double r202886 = r202884 * r202885;
double r202887 = z;
double r202888 = t;
double r202889 = r202887 * r202888;
double r202890 = 16.0;
double r202891 = r202889 / r202890;
double r202892 = r202886 + r202891;
double r202893 = a;
double r202894 = b;
double r202895 = r202893 * r202894;
double r202896 = 4.0;
double r202897 = r202895 / r202896;
double r202898 = r202892 - r202897;
double r202899 = c;
double r202900 = r202898 + r202899;
return r202900;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r202901 = x;
double r202902 = y;
double r202903 = r202901 * r202902;
double r202904 = z;
double r202905 = t;
double r202906 = r202904 * r202905;
double r202907 = 16.0;
double r202908 = r202906 / r202907;
double r202909 = r202903 + r202908;
double r202910 = a;
double r202911 = b;
double r202912 = r202910 * r202911;
double r202913 = 4.0;
double r202914 = r202912 / r202913;
double r202915 = r202909 - r202914;
double r202916 = c;
double r202917 = r202915 + r202916;
return r202917;
}



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