\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\frac{\left(x \cdot \frac{{\left(\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{1}\right)}^{1}}{\sqrt{\sqrt{e^{\mathsf{fma}\left(y, \log \left(\frac{1}{z}\right), \mathsf{fma}\left(\log \left(\frac{1}{a}\right), t, b\right)\right)}}} \cdot \sqrt{\sqrt{e^{\mathsf{fma}\left(y, \log \left(\frac{1}{z}\right), \mathsf{fma}\left(\log \left(\frac{1}{a}\right), t, b\right)\right)}}}}\right) \cdot \frac{{\left(\frac{\sqrt[3]{1}}{a}\right)}^{1}}{\sqrt{e^{\mathsf{fma}\left(y, \log \left(\frac{1}{z}\right), \mathsf{fma}\left(\log \left(\frac{1}{a}\right), t, b\right)\right)}}}}{y}double f(double x, double y, double z, double t, double a, double b) {
double r99657 = x;
double r99658 = y;
double r99659 = z;
double r99660 = log(r99659);
double r99661 = r99658 * r99660;
double r99662 = t;
double r99663 = 1.0;
double r99664 = r99662 - r99663;
double r99665 = a;
double r99666 = log(r99665);
double r99667 = r99664 * r99666;
double r99668 = r99661 + r99667;
double r99669 = b;
double r99670 = r99668 - r99669;
double r99671 = exp(r99670);
double r99672 = r99657 * r99671;
double r99673 = r99672 / r99658;
return r99673;
}
double f(double x, double y, double z, double t, double a, double b) {
double r99674 = x;
double r99675 = 1.0;
double r99676 = cbrt(r99675);
double r99677 = r99676 * r99676;
double r99678 = r99677 / r99675;
double r99679 = 1.0;
double r99680 = pow(r99678, r99679);
double r99681 = y;
double r99682 = z;
double r99683 = r99675 / r99682;
double r99684 = log(r99683);
double r99685 = a;
double r99686 = r99675 / r99685;
double r99687 = log(r99686);
double r99688 = t;
double r99689 = b;
double r99690 = fma(r99687, r99688, r99689);
double r99691 = fma(r99681, r99684, r99690);
double r99692 = exp(r99691);
double r99693 = sqrt(r99692);
double r99694 = sqrt(r99693);
double r99695 = r99694 * r99694;
double r99696 = r99680 / r99695;
double r99697 = r99674 * r99696;
double r99698 = r99676 / r99685;
double r99699 = pow(r99698, r99679);
double r99700 = r99699 / r99693;
double r99701 = r99697 * r99700;
double r99702 = r99701 / r99681;
return r99702;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t



Bits error versus a



Bits error versus b
Initial program 1.9
Taylor expanded around inf 2.0
Simplified1.2
rmApplied add-sqr-sqrt1.2
Applied *-un-lft-identity1.2
Applied add-cube-cbrt1.2
Applied times-frac1.2
Applied unpow-prod-down1.2
Applied times-frac1.2
Applied associate-*r*1.2
rmApplied add-sqr-sqrt1.2
Applied sqrt-prod1.2
Final simplification1.2
herbie shell --seed 2020062 +o rules:numerics
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2"
:precision binary64
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1) (log a))) b))) y))