Average Error: 0.8 → 0.2
Time: 3.9s
Precision: 64
\[\left(\sqrt{\left(\frac{x}{\left(1\right)}\right)}\right) - \left(\sqrt{x}\right)\]
↓
\[\frac{1}{\sqrt{x + 1} + \sqrt{x}}\]
\left(\sqrt{\left(\frac{x}{\left(1\right)}\right)}\right) - \left(\sqrt{x}\right)
↓
\frac{1}{\sqrt{x + 1} + \sqrt{x}}
double f(double x) {
        double r1658814 = x;
        double r1658815 = 1.0;
        double r1658816 = /* ERROR: no posit support in C */;
        double r1658817 = r1658814 + r1658816;
        double r1658818 = sqrt(r1658817);
        double r1658819 = sqrt(r1658814);
        double r1658820 = r1658818 - r1658819;
        return r1658820;
}
↓
double f(double x) {
        double r1658821 = 1.0;
        double r1658822 = x;
        double r1658823 = r1658822 + r1658821;
        double r1658824 = sqrt(r1658823);
        double r1658825 = sqrt(r1658822);
        double r1658826 = r1658824 + r1658825;
        double r1658827 = r1658821 / r1658826;
        return r1658827;
}

Error

Bits error versus x

Derivation

  1. Initial program 0.8

    \[\left(\sqrt{\left(\frac{x}{\left(1\right)}\right)}\right) - \left(\sqrt{x}\right)\]
  2. Using strategy rm

  3. Applied p16-flip--0.6

    \[\leadsto \color{blue}{\frac{\left(\left(\left(\sqrt{\left(\frac{x}{\left(1\right)}\right)}\right) \cdot \left(\sqrt{\left(\frac{x}{\left(1\right)}\right)}\right)\right) - \left(\left(\sqrt{x}\right) \cdot \left(\sqrt{x}\right)\right)\right)}{\left(\frac{\left(\sqrt{\left(\frac{x}{\left(1\right)}\right)}\right)}{\left(\sqrt{x}\right)}\right)}}\]

  4. Simplified0.2

    \[\leadsto \frac{\color{blue}{\left(1\right)}}{\left(\frac{\left(\sqrt{\left(\frac{x}{\left(1\right)}\right)}\right)}{\left(\sqrt{x}\right)}\right)}\]

  5. Final simplification0.2

    \[\leadsto \frac{1}{\sqrt{x + 1} + \sqrt{x}}\]

Reproduce

herbie shell --seed 2019125 +o rules:numerics
(FPCore (x)
  :name "2sqrt (example 3.1)"
  (-.p16 (sqrt.p16 (+.p16 x (real->posit16 1))) (sqrt.p16 x)))