\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\frac{x \cdot \left({\left(\sqrt{{e}^{\left(\sqrt[3]{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b} \cdot \sqrt[3]{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}\right)}}\right)}^{\left(\sqrt[3]{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}\right)} \cdot {\left(\sqrt{{e}^{\left(\sqrt[3]{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b} \cdot \sqrt[3]{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}\right)}}\right)}^{\left(\sqrt[3]{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}\right)}\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
(sqrt
(pow
E
(*
(cbrt (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))
(cbrt (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)))))
(cbrt (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)))
(pow
(sqrt
(pow
E
(*
(cbrt (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))
(cbrt (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)))))
(cbrt (- (+ (* 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 (((double) (x * ((double) exp(((double) (((double) (((double) (y * ((double) log(z)))) + ((double) (((double) (t - 1.0)) * ((double) log(a)))))) - b)))))) / y);
}
double code(double x, double y, double z, double t, double a, double b) {
return (((double) (x * ((double) (((double) pow(((double) sqrt(((double) pow(((double) M_E), ((double) (((double) cbrt(((double) (((double) (((double) (y * ((double) log(z)))) + ((double) (((double) (t - 1.0)) * ((double) log(a)))))) - b)))) * ((double) cbrt(((double) (((double) (((double) (y * ((double) log(z)))) + ((double) (((double) (t - 1.0)) * ((double) log(a)))))) - b)))))))))), ((double) cbrt(((double) (((double) (((double) (y * ((double) log(z)))) + ((double) (((double) (t - 1.0)) * ((double) log(a)))))) - b)))))) * ((double) pow(((double) sqrt(((double) pow(((double) M_E), ((double) (((double) cbrt(((double) (((double) (((double) (y * ((double) log(z)))) + ((double) (((double) (t - 1.0)) * ((double) log(a)))))) - b)))) * ((double) cbrt(((double) (((double) (((double) (y * ((double) log(z)))) + ((double) (((double) (t - 1.0)) * ((double) log(a)))))) - b)))))))))), ((double) cbrt(((double) (((double) (((double) (y * ((double) log(z)))) + ((double) (((double) (t - 1.0)) * ((double) log(a)))))) - b)))))))))) / 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 | 1.9 |
|---|---|
| Target | 11.1 |
| Herbie | 2.2 |
Initial program 1.9
rmApplied *-un-lft-identity_binary641.9
Applied exp-prod_binary641.9
Simplified1.9
rmApplied add-cube-cbrt_binary642.2
Applied pow-unpow_binary642.2
rmApplied add-sqr-sqrt_binary642.2
Applied unpow-prod-down_binary642.2
Final simplification2.2
herbie shell --seed 2020210
(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))