\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 r3145809 = x;
double r3145810 = y;
double r3145811 = r3145809 * r3145810;
double r3145812 = z;
double r3145813 = r3145811 + r3145812;
double r3145814 = r3145813 * r3145810;
double r3145815 = 27464.7644705;
double r3145816 = r3145814 + r3145815;
double r3145817 = r3145816 * r3145810;
double r3145818 = 230661.510616;
double r3145819 = r3145817 + r3145818;
double r3145820 = r3145819 * r3145810;
double r3145821 = t;
double r3145822 = r3145820 + r3145821;
double r3145823 = a;
double r3145824 = r3145810 + r3145823;
double r3145825 = r3145824 * r3145810;
double r3145826 = b;
double r3145827 = r3145825 + r3145826;
double r3145828 = r3145827 * r3145810;
double r3145829 = c;
double r3145830 = r3145828 + r3145829;
double r3145831 = r3145830 * r3145810;
double r3145832 = i;
double r3145833 = r3145831 + r3145832;
double r3145834 = r3145822 / r3145833;
return r3145834;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r3145835 = t;
double r3145836 = y;
double r3145837 = z;
double r3145838 = x;
double r3145839 = r3145838 * r3145836;
double r3145840 = r3145837 + r3145839;
double r3145841 = r3145836 * r3145840;
double r3145842 = 27464.7644705;
double r3145843 = r3145841 + r3145842;
double r3145844 = r3145836 * r3145843;
double r3145845 = 230661.510616;
double r3145846 = r3145844 + r3145845;
double r3145847 = r3145846 * r3145836;
double r3145848 = r3145835 + r3145847;
double r3145849 = 1.0;
double r3145850 = i;
double r3145851 = a;
double r3145852 = r3145851 + r3145836;
double r3145853 = r3145852 * r3145836;
double r3145854 = b;
double r3145855 = r3145853 + r3145854;
double r3145856 = r3145855 * r3145836;
double r3145857 = c;
double r3145858 = r3145856 + r3145857;
double r3145859 = r3145836 * r3145858;
double r3145860 = r3145850 + r3145859;
double r3145861 = r3145849 / r3145860;
double r3145862 = r3145848 * r3145861;
return r3145862;
}



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