Octave 3.8, jcobi/2

Average Error: 23.9 → 12.2

Time: 28.5s

Precision: 64

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

Math

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

↓

\[\begin{array}{l} \mathbf{if}\;\alpha \le 3.068220085942097902687466146641291923906 \cdot 10^{110}:\\ \;\;\;\;\frac{\left(\alpha + \beta\right) \cdot \frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} + 1}{2}\\ \mathbf{elif}\;\alpha \le 3.104548159948919354427366616967038689222 \cdot 10^{193} \lor \neg \left(\alpha \le 2.571518739036188604595402311839971481899 \cdot 10^{244}\right):\\ \;\;\;\;\frac{\frac{8}{{\alpha}^{3}} + \left(\frac{2}{\alpha} - \frac{4}{\alpha \cdot \alpha}\right)}{2}\\ \mathbf{else}:\\ \;\;\;\;\frac{\sqrt[3]{{\left(\frac{\alpha + \beta}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} \cdot \left(\sqrt[3]{\left(\alpha + \beta\right) + 2 \cdot i} \cdot \sqrt[3]{\left(\alpha + \beta\right) + 2 \cdot i}\right)} \cdot \frac{\frac{\beta - \alpha}{\sqrt[3]{\left(\alpha + \beta\right) + 2 \cdot i}}}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}} + 1\right)}^{3}}}{2}\\ \end{array}\]

C

double f(double alpha, double beta, double i) {
        double r96801 = alpha;
        double r96802 = beta;
        double r96803 = r96801 + r96802;
        double r96804 = r96802 - r96801;
        double r96805 = r96803 * r96804;
        double r96806 = 2.0;
        double r96807 = i;
        double r96808 = r96806 * r96807;
        double r96809 = r96803 + r96808;
        double r96810 = r96805 / r96809;
        double r96811 = r96809 + r96806;
        double r96812 = r96810 / r96811;
        double r96813 = 1.0;
        double r96814 = r96812 + r96813;
        double r96815 = r96814 / r96806;
        return r96815;
}

↓

double f(double alpha, double beta, double i) {
        double r96816 = alpha;
        double r96817 = 3.068220085942098e+110;
        bool r96818 = r96816 <= r96817;
        double r96819 = beta;
        double r96820 = r96816 + r96819;
        double r96821 = r96819 - r96816;
        double r96822 = 2.0;
        double r96823 = i;
        double r96824 = r96822 * r96823;
        double r96825 = r96820 + r96824;
        double r96826 = r96821 / r96825;
        double r96827 = r96825 + r96822;
        double r96828 = r96826 / r96827;
        double r96829 = r96820 * r96828;
        double r96830 = 1.0;
        double r96831 = r96829 + r96830;
        double r96832 = r96831 / r96822;
        double r96833 = 3.1045481599489194e+193;
        bool r96834 = r96816 <= r96833;
        double r96835 = 2.5715187390361886e+244;
        bool r96836 = r96816 <= r96835;
        double r96837 = !r96836;
        bool r96838 = r96834 || r96837;
        double r96839 = 8.0;
        double r96840 = 3.0;
        double r96841 = pow(r96816, r96840);
        double r96842 = r96839 / r96841;
        double r96843 = r96822 / r96816;
        double r96844 = 4.0;
        double r96845 = r96816 * r96816;
        double r96846 = r96844 / r96845;
        double r96847 = r96843 - r96846;
        double r96848 = r96842 + r96847;
        double r96849 = r96848 / r96822;
        double r96850 = sqrt(r96827);
        double r96851 = cbrt(r96825);
        double r96852 = r96851 * r96851;
        double r96853 = r96850 * r96852;
        double r96854 = r96820 / r96853;
        double r96855 = r96821 / r96851;
        double r96856 = r96855 / r96850;
        double r96857 = r96854 * r96856;
        double r96858 = r96857 + r96830;
        double r96859 = pow(r96858, r96840);
        double r96860 = cbrt(r96859);
        double r96861 = r96860 / r96822;
        double r96862 = r96838 ? r96849 : r96861;
        double r96863 = r96818 ? r96832 : r96862;
        return r96863;
}

Error

Bits error versus alpha

Bits error versus beta

Bits error versus i

Derivation

1. Split input into 3 regimes

2. if alpha < 3.068220085942098e+110

   1. Initial program 14.2

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

   3. Applied *-un-lft-identity 14.2

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

   4. Applied *-un-lft-identity 14.2

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

   5. Applied times-frac 3.5

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

   6. Applied times-frac 3.5

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

   7. Simplified 3.5

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

   9. Applied add-cbrt-cube 3.5

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

   10. Simplified 3.5

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

   12. Applied rem-cbrt-cube 3.5

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

   if 3.068220085942098e+110 < alpha < 3.1045481599489194e+193 or 2.5715187390361886e+244 < alpha

   1. Initial program 56.6

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

   3. Applied *-un-lft-identity 56.6

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

   4. Applied *-un-lft-identity 56.6

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

   5. Applied times-frac 42.4

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

   6. Applied times-frac 42.4

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

   7. Simplified 42.4

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

   9. Applied add-cbrt-cube 42.4

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

   10. Simplified 42.4

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

   11. Taylor expanded around inf 41.5

       \[\leadsto \frac{\color{blue}{\left(2 \cdot \frac{1}{\alpha} + 8 \cdot \frac{1}{{\alpha}^{3}}\right) - 4 \cdot \frac{1}{{\alpha}^{2}}}}{2}\]

   12. Simplified 41.5

       \[\leadsto \frac{\color{blue}{\frac{8}{{\alpha}^{3}} + \left(\frac{2}{\alpha} - \frac{4}{\alpha \cdot \alpha}\right)}}{2}\]

   if 3.1045481599489194e+193 < alpha < 2.5715187390361886e+244

   1. Initial program 64.0

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

   3. Applied *-un-lft-identity 64.0

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

   4. Applied *-un-lft-identity 64.0

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

   5. Applied times-frac 44.8

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

   6. Applied times-frac 44.8

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

   7. Simplified 44.8

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

   9. Applied add-cbrt-cube 44.8

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

   10. Simplified 44.8

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

   12. Applied add-sqr-sqrt 45.0

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

   13. Applied add-cube-cbrt 45.2

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

   14. Applied *-un-lft-identity 45.2

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

   15. Applied times-frac 45.1

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

   16. Applied times-frac 47.2

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

   17. Applied associate-*r* 47.2

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

   18. Simplified 47.4

       \[\leadsto \frac{\sqrt[3]{{\left(\color{blue}{\frac{\alpha + \beta}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} \cdot \left(\sqrt[3]{\left(\alpha + \beta\right) + 2 \cdot i} \cdot \sqrt[3]{\left(\alpha + \beta\right) + 2 \cdot i}\right)}} \cdot \frac{\frac{\beta - \alpha}{\sqrt[3]{\left(\alpha + \beta\right) + 2 \cdot i}}}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}} + 1\right)}^{3}}}{2}\]
3. Recombined 3 regimes into one program.

4. Final simplification 12.2

   \[\leadsto \begin{array}{l} \mathbf{if}\;\alpha \le 3.068220085942097902687466146641291923906 \cdot 10^{110}:\\ \;\;\;\;\frac{\left(\alpha + \beta\right) \cdot \frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} + 1}{2}\\ \mathbf{elif}\;\alpha \le 3.104548159948919354427366616967038689222 \cdot 10^{193} \lor \neg \left(\alpha \le 2.571518739036188604595402311839971481899 \cdot 10^{244}\right):\\ \;\;\;\;\frac{\frac{8}{{\alpha}^{3}} + \left(\frac{2}{\alpha} - \frac{4}{\alpha \cdot \alpha}\right)}{2}\\ \mathbf{else}:\\ \;\;\;\;\frac{\sqrt[3]{{\left(\frac{\alpha + \beta}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} \cdot \left(\sqrt[3]{\left(\alpha + \beta\right) + 2 \cdot i} \cdot \sqrt[3]{\left(\alpha + \beta\right) + 2 \cdot i}\right)} \cdot \frac{\frac{\beta - \alpha}{\sqrt[3]{\left(\alpha + \beta\right) + 2 \cdot i}}}{\sqrt{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2}} + 1\right)}^{3}}}{2}\\ \end{array}\]

Reproduce

herbie shell --seed 2019347 
(FPCore (alpha beta i)
  :name "Octave 3.8, jcobi/2"
  :precision binary64
  :pre (and (> alpha -1) (> beta -1) (> i 0.0))
  (/ (+ (/ (/ (* (+ alpha beta) (- beta alpha)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) 2)) 1) 2))