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




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 15.7 |
|---|---|
| Target | 8.0 |
| Herbie | 0.2 |
Initial program Error: 15.7 bits
rmApplied add-cube-cbrtError: 15.7 bits
Applied add-cube-cbrtError: 15.7 bits
Applied times-fracError: 15.7 bits
Applied log-prodError: 3.8 bits
SimplifiedError: 0.2 bits
Final simplificationError: 0.2 bits
herbie shell --seed 2020204
(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))