x \cdot \log \left(\frac{x}{y}\right) - z\left(\left(\left(2 \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) \cdot 1\right) \cdot x + x \cdot \log \left(\frac{\sqrt[3]{x}}{\sqrt[3]{y}}\right)\right) - zdouble f(double x, double y, double z) {
double r449883 = x;
double r449884 = y;
double r449885 = r449883 / r449884;
double r449886 = log(r449885);
double r449887 = r449883 * r449886;
double r449888 = z;
double r449889 = r449887 - r449888;
return r449889;
}
double f(double x, double y, double z) {
double r449890 = 2.0;
double r449891 = x;
double r449892 = cbrt(r449891);
double r449893 = cbrt(r449892);
double r449894 = r449893 * r449893;
double r449895 = r449894 * r449893;
double r449896 = y;
double r449897 = cbrt(r449896);
double r449898 = r449895 / r449897;
double r449899 = log(r449898);
double r449900 = r449890 * r449899;
double r449901 = 1.0;
double r449902 = r449900 * r449901;
double r449903 = r449902 * r449891;
double r449904 = r449892 / r449897;
double r449905 = log(r449904);
double r449906 = r449891 * r449905;
double r449907 = r449903 + r449906;
double r449908 = z;
double r449909 = r449907 - r449908;
return r449909;
}




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))