\frac{\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.7644705\right) \cdot y + 230661.510616\right) \cdot y + t}{\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}\left(t + \left(y \cdot \left(y \cdot \left(z + x \cdot y\right) + 27464.7644705\right) + 230661.510616\right) \cdot y\right) \cdot \frac{1}{i + y \cdot \left(\left(\left(a + y\right) \cdot y + b\right) \cdot y + c\right)}double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r13201204 = x;
double r13201205 = y;
double r13201206 = r13201204 * r13201205;
double r13201207 = z;
double r13201208 = r13201206 + r13201207;
double r13201209 = r13201208 * r13201205;
double r13201210 = 27464.7644705;
double r13201211 = r13201209 + r13201210;
double r13201212 = r13201211 * r13201205;
double r13201213 = 230661.510616;
double r13201214 = r13201212 + r13201213;
double r13201215 = r13201214 * r13201205;
double r13201216 = t;
double r13201217 = r13201215 + r13201216;
double r13201218 = a;
double r13201219 = r13201205 + r13201218;
double r13201220 = r13201219 * r13201205;
double r13201221 = b;
double r13201222 = r13201220 + r13201221;
double r13201223 = r13201222 * r13201205;
double r13201224 = c;
double r13201225 = r13201223 + r13201224;
double r13201226 = r13201225 * r13201205;
double r13201227 = i;
double r13201228 = r13201226 + r13201227;
double r13201229 = r13201217 / r13201228;
return r13201229;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r13201230 = t;
double r13201231 = y;
double r13201232 = z;
double r13201233 = x;
double r13201234 = r13201233 * r13201231;
double r13201235 = r13201232 + r13201234;
double r13201236 = r13201231 * r13201235;
double r13201237 = 27464.7644705;
double r13201238 = r13201236 + r13201237;
double r13201239 = r13201231 * r13201238;
double r13201240 = 230661.510616;
double r13201241 = r13201239 + r13201240;
double r13201242 = r13201241 * r13201231;
double r13201243 = r13201230 + r13201242;
double r13201244 = 1.0;
double r13201245 = i;
double r13201246 = a;
double r13201247 = r13201246 + r13201231;
double r13201248 = r13201247 * r13201231;
double r13201249 = b;
double r13201250 = r13201248 + r13201249;
double r13201251 = r13201250 * r13201231;
double r13201252 = c;
double r13201253 = r13201251 + r13201252;
double r13201254 = r13201231 * r13201253;
double r13201255 = r13201245 + r13201254;
double r13201256 = r13201244 / r13201255;
double r13201257 = r13201243 * r13201256;
return r13201257;
}



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.8
rmApplied div-inv28.9
Final simplification28.9
herbie shell --seed 2019104
(FPCore (x y z t a b c i)
:name "Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2"
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))