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 r27906399 = x;
double r27906400 = y;
double r27906401 = r27906399 / r27906400;
double r27906402 = log(r27906401);
double r27906403 = r27906399 * r27906402;
double r27906404 = z;
double r27906405 = r27906403 - r27906404;
return r27906405;
}
double f(double x, double y, double z) {
double r27906406 = x;
double r27906407 = cbrt(r27906406);
double r27906408 = y;
double r27906409 = cbrt(r27906408);
double r27906410 = r27906407 / r27906409;
double r27906411 = log(r27906410);
double r27906412 = r27906411 + r27906411;
double r27906413 = cbrt(r27906409);
double r27906414 = r27906413 * r27906413;
double r27906415 = r27906413 * r27906414;
double r27906416 = r27906407 / r27906415;
double r27906417 = log(r27906416);
double r27906418 = r27906412 + r27906417;
double r27906419 = r27906418 * r27906406;
double r27906420 = z;
double r27906421 = r27906419 - r27906420;
return r27906421;
}




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