x \cdot \log \left(\frac{x}{y}\right) - zx \cdot \left(\log \left(\frac{1}{\sqrt[3]{y} \cdot \sqrt[3]{y}}\right) + \left(\log \left(\frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{\sqrt[3]{\sqrt[3]{y} \cdot \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 r436516 = x;
double r436517 = y;
double r436518 = r436516 / r436517;
double r436519 = log(r436518);
double r436520 = r436516 * r436519;
double r436521 = z;
double r436522 = r436520 - r436521;
return r436522;
}
double f(double x, double y, double z) {
double r436523 = x;
double r436524 = 1.0;
double r436525 = y;
double r436526 = cbrt(r436525);
double r436527 = r436526 * r436526;
double r436528 = r436524 / r436527;
double r436529 = log(r436528);
double r436530 = cbrt(r436523);
double r436531 = r436530 * r436530;
double r436532 = cbrt(r436527);
double r436533 = r436531 / r436532;
double r436534 = log(r436533);
double r436535 = cbrt(r436526);
double r436536 = r436530 / r436535;
double r436537 = log(r436536);
double r436538 = r436534 + r436537;
double r436539 = r436529 + r436538;
double r436540 = r436523 * r436539;
double r436541 = z;
double r436542 = r436540 - r436541;
return r436542;
}




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 cbrt-prod4.5
Applied add-cube-cbrt4.5
Applied times-frac4.5
Applied log-prod0.3
Final simplification0.3
herbie shell --seed 2020089 +o rules:numerics
(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))