\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b}}{y}\sqrt[3]{\frac{x \cdot {e}^{\left(\left(\log a \cdot \left(t - 1.0\right) + \log z \cdot y\right) - b\right)}}{y}} \cdot \left(\sqrt[3]{\frac{x \cdot e^{\left(\log a \cdot \left(t - 1.0\right) + \log z \cdot y\right) - b}}{y}} \cdot \sqrt[3]{\frac{x \cdot e^{\left(\log a \cdot \left(t - 1.0\right) + \log z \cdot y\right) - b}}{y}}\right)double f(double x, double y, double z, double t, double a, double b) {
double r18918975 = x;
double r18918976 = y;
double r18918977 = z;
double r18918978 = log(r18918977);
double r18918979 = r18918976 * r18918978;
double r18918980 = t;
double r18918981 = 1.0;
double r18918982 = r18918980 - r18918981;
double r18918983 = a;
double r18918984 = log(r18918983);
double r18918985 = r18918982 * r18918984;
double r18918986 = r18918979 + r18918985;
double r18918987 = b;
double r18918988 = r18918986 - r18918987;
double r18918989 = exp(r18918988);
double r18918990 = r18918975 * r18918989;
double r18918991 = r18918990 / r18918976;
return r18918991;
}
double f(double x, double y, double z, double t, double a, double b) {
double r18918992 = x;
double r18918993 = exp(1.0);
double r18918994 = a;
double r18918995 = log(r18918994);
double r18918996 = t;
double r18918997 = 1.0;
double r18918998 = r18918996 - r18918997;
double r18918999 = r18918995 * r18918998;
double r18919000 = z;
double r18919001 = log(r18919000);
double r18919002 = y;
double r18919003 = r18919001 * r18919002;
double r18919004 = r18918999 + r18919003;
double r18919005 = b;
double r18919006 = r18919004 - r18919005;
double r18919007 = pow(r18918993, r18919006);
double r18919008 = r18918992 * r18919007;
double r18919009 = r18919008 / r18919002;
double r18919010 = cbrt(r18919009);
double r18919011 = exp(r18919006);
double r18919012 = r18918992 * r18919011;
double r18919013 = r18919012 / r18919002;
double r18919014 = cbrt(r18919013);
double r18919015 = r18919014 * r18919014;
double r18919016 = r18919010 * r18919015;
return r18919016;
}




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.9 |
| Herbie | 1.8 |
Initial program 1.8
rmApplied add-cube-cbrt1.8
rmApplied *-un-lft-identity1.8
Applied exp-prod1.8
Simplified1.8
Final simplification1.8
herbie shell --seed 2019163 +o rules:numerics
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2, A"
:herbie-target
(if (< t -0.8845848504127471) (/ (* x (/ (pow a (- t 1.0)) y)) (- (+ b 1) (* y (log z)))) (if (< t 852031.2288374073) (/ (* (/ x y) (pow a (- t 1.0))) (exp (- b (* (log z) y)))) (/ (* x (/ (pow a (- t 1.0)) y)) (- (+ b 1) (* y (log z))))))
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y))