\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{\frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\sqrt[3]{e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}} \cdot \sqrt[3]{e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}}} \cdot \frac{\sqrt[3]{x}}{\frac{\left(\sqrt[3]{\sqrt[3]{y}} \cdot \sqrt[3]{\sqrt[3]{y}}\right) \cdot \sqrt[3]{\sqrt[3]{y}}}{\sqrt[3]{{e}^{\left(\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b\right)}}}}double f(double x, double y, double z, double t, double a, double b) {
double r24444061 = x;
double r24444062 = y;
double r24444063 = z;
double r24444064 = log(r24444063);
double r24444065 = r24444062 * r24444064;
double r24444066 = t;
double r24444067 = 1.0;
double r24444068 = r24444066 - r24444067;
double r24444069 = a;
double r24444070 = log(r24444069);
double r24444071 = r24444068 * r24444070;
double r24444072 = r24444065 + r24444071;
double r24444073 = b;
double r24444074 = r24444072 - r24444073;
double r24444075 = exp(r24444074);
double r24444076 = r24444061 * r24444075;
double r24444077 = r24444076 / r24444062;
return r24444077;
}
double f(double x, double y, double z, double t, double a, double b) {
double r24444078 = x;
double r24444079 = cbrt(r24444078);
double r24444080 = r24444079 * r24444079;
double r24444081 = y;
double r24444082 = cbrt(r24444081);
double r24444083 = r24444082 * r24444082;
double r24444084 = z;
double r24444085 = log(r24444084);
double r24444086 = r24444081 * r24444085;
double r24444087 = t;
double r24444088 = 1.0;
double r24444089 = r24444087 - r24444088;
double r24444090 = a;
double r24444091 = log(r24444090);
double r24444092 = r24444089 * r24444091;
double r24444093 = r24444086 + r24444092;
double r24444094 = b;
double r24444095 = r24444093 - r24444094;
double r24444096 = exp(r24444095);
double r24444097 = cbrt(r24444096);
double r24444098 = r24444097 * r24444097;
double r24444099 = r24444083 / r24444098;
double r24444100 = r24444080 / r24444099;
double r24444101 = cbrt(r24444082);
double r24444102 = r24444101 * r24444101;
double r24444103 = r24444102 * r24444101;
double r24444104 = exp(1.0);
double r24444105 = pow(r24444104, r24444095);
double r24444106 = cbrt(r24444105);
double r24444107 = r24444103 / r24444106;
double r24444108 = r24444079 / r24444107;
double r24444109 = r24444100 * r24444108;
return r24444109;
}




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.5 |
| Herbie | 1.0 |
Initial program 2.0
rmApplied associate-/l*1.9
rmApplied add-cube-cbrt2.0
Applied add-cube-cbrt2.0
Applied times-frac2.0
Applied add-cube-cbrt2.0
Applied times-frac1.0
rmApplied add-cube-cbrt1.0
rmApplied *-un-lft-identity1.0
Applied exp-prod1.0
Simplified1.0
Final simplification1.0
herbie shell --seed 2019179
(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.0) (* 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.0) (* y (log z))))))
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y))