\[\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2.0} + 1.0}{2.0}\]
Test:
Octave 3.8, jcobi/1
Bits:
128 bits
Bits error versus alpha
Bits error versus beta
Time: 21.7 s
Input Error: 16.4
Output Error: 0.2
Log:
Profile: 🕒
\(\begin{cases} \frac{2.0 + \frac{\frac{8.0}{\alpha}}{\alpha}}{2.0 \cdot \alpha} + \left(\frac{\frac{\beta}{2.0}}{\alpha + \left(2.0 + \beta\right)} - \frac{\frac{4.0}{\alpha \cdot \alpha}}{2.0}\right) & \text{when } \frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2.0} \le -0.9999999999644088 \\ \frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \log \left(e^{\frac{\alpha}{\left(\alpha + \beta\right) + 2.0} - 1.0}\right)}{2.0} & \text{otherwise} \end{cases}\)

    if (/ (- beta alpha) (+ (+ alpha beta) 2.0)) < -0.9999999999644088

    1. Started with
      \[\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2.0} + 1.0}{2.0}\]
      60.3
    2. Using strategy rm
      60.3
    3. Applied div-sub to get
      \[\frac{\color{red}{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2.0}} + 1.0}{2.0} \leadsto \frac{\color{blue}{\left(\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \frac{\alpha}{\left(\alpha + \beta\right) + 2.0}\right)} + 1.0}{2.0}\]
      60.3
    4. Applied associate-+l- to get
      \[\frac{\color{red}{\left(\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \frac{\alpha}{\left(\alpha + \beta\right) + 2.0}\right) + 1.0}}{2.0} \leadsto \frac{\color{blue}{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2.0} - 1.0\right)}}{2.0}\]
      49.1
    5. Applied taylor to get
      \[\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2.0} - 1.0\right)}{2.0} \leadsto \frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \left(4.0 \cdot \frac{1}{{\alpha}^2} - \left(8.0 \cdot \frac{1}{{\alpha}^{3}} + 2.0 \cdot \frac{1}{\alpha}\right)\right)}{2.0}\]
      0.0
    6. Taylor expanded around inf to get
      \[\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \color{red}{\left(4.0 \cdot \frac{1}{{\alpha}^2} - \left(8.0 \cdot \frac{1}{{\alpha}^{3}} + 2.0 \cdot \frac{1}{\alpha}\right)\right)}}{2.0} \leadsto \frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \color{blue}{\left(4.0 \cdot \frac{1}{{\alpha}^2} - \left(8.0 \cdot \frac{1}{{\alpha}^{3}} + 2.0 \cdot \frac{1}{\alpha}\right)\right)}}{2.0}\]
      0.0
    7. Applied simplify to get
      \[\color{red}{\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \left(4.0 \cdot \frac{1}{{\alpha}^2} - \left(8.0 \cdot \frac{1}{{\alpha}^{3}} + 2.0 \cdot \frac{1}{\alpha}\right)\right)}{2.0}} \leadsto \color{blue}{\frac{2.0 + \frac{\frac{8.0}{\alpha}}{\alpha}}{2.0 \cdot \alpha} + \left(\frac{\frac{\beta}{2.0}}{\alpha + \left(2.0 + \beta\right)} - \frac{\frac{4.0}{\alpha \cdot \alpha}}{2.0}\right)}\]
      0.1

    if -0.9999999999644088 < (/ (- beta alpha) (+ (+ alpha beta) 2.0))

    1. Started with
      \[\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2.0} + 1.0}{2.0}\]
      0.3
    2. Using strategy rm
      0.3
    3. Applied div-sub to get
      \[\frac{\color{red}{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2.0}} + 1.0}{2.0} \leadsto \frac{\color{blue}{\left(\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \frac{\alpha}{\left(\alpha + \beta\right) + 2.0}\right)} + 1.0}{2.0}\]
      0.2
    4. Applied associate-+l- to get
      \[\frac{\color{red}{\left(\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \frac{\alpha}{\left(\alpha + \beta\right) + 2.0}\right) + 1.0}}{2.0} \leadsto \frac{\color{blue}{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \left(\frac{\alpha}{\left(\alpha + \beta\right) + 2.0} - 1.0\right)}}{2.0}\]
      0.2
    5. Using strategy rm
      0.2
    6. Applied add-log-exp to get
      \[\frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \color{red}{\left(\frac{\alpha}{\left(\alpha + \beta\right) + 2.0} - 1.0\right)}}{2.0} \leadsto \frac{\frac{\beta}{\left(\alpha + \beta\right) + 2.0} - \color{blue}{\log \left(e^{\frac{\alpha}{\left(\alpha + \beta\right) + 2.0} - 1.0}\right)}}{2.0}\]
      0.3

  1. Removed slow pow expressions

Original test:


(lambda ((alpha default) (beta default))
  #:name "Octave 3.8, jcobi/1"
  (/ (+ (/ (- beta alpha) (+ (+ alpha beta) 2.0)) 1.0) 2.0))