x \cdot \log \left(\frac{x}{y}\right) - zx \cdot \left(2 \cdot \log \left(\frac{\sqrt[3]{x}}{\sqrt[3]{y}}\right) + \log \left(\frac{\left(\sqrt[3]{\sqrt[3]{x}} \cdot \sqrt[3]{\sqrt[3]{x}}\right) \cdot \sqrt[3]{\sqrt[3]{x}}}{\sqrt[3]{y}}\right)\right) - zdouble code(double x, double y, double z) {
return ((x * log((x / y))) - z);
}
double code(double x, double y, double z) {
return ((x * ((2.0 * log((cbrt(x) / cbrt(y)))) + log((((cbrt(cbrt(x)) * cbrt(cbrt(x))) * cbrt(cbrt(x))) / cbrt(y))))) - z);
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 15.6 |
|---|---|
| Target | 8.0 |
| Herbie | 0.2 |
Initial program 15.6
rmApplied add-cube-cbrt15.6
Applied add-cube-cbrt15.6
Applied times-frac15.6
Applied log-prod3.8
Simplified0.2
rmApplied add-cube-cbrt0.2
Final simplification0.2
herbie shell --seed 2020100
(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))