\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 r11711479 = x;
double r11711480 = y;
double r11711481 = z;
double r11711482 = log(r11711481);
double r11711483 = r11711480 * r11711482;
double r11711484 = t;
double r11711485 = 1.0;
double r11711486 = r11711484 - r11711485;
double r11711487 = a;
double r11711488 = log(r11711487);
double r11711489 = r11711486 * r11711488;
double r11711490 = r11711483 + r11711489;
double r11711491 = b;
double r11711492 = r11711490 - r11711491;
double r11711493 = exp(r11711492);
double r11711494 = r11711479 * r11711493;
double r11711495 = r11711494 / r11711480;
return r11711495;
}
double f(double x, double y, double z, double t, double a, double b) {
double r11711496 = x;
double r11711497 = a;
double r11711498 = log(r11711497);
double r11711499 = t;
double r11711500 = 1.0;
double r11711501 = r11711499 - r11711500;
double r11711502 = r11711498 * r11711501;
double r11711503 = z;
double r11711504 = log(r11711503);
double r11711505 = y;
double r11711506 = r11711504 * r11711505;
double r11711507 = r11711502 + r11711506;
double r11711508 = b;
double r11711509 = r11711507 - r11711508;
double r11711510 = exp(r11711509);
double r11711511 = r11711496 * r11711510;
double r11711512 = r11711511 / r11711505;
double r11711513 = cbrt(r11711512);
double r11711514 = r11711513 * r11711513;
double r11711515 = r11711514 * r11711513;
double r11711516 = cbrt(r11711515);
double r11711517 = r11711516 * r11711513;
double r11711518 = r11711516 * r11711517;
return r11711518;
}



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 +o rules:numerics
(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))