\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 \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)}}double f(double x, double y, double z, double t, double a, double b, double c) {
double r409058 = x;
double r409059 = y;
double r409060 = 2.0;
double r409061 = z;
double r409062 = t;
double r409063 = a;
double r409064 = r409062 + r409063;
double r409065 = sqrt(r409064);
double r409066 = r409061 * r409065;
double r409067 = r409066 / r409062;
double r409068 = b;
double r409069 = c;
double r409070 = r409068 - r409069;
double r409071 = 5.0;
double r409072 = 6.0;
double r409073 = r409071 / r409072;
double r409074 = r409063 + r409073;
double r409075 = 3.0;
double r409076 = r409062 * r409075;
double r409077 = r409060 / r409076;
double r409078 = r409074 - r409077;
double r409079 = r409070 * r409078;
double r409080 = r409067 - r409079;
double r409081 = r409060 * r409080;
double r409082 = exp(r409081);
double r409083 = r409059 * r409082;
double r409084 = r409058 + r409083;
double r409085 = r409058 / r409084;
return r409085;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r409086 = x;
double r409087 = y;
double r409088 = 2.0;
double r409089 = z;
double r409090 = t;
double r409091 = cbrt(r409090);
double r409092 = r409091 * r409091;
double r409093 = r409089 / r409092;
double r409094 = a;
double r409095 = r409090 + r409094;
double r409096 = sqrt(r409095);
double r409097 = r409096 / r409091;
double r409098 = r409093 * r409097;
double r409099 = b;
double r409100 = c;
double r409101 = r409099 - r409100;
double r409102 = 5.0;
double r409103 = 6.0;
double r409104 = r409102 / r409103;
double r409105 = r409094 + r409104;
double r409106 = 3.0;
double r409107 = r409090 * r409106;
double r409108 = r409088 / r409107;
double r409109 = r409105 - r409108;
double r409110 = r409101 * r409109;
double r409111 = r409098 - r409110;
double r409112 = r409088 * r409111;
double r409113 = exp(r409112);
double r409114 = r409087 * r409113;
double r409115 = r409086 + r409114;
double r409116 = r409086 / r409115;
return r409116;
}




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.6 |
|---|---|
| Target | 2.7 |
| Herbie | 2.6 |
Initial program 3.6
rmApplied add-cube-cbrt3.6
Applied times-frac2.6
Final simplification2.6
herbie shell --seed 2020056
(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)))))))))))