\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 r189103 = x;
double r189104 = y;
double r189105 = r189103 * r189104;
double r189106 = z;
double r189107 = t;
double r189108 = r189106 * r189107;
double r189109 = 16.0;
double r189110 = r189108 / r189109;
double r189111 = r189105 + r189110;
double r189112 = a;
double r189113 = b;
double r189114 = r189112 * r189113;
double r189115 = 4.0;
double r189116 = r189114 / r189115;
double r189117 = r189111 - r189116;
double r189118 = c;
double r189119 = r189117 + r189118;
return r189119;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r189120 = x;
double r189121 = y;
double r189122 = r189120 * r189121;
double r189123 = z;
double r189124 = t;
double r189125 = r189123 * r189124;
double r189126 = 16.0;
double r189127 = r189125 / r189126;
double r189128 = r189122 + r189127;
double r189129 = a;
double r189130 = b;
double r189131 = r189129 * r189130;
double r189132 = 4.0;
double r189133 = r189131 / r189132;
double r189134 = r189128 - r189133;
double r189135 = c;
double r189136 = r189134 + r189135;
return r189136;
}



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