\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 code(double x, double y, double z, double t, double a) {
return (((log((x + y)) + log(z)) - t) + ((a - 0.5) * log(t)));
}
double code(double x, double y, double z, double t, double a) {
return (((log((x + y)) + log(z)) - t) + (((a - 0.5) * (2.0 * log(cbrt(t)))) + ((a - 0.5) * log(pow(t, 0.3333333333333333)))));
}




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 2020102 +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))))