100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}\begin{array}{l}
\mathbf{if}\;i \le -7.03094120304917353 \cdot 10^{-9}:\\
\;\;\;\;100 \cdot \frac{{\left({\left(1 + \frac{i}{n}\right)}^{n}\right)}^{3} - {1}^{3}}{\frac{\mathsf{fma}\left(1, {\left(1 + \frac{i}{n}\right)}^{n} + 1, {\left(1 + \frac{i}{n}\right)}^{\left(2 \cdot n\right)}\right) \cdot i}{n}}\\
\mathbf{elif}\;i \le 0.1391033594988507:\\
\;\;\;\;100 \cdot \left(\mathsf{fma}\left(-0.5 \cdot i, \log 1, \mathsf{fma}\left(i, 0.5, \mathsf{fma}\left(\frac{\log 1}{i}, n, 1\right)\right)\right) \cdot n\right)\\
\mathbf{elif}\;i \le 8.5994276340506012 \cdot 10^{168}:\\
\;\;\;\;\frac{100 \cdot \left({\left(1 + \frac{i}{n}\right)}^{n} - 1\right)}{\frac{i}{n}}\\
\mathbf{elif}\;i \le 3.601044061012584 \cdot 10^{286}:\\
\;\;\;\;100 \cdot \frac{\mathsf{fma}\left(1, i, \mathsf{fma}\left(\log 1, n, 1\right)\right) - 1}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;\left(100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{i}\right) \cdot n\\
\end{array}double code(double i, double n) {
return ((double) (100.0 * ((double) (((double) (((double) pow(((double) (1.0 + ((double) (i / n)))), n)) - 1.0)) / ((double) (i / n))))));
}
double code(double i, double n) {
double VAR;
if ((i <= -7.0309412030491735e-09)) {
VAR = ((double) (100.0 * ((double) (((double) (((double) pow(((double) pow(((double) (1.0 + ((double) (i / n)))), n)), 3.0)) - ((double) pow(1.0, 3.0)))) / ((double) (((double) (((double) fma(1.0, ((double) (((double) pow(((double) (1.0 + ((double) (i / n)))), n)) + 1.0)), ((double) pow(((double) (1.0 + ((double) (i / n)))), ((double) (2.0 * n)))))) * i)) / n))))));
} else {
double VAR_1;
if ((i <= 0.13910335949885066)) {
VAR_1 = ((double) (100.0 * ((double) (((double) fma(((double) -(((double) (0.5 * i)))), ((double) log(1.0)), ((double) fma(i, 0.5, ((double) fma(((double) (((double) log(1.0)) / i)), n, 1.0)))))) * n))));
} else {
double VAR_2;
if ((i <= 8.599427634050601e+168)) {
VAR_2 = ((double) (((double) (100.0 * ((double) (((double) pow(((double) (1.0 + ((double) (i / n)))), n)) - 1.0)))) / ((double) (i / n))));
} else {
double VAR_3;
if ((i <= 3.601044061012584e+286)) {
VAR_3 = ((double) (100.0 * ((double) (((double) (((double) fma(1.0, i, ((double) fma(((double) log(1.0)), n, 1.0)))) - 1.0)) / ((double) (i / n))))));
} else {
VAR_3 = ((double) (((double) (100.0 * ((double) (((double) (((double) pow(((double) (1.0 + ((double) (i / n)))), n)) - 1.0)) / i)))) * n));
}
VAR_2 = VAR_3;
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus i




Bits error versus n
Results
| Original | 47.3 |
|---|---|
| Target | 47.3 |
| Herbie | 16.6 |
if i < -7.0309412030491735e-09Initial program 28.0
rmApplied flip3--28.0
Applied associate-/l/28.0
Simplified28.0
if -7.0309412030491735e-09 < i < 0.13910335949885066Initial program 58.2
Taylor expanded around 0 25.5
Simplified25.5
rmApplied associate-/r/8.8
Taylor expanded around 0 8.8
Simplified8.8
if 0.13910335949885066 < i < 8.599427634050601e+168Initial program 30.8
rmApplied associate-*r/30.8
if 8.599427634050601e+168 < i < 3.601044061012584e+286Initial program 32.2
Taylor expanded around 0 35.0
Simplified35.0
if 3.601044061012584e+286 < i Initial program 27.7
rmApplied associate-/r/27.8
Applied associate-*r*27.8
Final simplification16.6
herbie shell --seed 2020121 +o rules:numerics
(FPCore (i n)
:name "Compound Interest"
:precision binary64
:herbie-target
(* 100 (/ (- (exp (* n (if (== (+ 1 (/ i n)) 1) (/ i n) (/ (* (/ i n) (log (+ 1 (/ i n)))) (- (+ (/ i n) 1) 1))))) 1) (/ i n)))
(* 100 (/ (- (pow (+ 1 (/ i n)) n) 1) (/ i n))))