\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 \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{0.6666666666666666}{t} - \left(a + \frac{5}{6}\right)\right)\right)}}double code(double x, double y, double z, double t, double a, double b, double c) {
return (x / ((double) (x + ((double) (y * ((double) exp(((double) (2.0 * ((double) ((((double) (z * ((double) sqrt(((double) (t + a)))))) / t) - ((double) (((double) (b - c)) * ((double) (((double) (a + (5.0 / 6.0))) - (2.0 / ((double) (t * 3.0))))))))))))))))));
}
double code(double x, double y, double z, double t, double a, double b, double c) {
return (x / ((double) (x + ((double) (y * ((double) exp(((double) (2.0 * ((double) ((((double) (z * ((double) sqrt(((double) (t + a)))))) / t) + ((double) (((double) (b - c)) * ((double) ((0.6666666666666666 / t) - ((double) (a + (5.0 / 6.0))))))))))))))))));
}




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 | 4.0 |
|---|---|
| Target | 3.1 |
| Herbie | 4.0 |
Initial program 4.0
Taylor expanded around 0 4.0
Final simplification4.0
herbie shell --seed 2020196
(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.0 (- (+ (* a c) (* 0.8333333333333334 c)) (* a b))))))) (if (< t 5.196588770651547e-123) (/ x (+ x (* y (exp (* 2.0 (/ (- (* (* z (sqrt (+ t a))) (* (* 3.0 t) (- a (/ 5.0 6.0)))) (* (- (* (+ (/ 5.0 6.0) a) (* 3.0 t)) 2.0) (* (- a (/ 5.0 6.0)) (* (- b c) t)))) (* (* (* t t) 3.0) (- a (/ 5.0 6.0))))))))) (/ x (+ x (* y (exp (* 2.0 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))))))))))
(/ x (+ x (* y (exp (* 2.0 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))))))))))