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 r26434384 = x;
double r26434385 = y;
double r26434386 = r26434384 / r26434385;
double r26434387 = log(r26434386);
double r26434388 = r26434384 * r26434387;
double r26434389 = z;
double r26434390 = r26434388 - r26434389;
return r26434390;
}
double f(double x, double y, double z) {
double r26434391 = x;
double r26434392 = cbrt(r26434391);
double r26434393 = y;
double r26434394 = cbrt(r26434393);
double r26434395 = r26434392 / r26434394;
double r26434396 = log(r26434395);
double r26434397 = r26434396 + r26434396;
double r26434398 = cbrt(r26434394);
double r26434399 = r26434398 * r26434398;
double r26434400 = r26434398 * r26434399;
double r26434401 = r26434392 / r26434400;
double r26434402 = log(r26434401);
double r26434403 = r26434397 + r26434402;
double r26434404 = r26434403 * r26434391;
double r26434405 = z;
double r26434406 = r26434404 - r26434405;
return r26434406;
}




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
Simplified3.6
rmApplied log-prod0.2
rmApplied add-cube-cbrt0.2
Final simplification0.2
herbie shell --seed 2019171
(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))