\left(x \cdot \log y - z\right) - y
\mathsf{fma}\left(x, \log y, -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, log(y), -z) - y);
}



Bits error versus x



Bits error versus y



Bits error versus z
Results
Initial program 0.1
rmApplied fma-neg0.1
Final simplification0.1
herbie shell --seed 2020056 +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))