\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 r59441 = x;
double r59442 = y;
double r59443 = r59441 * r59442;
double r59444 = z;
double r59445 = r59443 + r59444;
double r59446 = r59445 * r59442;
double r59447 = 27464.7644705;
double r59448 = r59446 + r59447;
double r59449 = r59448 * r59442;
double r59450 = 230661.510616;
double r59451 = r59449 + r59450;
double r59452 = r59451 * r59442;
double r59453 = t;
double r59454 = r59452 + r59453;
double r59455 = a;
double r59456 = r59442 + r59455;
double r59457 = r59456 * r59442;
double r59458 = b;
double r59459 = r59457 + r59458;
double r59460 = r59459 * r59442;
double r59461 = c;
double r59462 = r59460 + r59461;
double r59463 = r59462 * r59442;
double r59464 = i;
double r59465 = r59463 + r59464;
double r59466 = r59454 / r59465;
return r59466;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r59467 = x;
double r59468 = y;
double r59469 = r59467 * r59468;
double r59470 = z;
double r59471 = r59469 + r59470;
double r59472 = r59471 * r59468;
double r59473 = 27464.7644705;
double r59474 = r59472 + r59473;
double r59475 = r59474 * r59468;
double r59476 = 230661.510616;
double r59477 = r59475 + r59476;
double r59478 = r59477 * r59468;
double r59479 = t;
double r59480 = r59478 + r59479;
double r59481 = a;
double r59482 = r59468 + r59481;
double r59483 = r59482 * r59468;
double r59484 = b;
double r59485 = r59483 + r59484;
double r59486 = r59485 * r59468;
double r59487 = c;
double r59488 = r59486 + r59487;
double r59489 = r59488 * r59468;
double r59490 = i;
double r59491 = r59489 + r59490;
double r59492 = r59480 / r59491;
return r59492;
}



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 29.0
Final simplification29.0
herbie shell --seed 2019305
(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)))