\left(\left(\log \left(x + y\right) + \log z\right) - t\right) + \left(a - 0.5\right) \cdot \log t
\left(\left(\left(\log \left(x + y\right) + \log \left(\sqrt[3]{z} \cdot \sqrt[3]{z}\right)\right) + \log \left(\sqrt[3]{z}\right)\right) - t\right) + \left(a - 0.5\right) \cdot \log tdouble f(double x, double y, double z, double t, double a) {
double r352419 = x;
double r352420 = y;
double r352421 = r352419 + r352420;
double r352422 = log(r352421);
double r352423 = z;
double r352424 = log(r352423);
double r352425 = r352422 + r352424;
double r352426 = t;
double r352427 = r352425 - r352426;
double r352428 = a;
double r352429 = 0.5;
double r352430 = r352428 - r352429;
double r352431 = log(r352426);
double r352432 = r352430 * r352431;
double r352433 = r352427 + r352432;
return r352433;
}
double f(double x, double y, double z, double t, double a) {
double r352434 = x;
double r352435 = y;
double r352436 = r352434 + r352435;
double r352437 = log(r352436);
double r352438 = z;
double r352439 = cbrt(r352438);
double r352440 = r352439 * r352439;
double r352441 = log(r352440);
double r352442 = r352437 + r352441;
double r352443 = log(r352439);
double r352444 = r352442 + r352443;
double r352445 = t;
double r352446 = r352444 - r352445;
double r352447 = a;
double r352448 = 0.5;
double r352449 = r352447 - r352448;
double r352450 = log(r352445);
double r352451 = r352449 * r352450;
double r352452 = r352446 + r352451;
return r352452;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 0.2 |
|---|---|
| Target | 0.3 |
| Herbie | 0.3 |
Initial program 0.2
rmApplied add-cube-cbrt0.2
Applied log-prod0.3
Applied associate-+r+0.3
Final simplification0.3
herbie shell --seed 2020056
(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))))