\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(\left(x \cdot y + z\right) \cdot y + 27464.764470499998\right) \cdot y + 230661.510616000014\right) \cdot y + t\right) \cdot \left(\sqrt[3]{\frac{1}{\left({y}^{2} \cdot \left(y + a\right) + \left(y \cdot b + c\right)\right) \cdot y + i}} \cdot \sqrt[3]{\frac{1}{\left({y}^{2} \cdot \left(y + a\right) + \left(y \cdot b + c\right)\right) \cdot y + i}}\right)\right) \cdot \sqrt[3]{\frac{1}{\left({y}^{2} \cdot \left(y + a\right) + \left(y \cdot b + c\right)\right) \cdot y + i}}double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return (((double) (((double) (((double) (((double) (((double) (((double) (((double) (((double) (x * y)) + z)) * y)) + 27464.7644705)) * y)) + 230661.510616)) * y)) + t)) / ((double) (((double) (((double) (((double) (((double) (((double) (((double) (y + a)) * y)) + b)) * y)) + c)) * y)) + i)));
}
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return ((double) (((double) (((double) (((double) (((double) (((double) (((double) (((double) (((double) (((double) (x * y)) + z)) * y)) + 27464.7644705)) * y)) + 230661.510616)) * y)) + t)) * ((double) (((double) cbrt((1.0 / ((double) (((double) (((double) (((double) (((double) pow(y, 2.0)) * ((double) (y + a)))) + ((double) (((double) (y * b)) + c)))) * y)) + i))))) * ((double) cbrt((1.0 / ((double) (((double) (((double) (((double) (((double) pow(y, 2.0)) * ((double) (y + a)))) + ((double) (((double) (y * b)) + c)))) * y)) + i))))))))) * ((double) cbrt((1.0 / ((double) (((double) (((double) (((double) (((double) pow(y, 2.0)) * ((double) (y + a)))) + ((double) (((double) (y * b)) + c)))) * y)) + i)))))));
}



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
Taylor expanded around inf 29.4
rmApplied div-inv29.4
Simplified29.4
rmApplied add-cube-cbrt29.9
Applied associate-*r*29.9
Final simplification29.9
herbie shell --seed 2020182
(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)))