\left(\left(\log \left(x + y\right) + \log z\right) - t\right) + \left(a - 0.5\right) \cdot \log t
\left(\left(\log \left(\sqrt[3]{x + y} \cdot \sqrt[3]{x + y}\right) + \left(\log \left(\sqrt[3]{x + y}\right) + \log 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 r356108 = x;
double r356109 = y;
double r356110 = r356108 + r356109;
double r356111 = log(r356110);
double r356112 = z;
double r356113 = log(r356112);
double r356114 = r356111 + r356113;
double r356115 = t;
double r356116 = r356114 - r356115;
double r356117 = a;
double r356118 = 0.5;
double r356119 = r356117 - r356118;
double r356120 = log(r356115);
double r356121 = r356119 * r356120;
double r356122 = r356116 + r356121;
return r356122;
}
double f(double x, double y, double z, double t, double a) {
double r356123 = x;
double r356124 = y;
double r356125 = r356123 + r356124;
double r356126 = cbrt(r356125);
double r356127 = r356126 * r356126;
double r356128 = log(r356127);
double r356129 = log(r356126);
double r356130 = z;
double r356131 = log(r356130);
double r356132 = r356129 + r356131;
double r356133 = r356128 + r356132;
double r356134 = t;
double r356135 = r356133 - r356134;
double r356136 = a;
double r356137 = 0.5;
double r356138 = r356136 - r356137;
double r356139 = log(r356134);
double r356140 = r356138 * r356139;
double r356141 = r356135 + r356140;
return r356141;
}




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 associate-+l+0.3
Final simplification0.3
herbie shell --seed 2020018 +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))))