\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 r186715 = x;
double r186716 = y;
double r186717 = r186715 * r186716;
double r186718 = z;
double r186719 = t;
double r186720 = r186718 * r186719;
double r186721 = 16.0;
double r186722 = r186720 / r186721;
double r186723 = r186717 + r186722;
double r186724 = a;
double r186725 = b;
double r186726 = r186724 * r186725;
double r186727 = 4.0;
double r186728 = r186726 / r186727;
double r186729 = r186723 - r186728;
double r186730 = c;
double r186731 = r186729 + r186730;
return r186731;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r186732 = x;
double r186733 = y;
double r186734 = r186732 * r186733;
double r186735 = z;
double r186736 = t;
double r186737 = r186735 * r186736;
double r186738 = 16.0;
double r186739 = r186737 / r186738;
double r186740 = r186734 + r186739;
double r186741 = a;
double r186742 = b;
double r186743 = r186741 * r186742;
double r186744 = 4.0;
double r186745 = r186743 / r186744;
double r186746 = r186740 - r186745;
double r186747 = c;
double r186748 = r186746 + r186747;
return r186748;
}



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