\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.25161200547673655 \cdot 10^{-9}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(\left(a \cdot c + 0.83333333333333337 \cdot c\right) - a \cdot b\right)}}\\
\mathbf{elif}\;t \le 2.6093485179308713 \cdot 10^{-9}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \frac{\left(\frac{z}{\sqrt[3]{t} \cdot \sqrt[3]{t}} \cdot \sqrt{t + a}\right) \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right) - \sqrt[3]{t} \cdot \left(\left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right) - \left(a - \frac{5}{6}\right) \cdot 2\right)\right)}{\sqrt[3]{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 \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \log \left(e^{\frac{2}{t \cdot 3}}\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 temp;
if ((t <= -4.2516120054767366e-09)) {
temp = (x / (x + (y * exp((2.0 * (((a * c) + (0.8333333333333334 * c)) - (a * b)))))));
} else {
double temp_1;
if ((t <= 2.6093485179308713e-09)) {
temp_1 = (x / (x + (y * exp((2.0 * (((((z / (cbrt(t) * cbrt(t))) * sqrt((t + a))) * ((a - (5.0 / 6.0)) * (t * 3.0))) - (cbrt(t) * ((b - c) * (((a + (5.0 / 6.0)) * ((a - (5.0 / 6.0)) * (t * 3.0))) - ((a - (5.0 / 6.0)) * 2.0))))) / (cbrt(t) * ((a - (5.0 / 6.0)) * (t * 3.0)))))))));
} else {
temp_1 = (x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) - ((b - c) * ((a + (5.0 / 6.0)) - log(exp((2.0 / (t * 3.0))))))))))));
}
temp = temp_1;
}
return temp;
}



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.2516120054767366e-09Initial program 3.7
Taylor expanded around inf 5.7
if -4.2516120054767366e-09 < t < 2.6093485179308713e-09Initial program 4.8
rmApplied add-cube-cbrt4.8
Applied times-frac4.8
rmApplied flip-+8.1
Applied frac-sub8.1
Applied associate-*r/8.1
Applied associate-*r/8.1
Applied frac-sub7.4
rmApplied difference-of-squares7.4
Applied associate-*l*4.9
if 2.6093485179308713e-09 < t Initial program 2.6
rmApplied add-log-exp2.9
Final simplification4.0
herbie shell --seed 2020053
(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)))))))))))