\sqrt[3]{x + 1} - \sqrt[3]{x}\begin{array}{l}
\mathbf{if}\;x \leq -30244119.16774476:\\
\;\;\;\;0.3333333333333333 \cdot \sqrt[3]{\frac{1}{x \cdot x}}\\
\mathbf{elif}\;x \leq 0.02398425106893859:\\
\;\;\;\;\sqrt[3]{x + 1} - \sqrt[3]{\sqrt[3]{x} \cdot \left(\sqrt[3]{x} \cdot \sqrt[3]{x}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{{x}^{0.6666666666666666} + \sqrt[3]{x + 1} \cdot \left(\sqrt[3]{x + 1} + \sqrt[3]{x}\right)}\\
\end{array}(FPCore (x) :precision binary64 (- (cbrt (+ x 1.0)) (cbrt x)))
(FPCore (x)
:precision binary64
(if (<= x -30244119.16774476)
(* 0.3333333333333333 (cbrt (/ 1.0 (* x x))))
(if (<= x 0.02398425106893859)
(- (cbrt (+ x 1.0)) (cbrt (* (cbrt x) (* (cbrt x) (cbrt x)))))
(/
1.0
(+
(pow x 0.6666666666666666)
(* (cbrt (+ x 1.0)) (+ (cbrt (+ x 1.0)) (cbrt x))))))))double code(double x) {
return cbrt(x + 1.0) - cbrt(x);
}
double code(double x) {
double tmp;
if (x <= -30244119.16774476) {
tmp = 0.3333333333333333 * cbrt(1.0 / (x * x));
} else if (x <= 0.02398425106893859) {
tmp = cbrt(x + 1.0) - cbrt(cbrt(x) * (cbrt(x) * cbrt(x)));
} else {
tmp = 1.0 / (pow(x, 0.6666666666666666) + (cbrt(x + 1.0) * (cbrt(x + 1.0) + cbrt(x))));
}
return tmp;
}



Bits error versus x
Results
if x < -30244119.1677447595Initial program 60.9
Taylor expanded around inf 34.2
Simplified32.1
if -30244119.1677447595 < x < 0.023984251068938589Initial program 0.2
rmApplied add-cbrt-cube_binary64_4550.2
if 0.023984251068938589 < x Initial program 59.4
rmApplied flip3--_binary64_42359.3
Simplified1.0
Simplified4.5
Final simplification9.4
herbie shell --seed 2021098
(FPCore (x)
:name "2cbrt (problem 3.3.4)"
:precision binary64
(- (cbrt (+ x 1.0)) (cbrt x)))