\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}\;t \le -4.7682077658519153 \cdot 10^{-188}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \mathsf{fma}\left(\frac{z}{1}, \frac{\sqrt{t + a}}{t}, -\mathsf{fma}\left(0.66666666666666663, \frac{c}{t}, 0.83333333333333337 \cdot b - 0.66666666666666663 \cdot \frac{b}{t}\right)\right)}}\\
\mathbf{elif}\;t \le 2.56349525650802909 \cdot 10^{-235}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \frac{\left(z \cdot \sqrt{t + a}\right) \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right) - t \cdot \left(\left(b - c\right) \cdot \left(\left(a \cdot a - \frac{5}{6} \cdot \frac{5}{6}\right) \cdot \left(t \cdot 3\right) - \left(a - \frac{5}{6}\right) \cdot 2\right)\right)}{t \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right)}}}\\
\mathbf{else}:\\
\;\;\;\;\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)}}\\
\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 ((t <= -4.7682077658519153e-188)) {
VAR = (x / (x + (y * exp((2.0 * fma((z / 1.0), (sqrt((t + a)) / t), -fma(0.6666666666666666, (c / t), ((0.8333333333333334 * b) - (0.6666666666666666 * (b / t))))))))));
} else {
double VAR_1;
if ((t <= 2.563495256508029e-235)) {
VAR_1 = (x / (x + (y * exp((2.0 * ((((z * sqrt((t + a))) * ((a - (5.0 / 6.0)) * (t * 3.0))) - (t * ((b - c) * ((((a * a) - ((5.0 / 6.0) * (5.0 / 6.0))) * (t * 3.0)) - ((a - (5.0 / 6.0)) * 2.0))))) / (t * ((a - (5.0 / 6.0)) * (t * 3.0)))))))));
} else {
VAR_1 = (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))))))))));
}
VAR = VAR_1;
}
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 t < -4.7682077658519153e-188Initial program 3.9
rmApplied *-un-lft-identity3.9
Applied times-frac3.0
Applied fma-neg1.6
Taylor expanded around 0 9.4
Simplified9.4
if -4.7682077658519153e-188 < t < 2.563495256508029e-235Initial program 9.2
rmApplied flip-+12.0
Applied frac-sub12.0
Applied associate-*r/12.0
Applied frac-sub8.5
if 2.563495256508029e-235 < t Initial program 2.6
rmApplied *-un-lft-identity2.6
Applied times-frac1.2
Applied fma-neg0.6
Final simplification3.2
herbie shell --seed 2020106 +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)))))))))))