Average Error: 53.7 → 35.8
Time: 6.3s
Precision: binary64
\[\alpha > -1 \land \beta > -1 \land i > 1\]
\[\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}\]
↓
\[\begin{array}{l} \mathbf{if}\;\alpha \leq 5.8969926040431975 \cdot 10^{+212}:\\ \;\;\;\;\left(\frac{\sqrt{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}}{\alpha + \left(\beta + \left(i \cdot 2 + \sqrt{1}\right)\right)} \cdot \left(\frac{i}{\alpha + \left(\beta + i \cdot 2\right)} \cdot \frac{i + \left(\alpha + \beta\right)}{\alpha + \left(\beta + i \cdot 2\right)}\right)\right) \cdot \frac{\sqrt{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}}{\left(\alpha + \left(\beta + i \cdot 2\right)\right) - \sqrt{1}}\\ \mathbf{else}:\\ \;\;\;\;0\\ \end{array}\]
\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}
↓
\begin{array}{l}
\mathbf{if}\;\alpha \leq 5.8969926040431975 \cdot 10^{+212}:\\
\;\;\;\;\left(\frac{\sqrt{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}}{\alpha + \left(\beta + \left(i \cdot 2 + \sqrt{1}\right)\right)} \cdot \left(\frac{i}{\alpha + \left(\beta + i \cdot 2\right)} \cdot \frac{i + \left(\alpha + \beta\right)}{\alpha + \left(\beta + i \cdot 2\right)}\right)\right) \cdot \frac{\sqrt{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}}{\left(\alpha + \left(\beta + i \cdot 2\right)\right) - \sqrt{1}}\\

\mathbf{else}:\\
\;\;\;\;0\\

\end{array}
(FPCore (alpha beta i)
 :precision binary64
 (/
  (/
   (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i))))
   (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i))))
  (- (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i))) 1.0)))
↓
(FPCore (alpha beta i)
 :precision binary64
 (if (<= alpha 5.8969926040431975e+212)
   (*
    (*
     (/
      (sqrt (+ (* i (+ i (+ alpha beta))) (* alpha beta)))
      (+ alpha (+ beta (+ (* i 2.0) (sqrt 1.0)))))
     (*
      (/ i (+ alpha (+ beta (* i 2.0))))
      (/ (+ i (+ alpha beta)) (+ alpha (+ beta (* i 2.0))))))
    (/
     (sqrt (+ (* i (+ i (+ alpha beta))) (* alpha beta)))
     (- (+ alpha (+ beta (* i 2.0))) (sqrt 1.0))))
   0.0))
double code(double alpha, double beta, double i) {
	return ((((double) (((double) (i * ((double) (((double) (alpha + beta)) + i)))) * ((double) (((double) (beta * alpha)) + ((double) (i * ((double) (((double) (alpha + beta)) + i)))))))) / ((double) (((double) (((double) (alpha + beta)) + ((double) (2.0 * i)))) * ((double) (((double) (alpha + beta)) + ((double) (2.0 * i))))))) / ((double) (((double) (((double) (((double) (alpha + beta)) + ((double) (2.0 * i)))) * ((double) (((double) (alpha + beta)) + ((double) (2.0 * i)))))) - 1.0)));
}
↓
double code(double alpha, double beta, double i) {
	double VAR;
	if ((alpha <= 5.8969926040431975e+212)) {
		VAR = ((double) (((double) ((((double) sqrt(((double) (((double) (i * ((double) (i + ((double) (alpha + beta)))))) + ((double) (alpha * beta)))))) / ((double) (alpha + ((double) (beta + ((double) (((double) (i * 2.0)) + ((double) sqrt(1.0))))))))) * ((double) ((i / ((double) (alpha + ((double) (beta + ((double) (i * 2.0))))))) * (((double) (i + ((double) (alpha + beta)))) / ((double) (alpha + ((double) (beta + ((double) (i * 2.0))))))))))) * (((double) sqrt(((double) (((double) (i * ((double) (i + ((double) (alpha + beta)))))) + ((double) (alpha * beta)))))) / ((double) (((double) (alpha + ((double) (beta + ((double) (i * 2.0)))))) - ((double) sqrt(1.0)))))));
	} else {
		VAR = 0.0;
	}
	return VAR;
}

Error

Bits error versus alpha

Bits error versus beta

Bits error versus i

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Split input into 2 regimes

  2. if alpha < 5.8969926040431975e212

    1. Initial program 52.6

      \[\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}\]

    2. Simplified47.4

      \[\leadsto \color{blue}{\left(\frac{i}{\alpha + \left(\beta + i \cdot 2\right)} \cdot \left(i + \left(\alpha + \beta\right)\right)\right) \cdot \frac{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}{\left(\alpha + \left(\beta + i \cdot 2\right)\right) \cdot \left(\left(\alpha + \left(\beta + i \cdot 2\right)\right) \cdot \left(\alpha + \left(\beta + i \cdot 2\right)\right) - 1\right)}}\]
    3. Using strategy rm

    4. Applied *-un-lft-identity47.4

      \[\leadsto \left(\frac{i}{\alpha + \left(\beta + i \cdot 2\right)} \cdot \left(i + \left(\alpha + \beta\right)\right)\right) \cdot \frac{\color{blue}{1 \cdot \left(i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta\right)}}{\left(\alpha + \left(\beta + i \cdot 2\right)\right) \cdot \left(\left(\alpha + \left(\beta + i \cdot 2\right)\right) \cdot \left(\alpha + \left(\beta + i \cdot 2\right)\right) - 1\right)}\]

    5. Applied times-frac37.4

      \[\leadsto \left(\frac{i}{\alpha + \left(\beta + i \cdot 2\right)} \cdot \left(i + \left(\alpha + \beta\right)\right)\right) \cdot \color{blue}{\left(\frac{1}{\alpha + \left(\beta + i \cdot 2\right)} \cdot \frac{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}{\left(\alpha + \left(\beta + i \cdot 2\right)\right) \cdot \left(\alpha + \left(\beta + i \cdot 2\right)\right) - 1}\right)}\]

    6. Applied associate-*r*37.4

      \[\leadsto \color{blue}{\left(\left(\frac{i}{\alpha + \left(\beta + i \cdot 2\right)} \cdot \left(i + \left(\alpha + \beta\right)\right)\right) \cdot \frac{1}{\alpha + \left(\beta + i \cdot 2\right)}\right) \cdot \frac{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}{\left(\alpha + \left(\beta + i \cdot 2\right)\right) \cdot \left(\alpha + \left(\beta + i \cdot 2\right)\right) - 1}}\]

    7. Simplified37.4

      \[\leadsto \color{blue}{\left(\frac{i}{\alpha + \left(\beta + i \cdot 2\right)} \cdot \frac{i + \left(\alpha + \beta\right)}{\alpha + \left(\beta + i \cdot 2\right)}\right)} \cdot \frac{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}{\left(\alpha + \left(\beta + i \cdot 2\right)\right) \cdot \left(\alpha + \left(\beta + i \cdot 2\right)\right) - 1}\]
    8. Using strategy rm

    9. Applied add-sqr-sqrt37.4

      \[\leadsto \left(\frac{i}{\alpha + \left(\beta + i \cdot 2\right)} \cdot \frac{i + \left(\alpha + \beta\right)}{\alpha + \left(\beta + i \cdot 2\right)}\right) \cdot \frac{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}{\left(\alpha + \left(\beta + i \cdot 2\right)\right) \cdot \left(\alpha + \left(\beta + i \cdot 2\right)\right) - \color{blue}{\sqrt{1} \cdot \sqrt{1}}}\]

    10. Applied difference-of-squares37.4

      \[\leadsto \left(\frac{i}{\alpha + \left(\beta + i \cdot 2\right)} \cdot \frac{i + \left(\alpha + \beta\right)}{\alpha + \left(\beta + i \cdot 2\right)}\right) \cdot \frac{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}{\color{blue}{\left(\left(\alpha + \left(\beta + i \cdot 2\right)\right) + \sqrt{1}\right) \cdot \left(\left(\alpha + \left(\beta + i \cdot 2\right)\right) - \sqrt{1}\right)}}\]

    11. Applied add-sqr-sqrt37.4

      \[\leadsto \left(\frac{i}{\alpha + \left(\beta + i \cdot 2\right)} \cdot \frac{i + \left(\alpha + \beta\right)}{\alpha + \left(\beta + i \cdot 2\right)}\right) \cdot \frac{\color{blue}{\sqrt{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta} \cdot \sqrt{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}}}{\left(\left(\alpha + \left(\beta + i \cdot 2\right)\right) + \sqrt{1}\right) \cdot \left(\left(\alpha + \left(\beta + i \cdot 2\right)\right) - \sqrt{1}\right)}\]

    12. Applied times-frac34.9

      \[\leadsto \left(\frac{i}{\alpha + \left(\beta + i \cdot 2\right)} \cdot \frac{i + \left(\alpha + \beta\right)}{\alpha + \left(\beta + i \cdot 2\right)}\right) \cdot \color{blue}{\left(\frac{\sqrt{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}}{\left(\alpha + \left(\beta + i \cdot 2\right)\right) + \sqrt{1}} \cdot \frac{\sqrt{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}}{\left(\alpha + \left(\beta + i \cdot 2\right)\right) - \sqrt{1}}\right)}\]

    13. Applied associate-*r*34.9

      \[\leadsto \color{blue}{\left(\left(\frac{i}{\alpha + \left(\beta + i \cdot 2\right)} \cdot \frac{i + \left(\alpha + \beta\right)}{\alpha + \left(\beta + i \cdot 2\right)}\right) \cdot \frac{\sqrt{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}}{\left(\alpha + \left(\beta + i \cdot 2\right)\right) + \sqrt{1}}\right) \cdot \frac{\sqrt{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}}{\left(\alpha + \left(\beta + i \cdot 2\right)\right) - \sqrt{1}}}\]

    14. Simplified34.9

      \[\leadsto \color{blue}{\left(\frac{\sqrt{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}}{\alpha + \left(\beta + \left(i \cdot 2 + \sqrt{1}\right)\right)} \cdot \left(\frac{i}{\alpha + \left(\beta + i \cdot 2\right)} \cdot \frac{i + \left(\alpha + \beta\right)}{\alpha + \left(\beta + i \cdot 2\right)}\right)\right)} \cdot \frac{\sqrt{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}}{\left(\alpha + \left(\beta + i \cdot 2\right)\right) - \sqrt{1}}\]

    if 5.8969926040431975e212 < alpha

    1. Initial program 64.0

      \[\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}\]

    2. Simplified56.7

      \[\leadsto \color{blue}{\left(\frac{i}{\alpha + \left(\beta + i \cdot 2\right)} \cdot \left(i + \left(\alpha + \beta\right)\right)\right) \cdot \frac{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}{\left(\alpha + \left(\beta + i \cdot 2\right)\right) \cdot \left(\left(\alpha + \left(\beta + i \cdot 2\right)\right) \cdot \left(\alpha + \left(\beta + i \cdot 2\right)\right) - 1\right)}}\]

    3. Taylor expanded around inf 44.2

      \[\leadsto \color{blue}{0}\]
  3. Recombined 2 regimes into one program.

  4. Final simplification35.8

    \[\leadsto \begin{array}{l} \mathbf{if}\;\alpha \leq 5.8969926040431975 \cdot 10^{+212}:\\ \;\;\;\;\left(\frac{\sqrt{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}}{\alpha + \left(\beta + \left(i \cdot 2 + \sqrt{1}\right)\right)} \cdot \left(\frac{i}{\alpha + \left(\beta + i \cdot 2\right)} \cdot \frac{i + \left(\alpha + \beta\right)}{\alpha + \left(\beta + i \cdot 2\right)}\right)\right) \cdot \frac{\sqrt{i \cdot \left(i + \left(\alpha + \beta\right)\right) + \alpha \cdot \beta}}{\left(\alpha + \left(\beta + i \cdot 2\right)\right) - \sqrt{1}}\\ \mathbf{else}:\\ \;\;\;\;0\\ \end{array}\]

Reproduce

herbie shell --seed 2020198 
(FPCore (alpha beta i)
  :name "Octave 3.8, jcobi/4"
  :precision binary64
  :pre (and (> alpha -1.0) (> beta -1.0) (> i 1.0))
  (/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i)))) (- (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i))) 1.0)))