\frac{\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.764470499998\right) \cdot y + 230661.510616000014\right) \cdot y + t}{\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}\left(\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.764470499998\right) \cdot y + 230661.510616000014\right) \cdot y + t\right) \cdot \frac{1}{\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 r58812 = x;
double r58813 = y;
double r58814 = r58812 * r58813;
double r58815 = z;
double r58816 = r58814 + r58815;
double r58817 = r58816 * r58813;
double r58818 = 27464.7644705;
double r58819 = r58817 + r58818;
double r58820 = r58819 * r58813;
double r58821 = 230661.510616;
double r58822 = r58820 + r58821;
double r58823 = r58822 * r58813;
double r58824 = t;
double r58825 = r58823 + r58824;
double r58826 = a;
double r58827 = r58813 + r58826;
double r58828 = r58827 * r58813;
double r58829 = b;
double r58830 = r58828 + r58829;
double r58831 = r58830 * r58813;
double r58832 = c;
double r58833 = r58831 + r58832;
double r58834 = r58833 * r58813;
double r58835 = i;
double r58836 = r58834 + r58835;
double r58837 = r58825 / r58836;
return r58837;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r58838 = x;
double r58839 = y;
double r58840 = r58838 * r58839;
double r58841 = z;
double r58842 = r58840 + r58841;
double r58843 = r58842 * r58839;
double r58844 = 27464.7644705;
double r58845 = r58843 + r58844;
double r58846 = r58845 * r58839;
double r58847 = 230661.510616;
double r58848 = r58846 + r58847;
double r58849 = r58848 * r58839;
double r58850 = t;
double r58851 = r58849 + r58850;
double r58852 = 1.0;
double r58853 = a;
double r58854 = r58839 + r58853;
double r58855 = r58854 * r58839;
double r58856 = b;
double r58857 = r58855 + r58856;
double r58858 = r58857 * r58839;
double r58859 = c;
double r58860 = r58858 + r58859;
double r58861 = r58860 * r58839;
double r58862 = i;
double r58863 = r58861 + r58862;
double r58864 = r58852 / r58863;
double r58865 = r58851 * r58864;
return r58865;
}



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.2
rmApplied div-inv29.3
Final simplification29.3
herbie shell --seed 2020064
(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.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))