\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 r184996 = x;
double r184997 = y;
double r184998 = r184996 * r184997;
double r184999 = z;
double r185000 = t;
double r185001 = r184999 * r185000;
double r185002 = 16.0;
double r185003 = r185001 / r185002;
double r185004 = r184998 + r185003;
double r185005 = a;
double r185006 = b;
double r185007 = r185005 * r185006;
double r185008 = 4.0;
double r185009 = r185007 / r185008;
double r185010 = r185004 - r185009;
double r185011 = c;
double r185012 = r185010 + r185011;
return r185012;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r185013 = x;
double r185014 = y;
double r185015 = r185013 * r185014;
double r185016 = z;
double r185017 = t;
double r185018 = r185016 * r185017;
double r185019 = 16.0;
double r185020 = r185018 / r185019;
double r185021 = r185015 + r185020;
double r185022 = a;
double r185023 = b;
double r185024 = r185022 * r185023;
double r185025 = 4.0;
double r185026 = r185024 / r185025;
double r185027 = r185021 - r185026;
double r185028 = c;
double r185029 = r185027 + r185028;
return r185029;
}



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