
(FPCore (i n) :precision binary64 (* 100.0 (/ (- (pow (+ 1.0 (/ i n)) n) 1.0) (/ i n))))
double code(double i, double n) {
return 100.0 * ((pow((1.0 + (i / n)), n) - 1.0) / (i / n));
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
code = 100.0d0 * ((((1.0d0 + (i / n)) ** n) - 1.0d0) / (i / n))
end function
public static double code(double i, double n) {
return 100.0 * ((Math.pow((1.0 + (i / n)), n) - 1.0) / (i / n));
}
def code(i, n): return 100.0 * ((math.pow((1.0 + (i / n)), n) - 1.0) / (i / n))
function code(i, n) return Float64(100.0 * Float64(Float64((Float64(1.0 + Float64(i / n)) ^ n) - 1.0) / Float64(i / n))) end
function tmp = code(i, n) tmp = 100.0 * ((((1.0 + (i / n)) ^ n) - 1.0) / (i / n)); 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]
\begin{array}{l}
\\
100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 6 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (i n) :precision binary64 (* 100.0 (/ (- (pow (+ 1.0 (/ i n)) n) 1.0) (/ i n))))
double code(double i, double n) {
return 100.0 * ((pow((1.0 + (i / n)), n) - 1.0) / (i / n));
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
code = 100.0d0 * ((((1.0d0 + (i / n)) ** n) - 1.0d0) / (i / n))
end function
public static double code(double i, double n) {
return 100.0 * ((Math.pow((1.0 + (i / n)), n) - 1.0) / (i / n));
}
def code(i, n): return 100.0 * ((math.pow((1.0 + (i / n)), n) - 1.0) / (i / n))
function code(i, n) return Float64(100.0 * Float64(Float64((Float64(1.0 + Float64(i / n)) ^ n) - 1.0) / Float64(i / n))) end
function tmp = code(i, n) tmp = 100.0 * ((((1.0 + (i / n)) ^ n) - 1.0) / (i / n)); 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]
\begin{array}{l}
\\
100 \cdot \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}
\end{array}
(FPCore (i n) :precision binary64 (* 100.0 (/ n (/ i (expm1 i)))))
double code(double i, double n) {
return 100.0 * (n / (i / expm1(i)));
}
public static double code(double i, double n) {
return 100.0 * (n / (i / Math.expm1(i)));
}
def code(i, n): return 100.0 * (n / (i / math.expm1(i)))
function code(i, n) return Float64(100.0 * Float64(n / Float64(i / expm1(i)))) end
code[i_, n_] := N[(100.0 * N[(n / N[(i / N[(Exp[i] - 1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
100 \cdot \frac{n}{\frac{i}{\mathsf{expm1}\left(i\right)}}
\end{array}
Initial program 26.4%
associate-/r/26.6%
sub-neg26.6%
metadata-eval26.6%
Simplified26.6%
Taylor expanded in n around inf 33.2%
associate-/l*33.2%
expm1-def75.0%
Simplified75.0%
Final simplification75.0%
(FPCore (i n) :precision binary64 (* 100.0 (* n (+ 1.0 (* i (- 0.5 (/ 0.5 n)))))))
double code(double i, double n) {
return 100.0 * (n * (1.0 + (i * (0.5 - (0.5 / n)))));
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
code = 100.0d0 * (n * (1.0d0 + (i * (0.5d0 - (0.5d0 / n)))))
end function
public static double code(double i, double n) {
return 100.0 * (n * (1.0 + (i * (0.5 - (0.5 / n)))));
}
def code(i, n): return 100.0 * (n * (1.0 + (i * (0.5 - (0.5 / n)))))
function code(i, n) return Float64(100.0 * Float64(n * Float64(1.0 + Float64(i * Float64(0.5 - Float64(0.5 / n)))))) end
function tmp = code(i, n) tmp = 100.0 * (n * (1.0 + (i * (0.5 - (0.5 / n))))); end
code[i_, n_] := N[(100.0 * N[(n * N[(1.0 + N[(i * N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
100 \cdot \left(n \cdot \left(1 + i \cdot \left(0.5 - \frac{0.5}{n}\right)\right)\right)
\end{array}
Initial program 26.4%
associate-/r/26.6%
sub-neg26.6%
metadata-eval26.6%
Simplified26.6%
Taylor expanded in i around 0 53.4%
associate-*r/53.4%
metadata-eval53.4%
Simplified53.4%
Final simplification53.4%
(FPCore (i n) :precision binary64 (* 100.0 (+ n (* (* n i) (- 0.5 (/ 0.5 n))))))
double code(double i, double n) {
return 100.0 * (n + ((n * i) * (0.5 - (0.5 / n))));
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
code = 100.0d0 * (n + ((n * i) * (0.5d0 - (0.5d0 / n))))
end function
public static double code(double i, double n) {
return 100.0 * (n + ((n * i) * (0.5 - (0.5 / n))));
}
def code(i, n): return 100.0 * (n + ((n * i) * (0.5 - (0.5 / n))))
function code(i, n) return Float64(100.0 * Float64(n + Float64(Float64(n * i) * Float64(0.5 - Float64(0.5 / n))))) end
function tmp = code(i, n) tmp = 100.0 * (n + ((n * i) * (0.5 - (0.5 / n)))); end
code[i_, n_] := N[(100.0 * N[(n + N[(N[(n * i), $MachinePrecision] * N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
100 \cdot \left(n + \left(n \cdot i\right) \cdot \left(0.5 - \frac{0.5}{n}\right)\right)
\end{array}
Initial program 26.4%
associate-/r/26.6%
associate-*r*26.6%
*-commutative26.6%
associate-*r/26.6%
sub-neg26.6%
distribute-lft-in26.6%
metadata-eval26.6%
metadata-eval26.6%
metadata-eval26.6%
fma-def26.6%
metadata-eval26.6%
Simplified26.6%
Taylor expanded in i around 0 53.4%
distribute-lft-out53.4%
associate-*r*53.5%
associate-*r/53.5%
metadata-eval53.5%
Simplified53.5%
Final simplification53.5%
(FPCore (i n) :precision binary64 (* n (+ 100.0 (* i 50.0))))
double code(double i, double n) {
return n * (100.0 + (i * 50.0));
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
code = n * (100.0d0 + (i * 50.0d0))
end function
public static double code(double i, double n) {
return n * (100.0 + (i * 50.0));
}
def code(i, n): return n * (100.0 + (i * 50.0))
function code(i, n) return Float64(n * Float64(100.0 + Float64(i * 50.0))) end
function tmp = code(i, n) tmp = n * (100.0 + (i * 50.0)); end
code[i_, n_] := N[(n * N[(100.0 + N[(i * 50.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
n \cdot \left(100 + i \cdot 50\right)
\end{array}
Initial program 26.4%
associate-/r/26.6%
sub-neg26.6%
metadata-eval26.6%
Simplified26.6%
Taylor expanded in n around inf 33.2%
associate-/l*33.2%
expm1-def75.0%
Simplified75.0%
associate-*r/75.0%
*-commutative75.0%
clear-num75.1%
*-commutative75.1%
Applied egg-rr75.1%
Taylor expanded in i around 0 53.7%
+-commutative53.7%
associate-*r*53.7%
distribute-rgt-out53.7%
Simplified53.7%
Final simplification53.7%
(FPCore (i n) :precision binary64 (* i -50.0))
double code(double i, double n) {
return i * -50.0;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
code = i * (-50.0d0)
end function
public static double code(double i, double n) {
return i * -50.0;
}
def code(i, n): return i * -50.0
function code(i, n) return Float64(i * -50.0) end
function tmp = code(i, n) tmp = i * -50.0; end
code[i_, n_] := N[(i * -50.0), $MachinePrecision]
\begin{array}{l}
\\
i \cdot -50
\end{array}
Initial program 26.4%
associate-/r/26.6%
sub-neg26.6%
metadata-eval26.6%
Simplified26.6%
Taylor expanded in i around 0 53.4%
associate-*r/53.4%
metadata-eval53.4%
Simplified53.4%
Taylor expanded in n around 0 2.6%
*-commutative2.6%
Simplified2.6%
Final simplification2.6%
(FPCore (i n) :precision binary64 (* 100.0 n))
double code(double i, double n) {
return 100.0 * n;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
code = 100.0d0 * n
end function
public static double code(double i, double n) {
return 100.0 * n;
}
def code(i, n): return 100.0 * n
function code(i, n) return Float64(100.0 * n) end
function tmp = code(i, n) tmp = 100.0 * n; end
code[i_, n_] := N[(100.0 * n), $MachinePrecision]
\begin{array}{l}
\\
100 \cdot n
\end{array}
Initial program 26.4%
associate-/r/26.6%
associate-*r*26.6%
*-commutative26.6%
associate-*r/26.6%
sub-neg26.6%
distribute-lft-in26.6%
metadata-eval26.6%
metadata-eval26.6%
metadata-eval26.6%
fma-def26.6%
metadata-eval26.6%
Simplified26.6%
Taylor expanded in i around 0 48.8%
Final simplification48.8%
(FPCore (i n)
:precision binary64
(let* ((t_0 (+ 1.0 (/ i n))))
(*
100.0
(/
(-
(exp
(*
n
(if (== t_0 1.0)
(/ i n)
(/ (* (/ i n) (log t_0)) (- (+ (/ i n) 1.0) 1.0)))))
1.0)
(/ i n)))))
double code(double i, double n) {
double t_0 = 1.0 + (i / n);
double tmp;
if (t_0 == 1.0) {
tmp = i / n;
} else {
tmp = ((i / n) * log(t_0)) / (((i / n) + 1.0) - 1.0);
}
return 100.0 * ((exp((n * tmp)) - 1.0) / (i / n));
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: t_0
real(8) :: tmp
t_0 = 1.0d0 + (i / n)
if (t_0 == 1.0d0) then
tmp = i / n
else
tmp = ((i / n) * log(t_0)) / (((i / n) + 1.0d0) - 1.0d0)
end if
code = 100.0d0 * ((exp((n * tmp)) - 1.0d0) / (i / n))
end function
public static double code(double i, double n) {
double t_0 = 1.0 + (i / n);
double tmp;
if (t_0 == 1.0) {
tmp = i / n;
} else {
tmp = ((i / n) * Math.log(t_0)) / (((i / n) + 1.0) - 1.0);
}
return 100.0 * ((Math.exp((n * tmp)) - 1.0) / (i / n));
}
def code(i, n): t_0 = 1.0 + (i / n) tmp = 0 if t_0 == 1.0: tmp = i / n else: tmp = ((i / n) * math.log(t_0)) / (((i / n) + 1.0) - 1.0) return 100.0 * ((math.exp((n * tmp)) - 1.0) / (i / n))
function code(i, n) t_0 = Float64(1.0 + Float64(i / n)) tmp = 0.0 if (t_0 == 1.0) tmp = Float64(i / n); else tmp = Float64(Float64(Float64(i / n) * log(t_0)) / Float64(Float64(Float64(i / n) + 1.0) - 1.0)); end return Float64(100.0 * Float64(Float64(exp(Float64(n * tmp)) - 1.0) / Float64(i / n))) end
function tmp_2 = code(i, n) t_0 = 1.0 + (i / n); tmp = 0.0; if (t_0 == 1.0) tmp = i / n; else tmp = ((i / n) * log(t_0)) / (((i / n) + 1.0) - 1.0); end tmp_2 = 100.0 * ((exp((n * tmp)) - 1.0) / (i / n)); end
code[i_, n_] := Block[{t$95$0 = N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision]}, N[(100.0 * N[(N[(N[Exp[N[(n * If[Equal[t$95$0, 1.0], N[(i / n), $MachinePrecision], N[(N[(N[(i / n), $MachinePrecision] * N[Log[t$95$0], $MachinePrecision]), $MachinePrecision] / N[(N[(N[(i / n), $MachinePrecision] + 1.0), $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]], $MachinePrecision] - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 1 + \frac{i}{n}\\
100 \cdot \frac{e^{n \cdot \begin{array}{l}
\mathbf{if}\;t\_0 = 1:\\
\;\;\;\;\frac{i}{n}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{i}{n} \cdot \log t\_0}{\left(\frac{i}{n} + 1\right) - 1}\\
\end{array}} - 1}{\frac{i}{n}}
\end{array}
\end{array}
herbie shell --seed 2024033
(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))))