\left(x \cdot \log y - z\right) - y
\left(\mathsf{fma}\left(x, 2 \cdot \log \left(\sqrt[3]{y}\right), x \cdot \log \left(e^{\log \left({y}^{\frac{1}{3}}\right)}\right)\right) - z\right) - ydouble code(double x, double y, double z) {
return (((x * log(y)) - z) - y);
}
double code(double x, double y, double z) {
return ((fma(x, (2.0 * log(cbrt(y))), (x * log(exp(log(pow(y, 0.3333333333333333)))))) - z) - y);
}



Bits error versus x



Bits error versus y



Bits error versus z
Results
Initial program 0.1
rmApplied add-cube-cbrt0.1
Applied log-prod0.1
Applied distribute-lft-in0.1
Simplified0.1
rmApplied fma-def0.1
rmApplied add-exp-log0.1
Simplified0.1
Final simplification0.1
herbie shell --seed 2020049 +o rules:numerics
(FPCore (x y z)
:name "Statistics.Distribution.Poisson:$clogProbability from math-functions-0.1.5.2"
:precision binary64
(- (- (* x (log y)) z) y))