\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 r64854 = x;
double r64855 = y;
double r64856 = r64854 * r64855;
double r64857 = z;
double r64858 = r64856 + r64857;
double r64859 = r64858 * r64855;
double r64860 = 27464.7644705;
double r64861 = r64859 + r64860;
double r64862 = r64861 * r64855;
double r64863 = 230661.510616;
double r64864 = r64862 + r64863;
double r64865 = r64864 * r64855;
double r64866 = t;
double r64867 = r64865 + r64866;
double r64868 = a;
double r64869 = r64855 + r64868;
double r64870 = r64869 * r64855;
double r64871 = b;
double r64872 = r64870 + r64871;
double r64873 = r64872 * r64855;
double r64874 = c;
double r64875 = r64873 + r64874;
double r64876 = r64875 * r64855;
double r64877 = i;
double r64878 = r64876 + r64877;
double r64879 = r64867 / r64878;
return r64879;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r64880 = x;
double r64881 = y;
double r64882 = r64880 * r64881;
double r64883 = z;
double r64884 = r64882 + r64883;
double r64885 = r64884 * r64881;
double r64886 = 27464.7644705;
double r64887 = r64885 + r64886;
double r64888 = r64887 * r64881;
double r64889 = 230661.510616;
double r64890 = r64888 + r64889;
double r64891 = r64890 * r64881;
double r64892 = t;
double r64893 = r64891 + r64892;
double r64894 = a;
double r64895 = r64881 + r64894;
double r64896 = r64895 * r64881;
double r64897 = b;
double r64898 = r64896 + r64897;
double r64899 = r64898 * r64881;
double r64900 = c;
double r64901 = r64899 + r64900;
double r64902 = r64901 * r64881;
double r64903 = i;
double r64904 = r64902 + r64903;
double r64905 = r64893 / r64904;
return r64905;
}



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 +o rules:numerics
(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)))