x \cdot \log \left(\frac{x}{y}\right) - z\left(x \cdot \left(2 \cdot \log \left(\frac{\sqrt[3]{x}}{\sqrt[3]{y}}\right)\right) + x \cdot \log \left(\frac{\sqrt[3]{x}}{\sqrt[3]{\sqrt[3]{y}} \cdot \left(\sqrt[3]{\sqrt[3]{y}} \cdot \sqrt[3]{\sqrt[3]{y}}\right)}\right)\right) - z(FPCore (x y z) :precision binary64 (- (* x (log (/ x y))) z))
(FPCore (x y z)
:precision binary64
(-
(+
(* x (* 2.0 (log (/ (cbrt x) (cbrt y)))))
(*
x
(log
(/ (cbrt x) (* (cbrt (cbrt y)) (* (cbrt (cbrt y)) (cbrt (cbrt y))))))))
z))double 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)))) + (x * log(cbrt(x) / (cbrt(cbrt(y)) * (cbrt(cbrt(y)) * cbrt(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-cbrt_binary64_182915.6
Applied add-cube-cbrt_binary64_182915.6
Applied times-frac_binary64_185315.6
Applied log-prod_binary64_17763.5
Applied distribute-rgt-in_binary64_19013.5
Simplified0.2
Simplified0.2
rmApplied add-cube-cbrt_binary64_18290.2
Final simplification0.2
herbie shell --seed 2020231
(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))