| Alternative 1 | |
|---|---|
| Accuracy | 98.1% |
| Cost | 20232 |
(FPCore (x n) :precision binary64 (- (pow (+ x 1.0) (/ 1.0 n)) (pow x (/ 1.0 n))))
(FPCore (x n)
:precision binary64
(if (<= (/ 1.0 n) -5e-10)
(/ (exp (/ (log x) n)) (* n x))
(if (<= (/ 1.0 n) 5e-14)
(/ (log1p (/ 1.0 x)) n)
(-
(exp (* (* 0.3333333333333333 (/ (log1p x) n)) 3.0))
(pow x (/ 1.0 n))))))double code(double x, double n) {
return pow((x + 1.0), (1.0 / n)) - pow(x, (1.0 / n));
}
double code(double x, double n) {
double tmp;
if ((1.0 / n) <= -5e-10) {
tmp = exp((log(x) / n)) / (n * x);
} else if ((1.0 / n) <= 5e-14) {
tmp = log1p((1.0 / x)) / n;
} else {
tmp = exp(((0.3333333333333333 * (log1p(x) / n)) * 3.0)) - pow(x, (1.0 / n));
}
return tmp;
}
public static double code(double x, double n) {
return Math.pow((x + 1.0), (1.0 / n)) - Math.pow(x, (1.0 / n));
}
public static double code(double x, double n) {
double tmp;
if ((1.0 / n) <= -5e-10) {
tmp = Math.exp((Math.log(x) / n)) / (n * x);
} else if ((1.0 / n) <= 5e-14) {
tmp = Math.log1p((1.0 / x)) / n;
} else {
tmp = Math.exp(((0.3333333333333333 * (Math.log1p(x) / n)) * 3.0)) - Math.pow(x, (1.0 / n));
}
return tmp;
}
def code(x, n): return math.pow((x + 1.0), (1.0 / n)) - math.pow(x, (1.0 / n))
def code(x, n): tmp = 0 if (1.0 / n) <= -5e-10: tmp = math.exp((math.log(x) / n)) / (n * x) elif (1.0 / n) <= 5e-14: tmp = math.log1p((1.0 / x)) / n else: tmp = math.exp(((0.3333333333333333 * (math.log1p(x) / n)) * 3.0)) - math.pow(x, (1.0 / n)) return tmp
function code(x, n) return Float64((Float64(x + 1.0) ^ Float64(1.0 / n)) - (x ^ Float64(1.0 / n))) end
function code(x, n) tmp = 0.0 if (Float64(1.0 / n) <= -5e-10) tmp = Float64(exp(Float64(log(x) / n)) / Float64(n * x)); elseif (Float64(1.0 / n) <= 5e-14) tmp = Float64(log1p(Float64(1.0 / x)) / n); else tmp = Float64(exp(Float64(Float64(0.3333333333333333 * Float64(log1p(x) / n)) * 3.0)) - (x ^ Float64(1.0 / n))); end return tmp end
code[x_, n_] := N[(N[Power[N[(x + 1.0), $MachinePrecision], N[(1.0 / n), $MachinePrecision]], $MachinePrecision] - N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
code[x_, n_] := If[LessEqual[N[(1.0 / n), $MachinePrecision], -5e-10], N[(N[Exp[N[(N[Log[x], $MachinePrecision] / n), $MachinePrecision]], $MachinePrecision] / N[(n * x), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(1.0 / n), $MachinePrecision], 5e-14], N[(N[Log[1 + N[(1.0 / x), $MachinePrecision]], $MachinePrecision] / n), $MachinePrecision], N[(N[Exp[N[(N[(0.3333333333333333 * N[(N[Log[1 + x], $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision] * 3.0), $MachinePrecision]], $MachinePrecision] - N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]
{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}
\begin{array}{l}
\mathbf{if}\;\frac{1}{n} \leq -5 \cdot 10^{-10}:\\
\;\;\;\;\frac{e^{\frac{\log x}{n}}}{n \cdot x}\\
\mathbf{elif}\;\frac{1}{n} \leq 5 \cdot 10^{-14}:\\
\;\;\;\;\frac{\mathsf{log1p}\left(\frac{1}{x}\right)}{n}\\
\mathbf{else}:\\
\;\;\;\;e^{\left(0.3333333333333333 \cdot \frac{\mathsf{log1p}\left(x\right)}{n}\right) \cdot 3} - {x}^{\left(\frac{1}{n}\right)}\\
\end{array}
Results
if (/.f64 1 n) < -5.00000000000000031e-10Initial program 96.9%
Applied egg-rr96.8%
[Start]96.9 | \[ {\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}
\] |
|---|---|
add-cube-cbrt [=>]96.8 | \[ \color{blue}{\left(\sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)}} \cdot \sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)}}\right) \cdot \sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)}}} - {x}^{\left(\frac{1}{n}\right)}
\] |
pow3 [=>]96.8 | \[ \color{blue}{{\left(\sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)}}\right)}^{3}} - {x}^{\left(\frac{1}{n}\right)}
\] |
pow-to-exp [=>]96.8 | \[ \color{blue}{e^{\log \left(\sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)}}\right) \cdot 3}} - {x}^{\left(\frac{1}{n}\right)}
\] |
pow1/3 [=>]96.8 | \[ e^{\log \color{blue}{\left({\left({\left(x + 1\right)}^{\left(\frac{1}{n}\right)}\right)}^{0.3333333333333333}\right)} \cdot 3} - {x}^{\left(\frac{1}{n}\right)}
\] |
log-pow [=>]96.8 | \[ e^{\color{blue}{\left(0.3333333333333333 \cdot \log \left({\left(x + 1\right)}^{\left(\frac{1}{n}\right)}\right)\right)} \cdot 3} - {x}^{\left(\frac{1}{n}\right)}
\] |
pow-to-exp [=>]96.8 | \[ e^{\left(0.3333333333333333 \cdot \log \color{blue}{\left(e^{\log \left(x + 1\right) \cdot \frac{1}{n}}\right)}\right) \cdot 3} - {x}^{\left(\frac{1}{n}\right)}
\] |
add-log-exp [<=]96.8 | \[ e^{\left(0.3333333333333333 \cdot \color{blue}{\left(\log \left(x + 1\right) \cdot \frac{1}{n}\right)}\right) \cdot 3} - {x}^{\left(\frac{1}{n}\right)}
\] |
un-div-inv [=>]96.8 | \[ e^{\left(0.3333333333333333 \cdot \color{blue}{\frac{\log \left(x + 1\right)}{n}}\right) \cdot 3} - {x}^{\left(\frac{1}{n}\right)}
\] |
+-commutative [=>]96.8 | \[ e^{\left(0.3333333333333333 \cdot \frac{\log \color{blue}{\left(1 + x\right)}}{n}\right) \cdot 3} - {x}^{\left(\frac{1}{n}\right)}
\] |
log1p-udef [<=]96.8 | \[ e^{\left(0.3333333333333333 \cdot \frac{\color{blue}{\mathsf{log1p}\left(x\right)}}{n}\right) \cdot 3} - {x}^{\left(\frac{1}{n}\right)}
\] |
Taylor expanded in x around inf 97.5%
Simplified97.5%
[Start]97.5 | \[ \frac{e^{-1 \cdot \frac{\log \left(\frac{1}{x}\right)}{n}}}{n \cdot x}
\] |
|---|---|
log-rec [=>]97.5 | \[ \frac{e^{-1 \cdot \frac{\color{blue}{-\log x}}{n}}}{n \cdot x}
\] |
mul-1-neg [<=]97.5 | \[ \frac{e^{-1 \cdot \frac{\color{blue}{-1 \cdot \log x}}{n}}}{n \cdot x}
\] |
associate-*r/ [=>]97.5 | \[ \frac{e^{\color{blue}{\frac{-1 \cdot \left(-1 \cdot \log x\right)}{n}}}}{n \cdot x}
\] |
neg-mul-1 [<=]97.5 | \[ \frac{e^{\frac{\color{blue}{--1 \cdot \log x}}{n}}}{n \cdot x}
\] |
mul-1-neg [=>]97.5 | \[ \frac{e^{\frac{-\color{blue}{\left(-\log x\right)}}{n}}}{n \cdot x}
\] |
remove-double-neg [=>]97.5 | \[ \frac{e^{\frac{\color{blue}{\log x}}{n}}}{n \cdot x}
\] |
*-commutative [=>]97.5 | \[ \frac{e^{\frac{\log x}{n}}}{\color{blue}{x \cdot n}}
\] |
if -5.00000000000000031e-10 < (/.f64 1 n) < 5.0000000000000002e-14Initial program 30.1%
Taylor expanded in n around inf 77.5%
Simplified77.5%
[Start]77.5 | \[ \frac{\log \left(1 + x\right) - \log x}{n}
\] |
|---|---|
log1p-def [=>]77.5 | \[ \frac{\color{blue}{\mathsf{log1p}\left(x\right)} - \log x}{n}
\] |
Applied egg-rr77.6%
[Start]77.5 | \[ \frac{\mathsf{log1p}\left(x\right) - \log x}{n}
\] |
|---|---|
log1p-udef [=>]77.5 | \[ \frac{\color{blue}{\log \left(1 + x\right)} - \log x}{n}
\] |
diff-log [=>]77.6 | \[ \frac{\color{blue}{\log \left(\frac{1 + x}{x}\right)}}{n}
\] |
+-commutative [=>]77.6 | \[ \frac{\log \left(\frac{\color{blue}{x + 1}}{x}\right)}{n}
\] |
Taylor expanded in n around 0 77.6%
Simplified99.3%
[Start]77.6 | \[ \frac{\log \left(\frac{1 + x}{x}\right)}{n}
\] |
|---|---|
log-div [=>]77.5 | \[ \frac{\color{blue}{\log \left(1 + x\right) - \log x}}{n}
\] |
log1p-def [=>]77.5 | \[ \frac{\color{blue}{\mathsf{log1p}\left(x\right)} - \log x}{n}
\] |
log1p-def [<=]77.5 | \[ \frac{\color{blue}{\log \left(1 + x\right)} - \log x}{n}
\] |
log-div [<=]77.6 | \[ \frac{\color{blue}{\log \left(\frac{1 + x}{x}\right)}}{n}
\] |
+-commutative [<=]77.6 | \[ \frac{\log \left(\frac{\color{blue}{x + 1}}{x}\right)}{n}
\] |
*-lft-identity [<=]77.6 | \[ \frac{\log \left(\frac{\color{blue}{1 \cdot \left(x + 1\right)}}{x}\right)}{n}
\] |
associate-*l/ [<=]74.0 | \[ \frac{\log \color{blue}{\left(\frac{1}{x} \cdot \left(x + 1\right)\right)}}{n}
\] |
distribute-rgt-in [=>]74.0 | \[ \frac{\log \color{blue}{\left(x \cdot \frac{1}{x} + 1 \cdot \frac{1}{x}\right)}}{n}
\] |
*-lft-identity [=>]74.0 | \[ \frac{\log \left(x \cdot \frac{1}{x} + \color{blue}{\frac{1}{x}}\right)}{n}
\] |
rgt-mult-inverse [=>]77.6 | \[ \frac{\log \left(\color{blue}{1} + \frac{1}{x}\right)}{n}
\] |
log1p-def [=>]99.3 | \[ \frac{\color{blue}{\mathsf{log1p}\left(\frac{1}{x}\right)}}{n}
\] |
if 5.0000000000000002e-14 < (/.f64 1 n) Initial program 86.3%
Applied egg-rr91.3%
[Start]86.3 | \[ {\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}
\] |
|---|---|
add-cube-cbrt [=>]86.1 | \[ \color{blue}{\left(\sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)}} \cdot \sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)}}\right) \cdot \sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)}}} - {x}^{\left(\frac{1}{n}\right)}
\] |
pow3 [=>]86.1 | \[ \color{blue}{{\left(\sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)}}\right)}^{3}} - {x}^{\left(\frac{1}{n}\right)}
\] |
pow-to-exp [=>]86.1 | \[ \color{blue}{e^{\log \left(\sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)}}\right) \cdot 3}} - {x}^{\left(\frac{1}{n}\right)}
\] |
pow1/3 [=>]86.1 | \[ e^{\log \color{blue}{\left({\left({\left(x + 1\right)}^{\left(\frac{1}{n}\right)}\right)}^{0.3333333333333333}\right)} \cdot 3} - {x}^{\left(\frac{1}{n}\right)}
\] |
log-pow [=>]86.2 | \[ e^{\color{blue}{\left(0.3333333333333333 \cdot \log \left({\left(x + 1\right)}^{\left(\frac{1}{n}\right)}\right)\right)} \cdot 3} - {x}^{\left(\frac{1}{n}\right)}
\] |
pow-to-exp [=>]86.2 | \[ e^{\left(0.3333333333333333 \cdot \log \color{blue}{\left(e^{\log \left(x + 1\right) \cdot \frac{1}{n}}\right)}\right) \cdot 3} - {x}^{\left(\frac{1}{n}\right)}
\] |
add-log-exp [<=]86.2 | \[ e^{\left(0.3333333333333333 \cdot \color{blue}{\left(\log \left(x + 1\right) \cdot \frac{1}{n}\right)}\right) \cdot 3} - {x}^{\left(\frac{1}{n}\right)}
\] |
un-div-inv [=>]86.2 | \[ e^{\left(0.3333333333333333 \cdot \color{blue}{\frac{\log \left(x + 1\right)}{n}}\right) \cdot 3} - {x}^{\left(\frac{1}{n}\right)}
\] |
+-commutative [=>]86.2 | \[ e^{\left(0.3333333333333333 \cdot \frac{\log \color{blue}{\left(1 + x\right)}}{n}\right) \cdot 3} - {x}^{\left(\frac{1}{n}\right)}
\] |
log1p-udef [<=]91.3 | \[ e^{\left(0.3333333333333333 \cdot \frac{\color{blue}{\mathsf{log1p}\left(x\right)}}{n}\right) \cdot 3} - {x}^{\left(\frac{1}{n}\right)}
\] |
Final simplification98.1%
| Alternative 1 | |
|---|---|
| Accuracy | 98.1% |
| Cost | 20232 |
| Alternative 2 | |
|---|---|
| Accuracy | 97.5% |
| Cost | 14024 |
| Alternative 3 | |
|---|---|
| Accuracy | 97.4% |
| Cost | 13576 |
| Alternative 4 | |
|---|---|
| Accuracy | 97.5% |
| Cost | 13452 |
| Alternative 5 | |
|---|---|
| Accuracy | 97.8% |
| Cost | 8076 |
| Alternative 6 | |
|---|---|
| Accuracy | 97.8% |
| Cost | 7436 |
| Alternative 7 | |
|---|---|
| Accuracy | 97.4% |
| Cost | 7180 |
| Alternative 8 | |
|---|---|
| Accuracy | 97.6% |
| Cost | 7180 |
| Alternative 9 | |
|---|---|
| Accuracy | 89.0% |
| Cost | 6985 |
| Alternative 10 | |
|---|---|
| Accuracy | 74.9% |
| Cost | 6788 |
| Alternative 11 | |
|---|---|
| Accuracy | 54.5% |
| Cost | 585 |
| Alternative 12 | |
|---|---|
| Accuracy | 55.2% |
| Cost | 585 |
| Alternative 13 | |
|---|---|
| Accuracy | 55.2% |
| Cost | 584 |
| Alternative 14 | |
|---|---|
| Accuracy | 39.3% |
| Cost | 64 |
herbie shell --seed 2023147
(FPCore (x n)
:name "2nthrt (problem 3.4.6)"
:precision binary64
(- (pow (+ x 1.0) (/ 1.0 n)) (pow x (/ 1.0 n))))