\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 r206625 = x;
double r206626 = y;
double r206627 = r206625 * r206626;
double r206628 = z;
double r206629 = t;
double r206630 = r206628 * r206629;
double r206631 = 16.0;
double r206632 = r206630 / r206631;
double r206633 = r206627 + r206632;
double r206634 = a;
double r206635 = b;
double r206636 = r206634 * r206635;
double r206637 = 4.0;
double r206638 = r206636 / r206637;
double r206639 = r206633 - r206638;
double r206640 = c;
double r206641 = r206639 + r206640;
return r206641;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r206642 = x;
double r206643 = y;
double r206644 = r206642 * r206643;
double r206645 = z;
double r206646 = t;
double r206647 = r206645 * r206646;
double r206648 = 16.0;
double r206649 = r206647 / r206648;
double r206650 = r206644 + r206649;
double r206651 = a;
double r206652 = b;
double r206653 = r206651 * r206652;
double r206654 = 4.0;
double r206655 = r206653 / r206654;
double r206656 = r206650 - r206655;
double r206657 = c;
double r206658 = r206656 + r206657;
return r206658;
}



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