\frac{\frac{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1}\frac{\left(\alpha + 1\right) \cdot \left(\frac{\sqrt{1 + \beta}}{\left(\alpha + \beta\right) + 2} \cdot \frac{\sqrt{1 + \beta}}{\left(\alpha + \beta\right) + 2}\right)}{\alpha + \left(\beta + 3\right)}(FPCore (alpha beta) :precision binary64 (/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2.0 1.0))) (+ (+ alpha beta) (* 2.0 1.0))) (+ (+ (+ alpha beta) (* 2.0 1.0)) 1.0)))
(FPCore (alpha beta)
:precision binary64
(/
(*
(+ alpha 1.0)
(*
(/ (sqrt (+ 1.0 beta)) (+ (+ alpha beta) 2.0))
(/ (sqrt (+ 1.0 beta)) (+ (+ alpha beta) 2.0))))
(+ alpha (+ beta 3.0))))double code(double alpha, double beta) {
return (((((alpha + beta) + (beta * alpha)) + 1.0) / ((alpha + beta) + (2.0 * 1.0))) / ((alpha + beta) + (2.0 * 1.0))) / (((alpha + beta) + (2.0 * 1.0)) + 1.0);
}
double code(double alpha, double beta) {
return ((alpha + 1.0) * ((sqrt(1.0 + beta) / ((alpha + beta) + 2.0)) * (sqrt(1.0 + beta) / ((alpha + beta) + 2.0)))) / (alpha + (beta + 3.0));
}



Bits error versus alpha



Bits error versus beta
Results
Initial program 3.8
Simplified2.2
rmApplied add-sqr-sqrt_binary64_12012.3
Applied times-frac_binary64_12110.2
Final simplification0.2
herbie shell --seed 2020231
(FPCore (alpha beta)
:name "Octave 3.8, jcobi/3"
:precision binary64
:pre (and (> alpha -1.0) (> beta -1.0))
(/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2.0 1.0))) (+ (+ alpha beta) (* 2.0 1.0))) (+ (+ (+ alpha beta) (* 2.0 1.0)) 1.0)))