\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 r3153608 = x;
double r3153609 = y;
double r3153610 = r3153608 * r3153609;
double r3153611 = z;
double r3153612 = r3153610 + r3153611;
double r3153613 = r3153612 * r3153609;
double r3153614 = 27464.7644705;
double r3153615 = r3153613 + r3153614;
double r3153616 = r3153615 * r3153609;
double r3153617 = 230661.510616;
double r3153618 = r3153616 + r3153617;
double r3153619 = r3153618 * r3153609;
double r3153620 = t;
double r3153621 = r3153619 + r3153620;
double r3153622 = a;
double r3153623 = r3153609 + r3153622;
double r3153624 = r3153623 * r3153609;
double r3153625 = b;
double r3153626 = r3153624 + r3153625;
double r3153627 = r3153626 * r3153609;
double r3153628 = c;
double r3153629 = r3153627 + r3153628;
double r3153630 = r3153629 * r3153609;
double r3153631 = i;
double r3153632 = r3153630 + r3153631;
double r3153633 = r3153621 / r3153632;
return r3153633;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r3153634 = t;
double r3153635 = y;
double r3153636 = z;
double r3153637 = x;
double r3153638 = r3153637 * r3153635;
double r3153639 = r3153636 + r3153638;
double r3153640 = r3153635 * r3153639;
double r3153641 = 27464.7644705;
double r3153642 = r3153640 + r3153641;
double r3153643 = r3153635 * r3153642;
double r3153644 = 230661.510616;
double r3153645 = r3153643 + r3153644;
double r3153646 = r3153645 * r3153635;
double r3153647 = r3153634 + r3153646;
double r3153648 = 1.0;
double r3153649 = i;
double r3153650 = a;
double r3153651 = r3153650 + r3153635;
double r3153652 = r3153651 * r3153635;
double r3153653 = b;
double r3153654 = r3153652 + r3153653;
double r3153655 = r3153654 * r3153635;
double r3153656 = c;
double r3153657 = r3153655 + r3153656;
double r3153658 = r3153635 * r3153657;
double r3153659 = r3153649 + r3153658;
double r3153660 = r3153648 / r3153659;
double r3153661 = r3153647 * r3153660;
return r3153661;
}



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.2
rmApplied div-inv28.3
Final simplification28.3
herbie shell --seed 2019139
(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)))