\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 r157767 = x;
double r157768 = y;
double r157769 = r157767 * r157768;
double r157770 = z;
double r157771 = t;
double r157772 = r157770 * r157771;
double r157773 = 16.0;
double r157774 = r157772 / r157773;
double r157775 = r157769 + r157774;
double r157776 = a;
double r157777 = b;
double r157778 = r157776 * r157777;
double r157779 = 4.0;
double r157780 = r157778 / r157779;
double r157781 = r157775 - r157780;
double r157782 = c;
double r157783 = r157781 + r157782;
return r157783;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r157784 = x;
double r157785 = y;
double r157786 = r157784 * r157785;
double r157787 = z;
double r157788 = t;
double r157789 = r157787 * r157788;
double r157790 = 16.0;
double r157791 = r157789 / r157790;
double r157792 = r157786 + r157791;
double r157793 = a;
double r157794 = b;
double r157795 = r157793 * r157794;
double r157796 = 4.0;
double r157797 = r157795 / r157796;
double r157798 = r157792 - r157797;
double r157799 = c;
double r157800 = r157798 + r157799;
return r157800;
}



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