\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1.0\right) \cdot \log a\right) - b}}{y}\sqrt[3]{\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) \cdot \sqrt[3]{\frac{x \cdot e^{\left(\log a \cdot \left(t - 1.0\right) + \log z \cdot y\right) - b}}{y}}} \cdot \left(\sqrt[3]{\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) \cdot \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 r10805569 = x;
double r10805570 = y;
double r10805571 = z;
double r10805572 = log(r10805571);
double r10805573 = r10805570 * r10805572;
double r10805574 = t;
double r10805575 = 1.0;
double r10805576 = r10805574 - r10805575;
double r10805577 = a;
double r10805578 = log(r10805577);
double r10805579 = r10805576 * r10805578;
double r10805580 = r10805573 + r10805579;
double r10805581 = b;
double r10805582 = r10805580 - r10805581;
double r10805583 = exp(r10805582);
double r10805584 = r10805569 * r10805583;
double r10805585 = r10805584 / r10805570;
return r10805585;
}
double f(double x, double y, double z, double t, double a, double b) {
double r10805586 = x;
double r10805587 = a;
double r10805588 = log(r10805587);
double r10805589 = t;
double r10805590 = 1.0;
double r10805591 = r10805589 - r10805590;
double r10805592 = r10805588 * r10805591;
double r10805593 = z;
double r10805594 = log(r10805593);
double r10805595 = y;
double r10805596 = r10805594 * r10805595;
double r10805597 = r10805592 + r10805596;
double r10805598 = b;
double r10805599 = r10805597 - r10805598;
double r10805600 = exp(r10805599);
double r10805601 = r10805586 * r10805600;
double r10805602 = r10805601 / r10805595;
double r10805603 = cbrt(r10805602);
double r10805604 = r10805603 * r10805603;
double r10805605 = r10805604 * r10805603;
double r10805606 = cbrt(r10805605);
double r10805607 = r10805606 * r10805603;
double r10805608 = r10805606 * r10805607;
return r10805608;
}



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 add-cube-cbrt1.9
rmApplied add-cube-cbrt1.9
rmApplied add-cube-cbrt1.9
Final simplification1.9
herbie shell --seed 2019121
(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))