\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 r2521271 = x;
double r2521272 = y;
double r2521273 = r2521271 * r2521272;
double r2521274 = z;
double r2521275 = r2521273 + r2521274;
double r2521276 = r2521275 * r2521272;
double r2521277 = 27464.7644705;
double r2521278 = r2521276 + r2521277;
double r2521279 = r2521278 * r2521272;
double r2521280 = 230661.510616;
double r2521281 = r2521279 + r2521280;
double r2521282 = r2521281 * r2521272;
double r2521283 = t;
double r2521284 = r2521282 + r2521283;
double r2521285 = a;
double r2521286 = r2521272 + r2521285;
double r2521287 = r2521286 * r2521272;
double r2521288 = b;
double r2521289 = r2521287 + r2521288;
double r2521290 = r2521289 * r2521272;
double r2521291 = c;
double r2521292 = r2521290 + r2521291;
double r2521293 = r2521292 * r2521272;
double r2521294 = i;
double r2521295 = r2521293 + r2521294;
double r2521296 = r2521284 / r2521295;
return r2521296;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r2521297 = t;
double r2521298 = y;
double r2521299 = z;
double r2521300 = x;
double r2521301 = r2521300 * r2521298;
double r2521302 = r2521299 + r2521301;
double r2521303 = r2521298 * r2521302;
double r2521304 = 27464.7644705;
double r2521305 = r2521303 + r2521304;
double r2521306 = r2521298 * r2521305;
double r2521307 = 230661.510616;
double r2521308 = r2521306 + r2521307;
double r2521309 = r2521308 * r2521298;
double r2521310 = r2521297 + r2521309;
double r2521311 = 1.0;
double r2521312 = i;
double r2521313 = a;
double r2521314 = r2521313 + r2521298;
double r2521315 = r2521314 * r2521298;
double r2521316 = b;
double r2521317 = r2521315 + r2521316;
double r2521318 = r2521317 * r2521298;
double r2521319 = c;
double r2521320 = r2521318 + r2521319;
double r2521321 = r2521298 * r2521320;
double r2521322 = r2521312 + r2521321;
double r2521323 = r2521311 / r2521322;
double r2521324 = r2521310 * r2521323;
return r2521324;
}



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.5
rmApplied div-inv28.6
Final simplification28.6
herbie shell --seed 2019132
(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)))