\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b}}{y}\frac{\sqrt[3]{\left(\sqrt[3]{x} \cdot \sqrt[3]{x}\right) \cdot \sqrt[3]{x}}}{\frac{\sqrt[3]{y}}{\sqrt[3]{e^{\left(y \cdot \log z - 1.0 \cdot \log a\right) + \left(t \cdot \log a - b\right)}}}} \cdot \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.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}}}}double f(double x, double y, double z, double t, double a, double b) {
double r12633676 = x;
double r12633677 = y;
double r12633678 = z;
double r12633679 = log(r12633678);
double r12633680 = r12633677 * r12633679;
double r12633681 = t;
double r12633682 = 1.0;
double r12633683 = r12633681 - r12633682;
double r12633684 = a;
double r12633685 = log(r12633684);
double r12633686 = r12633683 * r12633685;
double r12633687 = r12633680 + r12633686;
double r12633688 = b;
double r12633689 = r12633687 - r12633688;
double r12633690 = exp(r12633689);
double r12633691 = r12633676 * r12633690;
double r12633692 = r12633691 / r12633677;
return r12633692;
}
double f(double x, double y, double z, double t, double a, double b) {
double r12633693 = x;
double r12633694 = cbrt(r12633693);
double r12633695 = r12633694 * r12633694;
double r12633696 = r12633695 * r12633694;
double r12633697 = cbrt(r12633696);
double r12633698 = y;
double r12633699 = cbrt(r12633698);
double r12633700 = z;
double r12633701 = log(r12633700);
double r12633702 = r12633698 * r12633701;
double r12633703 = 1.0;
double r12633704 = a;
double r12633705 = log(r12633704);
double r12633706 = r12633703 * r12633705;
double r12633707 = r12633702 - r12633706;
double r12633708 = t;
double r12633709 = r12633708 * r12633705;
double r12633710 = b;
double r12633711 = r12633709 - r12633710;
double r12633712 = r12633707 + r12633711;
double r12633713 = exp(r12633712);
double r12633714 = cbrt(r12633713);
double r12633715 = r12633699 / r12633714;
double r12633716 = r12633697 / r12633715;
double r12633717 = r12633699 * r12633699;
double r12633718 = r12633708 - r12633703;
double r12633719 = r12633718 * r12633705;
double r12633720 = r12633702 + r12633719;
double r12633721 = r12633720 - r12633710;
double r12633722 = exp(r12633721);
double r12633723 = cbrt(r12633722);
double r12633724 = r12633723 * r12633723;
double r12633725 = r12633717 / r12633724;
double r12633726 = r12633695 / r12633725;
double r12633727 = r12633716 * r12633726;
return r12633727;
}



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.9
rmApplied add-cube-cbrt1.9
Applied add-cube-cbrt1.9
Applied times-frac1.9
Applied add-cube-cbrt1.9
Applied times-frac1.0
Taylor expanded around inf 1.0
Simplified1.0
rmApplied add-cbrt-cube1.0
Final simplification1.0
herbie shell --seed 2019120
(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))