\frac{\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.7644704999984242022037506103515625\right) \cdot y + 230661.5106160000141244381666183471679688\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.7644704999984242022037506103515625\right) \cdot y + 230661.5106160000141244381666183471679688\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 r61193 = x;
double r61194 = y;
double r61195 = r61193 * r61194;
double r61196 = z;
double r61197 = r61195 + r61196;
double r61198 = r61197 * r61194;
double r61199 = 27464.7644705;
double r61200 = r61198 + r61199;
double r61201 = r61200 * r61194;
double r61202 = 230661.510616;
double r61203 = r61201 + r61202;
double r61204 = r61203 * r61194;
double r61205 = t;
double r61206 = r61204 + r61205;
double r61207 = a;
double r61208 = r61194 + r61207;
double r61209 = r61208 * r61194;
double r61210 = b;
double r61211 = r61209 + r61210;
double r61212 = r61211 * r61194;
double r61213 = c;
double r61214 = r61212 + r61213;
double r61215 = r61214 * r61194;
double r61216 = i;
double r61217 = r61215 + r61216;
double r61218 = r61206 / r61217;
return r61218;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r61219 = x;
double r61220 = y;
double r61221 = r61219 * r61220;
double r61222 = z;
double r61223 = r61221 + r61222;
double r61224 = r61223 * r61220;
double r61225 = 27464.7644705;
double r61226 = r61224 + r61225;
double r61227 = r61226 * r61220;
double r61228 = 230661.510616;
double r61229 = r61227 + r61228;
double r61230 = r61229 * r61220;
double r61231 = t;
double r61232 = r61230 + r61231;
double r61233 = a;
double r61234 = r61220 + r61233;
double r61235 = r61234 * r61220;
double r61236 = b;
double r61237 = r61235 + r61236;
double r61238 = r61237 * r61220;
double r61239 = c;
double r61240 = r61238 + r61239;
double r61241 = r61240 * r61220;
double r61242 = i;
double r61243 = r61241 + r61242;
double r61244 = r61232 / r61243;
return r61244;
}



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 2019208
(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.764470499998) y) 230661.510616000014) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))