x \cdot \log \left(\frac{x}{y}\right) - z\left(\left(\left(2 \cdot \log \left(\frac{\sqrt[3]{x}}{\left(\sqrt[3]{\sqrt[3]{y}} \cdot \sqrt[3]{\sqrt[3]{y}}\right) \cdot \sqrt[3]{\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 r455287 = x;
double r455288 = y;
double r455289 = r455287 / r455288;
double r455290 = log(r455289);
double r455291 = r455287 * r455290;
double r455292 = z;
double r455293 = r455291 - r455292;
return r455293;
}
double f(double x, double y, double z) {
double r455294 = 2.0;
double r455295 = x;
double r455296 = cbrt(r455295);
double r455297 = y;
double r455298 = cbrt(r455297);
double r455299 = cbrt(r455298);
double r455300 = r455299 * r455299;
double r455301 = r455300 * r455299;
double r455302 = r455296 / r455301;
double r455303 = log(r455302);
double r455304 = r455294 * r455303;
double r455305 = 1.0;
double r455306 = r455304 * r455305;
double r455307 = r455306 * r455295;
double r455308 = r455296 / r455298;
double r455309 = log(r455308);
double r455310 = r455295 * r455309;
double r455311 = r455307 + r455310;
double r455312 = z;
double r455313 = r455311 - r455312;
return r455313;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 14.7 |
|---|---|
| Target | 7.7 |
| Herbie | 0.2 |
Initial program 14.7
rmApplied add-cube-cbrt14.7
Applied add-cube-cbrt14.7
Applied times-frac14.7
Applied log-prod3.2
Applied distribute-lft-in3.2
Simplified0.2
rmApplied add-cube-cbrt0.2
Final simplification0.2
herbie shell --seed 2020020 +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))