\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 r4060993 = x;
double r4060994 = y;
double r4060995 = r4060993 * r4060994;
double r4060996 = z;
double r4060997 = r4060995 + r4060996;
double r4060998 = r4060997 * r4060994;
double r4060999 = 27464.7644705;
double r4061000 = r4060998 + r4060999;
double r4061001 = r4061000 * r4060994;
double r4061002 = 230661.510616;
double r4061003 = r4061001 + r4061002;
double r4061004 = r4061003 * r4060994;
double r4061005 = t;
double r4061006 = r4061004 + r4061005;
double r4061007 = a;
double r4061008 = r4060994 + r4061007;
double r4061009 = r4061008 * r4060994;
double r4061010 = b;
double r4061011 = r4061009 + r4061010;
double r4061012 = r4061011 * r4060994;
double r4061013 = c;
double r4061014 = r4061012 + r4061013;
double r4061015 = r4061014 * r4060994;
double r4061016 = i;
double r4061017 = r4061015 + r4061016;
double r4061018 = r4061006 / r4061017;
return r4061018;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r4061019 = t;
double r4061020 = y;
double r4061021 = z;
double r4061022 = x;
double r4061023 = r4061022 * r4061020;
double r4061024 = r4061021 + r4061023;
double r4061025 = r4061020 * r4061024;
double r4061026 = 27464.7644705;
double r4061027 = r4061025 + r4061026;
double r4061028 = r4061020 * r4061027;
double r4061029 = 230661.510616;
double r4061030 = r4061028 + r4061029;
double r4061031 = r4061030 * r4061020;
double r4061032 = r4061019 + r4061031;
double r4061033 = 1.0;
double r4061034 = i;
double r4061035 = a;
double r4061036 = r4061035 + r4061020;
double r4061037 = r4061036 * r4061020;
double r4061038 = b;
double r4061039 = r4061037 + r4061038;
double r4061040 = r4061039 * r4061020;
double r4061041 = c;
double r4061042 = r4061040 + r4061041;
double r4061043 = r4061020 * r4061042;
double r4061044 = r4061034 + r4061043;
double r4061045 = r4061033 / r4061044;
double r4061046 = r4061032 * r4061045;
return r4061046;
}



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