| Alternative 1 | |
|---|---|
| Accuracy | 94.0% |
| Cost | 29004 |

(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)) (t_1 (/ (+ t_0 -1.0) (/ i n))))
(if (<= t_1 -2000.0)
(* n (/ (+ (* t_0 100.0) -100.0) i))
(if (<= t_1 0.0)
(* 100.0 (/ (expm1 (* n (log1p (/ i n)))) (/ i n)))
(if (<= t_1 INFINITY)
(/ n (/ i (fma 100.0 t_0 -100.0)))
(* n 100.0))))))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);
double t_1 = (t_0 + -1.0) / (i / n);
double tmp;
if (t_1 <= -2000.0) {
tmp = n * (((t_0 * 100.0) + -100.0) / i);
} else if (t_1 <= 0.0) {
tmp = 100.0 * (expm1((n * log1p((i / n)))) / (i / n));
} else if (t_1 <= ((double) INFINITY)) {
tmp = n / (i / fma(100.0, t_0, -100.0));
} else {
tmp = n * 100.0;
}
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(1.0 + Float64(i / n)) ^ n t_1 = Float64(Float64(t_0 + -1.0) / Float64(i / n)) tmp = 0.0 if (t_1 <= -2000.0) tmp = Float64(n * Float64(Float64(Float64(t_0 * 100.0) + -100.0) / i)); elseif (t_1 <= 0.0) tmp = Float64(100.0 * Float64(expm1(Float64(n * log1p(Float64(i / n)))) / Float64(i / n))); elseif (t_1 <= Inf) tmp = Float64(n / Float64(i / fma(100.0, t_0, -100.0))); else tmp = Float64(n * 100.0); 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[Power[N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision], n], $MachinePrecision]}, Block[{t$95$1 = N[(N[(t$95$0 + -1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -2000.0], N[(n * N[(N[(N[(t$95$0 * 100.0), $MachinePrecision] + -100.0), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(100.0 * N[(N[(Exp[N[(n * N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]] - 1), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(n / N[(i / N[(100.0 * t$95$0 + -100.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(n * 100.0), $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}\\
t_1 := \frac{t_0 + -1}{\frac{i}{n}}\\
\mathbf{if}\;t_1 \leq -2000:\\
\;\;\;\;n \cdot \frac{t_0 \cdot 100 + -100}{i}\\
\mathbf{elif}\;t_1 \leq 0:\\
\;\;\;\;100 \cdot \frac{\mathsf{expm1}\left(n \cdot \mathsf{log1p}\left(\frac{i}{n}\right)\right)}{\frac{i}{n}}\\
\mathbf{elif}\;t_1 \leq \infty:\\
\;\;\;\;\frac{n}{\frac{i}{\mathsf{fma}\left(100, t_0, -100\right)}}\\
\mathbf{else}:\\
\;\;\;\;n \cdot 100\\
\end{array}
Herbie found 19 alternatives:
| Alternative | Accuracy | Speedup |
|---|
| Original | 28.6% |
|---|---|
| Target | 34.2% |
| Herbie | 95.8% |
if (/.f64 (-.f64 (pow.f64 (+.f64 1 (/.f64 i n)) n) 1) (/.f64 i n)) < -2e3Initial program 99.7%
Simplified100.0%
[Start]99.7 | \[ 100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}
\] |
|---|---|
associate-/r/ [=>]99.8 | \[ 100 \cdot \color{blue}{\left(\frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{i} \cdot n\right)}
\] |
associate-*r* [=>]100.0 | \[ \color{blue}{\left(100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{i}\right) \cdot n}
\] |
*-commutative [=>]100.0 | \[ \color{blue}{n \cdot \left(100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{i}\right)}
\] |
associate-*r/ [=>]100.0 | \[ n \cdot \color{blue}{\frac{100 \cdot \left({\left(1 + \frac{i}{n}\right)}^{n} - 1\right)}{i}}
\] |
sub-neg [=>]100.0 | \[ n \cdot \frac{100 \cdot \color{blue}{\left({\left(1 + \frac{i}{n}\right)}^{n} + \left(-1\right)\right)}}{i}
\] |
distribute-lft-in [=>]100.0 | \[ n \cdot \frac{\color{blue}{100 \cdot {\left(1 + \frac{i}{n}\right)}^{n} + 100 \cdot \left(-1\right)}}{i}
\] |
fma-def [=>]100.0 | \[ n \cdot \frac{\color{blue}{\mathsf{fma}\left(100, {\left(1 + \frac{i}{n}\right)}^{n}, 100 \cdot \left(-1\right)\right)}}{i}
\] |
metadata-eval [=>]100.0 | \[ n \cdot \frac{\mathsf{fma}\left(100, {\left(1 + \frac{i}{n}\right)}^{n}, 100 \cdot \color{blue}{-1}\right)}{i}
\] |
metadata-eval [=>]100.0 | \[ n \cdot \frac{\mathsf{fma}\left(100, {\left(1 + \frac{i}{n}\right)}^{n}, \color{blue}{-100}\right)}{i}
\] |
Applied egg-rr100.0%
[Start]100.0 | \[ n \cdot \frac{\mathsf{fma}\left(100, {\left(1 + \frac{i}{n}\right)}^{n}, -100\right)}{i}
\] |
|---|---|
fma-udef [=>]100.0 | \[ n \cdot \frac{\color{blue}{100 \cdot {\left(1 + \frac{i}{n}\right)}^{n} + -100}}{i}
\] |
if -2e3 < (/.f64 (-.f64 (pow.f64 (+.f64 1 (/.f64 i n)) n) 1) (/.f64 i n)) < 0.0Initial program 25.2%
Applied egg-rr99.7%
[Start]25.2 | \[ 100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}
\] |
|---|---|
add-log-exp [=>]25.2 | \[ 100 \cdot \frac{\color{blue}{\log \left(e^{{\left(1 + \frac{i}{n}\right)}^{n} - 1}\right)}}{\frac{i}{n}}
\] |
*-un-lft-identity [=>]25.2 | \[ 100 \cdot \frac{\log \color{blue}{\left(1 \cdot e^{{\left(1 + \frac{i}{n}\right)}^{n} - 1}\right)}}{\frac{i}{n}}
\] |
log-prod [=>]25.2 | \[ 100 \cdot \frac{\color{blue}{\log 1 + \log \left(e^{{\left(1 + \frac{i}{n}\right)}^{n} - 1}\right)}}{\frac{i}{n}}
\] |
metadata-eval [=>]25.2 | \[ 100 \cdot \frac{\color{blue}{0} + \log \left(e^{{\left(1 + \frac{i}{n}\right)}^{n} - 1}\right)}{\frac{i}{n}}
\] |
add-log-exp [<=]25.2 | \[ 100 \cdot \frac{0 + \color{blue}{\left({\left(1 + \frac{i}{n}\right)}^{n} - 1\right)}}{\frac{i}{n}}
\] |
pow-to-exp [=>]25.2 | \[ 100 \cdot \frac{0 + \left(\color{blue}{e^{\log \left(1 + \frac{i}{n}\right) \cdot n}} - 1\right)}{\frac{i}{n}}
\] |
expm1-def [=>]36.7 | \[ 100 \cdot \frac{0 + \color{blue}{\mathsf{expm1}\left(\log \left(1 + \frac{i}{n}\right) \cdot n\right)}}{\frac{i}{n}}
\] |
*-commutative [=>]36.7 | \[ 100 \cdot \frac{0 + \mathsf{expm1}\left(\color{blue}{n \cdot \log \left(1 + \frac{i}{n}\right)}\right)}{\frac{i}{n}}
\] |
log1p-udef [<=]99.7 | \[ 100 \cdot \frac{0 + \mathsf{expm1}\left(n \cdot \color{blue}{\mathsf{log1p}\left(\frac{i}{n}\right)}\right)}{\frac{i}{n}}
\] |
Simplified99.7%
[Start]99.7 | \[ 100 \cdot \frac{0 + \mathsf{expm1}\left(n \cdot \mathsf{log1p}\left(\frac{i}{n}\right)\right)}{\frac{i}{n}}
\] |
|---|---|
+-lft-identity [=>]99.7 | \[ 100 \cdot \frac{\color{blue}{\mathsf{expm1}\left(n \cdot \mathsf{log1p}\left(\frac{i}{n}\right)\right)}}{\frac{i}{n}}
\] |
if 0.0 < (/.f64 (-.f64 (pow.f64 (+.f64 1 (/.f64 i n)) n) 1) (/.f64 i n)) < +inf.0Initial program 99.6%
Simplified99.6%
[Start]99.6 | \[ 100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}
\] |
|---|---|
associate-/r/ [=>]99.4 | \[ 100 \cdot \color{blue}{\left(\frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{i} \cdot n\right)}
\] |
associate-*r* [=>]99.6 | \[ \color{blue}{\left(100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{i}\right) \cdot n}
\] |
*-commutative [=>]99.6 | \[ \color{blue}{n \cdot \left(100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{i}\right)}
\] |
associate-*r/ [=>]99.6 | \[ n \cdot \color{blue}{\frac{100 \cdot \left({\left(1 + \frac{i}{n}\right)}^{n} - 1\right)}{i}}
\] |
sub-neg [=>]99.6 | \[ n \cdot \frac{100 \cdot \color{blue}{\left({\left(1 + \frac{i}{n}\right)}^{n} + \left(-1\right)\right)}}{i}
\] |
distribute-lft-in [=>]99.6 | \[ n \cdot \frac{\color{blue}{100 \cdot {\left(1 + \frac{i}{n}\right)}^{n} + 100 \cdot \left(-1\right)}}{i}
\] |
fma-def [=>]99.6 | \[ n \cdot \frac{\color{blue}{\mathsf{fma}\left(100, {\left(1 + \frac{i}{n}\right)}^{n}, 100 \cdot \left(-1\right)\right)}}{i}
\] |
metadata-eval [=>]99.6 | \[ n \cdot \frac{\mathsf{fma}\left(100, {\left(1 + \frac{i}{n}\right)}^{n}, 100 \cdot \color{blue}{-1}\right)}{i}
\] |
metadata-eval [=>]99.6 | \[ n \cdot \frac{\mathsf{fma}\left(100, {\left(1 + \frac{i}{n}\right)}^{n}, \color{blue}{-100}\right)}{i}
\] |
Applied egg-rr99.7%
[Start]99.6 | \[ n \cdot \frac{\mathsf{fma}\left(100, {\left(1 + \frac{i}{n}\right)}^{n}, -100\right)}{i}
\] |
|---|---|
clear-num [=>]99.6 | \[ n \cdot \color{blue}{\frac{1}{\frac{i}{\mathsf{fma}\left(100, {\left(1 + \frac{i}{n}\right)}^{n}, -100\right)}}}
\] |
un-div-inv [=>]99.7 | \[ \color{blue}{\frac{n}{\frac{i}{\mathsf{fma}\left(100, {\left(1 + \frac{i}{n}\right)}^{n}, -100\right)}}}
\] |
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 1 (/.f64 i n)) n) 1) (/.f64 i n)) Initial program 0.0%
Taylor expanded in i around 0 71.3%
Simplified71.3%
[Start]71.3 | \[ 100 \cdot n
\] |
|---|---|
*-commutative [=>]71.3 | \[ \color{blue}{n \cdot 100}
\] |
Final simplification95.4%
| Alternative 1 | |
|---|---|
| Accuracy | 94.0% |
| Cost | 29004 |
| Alternative 2 | |
|---|---|
| Accuracy | 94.9% |
| Cost | 29004 |
| Alternative 3 | |
|---|---|
| Accuracy | 95.8% |
| Cost | 29004 |
| Alternative 4 | |
|---|---|
| Accuracy | 80.8% |
| Cost | 7632 |
| Alternative 5 | |
|---|---|
| Accuracy | 80.8% |
| Cost | 7632 |
| Alternative 6 | |
|---|---|
| Accuracy | 80.8% |
| Cost | 7632 |
| Alternative 7 | |
|---|---|
| Accuracy | 73.5% |
| Cost | 7376 |
| Alternative 8 | |
|---|---|
| Accuracy | 80.1% |
| Cost | 7113 |
| Alternative 9 | |
|---|---|
| Accuracy | 64.4% |
| Cost | 1992 |
| Alternative 10 | |
|---|---|
| Accuracy | 64.3% |
| Cost | 1097 |
| Alternative 11 | |
|---|---|
| Accuracy | 61.0% |
| Cost | 976 |
| Alternative 12 | |
|---|---|
| Accuracy | 61.4% |
| Cost | 968 |
| Alternative 13 | |
|---|---|
| Accuracy | 56.1% |
| Cost | 713 |
| Alternative 14 | |
|---|---|
| Accuracy | 61.9% |
| Cost | 713 |
| Alternative 15 | |
|---|---|
| Accuracy | 61.9% |
| Cost | 712 |
| Alternative 16 | |
|---|---|
| Accuracy | 53.9% |
| Cost | 452 |
| Alternative 17 | |
|---|---|
| Accuracy | 2.8% |
| Cost | 192 |
| Alternative 18 | |
|---|---|
| Accuracy | 48.6% |
| Cost | 192 |
herbie shell --seed 2023161
(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))))