\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 r511173 = x;
double r511174 = y;
double r511175 = z;
double r511176 = log(r511175);
double r511177 = r511174 * r511176;
double r511178 = t;
double r511179 = 1.0;
double r511180 = r511178 - r511179;
double r511181 = a;
double r511182 = log(r511181);
double r511183 = r511180 * r511182;
double r511184 = r511177 + r511183;
double r511185 = b;
double r511186 = r511184 - r511185;
double r511187 = exp(r511186);
double r511188 = r511173 * r511187;
double r511189 = r511188 / r511174;
return r511189;
}
double f(double x, double y, double z, double t, double a, double b) {
double r511190 = x;
double r511191 = z;
double r511192 = log(r511191);
double r511193 = y;
double r511194 = r511192 * r511193;
double r511195 = b;
double r511196 = a;
double r511197 = log(r511196);
double r511198 = t;
double r511199 = 1.0;
double r511200 = r511198 - r511199;
double r511201 = r511197 * r511200;
double r511202 = r511195 - r511201;
double r511203 = r511194 - r511202;
double r511204 = exp(r511203);
double r511205 = r511204 / r511193;
double r511206 = r511190 * r511205;
return r511206;
}




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))