\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 r228016 = x;
double r228017 = y;
double r228018 = r228016 * r228017;
double r228019 = z;
double r228020 = t;
double r228021 = r228019 * r228020;
double r228022 = 16.0;
double r228023 = r228021 / r228022;
double r228024 = r228018 + r228023;
double r228025 = a;
double r228026 = b;
double r228027 = r228025 * r228026;
double r228028 = 4.0;
double r228029 = r228027 / r228028;
double r228030 = r228024 - r228029;
double r228031 = c;
double r228032 = r228030 + r228031;
return r228032;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r228033 = x;
double r228034 = y;
double r228035 = r228033 * r228034;
double r228036 = z;
double r228037 = t;
double r228038 = r228036 * r228037;
double r228039 = 16.0;
double r228040 = r228038 / r228039;
double r228041 = r228035 + r228040;
double r228042 = a;
double r228043 = b;
double r228044 = r228042 * r228043;
double r228045 = 4.0;
double r228046 = r228044 / r228045;
double r228047 = r228041 - r228046;
double r228048 = c;
double r228049 = r228047 + r228048;
return r228049;
}



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