\left(\left(\left(x + y\right) + z\right) - z \cdot \log t\right) + \left(a - 0.5\right) \cdot b
\mathsf{fma}\left(b, a - 0.5, \mathsf{fma}\left(z, 1 - \log \left(\sqrt{t}\right), x + y\right) - z \cdot \log \left(\sqrt{t}\right)\right)double code(double x, double y, double z, double t, double a, double b) {
return ((double) (((double) (((double) (((double) (x + y)) + z)) - ((double) (z * ((double) log(t)))))) + ((double) (((double) (a - 0.5)) * b))));
}
double code(double x, double y, double z, double t, double a, double b) {
return ((double) fma(b, ((double) (a - 0.5)), ((double) (((double) fma(z, ((double) (1.0 - ((double) log(((double) sqrt(t)))))), ((double) (x + y)))) - ((double) (z * ((double) log(((double) sqrt(t))))))))));
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a




Bits error versus b
Results
| Original | 0.1 |
|---|---|
| Target | 0.3 |
| Herbie | 0.1 |
Initial program 0.1
Simplified0.1
rmApplied add-sqr-sqrt0.1
Applied log-prod0.1
Applied distribute-lft-in0.1
Applied associate--r+0.1
Simplified0.1
Final simplification0.1
herbie shell --seed 2020120 +o rules:numerics
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:logBeta from math-functions-0.1.5.2, A"
:precision binary64
:herbie-target
(+ (+ (+ x y) (/ (* (- 1 (pow (log t) 2)) z) (+ 1 (log t)))) (* (- a 0.5) b))
(+ (- (+ (+ x y) z) (* z (log t))) (* (- a 0.5) b)))