\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 r284578 = x;
double r284579 = y;
double r284580 = 2.0;
double r284581 = z;
double r284582 = t;
double r284583 = a;
double r284584 = r284582 + r284583;
double r284585 = sqrt(r284584);
double r284586 = r284581 * r284585;
double r284587 = r284586 / r284582;
double r284588 = b;
double r284589 = c;
double r284590 = r284588 - r284589;
double r284591 = 5.0;
double r284592 = 6.0;
double r284593 = r284591 / r284592;
double r284594 = r284583 + r284593;
double r284595 = 3.0;
double r284596 = r284582 * r284595;
double r284597 = r284580 / r284596;
double r284598 = r284594 - r284597;
double r284599 = r284590 * r284598;
double r284600 = r284587 - r284599;
double r284601 = r284580 * r284600;
double r284602 = exp(r284601);
double r284603 = r284579 * r284602;
double r284604 = r284578 + r284603;
double r284605 = r284578 / r284604;
return r284605;
}
double f(double x, double y, double z, double t, double a, double b, double c) {
double r284606 = x;
double r284607 = y;
double r284608 = 2.0;
double r284609 = z;
double r284610 = t;
double r284611 = cbrt(r284610);
double r284612 = r284611 * r284611;
double r284613 = r284609 / r284612;
double r284614 = a;
double r284615 = r284610 + r284614;
double r284616 = sqrt(r284615);
double r284617 = r284616 / r284611;
double r284618 = r284613 * r284617;
double r284619 = b;
double r284620 = c;
double r284621 = r284619 - r284620;
double r284622 = 5.0;
double r284623 = 6.0;
double r284624 = r284622 / r284623;
double r284625 = r284614 + r284624;
double r284626 = 3.0;
double r284627 = r284610 * r284626;
double r284628 = r284608 / r284627;
double r284629 = r284625 - r284628;
double r284630 = r284621 * r284629;
double r284631 = r284618 - r284630;
double r284632 = r284608 * r284631;
double r284633 = exp(r284632);
double r284634 = r284607 * r284633;
double r284635 = r284606 + r284634;
double r284636 = r284606 / r284635;
return r284636;
}




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.1 |
| Herbie | 2.6 |
Initial program 3.7
rmApplied add-cube-cbrt3.7
Applied times-frac2.6
Final simplification2.6
herbie shell --seed 2019209
(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.1183266448915811e-50) (/ x (+ x (* y (exp (* 2 (- (+ (* a c) (* 0.83333333333333337 c)) (* a b))))))) (if (< t 5.19658877065154709e-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)))))))))))