\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 r323931 = x;
double r323932 = y;
double r323933 = r323931 * r323932;
double r323934 = z;
double r323935 = t;
double r323936 = r323934 * r323935;
double r323937 = 16.0;
double r323938 = r323936 / r323937;
double r323939 = r323933 + r323938;
double r323940 = a;
double r323941 = b;
double r323942 = r323940 * r323941;
double r323943 = 4.0;
double r323944 = r323942 / r323943;
double r323945 = r323939 - r323944;
double r323946 = c;
double r323947 = r323945 + r323946;
return r323947;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r323948 = x;
double r323949 = y;
double r323950 = r323948 * r323949;
double r323951 = z;
double r323952 = t;
double r323953 = r323951 * r323952;
double r323954 = 16.0;
double r323955 = r323953 / r323954;
double r323956 = r323950 + r323955;
double r323957 = a;
double r323958 = b;
double r323959 = r323957 * r323958;
double r323960 = 4.0;
double r323961 = r323959 / r323960;
double r323962 = r323956 - r323961;
double r323963 = c;
double r323964 = r323962 + r323963;
return r323964;
}



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