\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}x \cdot \frac{e^{\log z \cdot y - \left(b - \log a \cdot \left(t - 1\right)\right)}}{y}double f(double x, double y, double z, double t, double a, double b) {
double r453028 = x;
double r453029 = y;
double r453030 = z;
double r453031 = log(r453030);
double r453032 = r453029 * r453031;
double r453033 = t;
double r453034 = 1.0;
double r453035 = r453033 - r453034;
double r453036 = a;
double r453037 = log(r453036);
double r453038 = r453035 * r453037;
double r453039 = r453032 + r453038;
double r453040 = b;
double r453041 = r453039 - r453040;
double r453042 = exp(r453041);
double r453043 = r453028 * r453042;
double r453044 = r453043 / r453029;
return r453044;
}
double f(double x, double y, double z, double t, double a, double b) {
double r453045 = x;
double r453046 = z;
double r453047 = log(r453046);
double r453048 = y;
double r453049 = r453047 * r453048;
double r453050 = b;
double r453051 = a;
double r453052 = log(r453051);
double r453053 = t;
double r453054 = 1.0;
double r453055 = r453053 - r453054;
double r453056 = r453052 * r453055;
double r453057 = r453050 - r453056;
double r453058 = r453049 - r453057;
double r453059 = exp(r453058);
double r453060 = r453059 / r453048;
double r453061 = r453045 * r453060;
return r453061;
}




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 | 11.2 |
| Herbie | 2.1 |
Initial program 2.0
rmApplied *-un-lft-identity2.0
Applied times-frac2.1
Simplified2.1
Simplified18.2
rmApplied add-exp-log18.9
Applied pow-exp18.9
Applied div-exp13.4
Applied add-exp-log13.4
Applied pow-exp13.4
Applied div-exp2.1
Final simplification2.1
herbie shell --seed 2019323 +o rules:numerics
(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))