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

    if alpha < 9.899003901161565e+151

    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.3
    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.8
    3. Using strategy rm
      35.8
    4. Applied add-sqr-sqrt 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)}}}{\color{red}{{\left(\left(\beta + \alpha\right) + 2 \cdot i\right)}^2 - 1.0}} \leadsto \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)}}}{\color{blue}{{\left(\sqrt{{\left(\left(\beta + \alpha\right) + 2 \cdot i\right)}^2 - 1.0}\right)}^2}}\]
      35.8
    5. Applied add-sqr-sqrt to get
      \[\frac{\color{red}{\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(\sqrt{{\left(\left(\beta + \alpha\right) + 2 \cdot i\right)}^2 - 1.0}\right)}^2} \leadsto \frac{\color{blue}{{\left(\sqrt{\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)}}}\right)}^2}}{{\left(\sqrt{{\left(\left(\beta + \alpha\right) + 2 \cdot i\right)}^2 - 1.0}\right)}^2}\]
      35.8
    6. Applied square-undiv to get
      \[\color{red}{\frac{{\left(\sqrt{\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)}}}\right)}^2}{{\left(\sqrt{{\left(\left(\beta + \alpha\right) + 2 \cdot i\right)}^2 - 1.0}\right)}^2}} \leadsto \color{blue}{{\left(\frac{\sqrt{\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)}}}}{\sqrt{{\left(\left(\beta + \alpha\right) + 2 \cdot i\right)}^2 - 1.0}}\right)}^2}\]
      35.8

    if 9.899003901161565e+151 < alpha

    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}}\]
      50.8
    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)))