x \cdot \log \left(\frac{x}{y}\right) - z\begin{array}{l}
\mathbf{if}\;y \le -5.435048613210959909118539310961078344525 \cdot 10^{-309}:\\
\;\;\;\;\left(\log \left(\frac{-1}{y}\right) - \log \left(\frac{-1}{x}\right)\right) \cdot x - z\\
\mathbf{else}:\\
\;\;\;\;\left(\left(2 \cdot \log \left(\frac{\sqrt[3]{x}}{\sqrt[3]{y}}\right)\right) \cdot x + x \cdot \log \left(\frac{\sqrt[3]{{x}^{\frac{2}{3}}} \cdot \sqrt[3]{\sqrt[3]{x}}}{\sqrt[3]{y}}\right)\right) - z\\
\end{array}double f(double x, double y, double z) {
double r306622 = x;
double r306623 = y;
double r306624 = r306622 / r306623;
double r306625 = log(r306624);
double r306626 = r306622 * r306625;
double r306627 = z;
double r306628 = r306626 - r306627;
return r306628;
}
double f(double x, double y, double z) {
double r306629 = y;
double r306630 = -5.43504861321096e-309;
bool r306631 = r306629 <= r306630;
double r306632 = -1.0;
double r306633 = r306632 / r306629;
double r306634 = log(r306633);
double r306635 = x;
double r306636 = r306632 / r306635;
double r306637 = log(r306636);
double r306638 = r306634 - r306637;
double r306639 = r306638 * r306635;
double r306640 = z;
double r306641 = r306639 - r306640;
double r306642 = 2.0;
double r306643 = cbrt(r306635);
double r306644 = cbrt(r306629);
double r306645 = r306643 / r306644;
double r306646 = log(r306645);
double r306647 = r306642 * r306646;
double r306648 = r306647 * r306635;
double r306649 = 0.6666666666666666;
double r306650 = pow(r306635, r306649);
double r306651 = cbrt(r306650);
double r306652 = cbrt(r306643);
double r306653 = r306651 * r306652;
double r306654 = r306653 / r306644;
double r306655 = log(r306654);
double r306656 = r306635 * r306655;
double r306657 = r306648 + r306656;
double r306658 = r306657 - r306640;
double r306659 = r306631 ? r306641 : r306658;
return r306659;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 15.8 |
|---|---|
| Target | 8.2 |
| Herbie | 0.3 |
if y < -5.43504861321096e-309Initial program 16.3
Taylor expanded around -inf 0.3
Simplified0.3
if -5.43504861321096e-309 < y Initial program 15.2
rmApplied add-cube-cbrt15.2
Applied add-cube-cbrt15.2
Applied times-frac15.2
Applied log-prod3.7
Applied distribute-lft-in3.7
Simplified0.2
rmApplied add-cube-cbrt0.2
Applied cbrt-prod0.2
Simplified0.2
Final simplification0.3
herbie shell --seed 2019323
(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))