| Alternative 1 | |
|---|---|
| Accuracy | 82.7% |
| Cost | 22408 |
(FPCore (i n) :precision binary64 (* 100.0 (/ (- (pow (+ 1.0 (/ i n)) n) 1.0) (/ i n))))
(FPCore (i n)
:precision binary64
(let* ((t_0 (+ (pow (+ 1.0 (/ i n)) n) -1.0)) (t_1 (/ t_0 (/ i n))))
(if (<= t_1 5e-173)
(/ 100.0 (/ (/ i n) (expm1 (* n (log1p (/ i n))))))
(if (<= t_1 20000000000.0)
(/ (* n t_0) (/ i 100.0))
(* 100.0 (/ n (+ 1.0 (* i (+ -0.5 (* i 0.08333333333333333))))))))))double code(double i, double n) {
return 100.0 * ((pow((1.0 + (i / n)), n) - 1.0) / (i / n));
}
double code(double i, double n) {
double t_0 = pow((1.0 + (i / n)), n) + -1.0;
double t_1 = t_0 / (i / n);
double tmp;
if (t_1 <= 5e-173) {
tmp = 100.0 / ((i / n) / expm1((n * log1p((i / n)))));
} else if (t_1 <= 20000000000.0) {
tmp = (n * t_0) / (i / 100.0);
} else {
tmp = 100.0 * (n / (1.0 + (i * (-0.5 + (i * 0.08333333333333333)))));
}
return tmp;
}
public static double code(double i, double n) {
return 100.0 * ((Math.pow((1.0 + (i / n)), n) - 1.0) / (i / n));
}
public static double code(double i, double n) {
double t_0 = Math.pow((1.0 + (i / n)), n) + -1.0;
double t_1 = t_0 / (i / n);
double tmp;
if (t_1 <= 5e-173) {
tmp = 100.0 / ((i / n) / Math.expm1((n * Math.log1p((i / n)))));
} else if (t_1 <= 20000000000.0) {
tmp = (n * t_0) / (i / 100.0);
} else {
tmp = 100.0 * (n / (1.0 + (i * (-0.5 + (i * 0.08333333333333333)))));
}
return tmp;
}
def code(i, n): return 100.0 * ((math.pow((1.0 + (i / n)), n) - 1.0) / (i / n))
def code(i, n): t_0 = math.pow((1.0 + (i / n)), n) + -1.0 t_1 = t_0 / (i / n) tmp = 0 if t_1 <= 5e-173: tmp = 100.0 / ((i / n) / math.expm1((n * math.log1p((i / n))))) elif t_1 <= 20000000000.0: tmp = (n * t_0) / (i / 100.0) else: tmp = 100.0 * (n / (1.0 + (i * (-0.5 + (i * 0.08333333333333333))))) return tmp
function code(i, n) return Float64(100.0 * Float64(Float64((Float64(1.0 + Float64(i / n)) ^ n) - 1.0) / Float64(i / n))) end
function code(i, n) t_0 = Float64((Float64(1.0 + Float64(i / n)) ^ n) + -1.0) t_1 = Float64(t_0 / Float64(i / n)) tmp = 0.0 if (t_1 <= 5e-173) tmp = Float64(100.0 / Float64(Float64(i / n) / expm1(Float64(n * log1p(Float64(i / n)))))); elseif (t_1 <= 20000000000.0) tmp = Float64(Float64(n * t_0) / Float64(i / 100.0)); else tmp = Float64(100.0 * Float64(n / Float64(1.0 + Float64(i * Float64(-0.5 + Float64(i * 0.08333333333333333)))))); end return tmp end
code[i_, n_] := N[(100.0 * N[(N[(N[Power[N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision], n], $MachinePrecision] - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[i_, n_] := Block[{t$95$0 = N[(N[Power[N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision], n], $MachinePrecision] + -1.0), $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 / N[(i / n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, 5e-173], N[(100.0 / N[(N[(i / n), $MachinePrecision] / N[(Exp[N[(n * N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 20000000000.0], N[(N[(n * t$95$0), $MachinePrecision] / N[(i / 100.0), $MachinePrecision]), $MachinePrecision], N[(100.0 * N[(n / N[(1.0 + N[(i * N[(-0.5 + N[(i * 0.08333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}
\begin{array}{l}
t_0 := {\left(1 + \frac{i}{n}\right)}^{n} + -1\\
t_1 := \frac{t_0}{\frac{i}{n}}\\
\mathbf{if}\;t_1 \leq 5 \cdot 10^{-173}:\\
\;\;\;\;\frac{100}{\frac{\frac{i}{n}}{\mathsf{expm1}\left(n \cdot \mathsf{log1p}\left(\frac{i}{n}\right)\right)}}\\
\mathbf{elif}\;t_1 \leq 20000000000:\\
\;\;\;\;\frac{n \cdot t_0}{\frac{i}{100}}\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \frac{n}{1 + i \cdot \left(-0.5 + i \cdot 0.08333333333333333\right)}\\
\end{array}
Results
| Original | 25.5% |
|---|---|
| Target | 25.7% |
| Herbie | 98.3% |
if (/.f64 (-.f64 (pow.f64 (+.f64 1 (/.f64 i n)) n) 1) (/.f64 i n)) < 5.0000000000000002e-173Initial program 28.6%
Applied egg-rr98.1%
[Start]28.6 | \[ 100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}
\] |
|---|---|
associate-*r/ [=>]28.6 | \[ \color{blue}{\frac{100 \cdot \left({\left(1 + \frac{i}{n}\right)}^{n} - 1\right)}{\frac{i}{n}}}
\] |
associate-/l* [=>]28.6 | \[ \color{blue}{\frac{100}{\frac{\frac{i}{n}}{{\left(1 + \frac{i}{n}\right)}^{n} - 1}}}
\] |
*-un-lft-identity [=>]28.6 | \[ \frac{100}{\frac{\frac{i}{n}}{\color{blue}{1 \cdot \left({\left(1 + \frac{i}{n}\right)}^{n} - 1\right)}}}
\] |
associate-/r* [=>]28.6 | \[ \frac{100}{\color{blue}{\frac{\frac{\frac{i}{n}}{1}}{{\left(1 + \frac{i}{n}\right)}^{n} - 1}}}
\] |
associate-/r* [<=]28.6 | \[ \frac{100}{\color{blue}{\frac{\frac{i}{n}}{1 \cdot \left({\left(1 + \frac{i}{n}\right)}^{n} - 1\right)}}}
\] |
*-un-lft-identity [<=]28.6 | \[ \frac{100}{\frac{\frac{i}{n}}{\color{blue}{{\left(1 + \frac{i}{n}\right)}^{n} - 1}}}
\] |
pow-to-exp [=>]27.7 | \[ \frac{100}{\frac{\frac{i}{n}}{\color{blue}{e^{\log \left(1 + \frac{i}{n}\right) \cdot n}} - 1}}
\] |
expm1-def [=>]39.0 | \[ \frac{100}{\frac{\frac{i}{n}}{\color{blue}{\mathsf{expm1}\left(\log \left(1 + \frac{i}{n}\right) \cdot n\right)}}}
\] |
*-commutative [=>]39.0 | \[ \frac{100}{\frac{\frac{i}{n}}{\mathsf{expm1}\left(\color{blue}{n \cdot \log \left(1 + \frac{i}{n}\right)}\right)}}
\] |
log1p-def [=>]98.1 | \[ \frac{100}{\frac{\frac{i}{n}}{\mathsf{expm1}\left(n \cdot \color{blue}{\mathsf{log1p}\left(\frac{i}{n}\right)}\right)}}
\] |
if 5.0000000000000002e-173 < (/.f64 (-.f64 (pow.f64 (+.f64 1 (/.f64 i n)) n) 1) (/.f64 i n)) < 2e10Initial program 97.0%
Simplified97.0%
[Start]97.0 | \[ 100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}
\] |
|---|---|
associate-/r/ [=>]97.0 | \[ 100 \cdot \color{blue}{\left(\frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{i} \cdot n\right)}
\] |
*-commutative [=>]97.0 | \[ 100 \cdot \color{blue}{\left(n \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{i}\right)}
\] |
*-rgt-identity [<=]97.0 | \[ 100 \cdot \left(\color{blue}{\left(n \cdot 1\right)} \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{i}\right)
\] |
associate-*l* [=>]97.0 | \[ 100 \cdot \color{blue}{\left(n \cdot \left(1 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{i}\right)\right)}
\] |
*-lft-identity [=>]97.0 | \[ 100 \cdot \left(n \cdot \color{blue}{\frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{i}}\right)
\] |
sub-neg [=>]97.0 | \[ 100 \cdot \left(n \cdot \frac{\color{blue}{{\left(1 + \frac{i}{n}\right)}^{n} + \left(-1\right)}}{i}\right)
\] |
metadata-eval [=>]97.0 | \[ 100 \cdot \left(n \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} + \color{blue}{-1}}{i}\right)
\] |
Applied egg-rr97.2%
[Start]97.0 | \[ 100 \cdot \left(n \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} + -1}{i}\right)
\] |
|---|---|
*-commutative [=>]97.0 | \[ \color{blue}{\left(n \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} + -1}{i}\right) \cdot 100}
\] |
associate-*r/ [=>]97.1 | \[ \color{blue}{\frac{n \cdot \left({\left(1 + \frac{i}{n}\right)}^{n} + -1\right)}{i}} \cdot 100
\] |
associate-*l/ [=>]97.2 | \[ \color{blue}{\frac{\left(n \cdot \left({\left(1 + \frac{i}{n}\right)}^{n} + -1\right)\right) \cdot 100}{i}}
\] |
associate-/l* [=>]97.2 | \[ \color{blue}{\frac{n \cdot \left({\left(1 + \frac{i}{n}\right)}^{n} + -1\right)}{\frac{i}{100}}}
\] |
if 2e10 < (/.f64 (-.f64 (pow.f64 (+.f64 1 (/.f64 i n)) n) 1) (/.f64 i n)) Initial program 0.6%
Taylor expanded in n around inf 2.2%
Simplified78.3%
[Start]2.2 | \[ 100 \cdot \frac{n \cdot \left(e^{i} - 1\right)}{i}
\] |
|---|---|
*-commutative [=>]2.2 | \[ \color{blue}{\frac{n \cdot \left(e^{i} - 1\right)}{i} \cdot 100}
\] |
associate-/l* [=>]2.2 | \[ \color{blue}{\frac{n}{\frac{i}{e^{i} - 1}}} \cdot 100
\] |
expm1-def [=>]78.3 | \[ \frac{n}{\frac{i}{\color{blue}{\mathsf{expm1}\left(i\right)}}} \cdot 100
\] |
Taylor expanded in i around 0 99.3%
Simplified99.3%
[Start]99.3 | \[ \frac{n}{1 + \left(0.08333333333333333 \cdot {i}^{2} + -0.5 \cdot i\right)} \cdot 100
\] |
|---|---|
+-commutative [=>]99.3 | \[ \frac{n}{1 + \color{blue}{\left(-0.5 \cdot i + 0.08333333333333333 \cdot {i}^{2}\right)}} \cdot 100
\] |
*-commutative [=>]99.3 | \[ \frac{n}{1 + \left(\color{blue}{i \cdot -0.5} + 0.08333333333333333 \cdot {i}^{2}\right)} \cdot 100
\] |
*-commutative [=>]99.3 | \[ \frac{n}{1 + \left(i \cdot -0.5 + \color{blue}{{i}^{2} \cdot 0.08333333333333333}\right)} \cdot 100
\] |
unpow2 [=>]99.3 | \[ \frac{n}{1 + \left(i \cdot -0.5 + \color{blue}{\left(i \cdot i\right)} \cdot 0.08333333333333333\right)} \cdot 100
\] |
associate-*l* [=>]99.3 | \[ \frac{n}{1 + \left(i \cdot -0.5 + \color{blue}{i \cdot \left(i \cdot 0.08333333333333333\right)}\right)} \cdot 100
\] |
distribute-lft-out [=>]99.3 | \[ \frac{n}{1 + \color{blue}{i \cdot \left(-0.5 + i \cdot 0.08333333333333333\right)}} \cdot 100
\] |
Final simplification98.3%
| Alternative 1 | |
|---|---|
| Accuracy | 82.7% |
| Cost | 22408 |
| Alternative 2 | |
|---|---|
| Accuracy | 82.7% |
| Cost | 21832 |
| Alternative 3 | |
|---|---|
| Accuracy | 82.7% |
| Cost | 21768 |
| Alternative 4 | |
|---|---|
| Accuracy | 82.7% |
| Cost | 21768 |
| Alternative 5 | |
|---|---|
| Accuracy | 97.7% |
| Cost | 21768 |
| Alternative 6 | |
|---|---|
| Accuracy | 82.2% |
| Cost | 7244 |
| Alternative 7 | |
|---|---|
| Accuracy | 82.3% |
| Cost | 7244 |
| Alternative 8 | |
|---|---|
| Accuracy | 82.2% |
| Cost | 7244 |
| Alternative 9 | |
|---|---|
| Accuracy | 82.2% |
| Cost | 7244 |
| Alternative 10 | |
|---|---|
| Accuracy | 69.8% |
| Cost | 968 |
| Alternative 11 | |
|---|---|
| Accuracy | 69.8% |
| Cost | 968 |
| Alternative 12 | |
|---|---|
| Accuracy | 69.8% |
| Cost | 841 |
| Alternative 13 | |
|---|---|
| Accuracy | 69.8% |
| Cost | 841 |
| Alternative 14 | |
|---|---|
| Accuracy | 64.0% |
| Cost | 713 |
| Alternative 15 | |
|---|---|
| Accuracy | 65.0% |
| Cost | 713 |
| Alternative 16 | |
|---|---|
| Accuracy | 68.9% |
| Cost | 713 |
| Alternative 17 | |
|---|---|
| Accuracy | 57.3% |
| Cost | 584 |
| Alternative 18 | |
|---|---|
| Accuracy | 3.0% |
| Cost | 192 |
| Alternative 19 | |
|---|---|
| Accuracy | 55.9% |
| Cost | 192 |
herbie shell --seed 2023135
(FPCore (i n)
:name "Compound Interest"
:precision binary64
:herbie-target
(* 100.0 (/ (- (exp (* n (if (== (+ 1.0 (/ i n)) 1.0) (/ i n) (/ (* (/ i n) (log (+ 1.0 (/ i n)))) (- (+ (/ i n) 1.0) 1.0))))) 1.0) (/ i n)))
(* 100.0 (/ (- (pow (+ 1.0 (/ i n)) n) 1.0) (/ i n))))