x \cdot \log \left(\frac{x}{y}\right) - z\left(\left(2 \cdot \log \left(\frac{\sqrt[3]{x}}{\sqrt[3]{y}}\right)\right) \cdot x + \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) \cdot x\right) - zdouble f(double x, double y, double z) {
double r538837 = x;
double r538838 = y;
double r538839 = r538837 / r538838;
double r538840 = log(r538839);
double r538841 = r538837 * r538840;
double r538842 = z;
double r538843 = r538841 - r538842;
return r538843;
}
double f(double x, double y, double z) {
double r538844 = 2.0;
double r538845 = x;
double r538846 = cbrt(r538845);
double r538847 = y;
double r538848 = cbrt(r538847);
double r538849 = r538846 / r538848;
double r538850 = log(r538849);
double r538851 = r538844 * r538850;
double r538852 = r538851 * r538845;
double r538853 = cbrt(r538846);
double r538854 = r538853 * r538853;
double r538855 = r538854 * r538853;
double r538856 = r538855 / r538848;
double r538857 = log(r538856);
double r538858 = r538857 * r538845;
double r538859 = r538852 + r538858;
double r538860 = z;
double r538861 = r538859 - r538860;
return r538861;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 15.0 |
|---|---|
| Target | 8.0 |
| Herbie | 0.2 |
Initial program 15.0
rmApplied add-cube-cbrt15.0
Applied add-cube-cbrt15.0
Applied times-frac15.0
Applied log-prod3.5
Applied distribute-lft-in3.5
Simplified0.2
Simplified0.2
rmApplied add-cube-cbrt0.2
Final simplification0.2
herbie shell --seed 2020047
(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))