\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(2 \cdot \log \left(\sqrt[3]{t}\right)\right) + \left(a - 0.5\right) \cdot \log \left({t}^{\frac{1}{3}}\right)\right)double f(double x, double y, double z, double t, double a) {
double r375595 = x;
double r375596 = y;
double r375597 = r375595 + r375596;
double r375598 = log(r375597);
double r375599 = z;
double r375600 = log(r375599);
double r375601 = r375598 + r375600;
double r375602 = t;
double r375603 = r375601 - r375602;
double r375604 = a;
double r375605 = 0.5;
double r375606 = r375604 - r375605;
double r375607 = log(r375602);
double r375608 = r375606 * r375607;
double r375609 = r375603 + r375608;
return r375609;
}
double f(double x, double y, double z, double t, double a) {
double r375610 = x;
double r375611 = y;
double r375612 = r375610 + r375611;
double r375613 = log(r375612);
double r375614 = z;
double r375615 = log(r375614);
double r375616 = r375613 + r375615;
double r375617 = t;
double r375618 = r375616 - r375617;
double r375619 = a;
double r375620 = 0.5;
double r375621 = r375619 - r375620;
double r375622 = 2.0;
double r375623 = cbrt(r375617);
double r375624 = log(r375623);
double r375625 = r375622 * r375624;
double r375626 = r375621 * r375625;
double r375627 = 0.3333333333333333;
double r375628 = pow(r375617, r375627);
double r375629 = log(r375628);
double r375630 = r375621 * r375629;
double r375631 = r375626 + r375630;
double r375632 = r375618 + r375631;
return r375632;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 0.3 |
|---|---|
| Target | 0.3 |
| Herbie | 0.3 |
Initial program 0.3
rmApplied add-cube-cbrt0.3
Applied log-prod0.3
Applied distribute-lft-in0.3
Simplified0.3
rmApplied pow1/30.3
Final simplification0.3
herbie shell --seed 2020064 +o rules:numerics
(FPCore (x y z t a)
:name "Numeric.SpecFunctions:logGammaL from math-functions-0.1.5.2"
:precision binary64
:herbie-target
(+ (log (+ x y)) (+ (- (log z) t) (* (- a 0.5) (log t))))
(+ (- (+ (log (+ x y)) (log z)) t) (* (- a 0.5) (log t))))