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



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t



Bits error versus a
Results
Initial program 0.2
Simplified0.2
rmApplied add-cube-cbrt0.2
Applied log-prod0.3
Applied associate-+l+0.3
Final simplification0.3
herbie shell --seed 2020100 +o rules:numerics
(FPCore (x y z t a)
:name "Numeric.SpecFunctions:logGammaL from math-functions-0.1.5.2"
:precision binary64
(+ (- (+ (log (+ x y)) (log z)) t) (* (- a 0.5) (log t))))