x \cdot \log \left(\frac{x}{y}\right) - z\left(\left(\left(2 \cdot \log \left(\frac{\sqrt[3]{x}}{\sqrt[3]{y}}\right)\right) \cdot 1\right) \cdot x + x \cdot \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 f(double x, double y, double z) {
double r368782 = x;
double r368783 = y;
double r368784 = r368782 / r368783;
double r368785 = log(r368784);
double r368786 = r368782 * r368785;
double r368787 = z;
double r368788 = r368786 - r368787;
return r368788;
}
double f(double x, double y, double z) {
double r368789 = 2.0;
double r368790 = x;
double r368791 = cbrt(r368790);
double r368792 = y;
double r368793 = cbrt(r368792);
double r368794 = r368791 / r368793;
double r368795 = log(r368794);
double r368796 = r368789 * r368795;
double r368797 = 1.0;
double r368798 = r368796 * r368797;
double r368799 = r368798 * r368790;
double r368800 = cbrt(r368791);
double r368801 = r368800 * r368800;
double r368802 = r368801 * r368800;
double r368803 = r368802 / r368793;
double r368804 = log(r368803);
double r368805 = r368790 * r368804;
double r368806 = r368799 + r368805;
double r368807 = z;
double r368808 = r368806 - r368807;
return r368808;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 15.3 |
|---|---|
| Target | 7.9 |
| Herbie | 0.2 |
Initial program 15.3
rmApplied add-cube-cbrt15.3
Applied add-cube-cbrt15.3
Applied times-frac15.3
Applied log-prod3.8
Applied distribute-lft-in3.8
Simplified0.2
rmApplied add-cube-cbrt0.2
Final simplification0.2
herbie shell --seed 2020062 +o rules:numerics
(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))