\left(\left(x \cdot \log y - y\right) - z\right) + \log t
\mathsf{fma}\left(\log y, x, \log t - \left(y + z\right)\right) + \mathsf{fma}\left(-\sqrt[3]{y + z}, \sqrt[3]{y + z} \cdot \sqrt[3]{y + z}, \sqrt[3]{y + z} \cdot \left(\sqrt[3]{y + z} \cdot \sqrt[3]{y + z}\right)\right)double code(double x, double y, double z, double t) {
return ((((x * log(y)) - y) - z) + log(t));
}
double code(double x, double y, double z, double t) {
return (fma(log(y), x, (log(t) - (y + z))) + fma(-cbrt((y + z)), (cbrt((y + z)) * cbrt((y + z))), (cbrt((y + z)) * (cbrt((y + z)) * cbrt((y + z))))));
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t
Results
Initial program 0.1
Simplified0.1
rmApplied fma-udef0.1
rmApplied add-cube-cbrt0.8
Applied pow10.8
Applied log-pow0.8
Applied prod-diff0.8
Applied associate-+r+0.8
Simplified0.1
Final simplification0.1
herbie shell --seed 2020075 +o rules:numerics
(FPCore (x y z t)
:name "Numeric.SpecFunctions:incompleteGamma from math-functions-0.1.5.2, A"
:precision binary64
(+ (- (- (* x (log y)) y) z) (log t)))