\frac{\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.764470499998\right) \cdot y + 230661.510616000014\right) \cdot y + t}{\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}\frac{\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.764470499998\right) \cdot y + 230661.510616000014\right) \cdot y + t}{\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r95171 = x;
double r95172 = y;
double r95173 = r95171 * r95172;
double r95174 = z;
double r95175 = r95173 + r95174;
double r95176 = r95175 * r95172;
double r95177 = 27464.7644705;
double r95178 = r95176 + r95177;
double r95179 = r95178 * r95172;
double r95180 = 230661.510616;
double r95181 = r95179 + r95180;
double r95182 = r95181 * r95172;
double r95183 = t;
double r95184 = r95182 + r95183;
double r95185 = a;
double r95186 = r95172 + r95185;
double r95187 = r95186 * r95172;
double r95188 = b;
double r95189 = r95187 + r95188;
double r95190 = r95189 * r95172;
double r95191 = c;
double r95192 = r95190 + r95191;
double r95193 = r95192 * r95172;
double r95194 = i;
double r95195 = r95193 + r95194;
double r95196 = r95184 / r95195;
return r95196;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r95197 = x;
double r95198 = y;
double r95199 = r95197 * r95198;
double r95200 = z;
double r95201 = r95199 + r95200;
double r95202 = r95201 * r95198;
double r95203 = 27464.7644705;
double r95204 = r95202 + r95203;
double r95205 = r95204 * r95198;
double r95206 = 230661.510616;
double r95207 = r95205 + r95206;
double r95208 = r95207 * r95198;
double r95209 = t;
double r95210 = r95208 + r95209;
double r95211 = a;
double r95212 = r95198 + r95211;
double r95213 = r95212 * r95198;
double r95214 = b;
double r95215 = r95213 + r95214;
double r95216 = r95215 * r95198;
double r95217 = c;
double r95218 = r95216 + r95217;
double r95219 = r95218 * r95198;
double r95220 = i;
double r95221 = r95219 + r95220;
double r95222 = r95210 / r95221;
return r95222;
}



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



Bits error versus i
Results
Initial program 28.7
Final simplification28.7
herbie shell --seed 2020027
(FPCore (x y z t a b c i)
:name "Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2"
:precision binary64
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))