\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 r241083 = x;
double r241084 = y;
double r241085 = r241083 * r241084;
double r241086 = z;
double r241087 = t;
double r241088 = r241086 * r241087;
double r241089 = 16.0;
double r241090 = r241088 / r241089;
double r241091 = r241085 + r241090;
double r241092 = a;
double r241093 = b;
double r241094 = r241092 * r241093;
double r241095 = 4.0;
double r241096 = r241094 / r241095;
double r241097 = r241091 - r241096;
double r241098 = c;
double r241099 = r241097 + r241098;
return r241099;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r241100 = x;
double r241101 = y;
double r241102 = r241100 * r241101;
double r241103 = z;
double r241104 = t;
double r241105 = r241103 * r241104;
double r241106 = 16.0;
double r241107 = r241105 / r241106;
double r241108 = r241102 + r241107;
double r241109 = a;
double r241110 = b;
double r241111 = r241109 * r241110;
double r241112 = 4.0;
double r241113 = r241111 / r241112;
double r241114 = r241108 - r241113;
double r241115 = c;
double r241116 = r241114 + r241115;
return r241116;
}



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