\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b}}{y}\frac{\sqrt[3]{x}}{\frac{\sqrt[3]{y}}{\sqrt{e^{\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b}}}} \cdot \frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{\frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\sqrt{\sqrt[3]{{e}^{\left(\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b\right)}} \cdot \left(\sqrt[3]{e^{\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b}} \cdot \sqrt[3]{e^{\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b}}\right)}}}double f(double x, double y, double z, double t, double a, double b) {
double r3169066 = x;
double r3169067 = y;
double r3169068 = z;
double r3169069 = log(r3169068);
double r3169070 = r3169067 * r3169069;
double r3169071 = t;
double r3169072 = 1.0;
double r3169073 = r3169071 - r3169072;
double r3169074 = a;
double r3169075 = log(r3169074);
double r3169076 = r3169073 * r3169075;
double r3169077 = r3169070 + r3169076;
double r3169078 = b;
double r3169079 = r3169077 - r3169078;
double r3169080 = exp(r3169079);
double r3169081 = r3169066 * r3169080;
double r3169082 = r3169081 / r3169067;
return r3169082;
}
double f(double x, double y, double z, double t, double a, double b) {
double r3169083 = x;
double r3169084 = cbrt(r3169083);
double r3169085 = y;
double r3169086 = cbrt(r3169085);
double r3169087 = z;
double r3169088 = log(r3169087);
double r3169089 = r3169085 * r3169088;
double r3169090 = t;
double r3169091 = 1.0;
double r3169092 = r3169090 - r3169091;
double r3169093 = a;
double r3169094 = log(r3169093);
double r3169095 = r3169092 * r3169094;
double r3169096 = r3169089 + r3169095;
double r3169097 = b;
double r3169098 = r3169096 - r3169097;
double r3169099 = exp(r3169098);
double r3169100 = sqrt(r3169099);
double r3169101 = r3169086 / r3169100;
double r3169102 = r3169084 / r3169101;
double r3169103 = r3169084 * r3169084;
double r3169104 = r3169086 * r3169086;
double r3169105 = exp(1.0);
double r3169106 = pow(r3169105, r3169098);
double r3169107 = cbrt(r3169106);
double r3169108 = cbrt(r3169099);
double r3169109 = r3169108 * r3169108;
double r3169110 = r3169107 * r3169109;
double r3169111 = sqrt(r3169110);
double r3169112 = r3169104 / r3169111;
double r3169113 = r3169103 / r3169112;
double r3169114 = r3169102 * r3169113;
return r3169114;
}



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
Initial program 1.9
rmApplied associate-/l*1.8
rmApplied add-sqr-sqrt1.8
Applied add-cube-cbrt1.8
Applied times-frac1.8
Applied add-cube-cbrt1.8
Applied times-frac0.9
rmApplied add-cube-cbrt0.9
rmApplied *-un-lft-identity0.9
Applied exp-prod0.9
Simplified0.9
Final simplification0.9
herbie shell --seed 2019134
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2"
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y))