\[\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.2 m
Input Error: 52.8
Output Error: 29.2
Log:
Profile: 🕒
\(\begin{cases} \frac{\frac{i \cdot \left(\beta + \left(i + \alpha\right)\right)}{{\left(\left(\beta + \alpha\right) + 2 \cdot i\right)}^2 \cdot \frac{1}{\alpha \cdot \beta + i \cdot \left(\beta + \left(i + \alpha\right)\right)}}}{{\left(\left(\beta + \alpha\right) + 2 \cdot i\right)}^2 - 1.0} & \text{when } \beta \le 4.731717689114216 \cdot 10^{+137} \\ 0 & \text{otherwise} \end{cases}\)

    if beta < 4.731717689114216e+137

    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.6
    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{\frac{i \cdot \left(\beta + \left(i + \alpha\right)\right)}{\frac{{\left(\left(\beta + \alpha\right) + 2 \cdot i\right)}^2}{\alpha \cdot \beta + i \cdot \left(\beta + \left(i + \alpha\right)\right)}}}{{\left(\left(\beta + \alpha\right) + 2 \cdot i\right)}^2 - 1.0}}\]
      35.4
    3. Using strategy rm
      35.4
    4. Applied div-inv to get
      \[\frac{\frac{i \cdot \left(\beta + \left(i + \alpha\right)\right)}{\color{red}{\frac{{\left(\left(\beta + \alpha\right) + 2 \cdot i\right)}^2}{\alpha \cdot \beta + i \cdot \left(\beta + \left(i + \alpha\right)\right)}}}}{{\left(\left(\beta + \alpha\right) + 2 \cdot i\right)}^2 - 1.0} \leadsto \frac{\frac{i \cdot \left(\beta + \left(i + \alpha\right)\right)}{\color{blue}{{\left(\left(\beta + \alpha\right) + 2 \cdot i\right)}^2 \cdot \frac{1}{\alpha \cdot \beta + i \cdot \left(\beta + \left(i + \alpha\right)\right)}}}}{{\left(\left(\beta + \alpha\right) + 2 \cdot i\right)}^2 - 1.0}\]
      35.4

    if 4.731717689114216e+137 < 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}\]
      63.0
    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{\frac{i \cdot \left(\beta + \left(i + \alpha\right)\right)}{\frac{{\left(\left(\beta + \alpha\right) + 2 \cdot i\right)}^2}{\alpha \cdot \beta + i \cdot \left(\beta + \left(i + \alpha\right)\right)}}}{{\left(\left(\beta + \alpha\right) + 2 \cdot i\right)}^2 - 1.0}}\]
      52.2
    3. Applied taylor to get
      \[\frac{\frac{i \cdot \left(\beta + \left(i + \alpha\right)\right)}{\frac{{\left(\left(\beta + \alpha\right) + 2 \cdot i\right)}^2}{\alpha \cdot \beta + i \cdot \left(\beta + \left(i + \alpha\right)\right)}}}{{\left(\left(\beta + \alpha\right) + 2 \cdot i\right)}^2 - 1.0} \leadsto 0\]
      0
    4. Taylor expanded around inf to get
      \[\color{red}{0} \leadsto \color{blue}{0}\]
      0

  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)))