\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\frac{\left(x \cdot {\left({\left(e^{\frac{1}{\sqrt{2}}}\right)}^{\left(\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b\right)}\right)}^{\left(\frac{1}{\sqrt{2}}\right)}\right) \cdot {e}^{\left(\frac{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}{2}\right)}}{y}(FPCore (x y z t a b) :precision binary64 (/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y))
(FPCore (x y z t a b)
:precision binary64
(/
(*
(*
x
(pow
(pow
(exp (/ 1.0 (sqrt 2.0)))
(- (+ (* y (log z)) (* (- t 1.0) (log a))) b))
(/ 1.0 (sqrt 2.0))))
(pow E (/ (- (+ (* y (log z)) (* (- t 1.0) (log a))) b) 2.0)))
y))double code(double x, double y, double z, double t, double a, double b) {
return (x * exp(((y * log(z)) + ((t - 1.0) * log(a))) - b)) / y;
}
double code(double x, double y, double z, double t, double a, double b) {
return ((x * pow(pow(exp(1.0 / sqrt(2.0)), (((y * log(z)) + ((t - 1.0) * log(a))) - b)), (1.0 / sqrt(2.0)))) * pow(((double) M_E), ((((y * log(z)) + ((t - 1.0) * log(a))) - b) / 2.0))) / y;
}




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 | 2.0 |
|---|---|
| Target | 10.9 |
| Herbie | 2.3 |
Initial program 2.0
rmApplied *-un-lft-identity_binary642.0
Applied exp-prod_binary642.0
Simplified2.0
rmApplied sqr-pow_binary642.0
Applied associate-*r*_binary642.0
rmApplied add-sqr-sqrt_binary642.2
Applied *-un-lft-identity_binary642.2
Applied times-frac_binary642.2
Applied pow-unpow_binary642.2
rmApplied div-inv_binary642.2
Applied pow-unpow_binary642.3
Simplified2.3
Final simplification2.3
herbie shell --seed 2020232
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2, A"
:precision binary64
:herbie-target
(if (< t -0.8845848504127471) (/ (* x (/ (pow a (- t 1.0)) y)) (- (+ b 1.0) (* y (log z)))) (if (< t 852031.2288374073) (/ (* (/ x y) (pow a (- t 1.0))) (exp (- b (* (log z) y)))) (/ (* x (/ (pow a (- t 1.0)) y)) (- (+ b 1.0) (* y (log z))))))
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y))