Average Error: 29.4 → 0.5
Time: 18.5s
Precision: 64
\[\sqrt[3]{x + 1} - \sqrt[3]{x}\]
\[\frac{1}{\mathsf{fma}\left(\sqrt[3]{x + 1}, \sqrt[3]{x + 1}, \left(\sqrt[3]{x + 1} + \sqrt[3]{x}\right) \cdot \sqrt[3]{x}\right)}\]
\sqrt[3]{x + 1} - \sqrt[3]{x}
\frac{1}{\mathsf{fma}\left(\sqrt[3]{x + 1}, \sqrt[3]{x + 1}, \left(\sqrt[3]{x + 1} + \sqrt[3]{x}\right) \cdot \sqrt[3]{x}\right)}
double f(double x) {
        double r2752443 = x;
        double r2752444 = 1.0;
        double r2752445 = r2752443 + r2752444;
        double r2752446 = cbrt(r2752445);
        double r2752447 = cbrt(r2752443);
        double r2752448 = r2752446 - r2752447;
        return r2752448;
}

double f(double x) {
        double r2752449 = 1.0;
        double r2752450 = x;
        double r2752451 = r2752450 + r2752449;
        double r2752452 = cbrt(r2752451);
        double r2752453 = cbrt(r2752450);
        double r2752454 = r2752452 + r2752453;
        double r2752455 = r2752454 * r2752453;
        double r2752456 = fma(r2752452, r2752452, r2752455);
        double r2752457 = r2752449 / r2752456;
        return r2752457;
}

Error

Bits error versus x

Derivation

  1. Initial program 29.4

    \[\sqrt[3]{x + 1} - \sqrt[3]{x}\]
  2. Using strategy rm
  3. Applied flip3--29.3

    \[\leadsto \color{blue}{\frac{{\left(\sqrt[3]{x + 1}\right)}^{3} - {\left(\sqrt[3]{x}\right)}^{3}}{\sqrt[3]{x + 1} \cdot \sqrt[3]{x + 1} + \left(\sqrt[3]{x} \cdot \sqrt[3]{x} + \sqrt[3]{x + 1} \cdot \sqrt[3]{x}\right)}}\]
  4. Simplified0.5

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

    \[\leadsto \frac{1 + 0}{\color{blue}{\mathsf{fma}\left(\sqrt[3]{1 + x}, \sqrt[3]{1 + x}, \sqrt[3]{x} \cdot \left(\sqrt[3]{1 + x} + \sqrt[3]{x}\right)\right)}}\]
  6. Final simplification0.5

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

Reproduce

herbie shell --seed 2019172 +o rules:numerics
(FPCore (x)
  :name "2cbrt (problem 3.3.4)"
  (- (cbrt (+ x 1.0)) (cbrt x)))