\left(\left(\log \left(x + y\right) + \log z\right) - t\right) + \left(a - 0.5\right) \cdot \log t
\left(\left(\log \left(x + y\right) + \log z\right) - t\right) + \left(\left(a - 0.5\right) \cdot \left(\log \left(\sqrt[3]{\sqrt[3]{t}} \cdot \sqrt[3]{\sqrt[3]{t}}\right) \cdot 2\right) + \mathsf{fma}\left(\log \left(\sqrt[3]{\sqrt[3]{t}}\right), 2, \log \left({t}^{\frac{1}{3}}\right)\right) \cdot \left(a - 0.5\right)\right)double f(double x, double y, double z, double t, double a) {
double r56700 = x;
double r56701 = y;
double r56702 = r56700 + r56701;
double r56703 = log(r56702);
double r56704 = z;
double r56705 = log(r56704);
double r56706 = r56703 + r56705;
double r56707 = t;
double r56708 = r56706 - r56707;
double r56709 = a;
double r56710 = 0.5;
double r56711 = r56709 - r56710;
double r56712 = log(r56707);
double r56713 = r56711 * r56712;
double r56714 = r56708 + r56713;
return r56714;
}
double f(double x, double y, double z, double t, double a) {
double r56715 = x;
double r56716 = y;
double r56717 = r56715 + r56716;
double r56718 = log(r56717);
double r56719 = z;
double r56720 = log(r56719);
double r56721 = r56718 + r56720;
double r56722 = t;
double r56723 = r56721 - r56722;
double r56724 = a;
double r56725 = 0.5;
double r56726 = r56724 - r56725;
double r56727 = cbrt(r56722);
double r56728 = cbrt(r56727);
double r56729 = r56728 * r56728;
double r56730 = log(r56729);
double r56731 = 2.0;
double r56732 = r56730 * r56731;
double r56733 = r56726 * r56732;
double r56734 = log(r56728);
double r56735 = 0.3333333333333333;
double r56736 = pow(r56722, r56735);
double r56737 = log(r56736);
double r56738 = fma(r56734, r56731, r56737);
double r56739 = r56738 * r56726;
double r56740 = r56733 + r56739;
double r56741 = r56723 + r56740;
return r56741;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t



Bits error versus a
Initial program 0.3
rmApplied add-cube-cbrt0.3
Applied log-prod0.3
Applied distribute-lft-in0.3
Simplified0.3
rmApplied add-cube-cbrt0.3
Applied log-prod0.3
Applied distribute-rgt-in0.3
Applied distribute-lft-in0.3
Applied associate-+l+0.3
Simplified0.3
Final simplification0.3
herbie shell --seed 2020001 +o rules:numerics
(FPCore (x y z t a)
:name "Numeric.SpecFunctions:logGammaL from math-functions-0.1.5.2"
:precision binary64
(+ (- (+ (log (+ x y)) (log z)) t) (* (- a 0.5) (log t))))