
(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 15 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
(let* ((t_0 (/ (- (pow (+ 1.0 (/ i n)) n) 1.0) (/ i n))))
(if (<= t_0 1e-302)
(/ (* 100.0 (expm1 (* (log1p (/ i n)) n))) (/ i n))
(if (<= t_0 INFINITY)
(- (* (* (/ (pow (+ (/ i n) 1.0) n) i) n) 100.0) (/ (* n 100.0) i))
(pow (fma (* 0.01 i) (- (/ 0.5 (* n n)) (/ 0.5 n)) (/ 0.01 n)) -1.0)))))
double code(double i, double n) {
double t_0 = (pow((1.0 + (i / n)), n) - 1.0) / (i / n);
double tmp;
if (t_0 <= 1e-302) {
tmp = (100.0 * expm1((log1p((i / n)) * n))) / (i / n);
} else if (t_0 <= ((double) INFINITY)) {
tmp = (((pow(((i / n) + 1.0), n) / i) * n) * 100.0) - ((n * 100.0) / i);
} else {
tmp = pow(fma((0.01 * i), ((0.5 / (n * n)) - (0.5 / n)), (0.01 / n)), -1.0);
}
return tmp;
}
function code(i, n) t_0 = Float64(Float64((Float64(1.0 + Float64(i / n)) ^ n) - 1.0) / Float64(i / n)) tmp = 0.0 if (t_0 <= 1e-302) tmp = Float64(Float64(100.0 * expm1(Float64(log1p(Float64(i / n)) * n))) / Float64(i / n)); elseif (t_0 <= Inf) tmp = Float64(Float64(Float64(Float64((Float64(Float64(i / n) + 1.0) ^ n) / i) * n) * 100.0) - Float64(Float64(n * 100.0) / i)); else tmp = fma(Float64(0.01 * i), Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), Float64(0.01 / n)) ^ -1.0; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(N[Power[N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision], n], $MachinePrecision] - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, 1e-302], N[(N[(100.0 * N[(Exp[N[(N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision] * n), $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, Infinity], N[(N[(N[(N[(N[Power[N[(N[(i / n), $MachinePrecision] + 1.0), $MachinePrecision], n], $MachinePrecision] / i), $MachinePrecision] * n), $MachinePrecision] * 100.0), $MachinePrecision] - N[(N[(n * 100.0), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], N[Power[N[(N[(0.01 * i), $MachinePrecision] * N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 / n), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}\\
\mathbf{if}\;t\_0 \leq 10^{-302}:\\
\;\;\;\;\frac{100 \cdot \mathsf{expm1}\left(\mathsf{log1p}\left(\frac{i}{n}\right) \cdot n\right)}{\frac{i}{n}}\\
\mathbf{elif}\;t\_0 \leq \infty:\\
\;\;\;\;\left(\frac{{\left(\frac{i}{n} + 1\right)}^{n}}{i} \cdot n\right) \cdot 100 - \frac{n \cdot 100}{i}\\
\mathbf{else}:\\
\;\;\;\;{\left(\mathsf{fma}\left(0.01 \cdot i, \frac{0.5}{n \cdot n} - \frac{0.5}{n}, \frac{0.01}{n}\right)\right)}^{-1}\\
\end{array}
\end{array}
if (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < 9.9999999999999996e-303Initial program 24.9%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6424.9
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6497.9
Applied rewrites97.9%
if 9.9999999999999996e-303 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < +inf.0Initial program 96.7%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
distribute-rgt-inN/A
lower-+.f64N/A
Applied rewrites97.0%
lift-*.f64N/A
lift-/.f64N/A
associate-*l/N/A
lower-/.f64N/A
lower-*.f6497.1
Applied rewrites97.1%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f640.0
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f640.0
Applied rewrites0.0%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
*-lft-identityN/A
lift-*.f64N/A
lift-expm1.f64N/A
lift-*.f64N/A
lift-log1p.f64N/A
pow-to-expN/A
lift-/.f64N/A
times-fracN/A
lower-*.f64N/A
metadata-evalN/A
lower-/.f64N/A
Applied rewrites0.0%
Taylor expanded in i around 0
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6499.7
Applied rewrites99.7%
Final simplification98.2%
(FPCore (i n)
:precision binary64
(let* ((t_0 (/ (- (pow (+ 1.0 (/ i n)) n) 1.0) (/ i n))))
(if (<= t_0 1e-302)
(* (/ n i) (* 100.0 (expm1 (* (log1p (/ i n)) n))))
(if (<= t_0 INFINITY)
(- (* (* (/ (pow (+ (/ i n) 1.0) n) i) n) 100.0) (/ (* n 100.0) i))
(pow (fma (* 0.01 i) (- (/ 0.5 (* n n)) (/ 0.5 n)) (/ 0.01 n)) -1.0)))))
double code(double i, double n) {
double t_0 = (pow((1.0 + (i / n)), n) - 1.0) / (i / n);
double tmp;
if (t_0 <= 1e-302) {
tmp = (n / i) * (100.0 * expm1((log1p((i / n)) * n)));
} else if (t_0 <= ((double) INFINITY)) {
tmp = (((pow(((i / n) + 1.0), n) / i) * n) * 100.0) - ((n * 100.0) / i);
} else {
tmp = pow(fma((0.01 * i), ((0.5 / (n * n)) - (0.5 / n)), (0.01 / n)), -1.0);
}
return tmp;
}
function code(i, n) t_0 = Float64(Float64((Float64(1.0 + Float64(i / n)) ^ n) - 1.0) / Float64(i / n)) tmp = 0.0 if (t_0 <= 1e-302) tmp = Float64(Float64(n / i) * Float64(100.0 * expm1(Float64(log1p(Float64(i / n)) * n)))); elseif (t_0 <= Inf) tmp = Float64(Float64(Float64(Float64((Float64(Float64(i / n) + 1.0) ^ n) / i) * n) * 100.0) - Float64(Float64(n * 100.0) / i)); else tmp = fma(Float64(0.01 * i), Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), Float64(0.01 / n)) ^ -1.0; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(N[Power[N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision], n], $MachinePrecision] - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, 1e-302], N[(N[(n / i), $MachinePrecision] * N[(100.0 * N[(Exp[N[(N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision] * n), $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, Infinity], N[(N[(N[(N[(N[Power[N[(N[(i / n), $MachinePrecision] + 1.0), $MachinePrecision], n], $MachinePrecision] / i), $MachinePrecision] * n), $MachinePrecision] * 100.0), $MachinePrecision] - N[(N[(n * 100.0), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], N[Power[N[(N[(0.01 * i), $MachinePrecision] * N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 / n), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}\\
\mathbf{if}\;t\_0 \leq 10^{-302}:\\
\;\;\;\;\frac{n}{i} \cdot \left(100 \cdot \mathsf{expm1}\left(\mathsf{log1p}\left(\frac{i}{n}\right) \cdot n\right)\right)\\
\mathbf{elif}\;t\_0 \leq \infty:\\
\;\;\;\;\left(\frac{{\left(\frac{i}{n} + 1\right)}^{n}}{i} \cdot n\right) \cdot 100 - \frac{n \cdot 100}{i}\\
\mathbf{else}:\\
\;\;\;\;{\left(\mathsf{fma}\left(0.01 \cdot i, \frac{0.5}{n \cdot n} - \frac{0.5}{n}, \frac{0.01}{n}\right)\right)}^{-1}\\
\end{array}
\end{array}
if (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < 9.9999999999999996e-303Initial program 24.9%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
clear-numN/A
associate-/r/N/A
lift-/.f64N/A
clear-numN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f6424.4
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6496.0
Applied rewrites96.0%
if 9.9999999999999996e-303 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < +inf.0Initial program 96.7%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
distribute-rgt-inN/A
lower-+.f64N/A
Applied rewrites97.0%
lift-*.f64N/A
lift-/.f64N/A
associate-*l/N/A
lower-/.f64N/A
lower-*.f6497.1
Applied rewrites97.1%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f640.0
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f640.0
Applied rewrites0.0%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
*-lft-identityN/A
lift-*.f64N/A
lift-expm1.f64N/A
lift-*.f64N/A
lift-log1p.f64N/A
pow-to-expN/A
lift-/.f64N/A
times-fracN/A
lower-*.f64N/A
metadata-evalN/A
lower-/.f64N/A
Applied rewrites0.0%
Taylor expanded in i around 0
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6499.7
Applied rewrites99.7%
Final simplification96.7%
(FPCore (i n)
:precision binary64
(let* ((t_0 (/ (- (pow (+ 1.0 (/ i n)) n) 1.0) (/ i n))))
(if (<= t_0 1e-302)
(* (* (/ n i) 100.0) (expm1 (* (log1p (/ i n)) n)))
(if (<= t_0 INFINITY)
(- (* (* (/ (pow (+ (/ i n) 1.0) n) i) n) 100.0) (/ (* n 100.0) i))
(pow (fma (* 0.01 i) (- (/ 0.5 (* n n)) (/ 0.5 n)) (/ 0.01 n)) -1.0)))))
double code(double i, double n) {
double t_0 = (pow((1.0 + (i / n)), n) - 1.0) / (i / n);
double tmp;
if (t_0 <= 1e-302) {
tmp = ((n / i) * 100.0) * expm1((log1p((i / n)) * n));
} else if (t_0 <= ((double) INFINITY)) {
tmp = (((pow(((i / n) + 1.0), n) / i) * n) * 100.0) - ((n * 100.0) / i);
} else {
tmp = pow(fma((0.01 * i), ((0.5 / (n * n)) - (0.5 / n)), (0.01 / n)), -1.0);
}
return tmp;
}
function code(i, n) t_0 = Float64(Float64((Float64(1.0 + Float64(i / n)) ^ n) - 1.0) / Float64(i / n)) tmp = 0.0 if (t_0 <= 1e-302) tmp = Float64(Float64(Float64(n / i) * 100.0) * expm1(Float64(log1p(Float64(i / n)) * n))); elseif (t_0 <= Inf) tmp = Float64(Float64(Float64(Float64((Float64(Float64(i / n) + 1.0) ^ n) / i) * n) * 100.0) - Float64(Float64(n * 100.0) / i)); else tmp = fma(Float64(0.01 * i), Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), Float64(0.01 / n)) ^ -1.0; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(N[Power[N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision], n], $MachinePrecision] - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, 1e-302], N[(N[(N[(n / i), $MachinePrecision] * 100.0), $MachinePrecision] * N[(Exp[N[(N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision] * n), $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, Infinity], N[(N[(N[(N[(N[Power[N[(N[(i / n), $MachinePrecision] + 1.0), $MachinePrecision], n], $MachinePrecision] / i), $MachinePrecision] * n), $MachinePrecision] * 100.0), $MachinePrecision] - N[(N[(n * 100.0), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], N[Power[N[(N[(0.01 * i), $MachinePrecision] * N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 / n), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{{\left(1 + \frac{i}{n}\right)}^{n} - 1}{\frac{i}{n}}\\
\mathbf{if}\;t\_0 \leq 10^{-302}:\\
\;\;\;\;\left(\frac{n}{i} \cdot 100\right) \cdot \mathsf{expm1}\left(\mathsf{log1p}\left(\frac{i}{n}\right) \cdot n\right)\\
\mathbf{elif}\;t\_0 \leq \infty:\\
\;\;\;\;\left(\frac{{\left(\frac{i}{n} + 1\right)}^{n}}{i} \cdot n\right) \cdot 100 - \frac{n \cdot 100}{i}\\
\mathbf{else}:\\
\;\;\;\;{\left(\mathsf{fma}\left(0.01 \cdot i, \frac{0.5}{n \cdot n} - \frac{0.5}{n}, \frac{0.01}{n}\right)\right)}^{-1}\\
\end{array}
\end{array}
if (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < 9.9999999999999996e-303Initial program 24.9%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
associate-/r/N/A
lift-/.f64N/A
clear-numN/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6424.4
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6495.6
Applied rewrites95.6%
if 9.9999999999999996e-303 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < +inf.0Initial program 96.7%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
distribute-rgt-inN/A
lower-+.f64N/A
Applied rewrites97.0%
lift-*.f64N/A
lift-/.f64N/A
associate-*l/N/A
lower-/.f64N/A
lower-*.f6497.1
Applied rewrites97.1%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f640.0
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f640.0
Applied rewrites0.0%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
*-lft-identityN/A
lift-*.f64N/A
lift-expm1.f64N/A
lift-*.f64N/A
lift-log1p.f64N/A
pow-to-expN/A
lift-/.f64N/A
times-fracN/A
lower-*.f64N/A
metadata-evalN/A
lower-/.f64N/A
Applied rewrites0.0%
Taylor expanded in i around 0
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6499.7
Applied rewrites99.7%
Final simplification96.4%
(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 0.0)
(* (* (/ (expm1 i) i) 100.0) n)
(if (<= t_1 INFINITY)
(* (* (/ n i) 100.0) t_0)
(pow (fma (* 0.01 i) (- (/ 0.5 (* n n)) (/ 0.5 n)) (/ 0.01 n)) -1.0)))))
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 <= 0.0) {
tmp = ((expm1(i) / i) * 100.0) * n;
} else if (t_1 <= ((double) INFINITY)) {
tmp = ((n / i) * 100.0) * t_0;
} else {
tmp = pow(fma((0.01 * i), ((0.5 / (n * n)) - (0.5 / n)), (0.01 / n)), -1.0);
}
return tmp;
}
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 <= 0.0) tmp = Float64(Float64(Float64(expm1(i) / i) * 100.0) * n); elseif (t_1 <= Inf) tmp = Float64(Float64(Float64(n / i) * 100.0) * t_0); else tmp = fma(Float64(0.01 * i), Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), Float64(0.01 / n)) ^ -1.0; end return tmp end
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, 0.0], N[(N[(N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(N[(N[(n / i), $MachinePrecision] * 100.0), $MachinePrecision] * t$95$0), $MachinePrecision], N[Power[N[(N[(0.01 * i), $MachinePrecision] * N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 / n), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision]]]]]
\begin{array}{l}
\\
\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 0:\\
\;\;\;\;\left(\frac{\mathsf{expm1}\left(i\right)}{i} \cdot 100\right) \cdot n\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;\left(\frac{n}{i} \cdot 100\right) \cdot t\_0\\
\mathbf{else}:\\
\;\;\;\;{\left(\mathsf{fma}\left(0.01 \cdot i, \frac{0.5}{n \cdot n} - \frac{0.5}{n}, \frac{0.01}{n}\right)\right)}^{-1}\\
\end{array}
\end{array}
if (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < 0.0Initial program 24.6%
Taylor expanded in n around inf
associate-/l*N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6475.8
Applied rewrites75.8%
if 0.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < +inf.0Initial program 95.7%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
associate-/r/N/A
lift-/.f64N/A
clear-numN/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6495.9
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6452.5
Applied rewrites52.5%
lift-expm1.f64N/A
lower--.f64N/A
lift-*.f64N/A
lift-log1p.f64N/A
exp-to-powN/A
+-commutativeN/A
lift-+.f64N/A
lift-pow.f6495.9
lift-+.f64N/A
+-commutativeN/A
lower-+.f6495.9
Applied rewrites95.9%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f640.0
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f640.0
Applied rewrites0.0%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
*-lft-identityN/A
lift-*.f64N/A
lift-expm1.f64N/A
lift-*.f64N/A
lift-log1p.f64N/A
pow-to-expN/A
lift-/.f64N/A
times-fracN/A
lower-*.f64N/A
metadata-evalN/A
lower-/.f64N/A
Applied rewrites0.0%
Taylor expanded in i around 0
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6499.7
Applied rewrites99.7%
Final simplification81.5%
(FPCore (i n) :precision binary64 (if (or (<= n -2.6e-209) (not (<= n 1.6e-7))) (* (* (/ (expm1 i) i) 100.0) n) (pow (fma (* 0.01 i) (- (/ 0.5 (* n n)) (/ 0.5 n)) (/ 0.01 n)) -1.0)))
double code(double i, double n) {
double tmp;
if ((n <= -2.6e-209) || !(n <= 1.6e-7)) {
tmp = ((expm1(i) / i) * 100.0) * n;
} else {
tmp = pow(fma((0.01 * i), ((0.5 / (n * n)) - (0.5 / n)), (0.01 / n)), -1.0);
}
return tmp;
}
function code(i, n) tmp = 0.0 if ((n <= -2.6e-209) || !(n <= 1.6e-7)) tmp = Float64(Float64(Float64(expm1(i) / i) * 100.0) * n); else tmp = fma(Float64(0.01 * i), Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), Float64(0.01 / n)) ^ -1.0; end return tmp end
code[i_, n_] := If[Or[LessEqual[n, -2.6e-209], N[Not[LessEqual[n, 1.6e-7]], $MachinePrecision]], N[(N[(N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], N[Power[N[(N[(0.01 * i), $MachinePrecision] * N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 / n), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -2.6 \cdot 10^{-209} \lor \neg \left(n \leq 1.6 \cdot 10^{-7}\right):\\
\;\;\;\;\left(\frac{\mathsf{expm1}\left(i\right)}{i} \cdot 100\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;{\left(\mathsf{fma}\left(0.01 \cdot i, \frac{0.5}{n \cdot n} - \frac{0.5}{n}, \frac{0.01}{n}\right)\right)}^{-1}\\
\end{array}
\end{array}
if n < -2.59999999999999984e-209 or 1.6e-7 < n Initial program 23.4%
Taylor expanded in n around inf
associate-/l*N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6485.6
Applied rewrites85.6%
if -2.59999999999999984e-209 < n < 1.6e-7Initial program 31.1%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6431.1
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6484.1
Applied rewrites84.1%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
*-lft-identityN/A
lift-*.f64N/A
lift-expm1.f64N/A
lift-*.f64N/A
lift-log1p.f64N/A
pow-to-expN/A
lift-/.f64N/A
times-fracN/A
lower-*.f64N/A
metadata-evalN/A
lower-/.f64N/A
Applied rewrites82.5%
Taylor expanded in i around 0
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6464.1
Applied rewrites64.1%
Final simplification79.6%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* (- 0.5 (/ 0.5 n)) n)))
(if (<= n -1.5e+95)
(/
(fma
(fma
(* 100.0 i)
(fma
(* i n)
(- (+ (/ 0.3333333333333333 (* n n)) 0.16666666666666666) (/ 0.5 n))
t_0)
(* 100.0 n))
i
0.0)
i)
(if (<= n 1.6e-7)
(pow (fma (* 0.01 i) (- (/ 0.5 (* n n)) (/ 0.5 n)) (/ 0.01 n)) -1.0)
(/ (fma (* 100.0 (fma t_0 i n)) i 0.0) i)))))
double code(double i, double n) {
double t_0 = (0.5 - (0.5 / n)) * n;
double tmp;
if (n <= -1.5e+95) {
tmp = fma(fma((100.0 * i), fma((i * n), (((0.3333333333333333 / (n * n)) + 0.16666666666666666) - (0.5 / n)), t_0), (100.0 * n)), i, 0.0) / i;
} else if (n <= 1.6e-7) {
tmp = pow(fma((0.01 * i), ((0.5 / (n * n)) - (0.5 / n)), (0.01 / n)), -1.0);
} else {
tmp = fma((100.0 * fma(t_0, i, n)), i, 0.0) / i;
}
return tmp;
}
function code(i, n) t_0 = Float64(Float64(0.5 - Float64(0.5 / n)) * n) tmp = 0.0 if (n <= -1.5e+95) tmp = Float64(fma(fma(Float64(100.0 * i), fma(Float64(i * n), Float64(Float64(Float64(0.3333333333333333 / Float64(n * n)) + 0.16666666666666666) - Float64(0.5 / n)), t_0), Float64(100.0 * n)), i, 0.0) / i); elseif (n <= 1.6e-7) tmp = fma(Float64(0.01 * i), Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), Float64(0.01 / n)) ^ -1.0; else tmp = Float64(fma(Float64(100.0 * fma(t_0, i, n)), i, 0.0) / i); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision]}, If[LessEqual[n, -1.5e+95], N[(N[(N[(N[(100.0 * i), $MachinePrecision] * N[(N[(i * n), $MachinePrecision] * N[(N[(N[(0.3333333333333333 / N[(n * n), $MachinePrecision]), $MachinePrecision] + 0.16666666666666666), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] + t$95$0), $MachinePrecision] + N[(100.0 * n), $MachinePrecision]), $MachinePrecision] * i + 0.0), $MachinePrecision] / i), $MachinePrecision], If[LessEqual[n, 1.6e-7], N[Power[N[(N[(0.01 * i), $MachinePrecision] * N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 / n), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision], N[(N[(N[(100.0 * N[(t$95$0 * i + n), $MachinePrecision]), $MachinePrecision] * i + 0.0), $MachinePrecision] / i), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(0.5 - \frac{0.5}{n}\right) \cdot n\\
\mathbf{if}\;n \leq -1.5 \cdot 10^{+95}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(100 \cdot i, \mathsf{fma}\left(i \cdot n, \left(\frac{0.3333333333333333}{n \cdot n} + 0.16666666666666666\right) - \frac{0.5}{n}, t\_0\right), 100 \cdot n\right), i, 0\right)}{i}\\
\mathbf{elif}\;n \leq 1.6 \cdot 10^{-7}:\\
\;\;\;\;{\left(\mathsf{fma}\left(0.01 \cdot i, \frac{0.5}{n \cdot n} - \frac{0.5}{n}, \frac{0.01}{n}\right)\right)}^{-1}\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{fma}\left(100 \cdot \mathsf{fma}\left(t\_0, i, n\right), i, 0\right)}{i}\\
\end{array}
\end{array}
if n < -1.49999999999999996e95Initial program 17.3%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
distribute-rgt-inN/A
lower-+.f64N/A
Applied rewrites17.4%
Taylor expanded in i around 0
distribute-rgt-outN/A
metadata-evalN/A
mul0-rgtN/A
lower-/.f642.1
Applied rewrites2.1%
Taylor expanded in i around 0
Applied rewrites2.1%
Taylor expanded in i around 0
lower-/.f64N/A
Applied rewrites61.5%
if -1.49999999999999996e95 < n < 1.6e-7Initial program 31.7%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6431.8
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6485.4
Applied rewrites85.4%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
*-lft-identityN/A
lift-*.f64N/A
lift-expm1.f64N/A
lift-*.f64N/A
lift-log1p.f64N/A
pow-to-expN/A
lift-/.f64N/A
times-fracN/A
lower-*.f64N/A
metadata-evalN/A
lower-/.f64N/A
Applied rewrites84.1%
Taylor expanded in i around 0
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6464.6
Applied rewrites64.6%
if 1.6e-7 < n Initial program 18.2%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
distribute-rgt-inN/A
lower-+.f64N/A
Applied rewrites18.5%
Taylor expanded in i around 0
lower-/.f64N/A
Applied rewrites84.0%
Final simplification68.9%
(FPCore (i n)
:precision binary64
(if (<= n -2.1e+97)
(* (fma (+ 50.0 (* i 16.666666666666668)) i 100.0) n)
(if (<= n 1.6e-7)
(pow (fma (* 0.01 i) (- (/ 0.5 (* n n)) (/ 0.5 n)) (/ 0.01 n)) -1.0)
(/ (fma (* 100.0 (fma (* (- 0.5 (/ 0.5 n)) n) i n)) i 0.0) i))))
double code(double i, double n) {
double tmp;
if (n <= -2.1e+97) {
tmp = fma((50.0 + (i * 16.666666666666668)), i, 100.0) * n;
} else if (n <= 1.6e-7) {
tmp = pow(fma((0.01 * i), ((0.5 / (n * n)) - (0.5 / n)), (0.01 / n)), -1.0);
} else {
tmp = fma((100.0 * fma(((0.5 - (0.5 / n)) * n), i, n)), i, 0.0) / i;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -2.1e+97) tmp = Float64(fma(Float64(50.0 + Float64(i * 16.666666666666668)), i, 100.0) * n); elseif (n <= 1.6e-7) tmp = fma(Float64(0.01 * i), Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), Float64(0.01 / n)) ^ -1.0; else tmp = Float64(fma(Float64(100.0 * fma(Float64(Float64(0.5 - Float64(0.5 / n)) * n), i, n)), i, 0.0) / i); end return tmp end
code[i_, n_] := If[LessEqual[n, -2.1e+97], N[(N[(N[(50.0 + N[(i * 16.666666666666668), $MachinePrecision]), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[n, 1.6e-7], N[Power[N[(N[(0.01 * i), $MachinePrecision] * N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 / n), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision], N[(N[(N[(100.0 * N[(N[(N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision] * i + n), $MachinePrecision]), $MachinePrecision] * i + 0.0), $MachinePrecision] / i), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -2.1 \cdot 10^{+97}:\\
\;\;\;\;\mathsf{fma}\left(50 + i \cdot 16.666666666666668, i, 100\right) \cdot n\\
\mathbf{elif}\;n \leq 1.6 \cdot 10^{-7}:\\
\;\;\;\;{\left(\mathsf{fma}\left(0.01 \cdot i, \frac{0.5}{n \cdot n} - \frac{0.5}{n}, \frac{0.01}{n}\right)\right)}^{-1}\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{fma}\left(100 \cdot \mathsf{fma}\left(\left(0.5 - \frac{0.5}{n}\right) \cdot n, i, n\right), i, 0\right)}{i}\\
\end{array}
\end{array}
if n < -2.10000000000000012e97Initial program 17.3%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6417.3
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6471.3
Applied rewrites71.3%
Taylor expanded in i around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites59.7%
Taylor expanded in n around inf
Applied rewrites59.7%
if -2.10000000000000012e97 < n < 1.6e-7Initial program 31.7%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6431.8
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6485.4
Applied rewrites85.4%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
*-lft-identityN/A
lift-*.f64N/A
lift-expm1.f64N/A
lift-*.f64N/A
lift-log1p.f64N/A
pow-to-expN/A
lift-/.f64N/A
times-fracN/A
lower-*.f64N/A
metadata-evalN/A
lower-/.f64N/A
Applied rewrites84.1%
Taylor expanded in i around 0
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6464.6
Applied rewrites64.6%
if 1.6e-7 < n Initial program 18.2%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
distribute-rgt-inN/A
lower-+.f64N/A
Applied rewrites18.5%
Taylor expanded in i around 0
lower-/.f64N/A
Applied rewrites84.0%
Final simplification68.6%
(FPCore (i n)
:precision binary64
(if (<= n -7.5e+66)
(* (fma (+ 50.0 (* i 16.666666666666668)) i 100.0) n)
(if (<= n 1.5e-7)
(/ (* 100.0 i) (/ i n))
(/ (fma (* 100.0 (fma (* (- 0.5 (/ 0.5 n)) n) i n)) i 0.0) i))))
double code(double i, double n) {
double tmp;
if (n <= -7.5e+66) {
tmp = fma((50.0 + (i * 16.666666666666668)), i, 100.0) * n;
} else if (n <= 1.5e-7) {
tmp = (100.0 * i) / (i / n);
} else {
tmp = fma((100.0 * fma(((0.5 - (0.5 / n)) * n), i, n)), i, 0.0) / i;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -7.5e+66) tmp = Float64(fma(Float64(50.0 + Float64(i * 16.666666666666668)), i, 100.0) * n); elseif (n <= 1.5e-7) tmp = Float64(Float64(100.0 * i) / Float64(i / n)); else tmp = Float64(fma(Float64(100.0 * fma(Float64(Float64(0.5 - Float64(0.5 / n)) * n), i, n)), i, 0.0) / i); end return tmp end
code[i_, n_] := If[LessEqual[n, -7.5e+66], N[(N[(N[(50.0 + N[(i * 16.666666666666668), $MachinePrecision]), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[n, 1.5e-7], N[(N[(100.0 * i), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision], N[(N[(N[(100.0 * N[(N[(N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision] * i + n), $MachinePrecision]), $MachinePrecision] * i + 0.0), $MachinePrecision] / i), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -7.5 \cdot 10^{+66}:\\
\;\;\;\;\mathsf{fma}\left(50 + i \cdot 16.666666666666668, i, 100\right) \cdot n\\
\mathbf{elif}\;n \leq 1.5 \cdot 10^{-7}:\\
\;\;\;\;\frac{100 \cdot i}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{fma}\left(100 \cdot \mathsf{fma}\left(\left(0.5 - \frac{0.5}{n}\right) \cdot n, i, n\right), i, 0\right)}{i}\\
\end{array}
\end{array}
if n < -7.50000000000000024e66Initial program 17.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6417.3
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6467.1
Applied rewrites67.1%
Taylor expanded in i around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites62.3%
Taylor expanded in n around inf
Applied rewrites62.3%
if -7.50000000000000024e66 < n < 1.4999999999999999e-7Initial program 32.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6432.4
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6487.7
Applied rewrites87.7%
Taylor expanded in i around 0
lower-*.f6460.2
Applied rewrites60.2%
if 1.4999999999999999e-7 < n Initial program 18.2%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
distribute-rgt-inN/A
lower-+.f64N/A
Applied rewrites18.5%
Taylor expanded in i around 0
lower-/.f64N/A
Applied rewrites84.0%
Final simplification66.7%
(FPCore (i n)
:precision binary64
(if (<= n -7.5e+66)
(* (fma (+ 50.0 (* i 16.666666666666668)) i 100.0) n)
(if (<= n 1.5e-7)
(/ (* 100.0 i) (/ i n))
(/ (* (* (fma (* (- 0.5 (/ 0.5 n)) i) 100.0 100.0) i) n) i))))
double code(double i, double n) {
double tmp;
if (n <= -7.5e+66) {
tmp = fma((50.0 + (i * 16.666666666666668)), i, 100.0) * n;
} else if (n <= 1.5e-7) {
tmp = (100.0 * i) / (i / n);
} else {
tmp = ((fma(((0.5 - (0.5 / n)) * i), 100.0, 100.0) * i) * n) / i;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -7.5e+66) tmp = Float64(fma(Float64(50.0 + Float64(i * 16.666666666666668)), i, 100.0) * n); elseif (n <= 1.5e-7) tmp = Float64(Float64(100.0 * i) / Float64(i / n)); else tmp = Float64(Float64(Float64(fma(Float64(Float64(0.5 - Float64(0.5 / n)) * i), 100.0, 100.0) * i) * n) / i); end return tmp end
code[i_, n_] := If[LessEqual[n, -7.5e+66], N[(N[(N[(50.0 + N[(i * 16.666666666666668), $MachinePrecision]), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[n, 1.5e-7], N[(N[(100.0 * i), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(N[(N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * i), $MachinePrecision] * 100.0 + 100.0), $MachinePrecision] * i), $MachinePrecision] * n), $MachinePrecision] / i), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -7.5 \cdot 10^{+66}:\\
\;\;\;\;\mathsf{fma}\left(50 + i \cdot 16.666666666666668, i, 100\right) \cdot n\\
\mathbf{elif}\;n \leq 1.5 \cdot 10^{-7}:\\
\;\;\;\;\frac{100 \cdot i}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(\mathsf{fma}\left(\left(0.5 - \frac{0.5}{n}\right) \cdot i, 100, 100\right) \cdot i\right) \cdot n}{i}\\
\end{array}
\end{array}
if n < -7.50000000000000024e66Initial program 17.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6417.3
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6467.1
Applied rewrites67.1%
Taylor expanded in i around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites62.3%
Taylor expanded in n around inf
Applied rewrites62.3%
if -7.50000000000000024e66 < n < 1.4999999999999999e-7Initial program 32.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6432.4
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6487.7
Applied rewrites87.7%
Taylor expanded in i around 0
lower-*.f6460.2
Applied rewrites60.2%
if 1.4999999999999999e-7 < n Initial program 18.2%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6418.2
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6464.3
Applied rewrites64.3%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6452.5
Applied rewrites52.5%
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
associate-*l/N/A
lower-/.f64N/A
lower-*.f6483.9
Applied rewrites83.9%
(FPCore (i n)
:precision binary64
(if (<= n -7.5e+66)
(* (fma (+ 50.0 (* i 16.666666666666668)) i 100.0) n)
(if (<= n 1.5e-7)
(/ (* 100.0 i) (/ i n))
(* (/ (* (fma (* (- 0.5 (/ 0.5 n)) i) 100.0 100.0) i) i) n))))
double code(double i, double n) {
double tmp;
if (n <= -7.5e+66) {
tmp = fma((50.0 + (i * 16.666666666666668)), i, 100.0) * n;
} else if (n <= 1.5e-7) {
tmp = (100.0 * i) / (i / n);
} else {
tmp = ((fma(((0.5 - (0.5 / n)) * i), 100.0, 100.0) * i) / i) * n;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -7.5e+66) tmp = Float64(fma(Float64(50.0 + Float64(i * 16.666666666666668)), i, 100.0) * n); elseif (n <= 1.5e-7) tmp = Float64(Float64(100.0 * i) / Float64(i / n)); else tmp = Float64(Float64(Float64(fma(Float64(Float64(0.5 - Float64(0.5 / n)) * i), 100.0, 100.0) * i) / i) * n); end return tmp end
code[i_, n_] := If[LessEqual[n, -7.5e+66], N[(N[(N[(50.0 + N[(i * 16.666666666666668), $MachinePrecision]), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[n, 1.5e-7], N[(N[(100.0 * i), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(N[(N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * i), $MachinePrecision] * 100.0 + 100.0), $MachinePrecision] * i), $MachinePrecision] / i), $MachinePrecision] * n), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -7.5 \cdot 10^{+66}:\\
\;\;\;\;\mathsf{fma}\left(50 + i \cdot 16.666666666666668, i, 100\right) \cdot n\\
\mathbf{elif}\;n \leq 1.5 \cdot 10^{-7}:\\
\;\;\;\;\frac{100 \cdot i}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\left(0.5 - \frac{0.5}{n}\right) \cdot i, 100, 100\right) \cdot i}{i} \cdot n\\
\end{array}
\end{array}
if n < -7.50000000000000024e66Initial program 17.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6417.3
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6467.1
Applied rewrites67.1%
Taylor expanded in i around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites62.3%
Taylor expanded in n around inf
Applied rewrites62.3%
if -7.50000000000000024e66 < n < 1.4999999999999999e-7Initial program 32.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6432.4
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6487.7
Applied rewrites87.7%
Taylor expanded in i around 0
lower-*.f6460.2
Applied rewrites60.2%
if 1.4999999999999999e-7 < n Initial program 18.2%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6418.2
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6464.3
Applied rewrites64.3%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6452.5
Applied rewrites52.5%
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f6481.2
Applied rewrites81.2%
(FPCore (i n) :precision binary64 (if (or (<= n -7.5e+66) (not (<= n 7e-8))) (* (fma (+ 50.0 (* i 16.666666666666668)) i 100.0) n) (/ (* 100.0 i) (/ i n))))
double code(double i, double n) {
double tmp;
if ((n <= -7.5e+66) || !(n <= 7e-8)) {
tmp = fma((50.0 + (i * 16.666666666666668)), i, 100.0) * n;
} else {
tmp = (100.0 * i) / (i / n);
}
return tmp;
}
function code(i, n) tmp = 0.0 if ((n <= -7.5e+66) || !(n <= 7e-8)) tmp = Float64(fma(Float64(50.0 + Float64(i * 16.666666666666668)), i, 100.0) * n); else tmp = Float64(Float64(100.0 * i) / Float64(i / n)); end return tmp end
code[i_, n_] := If[Or[LessEqual[n, -7.5e+66], N[Not[LessEqual[n, 7e-8]], $MachinePrecision]], N[(N[(N[(50.0 + N[(i * 16.666666666666668), $MachinePrecision]), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision], N[(N[(100.0 * i), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -7.5 \cdot 10^{+66} \lor \neg \left(n \leq 7 \cdot 10^{-8}\right):\\
\;\;\;\;\mathsf{fma}\left(50 + i \cdot 16.666666666666668, i, 100\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;\frac{100 \cdot i}{\frac{i}{n}}\\
\end{array}
\end{array}
if n < -7.50000000000000024e66 or 7.00000000000000048e-8 < n Initial program 17.8%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6417.8
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6465.6
Applied rewrites65.6%
Taylor expanded in i around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites72.0%
Taylor expanded in n around inf
Applied rewrites72.0%
if -7.50000000000000024e66 < n < 7.00000000000000048e-8Initial program 32.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6432.4
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6487.7
Applied rewrites87.7%
Taylor expanded in i around 0
lower-*.f6460.2
Applied rewrites60.2%
Final simplification65.7%
(FPCore (i n) :precision binary64 (if (or (<= n -4.2e-10) (not (<= n 1.45e-7))) (* (fma (+ 50.0 (* i 16.666666666666668)) i 100.0) n) (* (/ n i) (* 100.0 i))))
double code(double i, double n) {
double tmp;
if ((n <= -4.2e-10) || !(n <= 1.45e-7)) {
tmp = fma((50.0 + (i * 16.666666666666668)), i, 100.0) * n;
} else {
tmp = (n / i) * (100.0 * i);
}
return tmp;
}
function code(i, n) tmp = 0.0 if ((n <= -4.2e-10) || !(n <= 1.45e-7)) tmp = Float64(fma(Float64(50.0 + Float64(i * 16.666666666666668)), i, 100.0) * n); else tmp = Float64(Float64(n / i) * Float64(100.0 * i)); end return tmp end
code[i_, n_] := If[Or[LessEqual[n, -4.2e-10], N[Not[LessEqual[n, 1.45e-7]], $MachinePrecision]], N[(N[(N[(50.0 + N[(i * 16.666666666666668), $MachinePrecision]), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision], N[(N[(n / i), $MachinePrecision] * N[(100.0 * i), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -4.2 \cdot 10^{-10} \lor \neg \left(n \leq 1.45 \cdot 10^{-7}\right):\\
\;\;\;\;\mathsf{fma}\left(50 + i \cdot 16.666666666666668, i, 100\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;\frac{n}{i} \cdot \left(100 \cdot i\right)\\
\end{array}
\end{array}
if n < -4.2e-10 or 1.4499999999999999e-7 < n Initial program 22.3%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6422.3
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6466.2
Applied rewrites66.2%
Taylor expanded in i around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites69.0%
Taylor expanded in n around inf
Applied rewrites69.0%
if -4.2e-10 < n < 1.4499999999999999e-7Initial program 29.3%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
clear-numN/A
associate-/r/N/A
lift-/.f64N/A
clear-numN/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f6428.5
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6489.4
Applied rewrites89.4%
Taylor expanded in i around 0
lower-*.f6460.5
Applied rewrites60.5%
Final simplification65.1%
(FPCore (i n) :precision binary64 (if (<= i -7.2) 0.0 (* (fma (+ 50.0 (* i 16.666666666666668)) i 100.0) n)))
double code(double i, double n) {
double tmp;
if (i <= -7.2) {
tmp = 0.0;
} else {
tmp = fma((50.0 + (i * 16.666666666666668)), i, 100.0) * n;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (i <= -7.2) tmp = 0.0; else tmp = Float64(fma(Float64(50.0 + Float64(i * 16.666666666666668)), i, 100.0) * n); end return tmp end
code[i_, n_] := If[LessEqual[i, -7.2], 0.0, N[(N[(N[(50.0 + N[(i * 16.666666666666668), $MachinePrecision]), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -7.2:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(50 + i \cdot 16.666666666666668, i, 100\right) \cdot n\\
\end{array}
\end{array}
if i < -7.20000000000000018Initial program 61.9%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
distribute-rgt-inN/A
lower-+.f64N/A
Applied rewrites60.1%
Taylor expanded in i around 0
distribute-rgt-outN/A
metadata-evalN/A
mul0-rgtN/A
lower-/.f6436.0
Applied rewrites36.0%
Taylor expanded in i around 0
Applied rewrites36.0%
if -7.20000000000000018 < i Initial program 17.3%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6417.4
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6473.8
Applied rewrites73.8%
Taylor expanded in i around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites60.9%
Taylor expanded in n around inf
Applied rewrites69.5%
Final simplification63.4%
(FPCore (i n) :precision binary64 (if (<= i -7.2) 0.0 (if (<= i 1.8e+20) (* 100.0 n) 0.0)))
double code(double i, double n) {
double tmp;
if (i <= -7.2) {
tmp = 0.0;
} else if (i <= 1.8e+20) {
tmp = 100.0 * n;
} else {
tmp = 0.0;
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if (i <= (-7.2d0)) then
tmp = 0.0d0
else if (i <= 1.8d+20) then
tmp = 100.0d0 * n
else
tmp = 0.0d0
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (i <= -7.2) {
tmp = 0.0;
} else if (i <= 1.8e+20) {
tmp = 100.0 * n;
} else {
tmp = 0.0;
}
return tmp;
}
def code(i, n): tmp = 0 if i <= -7.2: tmp = 0.0 elif i <= 1.8e+20: tmp = 100.0 * n else: tmp = 0.0 return tmp
function code(i, n) tmp = 0.0 if (i <= -7.2) tmp = 0.0; elseif (i <= 1.8e+20) tmp = Float64(100.0 * n); else tmp = 0.0; end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (i <= -7.2) tmp = 0.0; elseif (i <= 1.8e+20) tmp = 100.0 * n; else tmp = 0.0; end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[i, -7.2], 0.0, If[LessEqual[i, 1.8e+20], N[(100.0 * n), $MachinePrecision], 0.0]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -7.2:\\
\;\;\;\;0\\
\mathbf{elif}\;i \leq 1.8 \cdot 10^{+20}:\\
\;\;\;\;100 \cdot n\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}
\end{array}
if i < -7.20000000000000018 or 1.8e20 < i Initial program 50.4%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
distribute-rgt-inN/A
lower-+.f64N/A
Applied rewrites47.6%
Taylor expanded in i around 0
distribute-rgt-outN/A
metadata-evalN/A
mul0-rgtN/A
lower-/.f6433.1
Applied rewrites33.1%
Taylor expanded in i around 0
Applied rewrites33.1%
if -7.20000000000000018 < i < 1.8e20Initial program 7.4%
Taylor expanded in i around 0
lower-*.f6483.2
Applied rewrites83.2%
Final simplification62.0%
(FPCore (i n) :precision binary64 0.0)
double code(double i, double n) {
return 0.0;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
code = 0.0d0
end function
public static double code(double i, double n) {
return 0.0;
}
def code(i, n): return 0.0
function code(i, n) return 0.0 end
function tmp = code(i, n) tmp = 0.0; end
code[i_, n_] := 0.0
\begin{array}{l}
\\
0
\end{array}
Initial program 25.5%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
distribute-rgt-inN/A
lower-+.f64N/A
Applied rewrites24.4%
Taylor expanded in i around 0
distribute-rgt-outN/A
metadata-evalN/A
mul0-rgtN/A
lower-/.f6418.5
Applied rewrites18.5%
Taylor expanded in i around 0
Applied rewrites18.5%
Final simplification18.5%
(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 2024298
(FPCore (i n)
:name "Compound Interest"
:precision binary64
:alt
(! :herbie-platform default (let ((lnbase (if (== (+ 1 (/ i n)) 1) (/ i n) (/ (* (/ i n) (log (+ 1 (/ i n)))) (- (+ (/ i n) 1) 1))))) (* 100 (/ (- (exp (* n lnbase)) 1) (/ i n)))))
(* 100.0 (/ (- (pow (+ 1.0 (/ i n)) n) 1.0) (/ i n))))