\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(\sqrt[3]{y}\right)\right) - z\right) - ydouble f(double x, double y, double z) {
double r21297 = x;
double r21298 = y;
double r21299 = log(r21298);
double r21300 = r21297 * r21299;
double r21301 = z;
double r21302 = r21300 - r21301;
double r21303 = r21302 - r21298;
return r21303;
}
double f(double x, double y, double z) {
double r21304 = x;
double r21305 = 2.0;
double r21306 = y;
double r21307 = cbrt(r21306);
double r21308 = log(r21307);
double r21309 = r21305 * r21308;
double r21310 = r21304 * r21308;
double r21311 = fma(r21304, r21309, r21310);
double r21312 = z;
double r21313 = r21311 - r21312;
double r21314 = r21313 - r21306;
return r21314;
}



Bits error versus x



Bits error versus y



Bits error versus z
Initial program 0.1
rmApplied add-cube-cbrt0.1
Applied log-prod0.1
Applied distribute-lft-in0.1
Simplified0.1
rmApplied fma-def0.1
Final simplification0.1
herbie shell --seed 2020036 +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))