\left(x \cdot \log y + z \cdot \log \left(1 - y\right)\right) - t
\left(x \cdot \log \left({\left(\sqrt[3]{y}\right)}^{1.6666666666666667} \cdot \sqrt[3]{\sqrt[3]{y}}\right) + \left(x \cdot \log \left({y}^{0.3333333333333333}\right) + z \cdot \left(\left(y \cdot y\right) \cdot \left(y \cdot -0.3333333333333333 + -0.5\right) - y\right)\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 (* (pow (cbrt y) 1.6666666666666667) (cbrt (cbrt y)))))
(+
(* x (log (pow y 0.3333333333333333)))
(* z (- (* (* y y) (+ (* y -0.3333333333333333) -0.5)) y))))
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(pow(cbrt(y), 1.6666666666666667) * cbrt(cbrt(y)))) + ((x * log(pow(y, 0.3333333333333333))) + (z * (((y * y) * ((y * -0.3333333333333333) + -0.5)) - y)))) - t;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 9.4 |
|---|---|
| Target | 0.3 |
| Herbie | 0.3 |
Initial program 9.4
Taylor expanded around 0 0.3
Simplified0.3
rmApplied add-cube-cbrt_binary640.3
Applied log-prod_binary640.3
Applied distribute-lft-in_binary640.3
Applied associate-+l+_binary640.3
rmApplied pow1/3_binary640.3
rmApplied add-cube-cbrt_binary640.3
Applied associate-*r*_binary640.3
Simplified0.3
Final simplification0.3
herbie shell --seed 2020224
(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))