\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 -7.929621432689752 \cdot 10^{-232} \lor \neg \left(t \le 3.2534287537046085 \cdot 10^{-299}\right):\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z}{\sqrt[3]{t} \cdot \sqrt[3]{t}} \cdot \frac{\sqrt{t + a}}{\sqrt[3]{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 \frac{\left(z \cdot \frac{\sqrt{t + a}}{\sqrt[3]{t}}\right) \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right) - \left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right) \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)}{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right) \cdot \left(\left(a - \frac{5}{6}\right) \cdot \left(t \cdot 3\right)\right)}}}\\
\end{array}double f(double x, double y, double z, double t, double a, double b, double c) {
double r436313 = x;
double r436314 = y;
double r436315 = 2.0;
double r436316 = z;
double r436317 = t;
double r436318 = a;
double r436319 = r436317 + r436318;
double r436320 = sqrt(r436319);
double r436321 = r436316 * r436320;
double r436322 = r436321 / r436317;
double r436323 = b;
double r436324 = c;
double r436325 = r436323 - r436324;
double r436326 = 5.0;
double r436327 = 6.0;
double r436328 = r436326 / r436327;
double r436329 = r436318 + r436328;
double r436330 = 3.0;
double r436331 = r436317 * r436330;
double r436332 = r436315 / r436331;
double r436333 = r436329 - r436332;
double r436334 = r436325 * r436333;
double r436335 = r436322 - r436334;
double r436336 = r436315 * r436335;
double r436337 = exp(r436336);
double r436338 = r436314 * r436337;
double r436339 = r436313 + r436338;
double r436340 = r436313 / r436339;
return r436340;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r436341 = t;
double r436342 = -7.929621432689752e-232;
bool r436343 = r436341 <= r436342;
double r436344 = 3.2534287537046085e-299;
bool r436345 = r436341 <= r436344;
double r436346 = !r436345;
bool r436347 = r436343 || r436346;
double r436348 = x;
double r436349 = y;
double r436350 = 2.0;
double r436351 = z;
double r436352 = cbrt(r436341);
double r436353 = r436352 * r436352;
double r436354 = r436351 / r436353;
double r436355 = a;
double r436356 = r436341 + r436355;
double r436357 = sqrt(r436356);
double r436358 = r436357 / r436352;
double r436359 = r436354 * r436358;
double r436360 = b;
double r436361 = c;
double r436362 = r436360 - r436361;
double r436363 = 5.0;
double r436364 = 6.0;
double r436365 = r436363 / r436364;
double r436366 = r436355 + r436365;
double r436367 = 3.0;
double r436368 = r436341 * r436367;
double r436369 = r436350 / r436368;
double r436370 = r436366 - r436369;
double r436371 = r436362 * r436370;
double r436372 = r436359 - r436371;
double r436373 = r436350 * r436372;
double r436374 = exp(r436373);
double r436375 = r436349 * r436374;
double r436376 = r436348 + r436375;
double r436377 = r436348 / r436376;
double r436378 = r436351 * r436358;
double r436379 = r436355 - r436365;
double r436380 = r436379 * r436368;
double r436381 = r436378 * r436380;
double r436382 = r436355 * r436355;
double r436383 = r436365 * r436365;
double r436384 = r436382 - r436383;
double r436385 = r436384 * r436368;
double r436386 = r436379 * r436350;
double r436387 = r436385 - r436386;
double r436388 = r436362 * r436387;
double r436389 = r436353 * r436388;
double r436390 = r436381 - r436389;
double r436391 = r436353 * r436380;
double r436392 = r436390 / r436391;
double r436393 = r436350 * r436392;
double r436394 = exp(r436393);
double r436395 = r436349 * r436394;
double r436396 = r436348 + r436395;
double r436397 = r436348 / r436396;
double r436398 = r436347 ? r436377 : r436397;
return r436398;
}




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
| Original | 3.9 |
|---|---|
| Target | 3.0 |
| Herbie | 2.6 |
if t < -7.929621432689752e-232 or 3.2534287537046085e-299 < t Initial program 3.3
rmApplied add-cube-cbrt3.3
Applied times-frac2.1
if -7.929621432689752e-232 < t < 3.2534287537046085e-299Initial program 12.0
rmApplied add-cube-cbrt12.0
Applied times-frac12.1
rmApplied flip-+15.1
Applied frac-sub15.1
Applied associate-*r/15.1
Applied associate-*l/15.0
Applied frac-sub8.9
Final simplification2.6
herbie shell --seed 2020035
(FPCore (x y z t a b c)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, I"
:precision binary64
:herbie-target
(if (< t -2.118326644891581e-50) (/ x (+ x (* y (exp (* 2 (- (+ (* a c) (* 0.8333333333333334 c)) (* a b))))))) (if (< t 5.196588770651547e-123) (/ x (+ x (* y (exp (* 2 (/ (- (* (* z (sqrt (+ t a))) (* (* 3 t) (- a (/ 5 6)))) (* (- (* (+ (/ 5 6) a) (* 3 t)) 2) (* (- a (/ 5 6)) (* (- b c) t)))) (* (* (* t t) 3) (- a (/ 5 6))))))))) (/ x (+ x (* y (exp (* 2 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5 6)) (/ 2 (* t 3))))))))))))
(/ x (+ x (* y (exp (* 2 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5 6)) (/ 2 (* t 3)))))))))))