
(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 16 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))
(t_1 (/ t_0 (/ i n)))
(t_2 (* (/ t_0 i) (* n 100.0))))
(if (<= t_1 (- INFINITY))
t_2
(if (<= t_1 1e-301)
(* (* (expm1 (* (log1p (/ i n)) n)) (/ 100.0 i)) n)
(if (<= t_1 INFINITY) t_2 (* 100.0 n))))))
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 t_2 = (t_0 / i) * (n * 100.0);
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = t_2;
} else if (t_1 <= 1e-301) {
tmp = (expm1((log1p((i / n)) * n)) * (100.0 / i)) * n;
} else if (t_1 <= ((double) INFINITY)) {
tmp = t_2;
} else {
tmp = 100.0 * n;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = Math.pow((1.0 + (i / n)), n) - 1.0;
double t_1 = t_0 / (i / n);
double t_2 = (t_0 / i) * (n * 100.0);
double tmp;
if (t_1 <= -Double.POSITIVE_INFINITY) {
tmp = t_2;
} else if (t_1 <= 1e-301) {
tmp = (Math.expm1((Math.log1p((i / n)) * n)) * (100.0 / i)) * n;
} else if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = t_2;
} else {
tmp = 100.0 * n;
}
return tmp;
}
def code(i, n): t_0 = math.pow((1.0 + (i / n)), n) - 1.0 t_1 = t_0 / (i / n) t_2 = (t_0 / i) * (n * 100.0) tmp = 0 if t_1 <= -math.inf: tmp = t_2 elif t_1 <= 1e-301: tmp = (math.expm1((math.log1p((i / n)) * n)) * (100.0 / i)) * n elif t_1 <= math.inf: tmp = t_2 else: tmp = 100.0 * n 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)) t_2 = Float64(Float64(t_0 / i) * Float64(n * 100.0)) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = t_2; elseif (t_1 <= 1e-301) tmp = Float64(Float64(expm1(Float64(log1p(Float64(i / n)) * n)) * Float64(100.0 / i)) * n); elseif (t_1 <= Inf) tmp = t_2; else tmp = Float64(100.0 * n); 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]}, Block[{t$95$2 = N[(N[(t$95$0 / i), $MachinePrecision] * N[(n * 100.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], t$95$2, If[LessEqual[t$95$1, 1e-301], N[(N[(N[(Exp[N[(N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision] * n), $MachinePrecision]] - 1), $MachinePrecision] * N[(100.0 / i), $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[t$95$1, Infinity], t$95$2, N[(100.0 * n), $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}}\\
t_2 := \frac{t\_0}{i} \cdot \left(n \cdot 100\right)\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq 10^{-301}:\\
\;\;\;\;\left(\mathsf{expm1}\left(\mathsf{log1p}\left(\frac{i}{n}\right) \cdot n\right) \cdot \frac{100}{i}\right) \cdot n\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;100 \cdot n\\
\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)) < -inf.0 or 1.00000000000000007e-301 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < +inf.0Initial program 99.7%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
associate-*l*N/A
lower-*.f64N/A
lower-/.f64N/A
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f64N/A
lower-*.f6443.5
Applied rewrites43.5%
lift-expm1.f64N/A
lift-*.f64N/A
lift-log1p.f64N/A
lift-/.f64N/A
pow-to-expN/A
lower--.f64N/A
lower-pow.f64N/A
lift-/.f64N/A
lower-+.f6499.7
Applied rewrites99.7%
if -inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < 1.00000000000000007e-301Initial program 20.8%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f64N/A
lower-/.f6497.9
Applied rewrites97.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%
Taylor expanded in i around 0
lower-*.f6477.0
Applied rewrites77.0%
(FPCore (i n)
:precision binary64
(if (<= n -5e-311)
(* (* (expm1 i) (/ 100.0 i)) n)
(if (<= n 7.5e-111)
(* (* n (/ (- (log i) (log n)) i)) (* n 100.0))
(* (/ (expm1 i) i) (* 100.0 n)))))
double code(double i, double n) {
double tmp;
if (n <= -5e-311) {
tmp = (expm1(i) * (100.0 / i)) * n;
} else if (n <= 7.5e-111) {
tmp = (n * ((log(i) - log(n)) / i)) * (n * 100.0);
} else {
tmp = (expm1(i) / i) * (100.0 * n);
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if (n <= -5e-311) {
tmp = (Math.expm1(i) * (100.0 / i)) * n;
} else if (n <= 7.5e-111) {
tmp = (n * ((Math.log(i) - Math.log(n)) / i)) * (n * 100.0);
} else {
tmp = (Math.expm1(i) / i) * (100.0 * n);
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -5e-311: tmp = (math.expm1(i) * (100.0 / i)) * n elif n <= 7.5e-111: tmp = (n * ((math.log(i) - math.log(n)) / i)) * (n * 100.0) else: tmp = (math.expm1(i) / i) * (100.0 * n) return tmp
function code(i, n) tmp = 0.0 if (n <= -5e-311) tmp = Float64(Float64(expm1(i) * Float64(100.0 / i)) * n); elseif (n <= 7.5e-111) tmp = Float64(Float64(n * Float64(Float64(log(i) - log(n)) / i)) * Float64(n * 100.0)); else tmp = Float64(Float64(expm1(i) / i) * Float64(100.0 * n)); end return tmp end
code[i_, n_] := If[LessEqual[n, -5e-311], N[(N[(N[(Exp[i] - 1), $MachinePrecision] * N[(100.0 / i), $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[n, 7.5e-111], N[(N[(n * N[(N[(N[Log[i], $MachinePrecision] - N[Log[n], $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision] * N[(n * 100.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision] * N[(100.0 * n), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -5 \cdot 10^{-311}:\\
\;\;\;\;\left(\mathsf{expm1}\left(i\right) \cdot \frac{100}{i}\right) \cdot n\\
\mathbf{elif}\;n \leq 7.5 \cdot 10^{-111}:\\
\;\;\;\;\left(n \cdot \frac{\log i - \log n}{i}\right) \cdot \left(n \cdot 100\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{expm1}\left(i\right)}{i} \cdot \left(100 \cdot n\right)\\
\end{array}
\end{array}
if n < -5.00000000000023e-311Initial program 29.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f64N/A
lower-/.f6476.6
Applied rewrites76.6%
Taylor expanded in n around inf
lower-expm1.f6477.5
Applied rewrites77.5%
if -5.00000000000023e-311 < n < 7.49999999999999965e-111Initial program 34.8%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
associate-*l*N/A
lower-*.f64N/A
lower-/.f64N/A
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f64N/A
lower-*.f6472.7
Applied rewrites72.7%
Taylor expanded in n around 0
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f64N/A
lower-log.f64N/A
lower-log.f6482.5
Applied rewrites82.5%
if 7.49999999999999965e-111 < n Initial program 14.8%
Taylor expanded in n around inf
lower-expm1.f6468.9
Applied rewrites68.9%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
associate-*l*N/A
lift-*.f64N/A
lower-*.f64N/A
lower-/.f6488.4
lift-*.f64N/A
*-commutativeN/A
lower-*.f6488.4
Applied rewrites88.4%
Final simplification83.0%
(FPCore (i n)
:precision binary64
(if (<= n -5e-311)
(* (* (expm1 i) (/ 100.0 i)) n)
(if (<= n 7.5e-111)
(* (* (* n (/ (- (log i) (log n)) i)) 100.0) n)
(* (/ (expm1 i) i) (* 100.0 n)))))
double code(double i, double n) {
double tmp;
if (n <= -5e-311) {
tmp = (expm1(i) * (100.0 / i)) * n;
} else if (n <= 7.5e-111) {
tmp = ((n * ((log(i) - log(n)) / i)) * 100.0) * n;
} else {
tmp = (expm1(i) / i) * (100.0 * n);
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if (n <= -5e-311) {
tmp = (Math.expm1(i) * (100.0 / i)) * n;
} else if (n <= 7.5e-111) {
tmp = ((n * ((Math.log(i) - Math.log(n)) / i)) * 100.0) * n;
} else {
tmp = (Math.expm1(i) / i) * (100.0 * n);
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -5e-311: tmp = (math.expm1(i) * (100.0 / i)) * n elif n <= 7.5e-111: tmp = ((n * ((math.log(i) - math.log(n)) / i)) * 100.0) * n else: tmp = (math.expm1(i) / i) * (100.0 * n) return tmp
function code(i, n) tmp = 0.0 if (n <= -5e-311) tmp = Float64(Float64(expm1(i) * Float64(100.0 / i)) * n); elseif (n <= 7.5e-111) tmp = Float64(Float64(Float64(n * Float64(Float64(log(i) - log(n)) / i)) * 100.0) * n); else tmp = Float64(Float64(expm1(i) / i) * Float64(100.0 * n)); end return tmp end
code[i_, n_] := If[LessEqual[n, -5e-311], N[(N[(N[(Exp[i] - 1), $MachinePrecision] * N[(100.0 / i), $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[n, 7.5e-111], N[(N[(N[(n * N[(N[(N[Log[i], $MachinePrecision] - N[Log[n], $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], N[(N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision] * N[(100.0 * n), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -5 \cdot 10^{-311}:\\
\;\;\;\;\left(\mathsf{expm1}\left(i\right) \cdot \frac{100}{i}\right) \cdot n\\
\mathbf{elif}\;n \leq 7.5 \cdot 10^{-111}:\\
\;\;\;\;\left(\left(n \cdot \frac{\log i - \log n}{i}\right) \cdot 100\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{expm1}\left(i\right)}{i} \cdot \left(100 \cdot n\right)\\
\end{array}
\end{array}
if n < -5.00000000000023e-311Initial program 29.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f64N/A
lower-/.f6476.6
Applied rewrites76.6%
Taylor expanded in n around inf
lower-expm1.f6477.5
Applied rewrites77.5%
if -5.00000000000023e-311 < n < 7.49999999999999965e-111Initial program 34.8%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f64N/A
lower-/.f6472.7
Applied rewrites72.7%
Taylor expanded in n around 0
*-commutativeN/A
lower-*.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f64N/A
lower-log.f64N/A
lower-log.f6482.5
Applied rewrites82.5%
if 7.49999999999999965e-111 < n Initial program 14.8%
Taylor expanded in n around inf
lower-expm1.f6468.9
Applied rewrites68.9%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
associate-*l*N/A
lift-*.f64N/A
lower-*.f64N/A
lower-/.f6488.4
lift-*.f64N/A
*-commutativeN/A
lower-*.f6488.4
Applied rewrites88.4%
Final simplification82.9%
(FPCore (i n)
:precision binary64
(if (<= n -5e-311)
(* (* (expm1 i) (/ 100.0 i)) n)
(if (<= n 7.5e-111)
(* (* 100.0 (* n n)) (/ (- (log i) (log n)) i))
(* (/ (expm1 i) i) (* 100.0 n)))))
double code(double i, double n) {
double tmp;
if (n <= -5e-311) {
tmp = (expm1(i) * (100.0 / i)) * n;
} else if (n <= 7.5e-111) {
tmp = (100.0 * (n * n)) * ((log(i) - log(n)) / i);
} else {
tmp = (expm1(i) / i) * (100.0 * n);
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if (n <= -5e-311) {
tmp = (Math.expm1(i) * (100.0 / i)) * n;
} else if (n <= 7.5e-111) {
tmp = (100.0 * (n * n)) * ((Math.log(i) - Math.log(n)) / i);
} else {
tmp = (Math.expm1(i) / i) * (100.0 * n);
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -5e-311: tmp = (math.expm1(i) * (100.0 / i)) * n elif n <= 7.5e-111: tmp = (100.0 * (n * n)) * ((math.log(i) - math.log(n)) / i) else: tmp = (math.expm1(i) / i) * (100.0 * n) return tmp
function code(i, n) tmp = 0.0 if (n <= -5e-311) tmp = Float64(Float64(expm1(i) * Float64(100.0 / i)) * n); elseif (n <= 7.5e-111) tmp = Float64(Float64(100.0 * Float64(n * n)) * Float64(Float64(log(i) - log(n)) / i)); else tmp = Float64(Float64(expm1(i) / i) * Float64(100.0 * n)); end return tmp end
code[i_, n_] := If[LessEqual[n, -5e-311], N[(N[(N[(Exp[i] - 1), $MachinePrecision] * N[(100.0 / i), $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[n, 7.5e-111], N[(N[(100.0 * N[(n * n), $MachinePrecision]), $MachinePrecision] * N[(N[(N[Log[i], $MachinePrecision] - N[Log[n], $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], N[(N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision] * N[(100.0 * n), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -5 \cdot 10^{-311}:\\
\;\;\;\;\left(\mathsf{expm1}\left(i\right) \cdot \frac{100}{i}\right) \cdot n\\
\mathbf{elif}\;n \leq 7.5 \cdot 10^{-111}:\\
\;\;\;\;\left(100 \cdot \left(n \cdot n\right)\right) \cdot \frac{\log i - \log n}{i}\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{expm1}\left(i\right)}{i} \cdot \left(100 \cdot n\right)\\
\end{array}
\end{array}
if n < -5.00000000000023e-311Initial program 29.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f64N/A
lower-/.f6476.6
Applied rewrites76.6%
Taylor expanded in n around inf
lower-expm1.f6477.5
Applied rewrites77.5%
if -5.00000000000023e-311 < n < 7.49999999999999965e-111Initial program 34.8%
Taylor expanded in n around 0
associate-/l*N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f64N/A
lower-log.f64N/A
lower-log.f6471.2
Applied rewrites71.2%
if 7.49999999999999965e-111 < n Initial program 14.8%
Taylor expanded in n around inf
lower-expm1.f6468.9
Applied rewrites68.9%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
associate-*l*N/A
lift-*.f64N/A
lower-*.f64N/A
lower-/.f6488.4
lift-*.f64N/A
*-commutativeN/A
lower-*.f6488.4
Applied rewrites88.4%
Final simplification81.2%
(FPCore (i n) :precision binary64 (if (or (<= n -7e-226) (not (<= n 4.5e-114))) (* (* (/ (expm1 i) i) 100.0) n) (* (/ (- 1.0 1.0) i) (* 100.0 n))))
double code(double i, double n) {
double tmp;
if ((n <= -7e-226) || !(n <= 4.5e-114)) {
tmp = ((expm1(i) / i) * 100.0) * n;
} else {
tmp = ((1.0 - 1.0) / i) * (100.0 * n);
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if ((n <= -7e-226) || !(n <= 4.5e-114)) {
tmp = ((Math.expm1(i) / i) * 100.0) * n;
} else {
tmp = ((1.0 - 1.0) / i) * (100.0 * n);
}
return tmp;
}
def code(i, n): tmp = 0 if (n <= -7e-226) or not (n <= 4.5e-114): tmp = ((math.expm1(i) / i) * 100.0) * n else: tmp = ((1.0 - 1.0) / i) * (100.0 * n) return tmp
function code(i, n) tmp = 0.0 if ((n <= -7e-226) || !(n <= 4.5e-114)) tmp = Float64(Float64(Float64(expm1(i) / i) * 100.0) * n); else tmp = Float64(Float64(Float64(1.0 - 1.0) / i) * Float64(100.0 * n)); end return tmp end
code[i_, n_] := If[Or[LessEqual[n, -7e-226], N[Not[LessEqual[n, 4.5e-114]], $MachinePrecision]], N[(N[(N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], N[(N[(N[(1.0 - 1.0), $MachinePrecision] / i), $MachinePrecision] * N[(100.0 * n), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -7 \cdot 10^{-226} \lor \neg \left(n \leq 4.5 \cdot 10^{-114}\right):\\
\;\;\;\;\left(\frac{\mathsf{expm1}\left(i\right)}{i} \cdot 100\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;\frac{1 - 1}{i} \cdot \left(100 \cdot n\right)\\
\end{array}
\end{array}
if n < -7e-226 or 4.49999999999999969e-114 < n Initial program 19.3%
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.f6483.3
Applied rewrites83.3%
if -7e-226 < n < 4.49999999999999969e-114Initial program 44.5%
Taylor expanded in i around 0
Applied rewrites67.9%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
associate-*l*N/A
lift-*.f64N/A
lower-*.f64N/A
lower-/.f6467.9
lift-*.f64N/A
*-commutativeN/A
lower-*.f6467.9
Applied rewrites67.9%
Final simplification80.5%
(FPCore (i n)
:precision binary64
(if (<= n -7.5e-226)
(* (* (expm1 i) (/ 100.0 i)) n)
(if (<= n 4.5e-114)
(* (/ (- 1.0 1.0) i) (* 100.0 n))
(* (/ (expm1 i) i) (* 100.0 n)))))
double code(double i, double n) {
double tmp;
if (n <= -7.5e-226) {
tmp = (expm1(i) * (100.0 / i)) * n;
} else if (n <= 4.5e-114) {
tmp = ((1.0 - 1.0) / i) * (100.0 * n);
} else {
tmp = (expm1(i) / i) * (100.0 * n);
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if (n <= -7.5e-226) {
tmp = (Math.expm1(i) * (100.0 / i)) * n;
} else if (n <= 4.5e-114) {
tmp = ((1.0 - 1.0) / i) * (100.0 * n);
} else {
tmp = (Math.expm1(i) / i) * (100.0 * n);
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -7.5e-226: tmp = (math.expm1(i) * (100.0 / i)) * n elif n <= 4.5e-114: tmp = ((1.0 - 1.0) / i) * (100.0 * n) else: tmp = (math.expm1(i) / i) * (100.0 * n) return tmp
function code(i, n) tmp = 0.0 if (n <= -7.5e-226) tmp = Float64(Float64(expm1(i) * Float64(100.0 / i)) * n); elseif (n <= 4.5e-114) tmp = Float64(Float64(Float64(1.0 - 1.0) / i) * Float64(100.0 * n)); else tmp = Float64(Float64(expm1(i) / i) * Float64(100.0 * n)); end return tmp end
code[i_, n_] := If[LessEqual[n, -7.5e-226], N[(N[(N[(Exp[i] - 1), $MachinePrecision] * N[(100.0 / i), $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[n, 4.5e-114], N[(N[(N[(1.0 - 1.0), $MachinePrecision] / i), $MachinePrecision] * N[(100.0 * n), $MachinePrecision]), $MachinePrecision], N[(N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision] * N[(100.0 * n), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -7.5 \cdot 10^{-226}:\\
\;\;\;\;\left(\mathsf{expm1}\left(i\right) \cdot \frac{100}{i}\right) \cdot n\\
\mathbf{elif}\;n \leq 4.5 \cdot 10^{-114}:\\
\;\;\;\;\frac{1 - 1}{i} \cdot \left(100 \cdot n\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{expm1}\left(i\right)}{i} \cdot \left(100 \cdot n\right)\\
\end{array}
\end{array}
if n < -7.50000000000000044e-226Initial program 24.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f64N/A
lower-/.f6474.9
Applied rewrites74.9%
Taylor expanded in n around inf
lower-expm1.f6477.8
Applied rewrites77.8%
if -7.50000000000000044e-226 < n < 4.49999999999999969e-114Initial program 44.5%
Taylor expanded in i around 0
Applied rewrites67.9%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
associate-*l*N/A
lift-*.f64N/A
lower-*.f64N/A
lower-/.f6467.9
lift-*.f64N/A
*-commutativeN/A
lower-*.f6467.9
Applied rewrites67.9%
if 4.49999999999999969e-114 < n Initial program 14.8%
Taylor expanded in n around inf
lower-expm1.f6468.9
Applied rewrites68.9%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
associate-*l*N/A
lift-*.f64N/A
lower-*.f64N/A
lower-/.f6488.4
lift-*.f64N/A
*-commutativeN/A
lower-*.f6488.4
Applied rewrites88.4%
Final simplification80.5%
(FPCore (i n)
:precision binary64
(if (<= n -7.5e-226)
(* (* (expm1 i) (/ 100.0 i)) n)
(if (<= n 4.5e-114)
(* (/ (- 1.0 1.0) i) (* 100.0 n))
(* (* (/ (expm1 i) i) 100.0) n))))
double code(double i, double n) {
double tmp;
if (n <= -7.5e-226) {
tmp = (expm1(i) * (100.0 / i)) * n;
} else if (n <= 4.5e-114) {
tmp = ((1.0 - 1.0) / i) * (100.0 * n);
} else {
tmp = ((expm1(i) / i) * 100.0) * n;
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if (n <= -7.5e-226) {
tmp = (Math.expm1(i) * (100.0 / i)) * n;
} else if (n <= 4.5e-114) {
tmp = ((1.0 - 1.0) / i) * (100.0 * n);
} else {
tmp = ((Math.expm1(i) / i) * 100.0) * n;
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -7.5e-226: tmp = (math.expm1(i) * (100.0 / i)) * n elif n <= 4.5e-114: tmp = ((1.0 - 1.0) / i) * (100.0 * n) else: tmp = ((math.expm1(i) / i) * 100.0) * n return tmp
function code(i, n) tmp = 0.0 if (n <= -7.5e-226) tmp = Float64(Float64(expm1(i) * Float64(100.0 / i)) * n); elseif (n <= 4.5e-114) tmp = Float64(Float64(Float64(1.0 - 1.0) / i) * Float64(100.0 * n)); else tmp = Float64(Float64(Float64(expm1(i) / i) * 100.0) * n); end return tmp end
code[i_, n_] := If[LessEqual[n, -7.5e-226], N[(N[(N[(Exp[i] - 1), $MachinePrecision] * N[(100.0 / i), $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[n, 4.5e-114], N[(N[(N[(1.0 - 1.0), $MachinePrecision] / i), $MachinePrecision] * N[(100.0 * n), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -7.5 \cdot 10^{-226}:\\
\;\;\;\;\left(\mathsf{expm1}\left(i\right) \cdot \frac{100}{i}\right) \cdot n\\
\mathbf{elif}\;n \leq 4.5 \cdot 10^{-114}:\\
\;\;\;\;\frac{1 - 1}{i} \cdot \left(100 \cdot n\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{\mathsf{expm1}\left(i\right)}{i} \cdot 100\right) \cdot n\\
\end{array}
\end{array}
if n < -7.50000000000000044e-226Initial program 24.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f64N/A
lower-/.f6474.9
Applied rewrites74.9%
Taylor expanded in n around inf
lower-expm1.f6477.8
Applied rewrites77.8%
if -7.50000000000000044e-226 < n < 4.49999999999999969e-114Initial program 44.5%
Taylor expanded in i around 0
Applied rewrites67.9%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
associate-*l*N/A
lift-*.f64N/A
lower-*.f64N/A
lower-/.f6467.9
lift-*.f64N/A
*-commutativeN/A
lower-*.f6467.9
Applied rewrites67.9%
if 4.49999999999999969e-114 < n Initial program 14.8%
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.f6488.3
Applied rewrites88.3%
Final simplification80.5%
(FPCore (i n)
:precision binary64
(if (or (<= n -1.05e-183) (not (<= n 4.5e-114)))
(*
100.0
(*
(fma (fma (fma 0.041666666666666664 i 0.16666666666666666) i 0.5) i 1.0)
n))
(* (/ (- 1.0 1.0) i) (* 100.0 n))))
double code(double i, double n) {
double tmp;
if ((n <= -1.05e-183) || !(n <= 4.5e-114)) {
tmp = 100.0 * (fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * n);
} else {
tmp = ((1.0 - 1.0) / i) * (100.0 * n);
}
return tmp;
}
function code(i, n) tmp = 0.0 if ((n <= -1.05e-183) || !(n <= 4.5e-114)) tmp = Float64(100.0 * Float64(fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * n)); else tmp = Float64(Float64(Float64(1.0 - 1.0) / i) * Float64(100.0 * n)); end return tmp end
code[i_, n_] := If[Or[LessEqual[n, -1.05e-183], N[Not[LessEqual[n, 4.5e-114]], $MachinePrecision]], N[(100.0 * N[(N[(N[(N[(0.041666666666666664 * i + 0.16666666666666666), $MachinePrecision] * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * n), $MachinePrecision]), $MachinePrecision], N[(N[(N[(1.0 - 1.0), $MachinePrecision] / i), $MachinePrecision] * N[(100.0 * n), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.05 \cdot 10^{-183} \lor \neg \left(n \leq 4.5 \cdot 10^{-114}\right):\\
\;\;\;\;100 \cdot \left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, i, 0.16666666666666666\right), i, 0.5\right), i, 1\right) \cdot n\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1 - 1}{i} \cdot \left(100 \cdot n\right)\\
\end{array}
\end{array}
if n < -1.0500000000000001e-183 or 4.49999999999999969e-114 < n Initial program 18.7%
Taylor expanded in n around inf
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6483.9
Applied rewrites83.9%
Taylor expanded in i around 0
Applied rewrites66.6%
if -1.0500000000000001e-183 < n < 4.49999999999999969e-114Initial program 44.4%
Taylor expanded in i around 0
Applied rewrites65.6%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
associate-*l*N/A
lift-*.f64N/A
lower-*.f64N/A
lower-/.f6465.6
lift-*.f64N/A
*-commutativeN/A
lower-*.f6465.6
Applied rewrites65.6%
Final simplification66.4%
(FPCore (i n) :precision binary64 (if (or (<= n -1.05e-183) (not (<= n 4.5e-114))) (* 100.0 (* (fma (fma 0.16666666666666666 i 0.5) i 1.0) n)) (* (/ (- 1.0 1.0) i) (* 100.0 n))))
double code(double i, double n) {
double tmp;
if ((n <= -1.05e-183) || !(n <= 4.5e-114)) {
tmp = 100.0 * (fma(fma(0.16666666666666666, i, 0.5), i, 1.0) * n);
} else {
tmp = ((1.0 - 1.0) / i) * (100.0 * n);
}
return tmp;
}
function code(i, n) tmp = 0.0 if ((n <= -1.05e-183) || !(n <= 4.5e-114)) tmp = Float64(100.0 * Float64(fma(fma(0.16666666666666666, i, 0.5), i, 1.0) * n)); else tmp = Float64(Float64(Float64(1.0 - 1.0) / i) * Float64(100.0 * n)); end return tmp end
code[i_, n_] := If[Or[LessEqual[n, -1.05e-183], N[Not[LessEqual[n, 4.5e-114]], $MachinePrecision]], N[(100.0 * N[(N[(N[(0.16666666666666666 * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * n), $MachinePrecision]), $MachinePrecision], N[(N[(N[(1.0 - 1.0), $MachinePrecision] / i), $MachinePrecision] * N[(100.0 * n), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.05 \cdot 10^{-183} \lor \neg \left(n \leq 4.5 \cdot 10^{-114}\right):\\
\;\;\;\;100 \cdot \left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, i, 0.5\right), i, 1\right) \cdot n\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1 - 1}{i} \cdot \left(100 \cdot n\right)\\
\end{array}
\end{array}
if n < -1.0500000000000001e-183 or 4.49999999999999969e-114 < n Initial program 18.7%
Taylor expanded in n around inf
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6483.9
Applied rewrites83.9%
Taylor expanded in i around 0
Applied rewrites64.3%
if -1.0500000000000001e-183 < n < 4.49999999999999969e-114Initial program 44.4%
Taylor expanded in i around 0
Applied rewrites65.6%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
associate-*l*N/A
lift-*.f64N/A
lower-*.f64N/A
lower-/.f6465.6
lift-*.f64N/A
*-commutativeN/A
lower-*.f6465.6
Applied rewrites65.6%
Final simplification64.6%
(FPCore (i n)
:precision binary64
(if (<= n -1.05e-183)
(* 100.0 (fma (fma (* n i) 0.16666666666666666 (* 0.5 n)) i n))
(if (<= n 4.5e-114)
(* (/ (- 1.0 1.0) i) (* 100.0 n))
(* 100.0 (* (fma (fma 0.16666666666666666 i 0.5) i 1.0) n)))))
double code(double i, double n) {
double tmp;
if (n <= -1.05e-183) {
tmp = 100.0 * fma(fma((n * i), 0.16666666666666666, (0.5 * n)), i, n);
} else if (n <= 4.5e-114) {
tmp = ((1.0 - 1.0) / i) * (100.0 * n);
} else {
tmp = 100.0 * (fma(fma(0.16666666666666666, i, 0.5), i, 1.0) * n);
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -1.05e-183) tmp = Float64(100.0 * fma(fma(Float64(n * i), 0.16666666666666666, Float64(0.5 * n)), i, n)); elseif (n <= 4.5e-114) tmp = Float64(Float64(Float64(1.0 - 1.0) / i) * Float64(100.0 * n)); else tmp = Float64(100.0 * Float64(fma(fma(0.16666666666666666, i, 0.5), i, 1.0) * n)); end return tmp end
code[i_, n_] := If[LessEqual[n, -1.05e-183], N[(100.0 * N[(N[(N[(n * i), $MachinePrecision] * 0.16666666666666666 + N[(0.5 * n), $MachinePrecision]), $MachinePrecision] * i + n), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 4.5e-114], N[(N[(N[(1.0 - 1.0), $MachinePrecision] / i), $MachinePrecision] * N[(100.0 * n), $MachinePrecision]), $MachinePrecision], N[(100.0 * N[(N[(N[(0.16666666666666666 * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * n), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.05 \cdot 10^{-183}:\\
\;\;\;\;100 \cdot \mathsf{fma}\left(\mathsf{fma}\left(n \cdot i, 0.16666666666666666, 0.5 \cdot n\right), i, n\right)\\
\mathbf{elif}\;n \leq 4.5 \cdot 10^{-114}:\\
\;\;\;\;\frac{1 - 1}{i} \cdot \left(100 \cdot n\right)\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, i, 0.5\right), i, 1\right) \cdot n\right)\\
\end{array}
\end{array}
if n < -1.0500000000000001e-183Initial program 23.4%
Taylor expanded in n around inf
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6478.8
Applied rewrites78.8%
Taylor expanded in i around 0
Applied rewrites60.2%
Taylor expanded in i around 0
Applied rewrites60.2%
if -1.0500000000000001e-183 < n < 4.49999999999999969e-114Initial program 44.4%
Taylor expanded in i around 0
Applied rewrites65.6%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
associate-*l*N/A
lift-*.f64N/A
lower-*.f64N/A
lower-/.f6465.6
lift-*.f64N/A
*-commutativeN/A
lower-*.f6465.6
Applied rewrites65.6%
if 4.49999999999999969e-114 < n Initial program 14.8%
Taylor expanded in n around inf
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6488.3
Applied rewrites88.3%
Taylor expanded in i around 0
Applied rewrites67.9%
(FPCore (i n) :precision binary64 (if (<= n -2.65e-226) (* 100.0 (* (fma (* 0.16666666666666666 i) i 1.0) n)) (if (<= n 2e-235) (* (/ 100.0 i) (* n i)) (* 100.0 (fma (* n i) 0.5 n)))))
double code(double i, double n) {
double tmp;
if (n <= -2.65e-226) {
tmp = 100.0 * (fma((0.16666666666666666 * i), i, 1.0) * n);
} else if (n <= 2e-235) {
tmp = (100.0 / i) * (n * i);
} else {
tmp = 100.0 * fma((n * i), 0.5, n);
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -2.65e-226) tmp = Float64(100.0 * Float64(fma(Float64(0.16666666666666666 * i), i, 1.0) * n)); elseif (n <= 2e-235) tmp = Float64(Float64(100.0 / i) * Float64(n * i)); else tmp = Float64(100.0 * fma(Float64(n * i), 0.5, n)); end return tmp end
code[i_, n_] := If[LessEqual[n, -2.65e-226], N[(100.0 * N[(N[(N[(0.16666666666666666 * i), $MachinePrecision] * i + 1.0), $MachinePrecision] * n), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 2e-235], N[(N[(100.0 / i), $MachinePrecision] * N[(n * i), $MachinePrecision]), $MachinePrecision], N[(100.0 * N[(N[(n * i), $MachinePrecision] * 0.5 + n), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -2.65 \cdot 10^{-226}:\\
\;\;\;\;100 \cdot \left(\mathsf{fma}\left(0.16666666666666666 \cdot i, i, 1\right) \cdot n\right)\\
\mathbf{elif}\;n \leq 2 \cdot 10^{-235}:\\
\;\;\;\;\frac{100}{i} \cdot \left(n \cdot i\right)\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \mathsf{fma}\left(n \cdot i, 0.5, n\right)\\
\end{array}
\end{array}
if n < -2.6500000000000002e-226Initial program 24.4%
Taylor expanded in n around inf
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6477.7
Applied rewrites77.7%
Taylor expanded in i around 0
Applied rewrites59.0%
Taylor expanded in i around inf
Applied rewrites58.3%
if -2.6500000000000002e-226 < n < 1.9999999999999999e-235Initial program 72.9%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f64N/A
lower-/.f6481.2
Applied rewrites81.2%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f643.3
Applied rewrites3.3%
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
lower-*.f6418.9
Applied rewrites18.9%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f6443.5
Applied rewrites43.5%
if 1.9999999999999999e-235 < n Initial program 16.1%
Taylor expanded in n around inf
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6477.1
Applied rewrites77.1%
Taylor expanded in i around 0
Applied rewrites60.7%
Taylor expanded in i around 0
Applied rewrites60.2%
(FPCore (i n) :precision binary64 (* 100.0 (* (fma (fma 0.16666666666666666 i 0.5) i 1.0) n)))
double code(double i, double n) {
return 100.0 * (fma(fma(0.16666666666666666, i, 0.5), i, 1.0) * n);
}
function code(i, n) return Float64(100.0 * Float64(fma(fma(0.16666666666666666, i, 0.5), i, 1.0) * n)) end
code[i_, n_] := N[(100.0 * N[(N[(N[(0.16666666666666666 * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * n), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
100 \cdot \left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, i, 0.5\right), i, 1\right) \cdot n\right)
\end{array}
Initial program 23.9%
Taylor expanded in n around inf
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6473.4
Applied rewrites73.4%
Taylor expanded in i around 0
Applied rewrites55.5%
(FPCore (i n) :precision binary64 (* 100.0 (* (fma (* 0.16666666666666666 i) i 1.0) n)))
double code(double i, double n) {
return 100.0 * (fma((0.16666666666666666 * i), i, 1.0) * n);
}
function code(i, n) return Float64(100.0 * Float64(fma(Float64(0.16666666666666666 * i), i, 1.0) * n)) end
code[i_, n_] := N[(100.0 * N[(N[(N[(0.16666666666666666 * i), $MachinePrecision] * i + 1.0), $MachinePrecision] * n), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
100 \cdot \left(\mathsf{fma}\left(0.16666666666666666 \cdot i, i, 1\right) \cdot n\right)
\end{array}
Initial program 23.9%
Taylor expanded in n around inf
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6473.4
Applied rewrites73.4%
Taylor expanded in i around 0
Applied rewrites55.5%
Taylor expanded in i around inf
Applied rewrites54.9%
(FPCore (i n) :precision binary64 (* 100.0 (fma (* n i) 0.5 n)))
double code(double i, double n) {
return 100.0 * fma((n * i), 0.5, n);
}
function code(i, n) return Float64(100.0 * fma(Float64(n * i), 0.5, n)) end
code[i_, n_] := N[(100.0 * N[(N[(n * i), $MachinePrecision] * 0.5 + n), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
100 \cdot \mathsf{fma}\left(n \cdot i, 0.5, n\right)
\end{array}
Initial program 23.9%
Taylor expanded in n around inf
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6473.4
Applied rewrites73.4%
Taylor expanded in i around 0
Applied rewrites55.5%
Taylor expanded in i around 0
Applied rewrites54.7%
(FPCore (i n) :precision binary64 (* (fma 50.0 i 100.0) n))
double code(double i, double n) {
return fma(50.0, i, 100.0) * n;
}
function code(i, n) return Float64(fma(50.0, i, 100.0) * n) end
code[i_, n_] := N[(N[(50.0 * i + 100.0), $MachinePrecision] * n), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(50, i, 100\right) \cdot n
\end{array}
Initial program 23.9%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f64N/A
lower-/.f6474.6
Applied rewrites74.6%
Taylor expanded in i around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6454.4
Applied rewrites54.4%
Taylor expanded in n around inf
Applied rewrites54.7%
(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 23.9%
Taylor expanded in i around 0
lower-*.f6446.6
Applied rewrites46.6%
(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 2024329
(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))))