\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 r235964 = x;
double r235965 = y;
double r235966 = r235964 * r235965;
double r235967 = z;
double r235968 = t;
double r235969 = r235967 * r235968;
double r235970 = 16.0;
double r235971 = r235969 / r235970;
double r235972 = r235966 + r235971;
double r235973 = a;
double r235974 = b;
double r235975 = r235973 * r235974;
double r235976 = 4.0;
double r235977 = r235975 / r235976;
double r235978 = r235972 - r235977;
double r235979 = c;
double r235980 = r235978 + r235979;
return r235980;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r235981 = x;
double r235982 = y;
double r235983 = r235981 * r235982;
double r235984 = z;
double r235985 = t;
double r235986 = r235984 * r235985;
double r235987 = 16.0;
double r235988 = r235986 / r235987;
double r235989 = r235983 + r235988;
double r235990 = a;
double r235991 = b;
double r235992 = r235990 * r235991;
double r235993 = 4.0;
double r235994 = r235992 / r235993;
double r235995 = r235989 - r235994;
double r235996 = c;
double r235997 = r235995 + r235996;
return r235997;
}



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