x \cdot \log \left(\frac{x}{y}\right) - zx \cdot \left(\log \left(\frac{1}{\sqrt[3]{y} \cdot \sqrt[3]{y}}\right) + \left(2 \cdot \log \left(\frac{\sqrt[3]{x}}{\sqrt[3]{\sqrt[3]{y}}}\right) + \log \left(\frac{\sqrt[3]{x}}{\sqrt[3]{\sqrt[3]{y}}}\right)\right)\right) - zdouble f(double x, double y, double z) {
double r559434 = x;
double r559435 = y;
double r559436 = r559434 / r559435;
double r559437 = log(r559436);
double r559438 = r559434 * r559437;
double r559439 = z;
double r559440 = r559438 - r559439;
return r559440;
}
double f(double x, double y, double z) {
double r559441 = x;
double r559442 = 1.0;
double r559443 = y;
double r559444 = cbrt(r559443);
double r559445 = r559444 * r559444;
double r559446 = r559442 / r559445;
double r559447 = log(r559446);
double r559448 = 2.0;
double r559449 = cbrt(r559441);
double r559450 = cbrt(r559444);
double r559451 = r559449 / r559450;
double r559452 = log(r559451);
double r559453 = r559448 * r559452;
double r559454 = r559453 + r559452;
double r559455 = r559447 + r559454;
double r559456 = r559441 * r559455;
double r559457 = z;
double r559458 = r559456 - r559457;
return r559458;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 15.3 |
|---|---|
| Target | 8.0 |
| Herbie | 0.3 |
Initial program 15.3
rmApplied add-cube-cbrt15.3
Applied *-un-lft-identity15.3
Applied times-frac15.3
Applied log-prod4.5
rmApplied add-cube-cbrt4.5
Applied add-cube-cbrt4.5
Applied times-frac4.5
Applied log-prod0.3
Simplified0.3
Final simplification0.3
herbie shell --seed 2020089
(FPCore (x y z)
:name "Numeric.SpecFunctions.Extra:bd0 from math-functions-0.1.5.2"
:precision binary64
:herbie-target
(if (< y 7.595077799083773e-308) (- (* x (log (/ x y))) z) (- (* x (- (log x) (log y))) z))
(- (* x (log (/ x y))) z))