\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 r219837 = x;
double r219838 = y;
double r219839 = r219837 * r219838;
double r219840 = z;
double r219841 = t;
double r219842 = r219840 * r219841;
double r219843 = 16.0;
double r219844 = r219842 / r219843;
double r219845 = r219839 + r219844;
double r219846 = a;
double r219847 = b;
double r219848 = r219846 * r219847;
double r219849 = 4.0;
double r219850 = r219848 / r219849;
double r219851 = r219845 - r219850;
double r219852 = c;
double r219853 = r219851 + r219852;
return r219853;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r219854 = x;
double r219855 = y;
double r219856 = r219854 * r219855;
double r219857 = z;
double r219858 = t;
double r219859 = r219857 * r219858;
double r219860 = 16.0;
double r219861 = r219859 / r219860;
double r219862 = r219856 + r219861;
double r219863 = a;
double r219864 = b;
double r219865 = r219863 * r219864;
double r219866 = 4.0;
double r219867 = r219865 / r219866;
double r219868 = r219862 - r219867;
double r219869 = c;
double r219870 = r219868 + r219869;
return r219870;
}



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