Average Error: 3.3 → 1.4
Time: 1.3m
Precision: 64
\[\alpha \gt \left(-1\right) \land \beta \gt \left(-1\right) \land i \gt \left(1\right)\]
\[\frac{\left(\frac{\left(\left(i \cdot \left(\frac{\left(\frac{\alpha}{\beta}\right)}{i}\right)\right) \cdot \left(\frac{\left(\beta \cdot \alpha\right)}{\left(i \cdot \left(\frac{\left(\frac{\alpha}{\beta}\right)}{i}\right)\right)}\right)\right)}{\left(\left(\frac{\left(\frac{\alpha}{\beta}\right)}{\left(\left(2\right) \cdot i\right)}\right) \cdot \left(\frac{\left(\frac{\alpha}{\beta}\right)}{\left(\left(2\right) \cdot i\right)}\right)\right)}\right)}{\left(\left(\left(\frac{\left(\frac{\alpha}{\beta}\right)}{\left(\left(2\right) \cdot i\right)}\right) \cdot \left(\frac{\left(\frac{\alpha}{\beta}\right)}{\left(\left(2\right) \cdot i\right)}\right)\right) - \left(1.0\right)\right)}\]
↓
\[\left(\frac{i}{1.0 + \left(\left(2 \cdot i + \alpha\right) + \beta\right)} \cdot \frac{\alpha + \left(\beta + i\right)}{\left(2 \cdot i + \alpha\right) + \beta}\right) \cdot \frac{1.0}{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1.0}{\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)} \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}\]
\frac{\left(\frac{\left(\left(i \cdot \left(\frac{\left(\frac{\alpha}{\beta}\right)}{i}\right)\right) \cdot \left(\frac{\left(\beta \cdot \alpha\right)}{\left(i \cdot \left(\frac{\left(\frac{\alpha}{\beta}\right)}{i}\right)\right)}\right)\right)}{\left(\left(\frac{\left(\frac{\alpha}{\beta}\right)}{\left(\left(2\right) \cdot i\right)}\right) \cdot \left(\frac{\left(\frac{\alpha}{\beta}\right)}{\left(\left(2\right) \cdot i\right)}\right)\right)}\right)}{\left(\left(\left(\frac{\left(\frac{\alpha}{\beta}\right)}{\left(\left(2\right) \cdot i\right)}\right) \cdot \left(\frac{\left(\frac{\alpha}{\beta}\right)}{\left(\left(2\right) \cdot i\right)}\right)\right) - \left(1.0\right)\right)}
↓
\left(\frac{i}{1.0 + \left(\left(2 \cdot i + \alpha\right) + \beta\right)} \cdot \frac{\alpha + \left(\beta + i\right)}{\left(2 \cdot i + \alpha\right) + \beta}\right) \cdot \frac{1.0}{\frac{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1.0}{\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)} \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}
double f(double alpha, double beta, double i) {
        double r4135534 = i;
        double r4135535 = alpha;
        double r4135536 = beta;
        double r4135537 = r4135535 + r4135536;
        double r4135538 = r4135537 + r4135534;
        double r4135539 = r4135534 * r4135538;
        double r4135540 = r4135536 * r4135535;
        double r4135541 = r4135540 + r4135539;
        double r4135542 = r4135539 * r4135541;
        double r4135543 = 2.0;
        double r4135544 = /* ERROR: no posit support in C */;
        double r4135545 = r4135544 * r4135534;
        double r4135546 = r4135537 + r4135545;
        double r4135547 = r4135546 * r4135546;
        double r4135548 = r4135542 / r4135547;
        double r4135549 = 1.0;
        double r4135550 = /* ERROR: no posit support in C */;
        double r4135551 = r4135547 - r4135550;
        double r4135552 = r4135548 / r4135551;
        return r4135552;
}
↓
double f(double alpha, double beta, double i) {
        double r4135553 = i;
        double r4135554 = 1.0;
        double r4135555 = 2.0;
        double r4135556 = r4135555 * r4135553;
        double r4135557 = alpha;
        double r4135558 = r4135556 + r4135557;
        double r4135559 = beta;
        double r4135560 = r4135558 + r4135559;
        double r4135561 = r4135554 + r4135560;
        double r4135562 = r4135553 / r4135561;
        double r4135563 = r4135559 + r4135553;
        double r4135564 = r4135557 + r4135563;
        double r4135565 = r4135564 / r4135560;
        double r4135566 = r4135562 * r4135565;
        double r4135567 = r4135557 + r4135559;
        double r4135568 = r4135567 + r4135556;
        double r4135569 = r4135568 - r4135554;
        double r4135570 = r4135559 * r4135557;
        double r4135571 = r4135567 + r4135553;
        double r4135572 = r4135553 * r4135571;
        double r4135573 = r4135570 + r4135572;
        double r4135574 = r4135569 / r4135573;
        double r4135575 = r4135574 * r4135568;
        double r4135576 = r4135554 / r4135575;
        double r4135577 = r4135566 * r4135576;
        return r4135577;
}

Error

Bits error versus alpha

Bits error versus beta

Bits error versus i

Derivation

  1. Initial program 3.3

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

  3. Applied p16-*-un-lft-identity3.3

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

  4. Applied difference-of-squares3.3

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

  5. Applied p16-times-frac1.7

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

  6. Applied p16-times-frac1.7

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

  8. Applied *p16-rgt-identity-expand1.7

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

  9. Applied distribute-lft1-in1.7

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

  10. Applied p16-*-un-lft-identity1.7

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

  11. Applied p16-times-frac1.5

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

  12. Applied p16-times-frac1.5

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

  13. Simplified1.4

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

  14. Simplified1.4

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

  16. Applied p16-*-un-lft-identity1.4

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

  17. Applied associate-/l*1.4

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

  19. Applied associate-/r/1.4

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

  20. Final simplification1.4

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

Reproduce

herbie shell --seed 2019158 
(FPCore (alpha beta i)
  :name "Octave 3.8, jcobi/4"
  :pre (and (>.p16 alpha (real->posit16 -1)) (>.p16 beta (real->posit16 -1)) (>.p16 i (real->posit16 1)))
  (/.p16 (/.p16 (*.p16 (*.p16 i (+.p16 (+.p16 alpha beta) i)) (+.p16 (*.p16 beta alpha) (*.p16 i (+.p16 (+.p16 alpha beta) i)))) (*.p16 (+.p16 (+.p16 alpha beta) (*.p16 (real->posit16 2) i)) (+.p16 (+.p16 alpha beta) (*.p16 (real->posit16 2) i)))) (-.p16 (*.p16 (+.p16 (+.p16 alpha beta) (*.p16 (real->posit16 2) i)) (+.p16 (+.p16 alpha beta) (*.p16 (real->posit16 2) i))) (real->posit16 1.0))))