\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\frac{1}{\frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\frac{{\left(\frac{1}{\sqrt{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)}}}}} \cdot \frac{x}{\frac{\sqrt[3]{y}}{\frac{{\left(\frac{1}{\sqrt{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)}}}}}double f(double x, double y, double z, double t, double a, double b) {
double r434296 = x;
double r434297 = y;
double r434298 = z;
double r434299 = log(r434298);
double r434300 = r434297 * r434299;
double r434301 = t;
double r434302 = 1.0;
double r434303 = r434301 - r434302;
double r434304 = a;
double r434305 = log(r434304);
double r434306 = r434303 * r434305;
double r434307 = r434300 + r434306;
double r434308 = b;
double r434309 = r434307 - r434308;
double r434310 = exp(r434309);
double r434311 = r434296 * r434310;
double r434312 = r434311 / r434297;
return r434312;
}
double f(double x, double y, double z, double t, double a, double b) {
double r434313 = 1.0;
double r434314 = y;
double r434315 = cbrt(r434314);
double r434316 = r434315 * r434315;
double r434317 = a;
double r434318 = sqrt(r434317);
double r434319 = r434313 / r434318;
double r434320 = 1.0;
double r434321 = pow(r434319, r434320);
double r434322 = z;
double r434323 = r434313 / r434322;
double r434324 = log(r434323);
double r434325 = r434313 / r434317;
double r434326 = log(r434325);
double r434327 = t;
double r434328 = b;
double r434329 = fma(r434326, r434327, r434328);
double r434330 = fma(r434314, r434324, r434329);
double r434331 = exp(r434330);
double r434332 = sqrt(r434331);
double r434333 = r434321 / r434332;
double r434334 = r434316 / r434333;
double r434335 = r434313 / r434334;
double r434336 = x;
double r434337 = r434315 / r434333;
double r434338 = r434336 / r434337;
double r434339 = r434335 * r434338;
return r434339;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a




Bits error versus b
| Original | 1.9 |
|---|---|
| Target | 11.7 |
| Herbie | 0.4 |
Initial program 1.9
Taylor expanded around inf 1.9
Simplified1.2
rmApplied associate-/l*1.3
rmApplied add-sqr-sqrt1.3
Applied add-sqr-sqrt1.3
Applied *-un-lft-identity1.3
Applied times-frac1.3
Applied unpow-prod-down1.3
Applied times-frac1.3
Applied add-cube-cbrt1.4
Applied times-frac1.4
Applied *-un-lft-identity1.4
Applied times-frac0.4
Final simplification0.4
herbie shell --seed 2020046 +o rules:numerics
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2, A"
:precision binary64
:herbie-target
(if (< t -0.8845848504127471) (/ (* x (/ (pow a (- t 1)) y)) (- (+ b 1) (* y (log z)))) (if (< t 852031.2288374073) (/ (* (/ x y) (pow a (- t 1))) (exp (- b (* (log z) y)))) (/ (* x (/ (pow a (- t 1)) y)) (- (+ b 1) (* y (log z))))))
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1) (log a))) b))) y))