\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(\mathsf{fma}\left(\mathsf{fma}\left(y + a, y, b\right), y, c\right) \cdot \sqrt[3]{y}\right) \cdot \sqrt[3]{y}\right) \cdot \sqrt[3]{y} + i}double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return (((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((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 (((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / ((((fma(fma((y + a), y, b), y, c) * cbrt(y)) * cbrt(y)) * cbrt(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 28.8
rmApplied add-cube-cbrt29.0
Applied associate-*r*29.0
Simplified29.0
Final simplification29.0
herbie shell --seed 2020056 +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.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))