\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 r79442 = x;
double r79443 = y;
double r79444 = r79442 * r79443;
double r79445 = z;
double r79446 = r79444 + r79445;
double r79447 = r79446 * r79443;
double r79448 = 27464.7644705;
double r79449 = r79447 + r79448;
double r79450 = r79449 * r79443;
double r79451 = 230661.510616;
double r79452 = r79450 + r79451;
double r79453 = r79452 * r79443;
double r79454 = t;
double r79455 = r79453 + r79454;
double r79456 = a;
double r79457 = r79443 + r79456;
double r79458 = r79457 * r79443;
double r79459 = b;
double r79460 = r79458 + r79459;
double r79461 = r79460 * r79443;
double r79462 = c;
double r79463 = r79461 + r79462;
double r79464 = r79463 * r79443;
double r79465 = i;
double r79466 = r79464 + r79465;
double r79467 = r79455 / r79466;
return r79467;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r79468 = x;
double r79469 = y;
double r79470 = r79468 * r79469;
double r79471 = z;
double r79472 = r79470 + r79471;
double r79473 = r79472 * r79469;
double r79474 = 27464.7644705;
double r79475 = r79473 + r79474;
double r79476 = r79475 * r79469;
double r79477 = 230661.510616;
double r79478 = r79476 + r79477;
double r79479 = r79478 * r79469;
double r79480 = t;
double r79481 = r79479 + r79480;
double r79482 = a;
double r79483 = r79469 + r79482;
double r79484 = r79483 * r79469;
double r79485 = b;
double r79486 = r79484 + r79485;
double r79487 = r79486 * r79469;
double r79488 = c;
double r79489 = r79487 + r79488;
double r79490 = r79489 * r79469;
double r79491 = i;
double r79492 = r79490 + r79491;
double r79493 = r79481 / r79492;
return r79493;
}



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