\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\begin{array}{l}
\mathbf{if}\;\left(t - 1\right) \cdot \log a \le -89857766.7366942018 \lor \neg \left(\left(t - 1\right) \cdot \log a \le 239.969958721699498\right):\\
\;\;\;\;\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot \left({a}^{t} \cdot e^{y \cdot \log z - b}\right)}{y \cdot {a}^{1}}\\
\end{array}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) {
double VAR;
if (((((t - 1.0) * log(a)) <= -89857766.7366942) || !(((t - 1.0) * log(a)) <= 239.9699587216995))) {
VAR = ((x * exp((((y * log(z)) + ((t - 1.0) * log(a))) - b))) / y);
} else {
VAR = ((x * (pow(a, t) * exp(((y * log(z)) - b)))) / (y * pow(a, 1.0)));
}
return VAR;
}




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.8 |
|---|---|
| Target | 10.8 |
| Herbie | 0.6 |
if (* (- t 1.0) (log a)) < -89857766.7366942 or 239.9699587216995 < (* (- t 1.0) (log a)) Initial program 0.3
if -89857766.7366942 < (* (- t 1.0) (log a)) < 239.9699587216995Initial program 4.4
rmApplied +-commutative4.4
Applied associate--l+4.4
Applied exp-sum4.4
Simplified3.0
rmApplied pow-sub3.0
Applied associate-*l/3.0
Applied associate-*r/2.9
Applied associate-/l/1.0
Final simplification0.6
herbie shell --seed 2020078
(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)) y)) (- (+ b 1) (* y (log z)))) (if (< t 852031.2288374073) (/ (* (/ x y) (pow a (- t 1))) (exp (- b (* (log z) y)))) (/ (* x (/ (pow a (- t 1)) y)) (- (+ b 1) (* y (log z))))))
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1) (log a))) b))) y))