\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\frac{\left(x \cdot {\left(\sqrt{e}\right)}^{\left(\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b\right)}\right) \cdot {\left(\sqrt{e}\right)}^{\left(\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b\right)}}{y}double f(double x, double y, double z, double t, double a, double b) {
double r271183 = x;
double r271184 = y;
double r271185 = z;
double r271186 = log(r271185);
double r271187 = r271184 * r271186;
double r271188 = t;
double r271189 = 1.0;
double r271190 = r271188 - r271189;
double r271191 = a;
double r271192 = log(r271191);
double r271193 = r271190 * r271192;
double r271194 = r271187 + r271193;
double r271195 = b;
double r271196 = r271194 - r271195;
double r271197 = exp(r271196);
double r271198 = r271183 * r271197;
double r271199 = r271198 / r271184;
return r271199;
}
double f(double x, double y, double z, double t, double a, double b) {
double r271200 = x;
double r271201 = exp(1.0);
double r271202 = sqrt(r271201);
double r271203 = y;
double r271204 = z;
double r271205 = log(r271204);
double r271206 = r271203 * r271205;
double r271207 = t;
double r271208 = 1.0;
double r271209 = r271207 - r271208;
double r271210 = a;
double r271211 = log(r271210);
double r271212 = r271209 * r271211;
double r271213 = r271206 + r271212;
double r271214 = b;
double r271215 = r271213 - r271214;
double r271216 = pow(r271202, r271215);
double r271217 = r271200 * r271216;
double r271218 = r271217 * r271216;
double r271219 = r271218 / r271203;
return r271219;
}




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 | 10.7 |
| Herbie | 2.0 |
Initial program 1.9
rmApplied *-un-lft-identity1.9
Applied exp-prod2.0
Simplified2.0
rmApplied add-sqr-sqrt2.2
Applied unpow-prod-down2.0
Applied associate-*r*2.0
Final simplification2.0
herbie shell --seed 2019322
(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))