| Alternative 1 | |
|---|---|
| Error | 11.7 |
| Cost | 7628 |
(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 (- (log i) (log n))) (t_1 (/ (* n 100.0) (/ i (expm1 i)))))
(if (<= n -5.4e-140)
t_1
(if (<= n -5e-310)
(/
(* 100.0 (expm1 (* n (- (log (/ -1.0 n)) (log (/ -1.0 i))))))
(/ i n))
(if (<= n 1.2e-92)
(/
(* 100.0 (+ (* n t_0) (* n (* n (fma 0.5 (pow t_0 2.0) (/ 1.0 i))))))
(/ i n))
t_1)))))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 = log(i) - log(n);
double t_1 = (n * 100.0) / (i / expm1(i));
double tmp;
if (n <= -5.4e-140) {
tmp = t_1;
} else if (n <= -5e-310) {
tmp = (100.0 * expm1((n * (log((-1.0 / n)) - log((-1.0 / i)))))) / (i / n);
} else if (n <= 1.2e-92) {
tmp = (100.0 * ((n * t_0) + (n * (n * fma(0.5, pow(t_0, 2.0), (1.0 / i)))))) / (i / n);
} else {
tmp = t_1;
}
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(log(i) - log(n)) t_1 = Float64(Float64(n * 100.0) / Float64(i / expm1(i))) tmp = 0.0 if (n <= -5.4e-140) tmp = t_1; elseif (n <= -5e-310) tmp = Float64(Float64(100.0 * expm1(Float64(n * Float64(log(Float64(-1.0 / n)) - log(Float64(-1.0 / i)))))) / Float64(i / n)); elseif (n <= 1.2e-92) tmp = Float64(Float64(100.0 * Float64(Float64(n * t_0) + Float64(n * Float64(n * fma(0.5, (t_0 ^ 2.0), Float64(1.0 / i)))))) / Float64(i / n)); else tmp = t_1; 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[Log[i], $MachinePrecision] - N[Log[n], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(n * 100.0), $MachinePrecision] / N[(i / N[(Exp[i] - 1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -5.4e-140], t$95$1, If[LessEqual[n, -5e-310], N[(N[(100.0 * N[(Exp[N[(n * N[(N[Log[N[(-1.0 / n), $MachinePrecision]], $MachinePrecision] - N[Log[N[(-1.0 / i), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.2e-92], N[(N[(100.0 * N[(N[(n * t$95$0), $MachinePrecision] + N[(n * N[(n * N[(0.5 * N[Power[t$95$0, 2.0], $MachinePrecision] + N[(1.0 / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}
\begin{array}{l}
t_0 := \log i - \log n\\
t_1 := \frac{n \cdot 100}{\frac{i}{\mathsf{expm1}\left(i\right)}}\\
\mathbf{if}\;n \leq -5.4 \cdot 10^{-140}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;n \leq -5 \cdot 10^{-310}:\\
\;\;\;\;\frac{100 \cdot \mathsf{expm1}\left(n \cdot \left(\log \left(\frac{-1}{n}\right) - \log \left(\frac{-1}{i}\right)\right)\right)}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 1.2 \cdot 10^{-92}:\\
\;\;\;\;\frac{100 \cdot \left(n \cdot t_0 + n \cdot \left(n \cdot \mathsf{fma}\left(0.5, {t_0}^{2}, \frac{1}{i}\right)\right)\right)}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
| Original | 47.7 |
|---|---|
| Target | 47.2 |
| Herbie | 12.0 |
if n < -5.4e-140 or 1.2000000000000001e-92 < n Initial program 51.8
Simplified51.6
[Start]51.8 | \[ 100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}
\] |
|---|---|
associate-/r/ [=>]51.6 | \[ 100 \cdot \color{blue}{\left(\frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{i} \cdot n\right)}
\] |
*-commutative [=>]51.6 | \[ 100 \cdot \color{blue}{\left(n \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{i}\right)}
\] |
*-rgt-identity [<=]51.6 | \[ 100 \cdot \left(\color{blue}{\left(n \cdot 1\right)} \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{i}\right)
\] |
associate-*l* [=>]51.6 | \[ 100 \cdot \color{blue}{\left(n \cdot \left(1 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{i}\right)\right)}
\] |
*-lft-identity [=>]51.6 | \[ 100 \cdot \left(n \cdot \color{blue}{\frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{i}}\right)
\] |
sub-neg [=>]51.6 | \[ 100 \cdot \left(n \cdot \frac{\color{blue}{{\left(1 + \frac{i}{n}\right)}^{n} + \left(-1\right)}}{i}\right)
\] |
metadata-eval [=>]51.6 | \[ 100 \cdot \left(n \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} + \color{blue}{-1}}{i}\right)
\] |
Taylor expanded in n around inf 48.5
Simplified9.9
[Start]48.5 | \[ 100 \cdot \left(n \cdot \frac{e^{i} - 1}{i}\right)
\] |
|---|---|
expm1-def [=>]9.9 | \[ 100 \cdot \left(n \cdot \frac{\color{blue}{\mathsf{expm1}\left(i\right)}}{i}\right)
\] |
Applied egg-rr10.0
if -5.4e-140 < n < -4.999999999999985e-310Initial program 23.4
Simplified23.4
[Start]23.4 | \[ 100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}
\] |
|---|---|
associate-*r/ [=>]23.4 | \[ \color{blue}{\frac{100 \cdot \left({\left(1 + \frac{i}{n}\right)}^{n} - 1\right)}{\frac{i}{n}}}
\] |
sub-neg [=>]23.4 | \[ \frac{100 \cdot \color{blue}{\left({\left(1 + \frac{i}{n}\right)}^{n} + \left(-1\right)\right)}}{\frac{i}{n}}
\] |
metadata-eval [=>]23.4 | \[ \frac{100 \cdot \left({\left(1 + \frac{i}{n}\right)}^{n} + \color{blue}{-1}\right)}{\frac{i}{n}}
\] |
Taylor expanded in i around -inf 23.8
Simplified15.9
[Start]23.8 | \[ \frac{100 \cdot \left(e^{n \cdot \left(-1 \cdot \log \left(\frac{-1}{i}\right) + \log \left(-\frac{1}{n}\right)\right)} - 1\right)}{\frac{i}{n}}
\] |
|---|---|
expm1-def [=>]15.9 | \[ \frac{100 \cdot \color{blue}{\mathsf{expm1}\left(n \cdot \left(-1 \cdot \log \left(\frac{-1}{i}\right) + \log \left(-\frac{1}{n}\right)\right)\right)}}{\frac{i}{n}}
\] |
+-commutative [=>]15.9 | \[ \frac{100 \cdot \mathsf{expm1}\left(n \cdot \color{blue}{\left(\log \left(-\frac{1}{n}\right) + -1 \cdot \log \left(\frac{-1}{i}\right)\right)}\right)}{\frac{i}{n}}
\] |
mul-1-neg [=>]15.9 | \[ \frac{100 \cdot \mathsf{expm1}\left(n \cdot \left(\log \left(-\frac{1}{n}\right) + \color{blue}{\left(-\log \left(\frac{-1}{i}\right)\right)}\right)\right)}{\frac{i}{n}}
\] |
unsub-neg [=>]15.9 | \[ \frac{100 \cdot \mathsf{expm1}\left(n \cdot \color{blue}{\left(\log \left(-\frac{1}{n}\right) - \log \left(\frac{-1}{i}\right)\right)}\right)}{\frac{i}{n}}
\] |
distribute-neg-frac [=>]15.9 | \[ \frac{100 \cdot \mathsf{expm1}\left(n \cdot \left(\log \color{blue}{\left(\frac{-1}{n}\right)} - \log \left(\frac{-1}{i}\right)\right)\right)}{\frac{i}{n}}
\] |
metadata-eval [=>]15.9 | \[ \frac{100 \cdot \mathsf{expm1}\left(n \cdot \left(\log \left(\frac{\color{blue}{-1}}{n}\right) - \log \left(\frac{-1}{i}\right)\right)\right)}{\frac{i}{n}}
\] |
if -4.999999999999985e-310 < n < 1.2000000000000001e-92Initial program 46.3
Simplified46.3
[Start]46.3 | \[ 100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}
\] |
|---|---|
associate-*r/ [=>]46.3 | \[ \color{blue}{\frac{100 \cdot \left({\left(1 + \frac{i}{n}\right)}^{n} - 1\right)}{\frac{i}{n}}}
\] |
sub-neg [=>]46.3 | \[ \frac{100 \cdot \color{blue}{\left({\left(1 + \frac{i}{n}\right)}^{n} + \left(-1\right)\right)}}{\frac{i}{n}}
\] |
metadata-eval [=>]46.3 | \[ \frac{100 \cdot \left({\left(1 + \frac{i}{n}\right)}^{n} + \color{blue}{-1}\right)}{\frac{i}{n}}
\] |
Taylor expanded in n around 0 19.9
Simplified19.7
[Start]19.9 | \[ \frac{100 \cdot \left({n}^{2} \cdot \left(0.5 \cdot {\left(-1 \cdot \log n + \log i\right)}^{2} + \frac{1}{i}\right)\right) + 100 \cdot \left(n \cdot \left(-1 \cdot \log n + \log i\right)\right)}{\frac{i}{n}}
\] |
|---|---|
distribute-lft-out [=>]19.9 | \[ \frac{\color{blue}{100 \cdot \left({n}^{2} \cdot \left(0.5 \cdot {\left(-1 \cdot \log n + \log i\right)}^{2} + \frac{1}{i}\right) + n \cdot \left(-1 \cdot \log n + \log i\right)\right)}}{\frac{i}{n}}
\] |
+-commutative [=>]19.9 | \[ \frac{100 \cdot \color{blue}{\left(n \cdot \left(-1 \cdot \log n + \log i\right) + {n}^{2} \cdot \left(0.5 \cdot {\left(-1 \cdot \log n + \log i\right)}^{2} + \frac{1}{i}\right)\right)}}{\frac{i}{n}}
\] |
+-commutative [=>]19.9 | \[ \frac{100 \cdot \left(n \cdot \color{blue}{\left(\log i + -1 \cdot \log n\right)} + {n}^{2} \cdot \left(0.5 \cdot {\left(-1 \cdot \log n + \log i\right)}^{2} + \frac{1}{i}\right)\right)}{\frac{i}{n}}
\] |
mul-1-neg [=>]19.9 | \[ \frac{100 \cdot \left(n \cdot \left(\log i + \color{blue}{\left(-\log n\right)}\right) + {n}^{2} \cdot \left(0.5 \cdot {\left(-1 \cdot \log n + \log i\right)}^{2} + \frac{1}{i}\right)\right)}{\frac{i}{n}}
\] |
unsub-neg [=>]19.9 | \[ \frac{100 \cdot \left(n \cdot \color{blue}{\left(\log i - \log n\right)} + {n}^{2} \cdot \left(0.5 \cdot {\left(-1 \cdot \log n + \log i\right)}^{2} + \frac{1}{i}\right)\right)}{\frac{i}{n}}
\] |
unpow2 [=>]19.9 | \[ \frac{100 \cdot \left(n \cdot \left(\log i - \log n\right) + \color{blue}{\left(n \cdot n\right)} \cdot \left(0.5 \cdot {\left(-1 \cdot \log n + \log i\right)}^{2} + \frac{1}{i}\right)\right)}{\frac{i}{n}}
\] |
associate-*l* [=>]19.7 | \[ \frac{100 \cdot \left(n \cdot \left(\log i - \log n\right) + \color{blue}{n \cdot \left(n \cdot \left(0.5 \cdot {\left(-1 \cdot \log n + \log i\right)}^{2} + \frac{1}{i}\right)\right)}\right)}{\frac{i}{n}}
\] |
fma-def [=>]19.7 | \[ \frac{100 \cdot \left(n \cdot \left(\log i - \log n\right) + n \cdot \left(n \cdot \color{blue}{\mathsf{fma}\left(0.5, {\left(-1 \cdot \log n + \log i\right)}^{2}, \frac{1}{i}\right)}\right)\right)}{\frac{i}{n}}
\] |
Final simplification12.0
| Alternative 1 | |
|---|---|
| Error | 11.7 |
| Cost | 7628 |
| Alternative 2 | |
|---|---|
| Error | 12.3 |
| Cost | 7244 |
| Alternative 3 | |
|---|---|
| Error | 12.4 |
| Cost | 7244 |
| Alternative 4 | |
|---|---|
| Error | 12.4 |
| Cost | 7244 |
| Alternative 5 | |
|---|---|
| Error | 19.6 |
| Cost | 969 |
| Alternative 6 | |
|---|---|
| Error | 19.6 |
| Cost | 969 |
| Alternative 7 | |
|---|---|
| Error | 23.0 |
| Cost | 841 |
| Alternative 8 | |
|---|---|
| Error | 23.0 |
| Cost | 841 |
| Alternative 9 | |
|---|---|
| Error | 21.3 |
| Cost | 841 |
| Alternative 10 | |
|---|---|
| Error | 19.6 |
| Cost | 841 |
| Alternative 11 | |
|---|---|
| Error | 22.9 |
| Cost | 840 |
| Alternative 12 | |
|---|---|
| Error | 23.3 |
| Cost | 713 |
| Alternative 13 | |
|---|---|
| Error | 23.0 |
| Cost | 713 |
| Alternative 14 | |
|---|---|
| Error | 62.1 |
| Cost | 192 |
| Alternative 15 | |
|---|---|
| Error | 28.4 |
| Cost | 192 |
herbie shell --seed 2023039
(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))))