\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 r3145924 = x;
double r3145925 = y;
double r3145926 = r3145924 * r3145925;
double r3145927 = z;
double r3145928 = r3145926 + r3145927;
double r3145929 = r3145928 * r3145925;
double r3145930 = 27464.7644705;
double r3145931 = r3145929 + r3145930;
double r3145932 = r3145931 * r3145925;
double r3145933 = 230661.510616;
double r3145934 = r3145932 + r3145933;
double r3145935 = r3145934 * r3145925;
double r3145936 = t;
double r3145937 = r3145935 + r3145936;
double r3145938 = a;
double r3145939 = r3145925 + r3145938;
double r3145940 = r3145939 * r3145925;
double r3145941 = b;
double r3145942 = r3145940 + r3145941;
double r3145943 = r3145942 * r3145925;
double r3145944 = c;
double r3145945 = r3145943 + r3145944;
double r3145946 = r3145945 * r3145925;
double r3145947 = i;
double r3145948 = r3145946 + r3145947;
double r3145949 = r3145937 / r3145948;
return r3145949;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r3145950 = t;
double r3145951 = y;
double r3145952 = z;
double r3145953 = x;
double r3145954 = r3145953 * r3145951;
double r3145955 = r3145952 + r3145954;
double r3145956 = r3145951 * r3145955;
double r3145957 = 27464.7644705;
double r3145958 = r3145956 + r3145957;
double r3145959 = r3145951 * r3145958;
double r3145960 = 230661.510616;
double r3145961 = r3145959 + r3145960;
double r3145962 = r3145961 * r3145951;
double r3145963 = r3145950 + r3145962;
double r3145964 = 1.0;
double r3145965 = i;
double r3145966 = a;
double r3145967 = r3145966 + r3145951;
double r3145968 = r3145967 * r3145951;
double r3145969 = b;
double r3145970 = r3145968 + r3145969;
double r3145971 = r3145970 * r3145951;
double r3145972 = c;
double r3145973 = r3145971 + r3145972;
double r3145974 = r3145951 * r3145973;
double r3145975 = r3145965 + r3145974;
double r3145976 = r3145964 / r3145975;
double r3145977 = r3145963 * r3145976;
return r3145977;
}



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)))