\frac{x \cdot e^{\left(y \cdot \log z + \left(t - 1\right) \cdot \log a\right) - b}}{y}\begin{array}{l}
\mathbf{if}\;a \le 4.0196862376713476 \cdot 10^{-20}:\\
\;\;\;\;\frac{x \cdot \frac{{\left(\frac{1}{a}\right)}^{1}}{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}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{\frac{{\left(\frac{1}{a}\right)}^{1}}{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}\\
\end{array}double f(double x, double y, double z, double t, double a, double b) {
double r444328 = x;
double r444329 = y;
double r444330 = z;
double r444331 = log(r444330);
double r444332 = r444329 * r444331;
double r444333 = t;
double r444334 = 1.0;
double r444335 = r444333 - r444334;
double r444336 = a;
double r444337 = log(r444336);
double r444338 = r444335 * r444337;
double r444339 = r444332 + r444338;
double r444340 = b;
double r444341 = r444339 - r444340;
double r444342 = exp(r444341);
double r444343 = r444328 * r444342;
double r444344 = r444343 / r444329;
return r444344;
}
double f(double x, double y, double z, double t, double a, double b) {
double r444345 = a;
double r444346 = 4.0196862376713476e-20;
bool r444347 = r444345 <= r444346;
double r444348 = x;
double r444349 = 1.0;
double r444350 = r444349 / r444345;
double r444351 = 1.0;
double r444352 = pow(r444350, r444351);
double r444353 = y;
double r444354 = z;
double r444355 = r444349 / r444354;
double r444356 = log(r444355);
double r444357 = log(r444350);
double r444358 = t;
double r444359 = b;
double r444360 = fma(r444357, r444358, r444359);
double r444361 = fma(r444353, r444356, r444360);
double r444362 = exp(r444361);
double r444363 = r444352 / r444362;
double r444364 = r444348 * r444363;
double r444365 = r444364 / r444353;
double r444366 = r444363 / r444353;
double r444367 = r444348 * r444366;
double r444368 = r444347 ? r444365 : r444367;
return r444368;
}




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.0 |
| Herbie | 0.1 |
if a < 4.0196862376713476e-20Initial program 0.7
Taylor expanded around inf 0.7
Simplified0.1
if 4.0196862376713476e-20 < a Initial program 2.9
Taylor expanded around inf 2.9
Simplified2.1
rmApplied *-un-lft-identity2.1
Applied times-frac0.1
Simplified0.1
Final simplification0.1
herbie shell --seed 2020039 +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))