\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{t \cdot 3}\right)\right)}}\begin{array}{l}
\mathbf{if}\;a \le 3.7163425386581422 \cdot 10^{33}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \mathsf{fma}\left(\frac{z}{1}, \frac{\sqrt{t + a}}{t}, -\left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{t \cdot 3}\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(\mathsf{fma}\left(a + \left(\frac{5}{6} - \frac{2}{t \cdot 3}\right), -\left(b - c\right), \frac{z \cdot \sqrt{t + a}}{t}\right) + \left(a + \left(\frac{5}{6} - \frac{2}{t \cdot 3}\right)\right) \cdot \left(\left(-\left(b - c\right)\right) + \left(b - c\right)\right)\right)}}\\
\end{array}double code(double x, double y, double z, double t, double a, double b, double c) {
return (x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0))))))))));
}
double code(double x, double y, double z, double t, double a, double b, double c) {
double VAR;
if ((a <= 3.716342538658142e+33)) {
VAR = (x / (x + (y * exp((2.0 * fma((z / 1.0), (sqrt((t + a)) / t), -((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0))))))))));
} else {
VAR = (x / (x + (y * exp((2.0 * (fma((a + ((5.0 / 6.0) - (2.0 / (t * 3.0)))), -(b - c), ((z * sqrt((t + a))) / t)) + ((a + ((5.0 / 6.0) - (2.0 / (t * 3.0)))) * (-(b - c) + (b - c)))))))));
}
return VAR;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t



Bits error versus a



Bits error versus b



Bits error versus c
Results
if a < 3.716342538658142e+33Initial program 2.6
rmApplied *-un-lft-identity2.6
Applied times-frac1.3
Applied fma-neg0.4
if 3.716342538658142e+33 < a Initial program 5.8
rmApplied log1p-expm1-u15.7
rmApplied add-sqr-sqrt39.5
Applied prod-diff56.4
Simplified49.4
Simplified3.5
Final simplification1.6
herbie shell --seed 2020100 +o rules:numerics
(FPCore (x y z t a b c)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2"
:precision binary64
(/ x (+ x (* y (exp (* 2 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5 6)) (/ 2 (* t 3)))))))))))