\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\frac{x \cdot \frac{{\left(\frac{1}{a}\right)}^{1}}{e^{y \cdot \log \left(\frac{1}{z}\right) + \left(\left(\log \left(\frac{1}{a}\right) \cdot \left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)\right) \cdot \sqrt[3]{t} + b\right)}}}{y}double f(double x, double y, double z, double t, double a, double b) {
double r475273 = x;
double r475274 = y;
double r475275 = z;
double r475276 = log(r475275);
double r475277 = r475274 * r475276;
double r475278 = t;
double r475279 = 1.0;
double r475280 = r475278 - r475279;
double r475281 = a;
double r475282 = log(r475281);
double r475283 = r475280 * r475282;
double r475284 = r475277 + r475283;
double r475285 = b;
double r475286 = r475284 - r475285;
double r475287 = exp(r475286);
double r475288 = r475273 * r475287;
double r475289 = r475288 / r475274;
return r475289;
}
double f(double x, double y, double z, double t, double a, double b) {
double r475290 = x;
double r475291 = 1.0;
double r475292 = a;
double r475293 = r475291 / r475292;
double r475294 = 1.0;
double r475295 = pow(r475293, r475294);
double r475296 = y;
double r475297 = z;
double r475298 = r475291 / r475297;
double r475299 = log(r475298);
double r475300 = r475296 * r475299;
double r475301 = log(r475293);
double r475302 = t;
double r475303 = cbrt(r475302);
double r475304 = r475303 * r475303;
double r475305 = r475301 * r475304;
double r475306 = r475305 * r475303;
double r475307 = b;
double r475308 = r475306 + r475307;
double r475309 = r475300 + r475308;
double r475310 = exp(r475309);
double r475311 = r475295 / r475310;
double r475312 = r475290 * r475311;
double r475313 = r475312 / r475296;
return r475313;
}




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.4 |
| Herbie | 1.2 |
Initial program 1.9
Taylor expanded around inf 1.9
Simplified1.2
rmApplied add-cube-cbrt1.2
Applied associate-*r*1.2
Final simplification1.2
herbie shell --seed 2020021
(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))