\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\frac{\frac{x}{\sqrt[3]{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[3]{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 \frac{{\left(\frac{\sqrt{1}}{a}\right)}^{1}}{\sqrt[3]{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 r81088 = x;
double r81089 = y;
double r81090 = z;
double r81091 = log(r81090);
double r81092 = r81089 * r81091;
double r81093 = t;
double r81094 = 1.0;
double r81095 = r81093 - r81094;
double r81096 = a;
double r81097 = log(r81096);
double r81098 = r81095 * r81097;
double r81099 = r81092 + r81098;
double r81100 = b;
double r81101 = r81099 - r81100;
double r81102 = exp(r81101);
double r81103 = r81088 * r81102;
double r81104 = r81103 / r81089;
return r81104;
}
double f(double x, double y, double z, double t, double a, double b) {
double r81105 = x;
double r81106 = y;
double r81107 = 1.0;
double r81108 = z;
double r81109 = r81107 / r81108;
double r81110 = log(r81109);
double r81111 = a;
double r81112 = r81107 / r81111;
double r81113 = log(r81112);
double r81114 = t;
double r81115 = b;
double r81116 = fma(r81113, r81114, r81115);
double r81117 = fma(r81106, r81110, r81116);
double r81118 = exp(r81117);
double r81119 = cbrt(r81118);
double r81120 = r81119 * r81119;
double r81121 = r81105 / r81120;
double r81122 = sqrt(r81107);
double r81123 = r81122 / r81111;
double r81124 = 1.0;
double r81125 = pow(r81123, r81124);
double r81126 = r81125 / r81119;
double r81127 = r81121 * r81126;
double r81128 = r81127 / r81106;
return r81128;
}



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.8
Taylor expanded around inf 1.8
Simplified1.1
rmApplied add-cube-cbrt1.1
Applied *-un-lft-identity1.1
Applied add-sqr-sqrt1.1
Applied times-frac1.1
Applied unpow-prod-down1.1
Applied times-frac1.1
Applied associate-*r*1.1
Simplified1.1
Final simplification1.1
herbie shell --seed 2020025 +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))