\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 r462083 = x;
double r462084 = y;
double r462085 = 2.0;
double r462086 = z;
double r462087 = t;
double r462088 = a;
double r462089 = r462087 + r462088;
double r462090 = sqrt(r462089);
double r462091 = r462086 * r462090;
double r462092 = r462091 / r462087;
double r462093 = b;
double r462094 = c;
double r462095 = r462093 - r462094;
double r462096 = 5.0;
double r462097 = 6.0;
double r462098 = r462096 / r462097;
double r462099 = r462088 + r462098;
double r462100 = 3.0;
double r462101 = r462087 * r462100;
double r462102 = r462085 / r462101;
double r462103 = r462099 - r462102;
double r462104 = r462095 * r462103;
double r462105 = r462092 - r462104;
double r462106 = r462085 * r462105;
double r462107 = exp(r462106);
double r462108 = r462084 * r462107;
double r462109 = r462083 + r462108;
double r462110 = r462083 / r462109;
return r462110;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r462111 = x;
double r462112 = y;
double r462113 = 2.0;
double r462114 = z;
double r462115 = t;
double r462116 = cbrt(r462115);
double r462117 = r462116 * r462116;
double r462118 = r462114 / r462117;
double r462119 = a;
double r462120 = r462115 + r462119;
double r462121 = sqrt(r462120);
double r462122 = r462121 / r462116;
double r462123 = r462118 * r462122;
double r462124 = b;
double r462125 = c;
double r462126 = r462124 - r462125;
double r462127 = 5.0;
double r462128 = 6.0;
double r462129 = r462127 / r462128;
double r462130 = r462119 + r462129;
double r462131 = 3.0;
double r462132 = r462115 * r462131;
double r462133 = r462113 / r462132;
double r462134 = r462130 - r462133;
double r462135 = r462126 * r462134;
double r462136 = r462123 - r462135;
double r462137 = r462113 * r462136;
double r462138 = exp(r462137);
double r462139 = r462112 * r462138;
double r462140 = r462111 + r462139;
double r462141 = r462111 / r462140;
return r462141;
}




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.7 |
|---|---|
| Target | 3.0 |
| Herbie | 2.5 |
Initial program 3.7
rmApplied add-cube-cbrt3.7
Applied times-frac2.5
Final simplification2.5
herbie shell --seed 2020081
(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)))))))))))