\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b}}{y}\frac{\left(\left(\sqrt[3]{e^{\left(\log a \cdot \left(t - 1.0\right) + \log z \cdot y\right) - b}} \cdot \sqrt[3]{\left(\sqrt[3]{e^{\left(\log a \cdot \left(t - 1.0\right) + \log z \cdot y\right) - b}} \cdot \sqrt[3]{e^{\left(\log a \cdot \left(t - 1.0\right) + \log z \cdot y\right) - b}}\right) \cdot \sqrt[3]{{e}^{\left(\left(\log a \cdot \left(t - 1.0\right) + \log z \cdot y\right) - b\right)}}}\right) \cdot \sqrt[3]{e^{\left(\log a \cdot \left(t - 1.0\right) + \log z \cdot y\right) - b}}\right) \cdot x}{y}double f(double x, double y, double z, double t, double a, double b) {
double r1206956 = x;
double r1206957 = y;
double r1206958 = z;
double r1206959 = log(r1206958);
double r1206960 = r1206957 * r1206959;
double r1206961 = t;
double r1206962 = 1.0;
double r1206963 = r1206961 - r1206962;
double r1206964 = a;
double r1206965 = log(r1206964);
double r1206966 = r1206963 * r1206965;
double r1206967 = r1206960 + r1206966;
double r1206968 = b;
double r1206969 = r1206967 - r1206968;
double r1206970 = exp(r1206969);
double r1206971 = r1206956 * r1206970;
double r1206972 = r1206971 / r1206957;
return r1206972;
}
double f(double x, double y, double z, double t, double a, double b) {
double r1206973 = a;
double r1206974 = log(r1206973);
double r1206975 = t;
double r1206976 = 1.0;
double r1206977 = r1206975 - r1206976;
double r1206978 = r1206974 * r1206977;
double r1206979 = z;
double r1206980 = log(r1206979);
double r1206981 = y;
double r1206982 = r1206980 * r1206981;
double r1206983 = r1206978 + r1206982;
double r1206984 = b;
double r1206985 = r1206983 - r1206984;
double r1206986 = exp(r1206985);
double r1206987 = cbrt(r1206986);
double r1206988 = r1206987 * r1206987;
double r1206989 = exp(1.0);
double r1206990 = pow(r1206989, r1206985);
double r1206991 = cbrt(r1206990);
double r1206992 = r1206988 * r1206991;
double r1206993 = cbrt(r1206992);
double r1206994 = r1206987 * r1206993;
double r1206995 = r1206994 * r1206987;
double r1206996 = x;
double r1206997 = r1206995 * r1206996;
double r1206998 = r1206997 / r1206981;
return r1206998;
}



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.8
rmApplied add-cube-cbrt1.8
rmApplied add-cube-cbrt1.8
rmApplied *-un-lft-identity1.8
Applied exp-prod1.8
Simplified1.8
Final simplification1.8
herbie shell --seed 2019154
(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))