\left(x \cdot \log y + z \cdot \log \left(1 - y\right)\right) - t
\left(x \cdot \log y + \left(\left(\left(y \cdot y\right) \cdot \left(y \cdot -0.3333333333333333 - 0.5\right)\right) \cdot z - y \cdot z\right)\right) - t
(FPCore (x y z t) :precision binary64 (- (+ (* x (log y)) (* z (log (- 1.0 y)))) t))
(FPCore (x y z t) :precision binary64 (- (+ (* x (log y)) (- (* (* (* y y) (- (* y -0.3333333333333333) 0.5)) z) (* y z))) t))
double code(double x, double y, double z, double t) {
return ((x * log(y)) + (z * log(1.0 - y))) - t;
}
double code(double x, double y, double z, double t) {
return ((x * log(y)) + ((((y * y) * ((y * -0.3333333333333333) - 0.5)) * z) - (y * z))) - t;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 9.9 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 9.9
Taylor expanded around 0 0.2
Simplified0.2
rmApplied sub-neg_binary640.2
Applied distribute-rgt-in_binary640.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2020232
(FPCore (x y z t)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, B"
:precision binary64
:herbie-target
(- (* (- z) (+ (+ (* 0.5 (* y y)) y) (* (/ 0.3333333333333333 (* 1.0 (* 1.0 1.0))) (* y (* y y))))) (- t (* x (log y))))
(- (+ (* x (log y)) (* z (log (- 1.0 y)))) t))