Octave 3.8, jcobi/4

Average Error: 53.8 → 36.7

Time: 20.7s

Precision: 64

\[\alpha \gt -1 \land \beta \gt -1 \land i \gt 1\]

Math

\[\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}\;i \le 1.332923833371911853732656743604909165254 \cdot 10^{154}:\\ \;\;\;\;\left(\sqrt{\frac{i}{\left(\sqrt{1} + \mathsf{fma}\left(2, i, \alpha + \beta\right)\right) \cdot \frac{\mathsf{fma}\left(2, i, \alpha + \beta\right)}{\sqrt{\mathsf{fma}\left(\beta, \alpha, i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}}}} \cdot \sqrt{\frac{i}{\left(\sqrt{1} + \mathsf{fma}\left(2, i, \alpha + \beta\right)\right) \cdot \frac{\mathsf{fma}\left(2, i, \alpha + \beta\right)}{\sqrt{\mathsf{fma}\left(\beta, \alpha, i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}}}}\right) \cdot \frac{\frac{\left(\alpha + \beta\right) + i}{\frac{\mathsf{fma}\left(2, i, \alpha + \beta\right)}{\sqrt{\mathsf{fma}\left(\beta, \alpha, i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}}}}{\mathsf{fma}\left(2, i, \alpha + \beta\right) - \sqrt{1}}\\ \mathbf{else}:\\ \;\;\;\;\frac{\frac{0}{\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}\\ \end{array}\]

C

double f(double alpha, double beta, double i) {
        double r65093 = i;
        double r65094 = alpha;
        double r65095 = beta;
        double r65096 = r65094 + r65095;
        double r65097 = r65096 + r65093;
        double r65098 = r65093 * r65097;
        double r65099 = r65095 * r65094;
        double r65100 = r65099 + r65098;
        double r65101 = r65098 * r65100;
        double r65102 = 2.0;
        double r65103 = r65102 * r65093;
        double r65104 = r65096 + r65103;
        double r65105 = r65104 * r65104;
        double r65106 = r65101 / r65105;
        double r65107 = 1.0;
        double r65108 = r65105 - r65107;
        double r65109 = r65106 / r65108;
        return r65109;
}

↓

double f(double alpha, double beta, double i) {
        double r65110 = i;
        double r65111 = 1.3329238333719119e+154;
        bool r65112 = r65110 <= r65111;
        double r65113 = 1.0;
        double r65114 = sqrt(r65113);
        double r65115 = 2.0;
        double r65116 = alpha;
        double r65117 = beta;
        double r65118 = r65116 + r65117;
        double r65119 = fma(r65115, r65110, r65118);
        double r65120 = r65114 + r65119;
        double r65121 = r65118 + r65110;
        double r65122 = r65110 * r65121;
        double r65123 = fma(r65117, r65116, r65122);
        double r65124 = sqrt(r65123);
        double r65125 = r65119 / r65124;
        double r65126 = r65120 * r65125;
        double r65127 = r65110 / r65126;
        double r65128 = sqrt(r65127);
        double r65129 = r65128 * r65128;
        double r65130 = r65121 / r65125;
        double r65131 = r65119 - r65114;
        double r65132 = r65130 / r65131;
        double r65133 = r65129 * r65132;
        double r65134 = 0.0;
        double r65135 = r65115 * r65110;
        double r65136 = r65118 + r65135;
        double r65137 = r65136 * r65136;
        double r65138 = r65134 / r65137;
        double r65139 = r65137 - r65113;
        double r65140 = r65138 / r65139;
        double r65141 = r65112 ? r65133 : r65140;
        return r65141;
}

Error

Bits error versus alpha

Bits error versus beta

Bits error versus i

Derivation

1. Split input into 2 regimes

2. if i < 1.3329238333719119e+154

   1. Initial program 43.7

      \[\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. Using strategy rm

   3. Applied associate-/l* 15.6

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

   4. Simplified 15.6

      \[\leadsto \frac{\frac{i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\color{blue}{\frac{\mathsf{fma}\left(2, i, \alpha + \beta\right) \cdot \mathsf{fma}\left(2, i, \alpha + \beta\right)}{\mathsf{fma}\left(\beta, \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) - 1}\]
   5. Using strategy rm

   6. Applied add-sqr-sqrt 15.6

      \[\leadsto \frac{\frac{i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\frac{\mathsf{fma}\left(2, i, \alpha + \beta\right) \cdot \mathsf{fma}\left(2, i, \alpha + \beta\right)}{\mathsf{fma}\left(\beta, \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) - \color{blue}{\sqrt{1} \cdot \sqrt{1}}}\]

   7. Applied difference-of-squares 15.6

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

   8. Applied add-sqr-sqrt 15.6

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

   9. Applied times-frac 15.6

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

   10. Applied times-frac 15.6

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

   11. Applied times-frac 11.7

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

   12. Simplified 11.7

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

   13. Simplified 11.7

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

   15. Applied add-sqr-sqrt 11.8

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

   if 1.3329238333719119e+154 < i

   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. Taylor expanded around 0 61.9

      \[\leadsto \frac{\frac{\color{blue}{0}}{\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}\]
3. Recombined 2 regimes into one program.

4. Final simplification 36.7

   \[\leadsto \begin{array}{l} \mathbf{if}\;i \le 1.332923833371911853732656743604909165254 \cdot 10^{154}:\\ \;\;\;\;\left(\sqrt{\frac{i}{\left(\sqrt{1} + \mathsf{fma}\left(2, i, \alpha + \beta\right)\right) \cdot \frac{\mathsf{fma}\left(2, i, \alpha + \beta\right)}{\sqrt{\mathsf{fma}\left(\beta, \alpha, i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}}}} \cdot \sqrt{\frac{i}{\left(\sqrt{1} + \mathsf{fma}\left(2, i, \alpha + \beta\right)\right) \cdot \frac{\mathsf{fma}\left(2, i, \alpha + \beta\right)}{\sqrt{\mathsf{fma}\left(\beta, \alpha, i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}}}}\right) \cdot \frac{\frac{\left(\alpha + \beta\right) + i}{\frac{\mathsf{fma}\left(2, i, \alpha + \beta\right)}{\sqrt{\mathsf{fma}\left(\beta, \alpha, i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}}}}{\mathsf{fma}\left(2, i, \alpha + \beta\right) - \sqrt{1}}\\ \mathbf{else}:\\ \;\;\;\;\frac{\frac{0}{\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}\\ \end{array}\]

Reproduce

herbie shell --seed 2019304 +o rules:numerics
(FPCore (alpha beta i)
  :name "Octave 3.8, jcobi/4"
  :precision binary64
  :pre (and (> alpha -1) (> beta -1) (> i 1))
  (/ (/ (* (* 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)))