\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)}}\frac{x}{x + y \cdot e^{2 \cdot \mathsf{fma}\left(\frac{z}{1}, \frac{\sqrt{t + a}}{t}, \left(-\mathsf{fma}\left(1, a + \frac{5}{6}, -\frac{2}{t \cdot 3}\right)\right) \cdot \left(b - c\right)\right)}}double f(double x, double y, double z, double t, double a, double b, double c) {
double r633971 = x;
double r633972 = y;
double r633973 = 2.0;
double r633974 = z;
double r633975 = t;
double r633976 = a;
double r633977 = r633975 + r633976;
double r633978 = sqrt(r633977);
double r633979 = r633974 * r633978;
double r633980 = r633979 / r633975;
double r633981 = b;
double r633982 = c;
double r633983 = r633981 - r633982;
double r633984 = 5.0;
double r633985 = 6.0;
double r633986 = r633984 / r633985;
double r633987 = r633976 + r633986;
double r633988 = 3.0;
double r633989 = r633975 * r633988;
double r633990 = r633973 / r633989;
double r633991 = r633987 - r633990;
double r633992 = r633983 * r633991;
double r633993 = r633980 - r633992;
double r633994 = r633973 * r633993;
double r633995 = exp(r633994);
double r633996 = r633972 * r633995;
double r633997 = r633971 + r633996;
double r633998 = r633971 / r633997;
return r633998;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r633999 = x;
double r634000 = y;
double r634001 = 2.0;
double r634002 = z;
double r634003 = 1.0;
double r634004 = r634002 / r634003;
double r634005 = t;
double r634006 = a;
double r634007 = r634005 + r634006;
double r634008 = sqrt(r634007);
double r634009 = r634008 / r634005;
double r634010 = 5.0;
double r634011 = 6.0;
double r634012 = r634010 / r634011;
double r634013 = r634006 + r634012;
double r634014 = 3.0;
double r634015 = r634005 * r634014;
double r634016 = r634001 / r634015;
double r634017 = -r634016;
double r634018 = fma(r634003, r634013, r634017);
double r634019 = -r634018;
double r634020 = b;
double r634021 = c;
double r634022 = r634020 - r634021;
double r634023 = r634019 * r634022;
double r634024 = fma(r634004, r634009, r634023);
double r634025 = r634001 * r634024;
double r634026 = exp(r634025);
double r634027 = r634000 * r634026;
double r634028 = r633999 + r634027;
double r634029 = r633999 / r634028;
return r634029;
}




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
| Original | 4.1 |
|---|---|
| Target | 3.1 |
| Herbie | 2.3 |
Initial program 4.1
rmApplied add-cbrt-cube4.1
Applied add-cbrt-cube7.1
Applied cbrt-unprod7.1
Applied add-cbrt-cube7.1
Applied cbrt-undiv7.3
Simplified7.3
rmApplied *-un-lft-identity7.3
Applied times-frac6.6
Applied fma-neg5.4
Simplified2.3
Final simplification2.3
herbie shell --seed 2019362 +o rules:numerics
(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)))))))))))