2isqrt (example 3.6)

Average Error: 20.2 → 0.4

Time: 17.3s

Precision: 64

Math

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

↓

\[\frac{\frac{1}{1 \cdot \mathsf{fma}\left(\sqrt{\sqrt[3]{x + 1} \cdot \sqrt[3]{x + 1}}, \sqrt{\sqrt[3]{x + 1}}, \sqrt{x}\right)}}{\sqrt{x + 1} \cdot \sqrt{x}}\]

C

double f(double x) {
        double r153166 = 1.0;
        double r153167 = x;
        double r153168 = sqrt(r153167);
        double r153169 = r153166 / r153168;
        double r153170 = r153167 + r153166;
        double r153171 = sqrt(r153170);
        double r153172 = r153166 / r153171;
        double r153173 = r153169 - r153172;
        return r153173;
}

↓

double f(double x) {
        double r153174 = 1.0;
        double r153175 = x;
        double r153176 = r153175 + r153174;
        double r153177 = cbrt(r153176);
        double r153178 = r153177 * r153177;
        double r153179 = sqrt(r153178);
        double r153180 = sqrt(r153177);
        double r153181 = sqrt(r153175);
        double r153182 = fma(r153179, r153180, r153181);
        double r153183 = r153174 * r153182;
        double r153184 = r153174 / r153183;
        double r153185 = sqrt(r153176);
        double r153186 = r153185 * r153181;
        double r153187 = r153184 / r153186;
        return r153187;
}

Error

Bits error versus x

Target

Original	20.2
Target	0.7
Herbie	0.4

\[\frac{1}{\left(x + 1\right) \cdot \sqrt{x} + x \cdot \sqrt{x + 1}}\]

Derivation

1. Initial program 20.2

   \[\frac{1}{\sqrt{x}} - \frac{1}{\sqrt{x + 1}}\]
2. Using strategy rm

3. Applied frac-sub 20.2

   \[\leadsto \color{blue}{\frac{1 \cdot \sqrt{x + 1} - \sqrt{x} \cdot 1}{\sqrt{x} \cdot \sqrt{x + 1}}}\]

4. Simplified 20.2

   \[\leadsto \frac{1 \cdot \sqrt{x + 1} - \sqrt{x} \cdot 1}{\color{blue}{\sqrt{x + 1} \cdot \sqrt{x}}}\]
5. Using strategy rm

6. Applied flip-- 20.0

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

7. Simplified 19.6

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

8. Simplified 19.6

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

9. Taylor expanded around 0 0.4

   \[\leadsto \frac{\frac{\color{blue}{1}}{1 \cdot \left(\sqrt{x + 1} + \sqrt{x}\right)}}{\sqrt{x + 1} \cdot \sqrt{x}}\]
10. Using strategy rm

11. Applied add-cube-cbrt 0.4

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

12. Applied sqrt-prod 0.4

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

13. Applied fma-def 0.4

    \[\leadsto \frac{\frac{1}{1 \cdot \color{blue}{\mathsf{fma}\left(\sqrt{\sqrt[3]{x + 1} \cdot \sqrt[3]{x + 1}}, \sqrt{\sqrt[3]{x + 1}}, \sqrt{x}\right)}}}{\sqrt{x + 1} \cdot \sqrt{x}}\]

14. Final simplification 0.4

    \[\leadsto \frac{\frac{1}{1 \cdot \mathsf{fma}\left(\sqrt{\sqrt[3]{x + 1} \cdot \sqrt[3]{x + 1}}, \sqrt{\sqrt[3]{x + 1}}, \sqrt{x}\right)}}{\sqrt{x + 1} \cdot \sqrt{x}}\]

Reproduce

herbie shell --seed 2019194 +o rules:numerics
(FPCore (x)
  :name "2isqrt (example 3.6)"

  :herbie-target
  (/ 1.0 (+ (* (+ x 1.0) (sqrt x)) (* x (sqrt (+ x 1.0)))))

  (- (/ 1.0 (sqrt x)) (/ 1.0 (sqrt (+ x 1.0)))))