\[\frac{\frac{\left(i \cdot \left(\left(\alpha + \beta\right) + i\right)\right) \cdot \left(\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1.0}\]
Test:
Octave 3.8, jcobi/4
Bits:
128 bits
Bits error versus alpha
Bits error versus beta
Bits error versus i
Time: 1.3 m
Input Error: 52.4
Output Error: 38.3
Log:
Profile: 🕒
\(\begin{cases} \frac{1}{\frac{{\left(\beta + (i * 2 + \alpha)_*\right)}^2 - 1.0}{\left(\frac{\left(\beta + \alpha\right) + i}{\beta + (i * 2 + \alpha)_*} \cdot (i * \left(\left(\beta + \alpha\right) + i\right) + \left(\alpha \cdot \beta\right))_*\right) \cdot \frac{i}{\beta + (i * 2 + \alpha)_*}}} & \text{when } \beta \le 5.551138825541697 \cdot 10^{+153} \\ {\left(\frac{\sqrt{\left(\frac{\left(\beta + \alpha\right) + i}{\beta + (i * 2 + \alpha)_*} \cdot (i * \left(\left(\beta + \alpha\right) + i\right) + \left(\alpha \cdot \beta\right))_*\right) \cdot \frac{i}{\beta + (i * 2 + \alpha)_*}}}{\left(\beta + (\left(\frac{1}{i}\right) * 2 + \left(\frac{1}{\alpha}\right))_*\right) - 0.5 \cdot \frac{1}{\beta}}\right)}^2 & \text{otherwise} \end{cases}\)

    if beta < 5.551138825541697e+153

    1. Started with
      \[\frac{\frac{\left(i \cdot \left(\left(\alpha + \beta\right) + i\right)\right) \cdot \left(\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1.0}\]
      50.4
    2. Applied simplify to get
      \[\color{red}{\frac{\frac{\left(i \cdot \left(\left(\alpha + \beta\right) + i\right)\right) \cdot \left(\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1.0}} \leadsto \color{blue}{\frac{\left(\frac{\left(\beta + \alpha\right) + i}{\beta + (i * 2 + \alpha)_*} \cdot (i * \left(\left(\beta + \alpha\right) + i\right) + \left(\alpha \cdot \beta\right))_*\right) \cdot \frac{i}{\beta + (i * 2 + \alpha)_*}}{{\left(\beta + (i * 2 + \alpha)_*\right)}^2 - 1.0}}\]
      35.9
    3. Using strategy rm
      35.9
    4. Applied clear-num to get
      \[\color{red}{\frac{\left(\frac{\left(\beta + \alpha\right) + i}{\beta + (i * 2 + \alpha)_*} \cdot (i * \left(\left(\beta + \alpha\right) + i\right) + \left(\alpha \cdot \beta\right))_*\right) \cdot \frac{i}{\beta + (i * 2 + \alpha)_*}}{{\left(\beta + (i * 2 + \alpha)_*\right)}^2 - 1.0}} \leadsto \color{blue}{\frac{1}{\frac{{\left(\beta + (i * 2 + \alpha)_*\right)}^2 - 1.0}{\left(\frac{\left(\beta + \alpha\right) + i}{\beta + (i * 2 + \alpha)_*} \cdot (i * \left(\left(\beta + \alpha\right) + i\right) + \left(\alpha \cdot \beta\right))_*\right) \cdot \frac{i}{\beta + (i * 2 + \alpha)_*}}}}\]
      35.9

    if 5.551138825541697e+153 < beta

    1. Started with
      \[\frac{\frac{\left(i \cdot \left(\left(\alpha + \beta\right) + i\right)\right) \cdot \left(\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1.0}\]
      62.4
    2. Applied simplify to get
      \[\color{red}{\frac{\frac{\left(i \cdot \left(\left(\alpha + \beta\right) + i\right)\right) \cdot \left(\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1.0}} \leadsto \color{blue}{\frac{\left(\frac{\left(\beta + \alpha\right) + i}{\beta + (i * 2 + \alpha)_*} \cdot (i * \left(\left(\beta + \alpha\right) + i\right) + \left(\alpha \cdot \beta\right))_*\right) \cdot \frac{i}{\beta + (i * 2 + \alpha)_*}}{{\left(\beta + (i * 2 + \alpha)_*\right)}^2 - 1.0}}\]
      55.2
    3. Using strategy rm
      55.2
    4. Applied add-sqr-sqrt to get
      \[\frac{\left(\frac{\left(\beta + \alpha\right) + i}{\beta + (i * 2 + \alpha)_*} \cdot (i * \left(\left(\beta + \alpha\right) + i\right) + \left(\alpha \cdot \beta\right))_*\right) \cdot \frac{i}{\beta + (i * 2 + \alpha)_*}}{\color{red}{{\left(\beta + (i * 2 + \alpha)_*\right)}^2 - 1.0}} \leadsto \frac{\left(\frac{\left(\beta + \alpha\right) + i}{\beta + (i * 2 + \alpha)_*} \cdot (i * \left(\left(\beta + \alpha\right) + i\right) + \left(\alpha \cdot \beta\right))_*\right) \cdot \frac{i}{\beta + (i * 2 + \alpha)_*}}{\color{blue}{{\left(\sqrt{{\left(\beta + (i * 2 + \alpha)_*\right)}^2 - 1.0}\right)}^2}}\]
      55.2
    5. Applied add-sqr-sqrt to get
      \[\frac{\color{red}{\left(\frac{\left(\beta + \alpha\right) + i}{\beta + (i * 2 + \alpha)_*} \cdot (i * \left(\left(\beta + \alpha\right) + i\right) + \left(\alpha \cdot \beta\right))_*\right) \cdot \frac{i}{\beta + (i * 2 + \alpha)_*}}}{{\left(\sqrt{{\left(\beta + (i * 2 + \alpha)_*\right)}^2 - 1.0}\right)}^2} \leadsto \frac{\color{blue}{{\left(\sqrt{\left(\frac{\left(\beta + \alpha\right) + i}{\beta + (i * 2 + \alpha)_*} \cdot (i * \left(\left(\beta + \alpha\right) + i\right) + \left(\alpha \cdot \beta\right))_*\right) \cdot \frac{i}{\beta + (i * 2 + \alpha)_*}}\right)}^2}}{{\left(\sqrt{{\left(\beta + (i * 2 + \alpha)_*\right)}^2 - 1.0}\right)}^2}\]
      55.2
    6. Applied square-undiv to get
      \[\color{red}{\frac{{\left(\sqrt{\left(\frac{\left(\beta + \alpha\right) + i}{\beta + (i * 2 + \alpha)_*} \cdot (i * \left(\left(\beta + \alpha\right) + i\right) + \left(\alpha \cdot \beta\right))_*\right) \cdot \frac{i}{\beta + (i * 2 + \alpha)_*}}\right)}^2}{{\left(\sqrt{{\left(\beta + (i * 2 + \alpha)_*\right)}^2 - 1.0}\right)}^2}} \leadsto \color{blue}{{\left(\frac{\sqrt{\left(\frac{\left(\beta + \alpha\right) + i}{\beta + (i * 2 + \alpha)_*} \cdot (i * \left(\left(\beta + \alpha\right) + i\right) + \left(\alpha \cdot \beta\right))_*\right) \cdot \frac{i}{\beta + (i * 2 + \alpha)_*}}}{\sqrt{{\left(\beta + (i * 2 + \alpha)_*\right)}^2 - 1.0}}\right)}^2}\]
      55.2
    7. Applied taylor to get
      \[{\left(\frac{\sqrt{\left(\frac{\left(\beta + \alpha\right) + i}{\beta + (i * 2 + \alpha)_*} \cdot (i * \left(\left(\beta + \alpha\right) + i\right) + \left(\alpha \cdot \beta\right))_*\right) \cdot \frac{i}{\beta + (i * 2 + \alpha)_*}}}{\sqrt{{\left(\beta + (i * 2 + \alpha)_*\right)}^2 - 1.0}}\right)}^2 \leadsto {\left(\frac{\sqrt{\left(\frac{\left(\beta + \alpha\right) + i}{\beta + (i * 2 + \alpha)_*} \cdot (i * \left(\left(\beta + \alpha\right) + i\right) + \left(\alpha \cdot \beta\right))_*\right) \cdot \frac{i}{\beta + (i * 2 + \alpha)_*}}}{\left(\beta + (\left(\frac{1}{i}\right) * 2 + \left(\frac{1}{\alpha}\right))_*\right) - 0.5 \cdot \frac{1}{\beta}}\right)}^2\]
      50.2
    8. Taylor expanded around inf to get
      \[{\left(\frac{\sqrt{\left(\frac{\left(\beta + \alpha\right) + i}{\beta + (i * 2 + \alpha)_*} \cdot (i * \left(\left(\beta + \alpha\right) + i\right) + \left(\alpha \cdot \beta\right))_*\right) \cdot \frac{i}{\beta + (i * 2 + \alpha)_*}}}{\color{red}{\left(\beta + (\left(\frac{1}{i}\right) * 2 + \left(\frac{1}{\alpha}\right))_*\right) - 0.5 \cdot \frac{1}{\beta}}}\right)}^2 \leadsto {\left(\frac{\sqrt{\left(\frac{\left(\beta + \alpha\right) + i}{\beta + (i * 2 + \alpha)_*} \cdot (i * \left(\left(\beta + \alpha\right) + i\right) + \left(\alpha \cdot \beta\right))_*\right) \cdot \frac{i}{\beta + (i * 2 + \alpha)_*}}}{\color{blue}{\left(\beta + (\left(\frac{1}{i}\right) * 2 + \left(\frac{1}{\alpha}\right))_*\right) - 0.5 \cdot \frac{1}{\beta}}}\right)}^2\]
      50.2

  1. Removed slow pow expressions

Original test:


(lambda ((alpha default) (beta default) (i default))
  #:name "Octave 3.8, jcobi/4"
  (/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)))