| Alternative 1 | |
|---|---|
| Accuracy | 99.4% |
| Cost | 39300 |

(FPCore (x) :precision binary64 (- (cbrt (+ x 1.0)) (cbrt x)))
(FPCore (x)
:precision binary64
(if (<= x 4e+54)
(/
1.0
(fma (cbrt x) (+ (cbrt (+ 1.0 x)) (cbrt x)) (cbrt (pow (+ 1.0 x) 2.0))))
(/ (* (cbrt (/ -1.0 (/ -1.0 x))) 0.3333333333333333) x)))double code(double x) {
return cbrt((x + 1.0)) - cbrt(x);
}
double code(double x) {
double tmp;
if (x <= 4e+54) {
tmp = 1.0 / fma(cbrt(x), (cbrt((1.0 + x)) + cbrt(x)), cbrt(pow((1.0 + x), 2.0)));
} else {
tmp = (cbrt((-1.0 / (-1.0 / x))) * 0.3333333333333333) / x;
}
return tmp;
}
function code(x) return Float64(cbrt(Float64(x + 1.0)) - cbrt(x)) end
function code(x) tmp = 0.0 if (x <= 4e+54) tmp = Float64(1.0 / fma(cbrt(x), Float64(cbrt(Float64(1.0 + x)) + cbrt(x)), cbrt((Float64(1.0 + x) ^ 2.0)))); else tmp = Float64(Float64(cbrt(Float64(-1.0 / Float64(-1.0 / x))) * 0.3333333333333333) / x); end return tmp end
code[x_] := N[(N[Power[N[(x + 1.0), $MachinePrecision], 1/3], $MachinePrecision] - N[Power[x, 1/3], $MachinePrecision]), $MachinePrecision]
code[x_] := If[LessEqual[x, 4e+54], N[(1.0 / N[(N[Power[x, 1/3], $MachinePrecision] * N[(N[Power[N[(1.0 + x), $MachinePrecision], 1/3], $MachinePrecision] + N[Power[x, 1/3], $MachinePrecision]), $MachinePrecision] + N[Power[N[Power[N[(1.0 + x), $MachinePrecision], 2.0], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[Power[N[(-1.0 / N[(-1.0 / x), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] * 0.3333333333333333), $MachinePrecision] / x), $MachinePrecision]]
\sqrt[3]{x + 1} - \sqrt[3]{x}
\begin{array}{l}
\mathbf{if}\;x \leq 4 \cdot 10^{+54}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(\sqrt[3]{x}, \sqrt[3]{1 + x} + \sqrt[3]{x}, \sqrt[3]{{\left(1 + x\right)}^{2}}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt[3]{\frac{-1}{\frac{-1}{x}}} \cdot 0.3333333333333333}{x}\\
\end{array}
Herbie found 15 alternatives:
| Alternative | Accuracy | Speedup |
|---|
if x < 4.0000000000000003e54Initial program 82.1%
Applied egg-rr85.0%
[Start]82.1% | \[ \sqrt[3]{x + 1} - \sqrt[3]{x}
\] |
|---|---|
flip3-- [=>]82.4% | \[ \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)}}
\] |
div-inv [=>]82.4% | \[ \color{blue}{\left({\left(\sqrt[3]{x + 1}\right)}^{3} - {\left(\sqrt[3]{x}\right)}^{3}\right) \cdot \frac{1}{\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)}}
\] |
rem-cube-cbrt [=>]82.8% | \[ \left(\color{blue}{\left(x + 1\right)} - {\left(\sqrt[3]{x}\right)}^{3}\right) \cdot \frac{1}{\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)}
\] |
rem-cube-cbrt [=>]84.9% | \[ \left(\left(x + 1\right) - \color{blue}{x}\right) \cdot \frac{1}{\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)}
\] |
cbrt-unprod [=>]85.0% | \[ \left(\left(x + 1\right) - x\right) \cdot \frac{1}{\color{blue}{\sqrt[3]{\left(x + 1\right) \cdot \left(x + 1\right)}} + \left(\sqrt[3]{x} \cdot \sqrt[3]{x} + \sqrt[3]{x + 1} \cdot \sqrt[3]{x}\right)}
\] |
pow2 [=>]85.0% | \[ \left(\left(x + 1\right) - x\right) \cdot \frac{1}{\sqrt[3]{\color{blue}{{\left(x + 1\right)}^{2}}} + \left(\sqrt[3]{x} \cdot \sqrt[3]{x} + \sqrt[3]{x + 1} \cdot \sqrt[3]{x}\right)}
\] |
distribute-rgt-out [=>]85.0% | \[ \left(\left(x + 1\right) - x\right) \cdot \frac{1}{\sqrt[3]{{\left(x + 1\right)}^{2}} + \color{blue}{\sqrt[3]{x} \cdot \left(\sqrt[3]{x} + \sqrt[3]{x + 1}\right)}}
\] |
+-commutative [<=]85.0% | \[ \left(\left(x + 1\right) - x\right) \cdot \frac{1}{\sqrt[3]{{\left(x + 1\right)}^{2}} + \sqrt[3]{x} \cdot \color{blue}{\left(\sqrt[3]{x + 1} + \sqrt[3]{x}\right)}}
\] |
Simplified99.7%
[Start]85.0% | \[ \left(\left(x + 1\right) - x\right) \cdot \frac{1}{\sqrt[3]{{\left(x + 1\right)}^{2}} + \sqrt[3]{x} \cdot \left(\sqrt[3]{x + 1} + \sqrt[3]{x}\right)}
\] |
|---|---|
associate-*r/ [=>]85.0% | \[ \color{blue}{\frac{\left(\left(x + 1\right) - x\right) \cdot 1}{\sqrt[3]{{\left(x + 1\right)}^{2}} + \sqrt[3]{x} \cdot \left(\sqrt[3]{x + 1} + \sqrt[3]{x}\right)}}
\] |
*-rgt-identity [=>]85.0% | \[ \frac{\color{blue}{\left(x + 1\right) - x}}{\sqrt[3]{{\left(x + 1\right)}^{2}} + \sqrt[3]{x} \cdot \left(\sqrt[3]{x + 1} + \sqrt[3]{x}\right)}
\] |
+-commutative [=>]85.0% | \[ \frac{\color{blue}{\left(1 + x\right)} - x}{\sqrt[3]{{\left(x + 1\right)}^{2}} + \sqrt[3]{x} \cdot \left(\sqrt[3]{x + 1} + \sqrt[3]{x}\right)}
\] |
associate--l+ [=>]99.6% | \[ \frac{\color{blue}{1 + \left(x - x\right)}}{\sqrt[3]{{\left(x + 1\right)}^{2}} + \sqrt[3]{x} \cdot \left(\sqrt[3]{x + 1} + \sqrt[3]{x}\right)}
\] |
+-inverses [=>]99.6% | \[ \frac{1 + \color{blue}{0}}{\sqrt[3]{{\left(x + 1\right)}^{2}} + \sqrt[3]{x} \cdot \left(\sqrt[3]{x + 1} + \sqrt[3]{x}\right)}
\] |
metadata-eval [=>]99.6% | \[ \frac{\color{blue}{1}}{\sqrt[3]{{\left(x + 1\right)}^{2}} + \sqrt[3]{x} \cdot \left(\sqrt[3]{x + 1} + \sqrt[3]{x}\right)}
\] |
+-commutative [=>]99.6% | \[ \frac{1}{\color{blue}{\sqrt[3]{x} \cdot \left(\sqrt[3]{x + 1} + \sqrt[3]{x}\right) + \sqrt[3]{{\left(x + 1\right)}^{2}}}}
\] |
fma-def [=>]99.7% | \[ \frac{1}{\color{blue}{\mathsf{fma}\left(\sqrt[3]{x}, \sqrt[3]{x + 1} + \sqrt[3]{x}, \sqrt[3]{{\left(x + 1\right)}^{2}}\right)}}
\] |
+-commutative [=>]99.7% | \[ \frac{1}{\mathsf{fma}\left(\sqrt[3]{x}, \sqrt[3]{\color{blue}{1 + x}} + \sqrt[3]{x}, \sqrt[3]{{\left(x + 1\right)}^{2}}\right)}
\] |
+-commutative [=>]99.7% | \[ \frac{1}{\mathsf{fma}\left(\sqrt[3]{x}, \sqrt[3]{1 + x} + \sqrt[3]{x}, \sqrt[3]{{\color{blue}{\left(1 + x\right)}}^{2}}\right)}
\] |
if 4.0000000000000003e54 < x Initial program 4.3%
Applied egg-rr4.3%
[Start]4.3% | \[ \sqrt[3]{x + 1} - \sqrt[3]{x}
\] |
|---|---|
add-sqr-sqrt [=>]3.7% | \[ \color{blue}{\sqrt{\sqrt[3]{x + 1}} \cdot \sqrt{\sqrt[3]{x + 1}}} - \sqrt[3]{x}
\] |
add-sqr-sqrt [=>]4.3% | \[ \sqrt{\sqrt[3]{x + 1}} \cdot \sqrt{\sqrt[3]{x + 1}} - \color{blue}{\sqrt{\sqrt[3]{x}} \cdot \sqrt{\sqrt[3]{x}}}
\] |
difference-of-squares [=>]4.3% | \[ \color{blue}{\left(\sqrt{\sqrt[3]{x + 1}} + \sqrt{\sqrt[3]{x}}\right) \cdot \left(\sqrt{\sqrt[3]{x + 1}} - \sqrt{\sqrt[3]{x}}\right)}
\] |
pow1/3 [=>]4.3% | \[ \left(\sqrt{\color{blue}{{\left(x + 1\right)}^{0.3333333333333333}}} + \sqrt{\sqrt[3]{x}}\right) \cdot \left(\sqrt{\sqrt[3]{x + 1}} - \sqrt{\sqrt[3]{x}}\right)
\] |
sqrt-pow1 [=>]4.3% | \[ \left(\color{blue}{{\left(x + 1\right)}^{\left(\frac{0.3333333333333333}{2}\right)}} + \sqrt{\sqrt[3]{x}}\right) \cdot \left(\sqrt{\sqrt[3]{x + 1}} - \sqrt{\sqrt[3]{x}}\right)
\] |
metadata-eval [=>]4.3% | \[ \left({\left(x + 1\right)}^{\color{blue}{0.16666666666666666}} + \sqrt{\sqrt[3]{x}}\right) \cdot \left(\sqrt{\sqrt[3]{x + 1}} - \sqrt{\sqrt[3]{x}}\right)
\] |
pow1/3 [=>]4.3% | \[ \left({\left(x + 1\right)}^{0.16666666666666666} + \sqrt{\color{blue}{{x}^{0.3333333333333333}}}\right) \cdot \left(\sqrt{\sqrt[3]{x + 1}} - \sqrt{\sqrt[3]{x}}\right)
\] |
sqrt-pow1 [=>]4.3% | \[ \left({\left(x + 1\right)}^{0.16666666666666666} + \color{blue}{{x}^{\left(\frac{0.3333333333333333}{2}\right)}}\right) \cdot \left(\sqrt{\sqrt[3]{x + 1}} - \sqrt{\sqrt[3]{x}}\right)
\] |
metadata-eval [=>]4.3% | \[ \left({\left(x + 1\right)}^{0.16666666666666666} + {x}^{\color{blue}{0.16666666666666666}}\right) \cdot \left(\sqrt{\sqrt[3]{x + 1}} - \sqrt{\sqrt[3]{x}}\right)
\] |
pow1/3 [=>]1.8% | \[ \left({\left(x + 1\right)}^{0.16666666666666666} + {x}^{0.16666666666666666}\right) \cdot \left(\sqrt{\color{blue}{{\left(x + 1\right)}^{0.3333333333333333}}} - \sqrt{\sqrt[3]{x}}\right)
\] |
sqrt-pow1 [=>]1.8% | \[ \left({\left(x + 1\right)}^{0.16666666666666666} + {x}^{0.16666666666666666}\right) \cdot \left(\color{blue}{{\left(x + 1\right)}^{\left(\frac{0.3333333333333333}{2}\right)}} - \sqrt{\sqrt[3]{x}}\right)
\] |
metadata-eval [=>]1.8% | \[ \left({\left(x + 1\right)}^{0.16666666666666666} + {x}^{0.16666666666666666}\right) \cdot \left({\left(x + 1\right)}^{\color{blue}{0.16666666666666666}} - \sqrt{\sqrt[3]{x}}\right)
\] |
pow1/3 [=>]4.2% | \[ \left({\left(x + 1\right)}^{0.16666666666666666} + {x}^{0.16666666666666666}\right) \cdot \left({\left(x + 1\right)}^{0.16666666666666666} - \sqrt{\color{blue}{{x}^{0.3333333333333333}}}\right)
\] |
sqrt-pow1 [=>]4.3% | \[ \left({\left(x + 1\right)}^{0.16666666666666666} + {x}^{0.16666666666666666}\right) \cdot \left({\left(x + 1\right)}^{0.16666666666666666} - \color{blue}{{x}^{\left(\frac{0.3333333333333333}{2}\right)}}\right)
\] |
metadata-eval [=>]4.3% | \[ \left({\left(x + 1\right)}^{0.16666666666666666} + {x}^{0.16666666666666666}\right) \cdot \left({\left(x + 1\right)}^{0.16666666666666666} - {x}^{\color{blue}{0.16666666666666666}}\right)
\] |
Taylor expanded in x around -inf 0.0%
Simplified99.1%
[Start]0.0% | \[ 0.3333333333333333 \cdot \frac{{\left(e^{0.16666666666666666 \cdot \left(-1 \cdot \log \left(\frac{-1}{x}\right) + \log -1\right)}\right)}^{2}}{x}
\] |
|---|---|
*-commutative [=>]0.0% | \[ \color{blue}{\frac{{\left(e^{0.16666666666666666 \cdot \left(-1 \cdot \log \left(\frac{-1}{x}\right) + \log -1\right)}\right)}^{2}}{x} \cdot 0.3333333333333333}
\] |
associate-*l/ [=>]0.0% | \[ \color{blue}{\frac{{\left(e^{0.16666666666666666 \cdot \left(-1 \cdot \log \left(\frac{-1}{x}\right) + \log -1\right)}\right)}^{2} \cdot 0.3333333333333333}{x}}
\] |
Applied egg-rr99.1%
[Start]99.1% | \[ \frac{\sqrt[3]{-1 \cdot {\left(\frac{-1}{x}\right)}^{-1}} \cdot 0.3333333333333333}{x}
\] |
|---|---|
unpow-1 [=>]99.1% | \[ \frac{\sqrt[3]{-1 \cdot \color{blue}{\frac{1}{\frac{-1}{x}}}} \cdot 0.3333333333333333}{x}
\] |
Final simplification99.5%
| Alternative 1 | |
|---|---|
| Accuracy | 99.4% |
| Cost | 39300 |
| Alternative 2 | |
|---|---|
| Accuracy | 99.2% |
| Cost | 58496 |
| Alternative 3 | |
|---|---|
| Accuracy | 99.1% |
| Cost | 52032 |
| Alternative 4 | |
|---|---|
| Accuracy | 99.4% |
| Cost | 40132 |
| Alternative 5 | |
|---|---|
| Accuracy | 99.2% |
| Cost | 39168 |
| Alternative 6 | |
|---|---|
| Accuracy | 99.5% |
| Cost | 32900 |
| Alternative 7 | |
|---|---|
| Accuracy | 99.1% |
| Cost | 26564 |
| Alternative 8 | |
|---|---|
| Accuracy | 99.1% |
| Cost | 26372 |
| Alternative 9 | |
|---|---|
| Accuracy | 99.1% |
| Cost | 26308 |
| Alternative 10 | |
|---|---|
| Accuracy | 98.1% |
| Cost | 7108 |
| Alternative 11 | |
|---|---|
| Accuracy | 98.1% |
| Cost | 6980 |
| Alternative 12 | |
|---|---|
| Accuracy | 97.6% |
| Cost | 6852 |
| Alternative 13 | |
|---|---|
| Accuracy | 52.3% |
| Cost | 6724 |
| Alternative 14 | |
|---|---|
| Accuracy | 3.6% |
| Cost | 64 |
| Alternative 15 | |
|---|---|
| Accuracy | 50.2% |
| Cost | 64 |
herbie shell --seed 2023263
(FPCore (x)
:name "2cbrt (problem 3.3.4)"
:precision binary64
(- (cbrt (+ x 1.0)) (cbrt x)))