x \cdot \log \left(\frac{x}{y}\right) - z\left(\left(\log \left(\frac{\sqrt[3]{x}}{\sqrt[3]{y}}\right) + \log \left(\frac{\sqrt[3]{x}}{\sqrt[3]{y}}\right)\right) + \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) \cdot x - zdouble f(double x, double y, double z) {
double r18806827 = x;
double r18806828 = y;
double r18806829 = r18806827 / r18806828;
double r18806830 = log(r18806829);
double r18806831 = r18806827 * r18806830;
double r18806832 = z;
double r18806833 = r18806831 - r18806832;
return r18806833;
}
double f(double x, double y, double z) {
double r18806834 = x;
double r18806835 = cbrt(r18806834);
double r18806836 = y;
double r18806837 = cbrt(r18806836);
double r18806838 = r18806835 / r18806837;
double r18806839 = log(r18806838);
double r18806840 = r18806839 + r18806839;
double r18806841 = cbrt(r18806837);
double r18806842 = r18806841 * r18806841;
double r18806843 = r18806841 * r18806842;
double r18806844 = r18806835 / r18806843;
double r18806845 = log(r18806844);
double r18806846 = r18806840 + r18806845;
double r18806847 = r18806846 * r18806834;
double r18806848 = z;
double r18806849 = r18806847 - r18806848;
return r18806849;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 15.4 |
|---|---|
| Target | 7.9 |
| Herbie | 0.2 |
Initial program 15.4
rmApplied add-cube-cbrt15.4
Applied add-cube-cbrt15.4
Applied times-frac15.4
Applied log-prod3.6
Simplified0.2
rmApplied add-cube-cbrt0.2
Final simplification0.2
herbie shell --seed 2019171 +o rules:numerics
(FPCore (x y z)
:name "Numeric.SpecFunctions.Extra:bd0 from math-functions-0.1.5.2"
:herbie-target
(if (< y 7.595077799083773e-308) (- (* x (log (/ x y))) z) (- (* x (- (log x) (log y))) z))
(- (* x (log (/ x y))) z))