Average Error: 3.8 → 2.7
Time: 2.7m
Precision: 64
\[\alpha \gt -1 \land \beta \gt -1\]
\[\frac{\frac{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1}\]
↓
\[\begin{array}{l} \mathbf{if}\;\alpha \le 9.316239619114071471332105808691819719617 \cdot 10^{188}:\\ \;\;\;\;\frac{\frac{\sqrt[3]{{\left(\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}{\left(\alpha + \beta\right) + 2 \cdot 1}\right)}^{3}}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1}\\ \mathbf{else}:\\ \;\;\;\;\frac{\frac{\sqrt[3]{0.1875 \cdot \alpha + \left(0.125 + 0.1875 \cdot \beta\right)}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1}\\ \end{array}\]
\frac{\frac{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1}
↓
\begin{array}{l}
\mathbf{if}\;\alpha \le 9.316239619114071471332105808691819719617 \cdot 10^{188}:\\
\;\;\;\;\frac{\frac{\sqrt[3]{{\left(\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}{\left(\alpha + \beta\right) + 2 \cdot 1}\right)}^{3}}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1}\\

\mathbf{else}:\\
\;\;\;\;\frac{\frac{\sqrt[3]{0.1875 \cdot \alpha + \left(0.125 + 0.1875 \cdot \beta\right)}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1}\\

\end{array}
double f(double alpha, double beta) {
        double r385435 = alpha;
        double r385436 = beta;
        double r385437 = r385435 + r385436;
        double r385438 = r385436 * r385435;
        double r385439 = r385437 + r385438;
        double r385440 = 1.0;
        double r385441 = r385439 + r385440;
        double r385442 = 2.0;
        double r385443 = r385442 * r385440;
        double r385444 = r385437 + r385443;
        double r385445 = r385441 / r385444;
        double r385446 = r385445 / r385444;
        double r385447 = r385444 + r385440;
        double r385448 = r385446 / r385447;
        return r385448;
}
↓
double f(double alpha, double beta) {
        double r385449 = alpha;
        double r385450 = 9.316239619114071e+188;
        bool r385451 = r385449 <= r385450;
        double r385452 = beta;
        double r385453 = r385449 + r385452;
        double r385454 = r385452 * r385449;
        double r385455 = r385453 + r385454;
        double r385456 = 1.0;
        double r385457 = r385455 + r385456;
        double r385458 = 2.0;
        double r385459 = r385458 * r385456;
        double r385460 = r385453 + r385459;
        double r385461 = r385457 / r385460;
        double r385462 = 3.0;
        double r385463 = pow(r385461, r385462);
        double r385464 = cbrt(r385463);
        double r385465 = r385464 / r385460;
        double r385466 = r385460 + r385456;
        double r385467 = r385465 / r385466;
        double r385468 = 0.1875;
        double r385469 = r385468 * r385449;
        double r385470 = 0.125;
        double r385471 = r385468 * r385452;
        double r385472 = r385470 + r385471;
        double r385473 = r385469 + r385472;
        double r385474 = cbrt(r385473);
        double r385475 = r385474 / r385460;
        double r385476 = r385475 / r385466;
        double r385477 = r385451 ? r385467 : r385476;
        return r385477;
}

Error

Bits error versus alpha

Bits error versus beta

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Split input into 2 regimes

  2. if alpha < 9.316239619114071e+188

    1. Initial program 1.7

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

    3. Applied add-cbrt-cube8.8

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

    4. Applied add-cbrt-cube20.4

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

    5. Applied cbrt-undiv20.4

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

    6. Simplified2.1

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

    if 9.316239619114071e+188 < alpha

    1. Initial program 18.3

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

    3. Applied add-cbrt-cube18.4

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

    4. Applied add-cbrt-cube64.0

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

    5. Applied cbrt-undiv64.0

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

    6. Simplified18.4

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

    7. Taylor expanded around 0 6.9

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

  4. Final simplification2.7

    \[\leadsto \begin{array}{l} \mathbf{if}\;\alpha \le 9.316239619114071471332105808691819719617 \cdot 10^{188}:\\ \;\;\;\;\frac{\frac{\sqrt[3]{{\left(\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}{\left(\alpha + \beta\right) + 2 \cdot 1}\right)}^{3}}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1}\\ \mathbf{else}:\\ \;\;\;\;\frac{\frac{\sqrt[3]{0.1875 \cdot \alpha + \left(0.125 + 0.1875 \cdot \beta\right)}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1}\\ \end{array}\]

Reproduce

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