x \cdot 0.5 + y \cdot \left(\left(1 - z\right) + \log z\right)
\left(\left(\mathsf{fma}\left(2, \log \left(\sqrt[3]{z}\right), 1 - z\right) \cdot y + \log \left({\left(\sqrt{z}\right)}^{\frac{1}{3}}\right) \cdot y\right) + \log \left({\left(\sqrt{z}\right)}^{\frac{1}{3}}\right) \cdot y\right) + x \cdot 0.5double f(double x, double y, double z) {
double r12020357 = x;
double r12020358 = 0.5;
double r12020359 = r12020357 * r12020358;
double r12020360 = y;
double r12020361 = 1.0;
double r12020362 = z;
double r12020363 = r12020361 - r12020362;
double r12020364 = log(r12020362);
double r12020365 = r12020363 + r12020364;
double r12020366 = r12020360 * r12020365;
double r12020367 = r12020359 + r12020366;
return r12020367;
}
double f(double x, double y, double z) {
double r12020368 = 2.0;
double r12020369 = z;
double r12020370 = cbrt(r12020369);
double r12020371 = log(r12020370);
double r12020372 = 1.0;
double r12020373 = r12020372 - r12020369;
double r12020374 = fma(r12020368, r12020371, r12020373);
double r12020375 = y;
double r12020376 = r12020374 * r12020375;
double r12020377 = sqrt(r12020369);
double r12020378 = 0.3333333333333333;
double r12020379 = pow(r12020377, r12020378);
double r12020380 = log(r12020379);
double r12020381 = r12020380 * r12020375;
double r12020382 = r12020376 + r12020381;
double r12020383 = r12020382 + r12020381;
double r12020384 = x;
double r12020385 = 0.5;
double r12020386 = r12020384 * r12020385;
double r12020387 = r12020383 + r12020386;
return r12020387;
}




Bits error versus x




Bits error versus y




Bits error versus z
| Original | 0.1 |
|---|---|
| Target | 0.1 |
| Herbie | 0.1 |
Initial program 0.1
rmApplied distribute-lft-in0.1
rmApplied add-cube-cbrt0.1
Applied log-prod0.1
Applied distribute-lft-in0.1
Applied associate-+r+0.1
Simplified0.1
rmApplied pow1/30.1
rmApplied add-sqr-sqrt0.1
Applied unpow-prod-down0.1
Applied log-prod0.1
Applied distribute-lft-in0.1
Applied associate-+r+0.1
Final simplification0.1
herbie shell --seed 2019179 +o rules:numerics
(FPCore (x y z)
:name "System.Random.MWC.Distributions:gamma from mwc-random-0.13.3.2"
:herbie-target
(- (+ y (* 0.5 x)) (* y (- z (log z))))
(+ (* x 0.5) (* y (+ (- 1.0 z) (log z)))))