\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 r158107 = x;
double r158108 = y;
double r158109 = r158107 * r158108;
double r158110 = z;
double r158111 = t;
double r158112 = r158110 * r158111;
double r158113 = 16.0;
double r158114 = r158112 / r158113;
double r158115 = r158109 + r158114;
double r158116 = a;
double r158117 = b;
double r158118 = r158116 * r158117;
double r158119 = 4.0;
double r158120 = r158118 / r158119;
double r158121 = r158115 - r158120;
double r158122 = c;
double r158123 = r158121 + r158122;
return r158123;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r158124 = x;
double r158125 = y;
double r158126 = r158124 * r158125;
double r158127 = z;
double r158128 = t;
double r158129 = r158127 * r158128;
double r158130 = 16.0;
double r158131 = r158129 / r158130;
double r158132 = r158126 + r158131;
double r158133 = a;
double r158134 = b;
double r158135 = r158133 * r158134;
double r158136 = 4.0;
double r158137 = r158135 / r158136;
double r158138 = r158132 - r158137;
double r158139 = c;
double r158140 = r158138 + r158139;
return r158140;
}



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