\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}\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}double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r49942 = x;
double r49943 = y;
double r49944 = r49942 * r49943;
double r49945 = z;
double r49946 = r49944 + r49945;
double r49947 = r49946 * r49943;
double r49948 = 27464.7644705;
double r49949 = r49947 + r49948;
double r49950 = r49949 * r49943;
double r49951 = 230661.510616;
double r49952 = r49950 + r49951;
double r49953 = r49952 * r49943;
double r49954 = t;
double r49955 = r49953 + r49954;
double r49956 = a;
double r49957 = r49943 + r49956;
double r49958 = r49957 * r49943;
double r49959 = b;
double r49960 = r49958 + r49959;
double r49961 = r49960 * r49943;
double r49962 = c;
double r49963 = r49961 + r49962;
double r49964 = r49963 * r49943;
double r49965 = i;
double r49966 = r49964 + r49965;
double r49967 = r49955 / r49966;
return r49967;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r49968 = x;
double r49969 = y;
double r49970 = r49968 * r49969;
double r49971 = z;
double r49972 = r49970 + r49971;
double r49973 = r49972 * r49969;
double r49974 = 27464.7644705;
double r49975 = r49973 + r49974;
double r49976 = r49975 * r49969;
double r49977 = 230661.510616;
double r49978 = r49976 + r49977;
double r49979 = r49978 * r49969;
double r49980 = t;
double r49981 = r49979 + r49980;
double r49982 = a;
double r49983 = r49969 + r49982;
double r49984 = r49983 * r49969;
double r49985 = b;
double r49986 = r49984 + r49985;
double r49987 = r49986 * r49969;
double r49988 = c;
double r49989 = r49987 + r49988;
double r49990 = r49989 * r49969;
double r49991 = i;
double r49992 = r49990 + r49991;
double r49993 = r49981 / r49992;
return r49993;
}



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.3
Final simplification29.3
herbie shell --seed 2020045
(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)))