
(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 0.0)
(* n (* (expm1 (* n (log1p (/ i n)))) (/ 100.0 i)))
(if (<= t_0 INFINITY) (* t_0 100.0) (* n 100.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 <= 0.0) {
tmp = n * (expm1((n * log1p((i / n)))) * (100.0 / i));
} else if (t_0 <= ((double) INFINITY)) {
tmp = t_0 * 100.0;
} else {
tmp = n * 100.0;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = (Math.pow((1.0 + (i / n)), n) + -1.0) / (i / n);
double tmp;
if (t_0 <= 0.0) {
tmp = n * (Math.expm1((n * Math.log1p((i / n)))) * (100.0 / i));
} else if (t_0 <= Double.POSITIVE_INFINITY) {
tmp = t_0 * 100.0;
} else {
tmp = n * 100.0;
}
return tmp;
}
def code(i, n): t_0 = (math.pow((1.0 + (i / n)), n) + -1.0) / (i / n) tmp = 0 if t_0 <= 0.0: tmp = n * (math.expm1((n * math.log1p((i / n)))) * (100.0 / i)) elif t_0 <= math.inf: tmp = t_0 * 100.0 else: tmp = n * 100.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 <= 0.0) tmp = Float64(n * Float64(expm1(Float64(n * log1p(Float64(i / n)))) * Float64(100.0 / i))); elseif (t_0 <= Inf) tmp = Float64(t_0 * 100.0); else tmp = Float64(n * 100.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, 0.0], N[(n * N[(N[(Exp[N[(n * N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]] - 1), $MachinePrecision] * N[(100.0 / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, Infinity], N[(t$95$0 * 100.0), $MachinePrecision], N[(n * 100.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 0:\\
\;\;\;\;n \cdot \left(\mathsf{expm1}\left(n \cdot \mathsf{log1p}\left(\frac{i}{n}\right)\right) \cdot \frac{100}{i}\right)\\
\mathbf{elif}\;t\_0 \leq \infty:\\
\;\;\;\;t\_0 \cdot 100\\
\mathbf{else}:\\
\;\;\;\;n \cdot 100\\
\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 20.3%
associate-/r/20.0%
associate-*r*20.0%
*-commutative20.0%
associate-*r/20.0%
sub-neg20.0%
distribute-lft-in20.0%
metadata-eval20.0%
metadata-eval20.0%
metadata-eval20.0%
fma-define20.0%
metadata-eval20.0%
Simplified20.0%
fma-undefine20.0%
metadata-eval20.0%
metadata-eval20.0%
distribute-lft-in20.0%
sub-neg20.0%
*-commutative20.0%
add-exp-log20.0%
expm1-define20.0%
log-pow33.7%
log1p-define98.3%
Applied egg-rr98.3%
associate-/l*98.3%
Applied egg-rr98.3%
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 92.8%
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%
associate-/r/1.9%
associate-*r*1.9%
*-commutative1.9%
associate-*r/1.9%
sub-neg1.9%
distribute-lft-in1.9%
metadata-eval1.9%
metadata-eval1.9%
metadata-eval1.9%
fma-define1.9%
metadata-eval1.9%
Simplified1.9%
Taylor expanded in i around 0 71.5%
*-commutative71.5%
Simplified71.5%
Final simplification93.0%
(FPCore (i n)
:precision binary64
(if (<= n -4.8e-248)
(* n (* (/ 100.0 i) (expm1 i)))
(if (<= n 1.05e-240)
(/ 0.0 (/ i n))
(if (<= n 1.05e+32)
(* 100.0 (/ i (/ i n)))
(* n (/ (* 100.0 (expm1 i)) i))))))
double code(double i, double n) {
double tmp;
if (n <= -4.8e-248) {
tmp = n * ((100.0 / i) * expm1(i));
} else if (n <= 1.05e-240) {
tmp = 0.0 / (i / n);
} else if (n <= 1.05e+32) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = n * ((100.0 * expm1(i)) / i);
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if (n <= -4.8e-248) {
tmp = n * ((100.0 / i) * Math.expm1(i));
} else if (n <= 1.05e-240) {
tmp = 0.0 / (i / n);
} else if (n <= 1.05e+32) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = n * ((100.0 * Math.expm1(i)) / i);
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -4.8e-248: tmp = n * ((100.0 / i) * math.expm1(i)) elif n <= 1.05e-240: tmp = 0.0 / (i / n) elif n <= 1.05e+32: tmp = 100.0 * (i / (i / n)) else: tmp = n * ((100.0 * math.expm1(i)) / i) return tmp
function code(i, n) tmp = 0.0 if (n <= -4.8e-248) tmp = Float64(n * Float64(Float64(100.0 / i) * expm1(i))); elseif (n <= 1.05e-240) tmp = Float64(0.0 / Float64(i / n)); elseif (n <= 1.05e+32) tmp = Float64(100.0 * Float64(i / Float64(i / n))); else tmp = Float64(n * Float64(Float64(100.0 * expm1(i)) / i)); end return tmp end
code[i_, n_] := If[LessEqual[n, -4.8e-248], N[(n * N[(N[(100.0 / i), $MachinePrecision] * N[(Exp[i] - 1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.05e-240], N[(0.0 / N[(i / n), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.05e+32], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(n * N[(N[(100.0 * N[(Exp[i] - 1), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -4.8 \cdot 10^{-248}:\\
\;\;\;\;n \cdot \left(\frac{100}{i} \cdot \mathsf{expm1}\left(i\right)\right)\\
\mathbf{elif}\;n \leq 1.05 \cdot 10^{-240}:\\
\;\;\;\;\frac{0}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 1.05 \cdot 10^{+32}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;n \cdot \frac{100 \cdot \mathsf{expm1}\left(i\right)}{i}\\
\end{array}
\end{array}
if n < -4.80000000000000006e-248Initial program 25.0%
associate-/r/25.2%
associate-*r*25.3%
*-commutative25.3%
associate-*r/25.3%
sub-neg25.3%
distribute-lft-in25.3%
metadata-eval25.3%
metadata-eval25.3%
metadata-eval25.3%
fma-define25.3%
metadata-eval25.3%
Simplified25.3%
fma-undefine25.3%
metadata-eval25.3%
metadata-eval25.3%
distribute-lft-in25.3%
sub-neg25.3%
*-commutative25.3%
add-exp-log25.3%
expm1-define25.3%
log-pow24.1%
log1p-define77.9%
Applied egg-rr77.9%
associate-/l*78.0%
Applied egg-rr78.0%
Taylor expanded in n around inf 32.3%
expm1-define83.1%
Simplified83.1%
if -4.80000000000000006e-248 < n < 1.04999999999999997e-240Initial program 51.1%
associate-*r/51.1%
sub-neg51.1%
distribute-rgt-in51.1%
metadata-eval51.1%
metadata-eval51.1%
Simplified51.1%
Taylor expanded in i around 0 86.1%
if 1.04999999999999997e-240 < n < 1.05e32Initial program 13.9%
Taylor expanded in i around 0 56.7%
if 1.05e32 < n Initial program 17.1%
associate-/r/17.6%
associate-*r*17.6%
*-commutative17.6%
associate-*r/17.6%
sub-neg17.6%
distribute-lft-in17.6%
metadata-eval17.6%
metadata-eval17.6%
metadata-eval17.6%
fma-define17.6%
metadata-eval17.6%
Simplified17.6%
Taylor expanded in n around inf 50.5%
sub-neg50.5%
metadata-eval50.5%
metadata-eval50.5%
distribute-lft-in50.4%
metadata-eval50.4%
sub-neg50.4%
expm1-define98.3%
Simplified98.3%
Final simplification81.7%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* n (* (/ 100.0 i) (expm1 i)))))
(if (<= n -4.4e-248)
t_0
(if (<= n 1.06e-240)
(/ 0.0 (/ i n))
(if (<= n 2.8) (* 100.0 (/ i (/ i n))) t_0)))))
double code(double i, double n) {
double t_0 = n * ((100.0 / i) * expm1(i));
double tmp;
if (n <= -4.4e-248) {
tmp = t_0;
} else if (n <= 1.06e-240) {
tmp = 0.0 / (i / n);
} else if (n <= 2.8) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = t_0;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = n * ((100.0 / i) * Math.expm1(i));
double tmp;
if (n <= -4.4e-248) {
tmp = t_0;
} else if (n <= 1.06e-240) {
tmp = 0.0 / (i / n);
} else if (n <= 2.8) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = n * ((100.0 / i) * math.expm1(i)) tmp = 0 if n <= -4.4e-248: tmp = t_0 elif n <= 1.06e-240: tmp = 0.0 / (i / n) elif n <= 2.8: tmp = 100.0 * (i / (i / n)) else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(n * Float64(Float64(100.0 / i) * expm1(i))) tmp = 0.0 if (n <= -4.4e-248) tmp = t_0; elseif (n <= 1.06e-240) tmp = Float64(0.0 / Float64(i / n)); elseif (n <= 2.8) tmp = Float64(100.0 * Float64(i / Float64(i / n))); else tmp = t_0; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(n * N[(N[(100.0 / i), $MachinePrecision] * N[(Exp[i] - 1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -4.4e-248], t$95$0, If[LessEqual[n, 1.06e-240], N[(0.0 / N[(i / n), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 2.8], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := n \cdot \left(\frac{100}{i} \cdot \mathsf{expm1}\left(i\right)\right)\\
\mathbf{if}\;n \leq -4.4 \cdot 10^{-248}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 1.06 \cdot 10^{-240}:\\
\;\;\;\;\frac{0}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 2.8:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -4.39999999999999999e-248 or 2.7999999999999998 < n Initial program 21.8%
associate-/r/22.1%
associate-*r*22.1%
*-commutative22.1%
associate-*r/22.2%
sub-neg22.2%
distribute-lft-in22.1%
metadata-eval22.1%
metadata-eval22.1%
metadata-eval22.1%
fma-define22.2%
metadata-eval22.2%
Simplified22.2%
fma-undefine22.1%
metadata-eval22.1%
metadata-eval22.1%
distribute-lft-in22.2%
sub-neg22.2%
*-commutative22.2%
add-exp-log22.2%
expm1-define22.2%
log-pow20.9%
log1p-define77.2%
Applied egg-rr77.2%
associate-/l*77.3%
Applied egg-rr77.3%
Taylor expanded in n around inf 37.5%
expm1-define88.6%
Simplified88.6%
if -4.39999999999999999e-248 < n < 1.06e-240Initial program 51.1%
associate-*r/51.1%
sub-neg51.1%
distribute-rgt-in51.1%
metadata-eval51.1%
metadata-eval51.1%
Simplified51.1%
Taylor expanded in i around 0 86.1%
if 1.06e-240 < n < 2.7999999999999998Initial program 15.1%
Taylor expanded in i around 0 52.1%
Final simplification81.7%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* 100.0 (/ (expm1 i) (/ i n)))))
(if (<= i -1e-13)
t_0
(if (<= i 9.2e-14)
(* n (+ 100.0 (* 100.0 (* i (- 0.5 (/ 0.5 n))))))
(if (<= i 1.26e+150) t_0 (/ 0.0 (/ i n)))))))
double code(double i, double n) {
double t_0 = 100.0 * (expm1(i) / (i / n));
double tmp;
if (i <= -1e-13) {
tmp = t_0;
} else if (i <= 9.2e-14) {
tmp = n * (100.0 + (100.0 * (i * (0.5 - (0.5 / n)))));
} else if (i <= 1.26e+150) {
tmp = t_0;
} else {
tmp = 0.0 / (i / n);
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = 100.0 * (Math.expm1(i) / (i / n));
double tmp;
if (i <= -1e-13) {
tmp = t_0;
} else if (i <= 9.2e-14) {
tmp = n * (100.0 + (100.0 * (i * (0.5 - (0.5 / n)))));
} else if (i <= 1.26e+150) {
tmp = t_0;
} else {
tmp = 0.0 / (i / n);
}
return tmp;
}
def code(i, n): t_0 = 100.0 * (math.expm1(i) / (i / n)) tmp = 0 if i <= -1e-13: tmp = t_0 elif i <= 9.2e-14: tmp = n * (100.0 + (100.0 * (i * (0.5 - (0.5 / n))))) elif i <= 1.26e+150: tmp = t_0 else: tmp = 0.0 / (i / n) return tmp
function code(i, n) t_0 = Float64(100.0 * Float64(expm1(i) / Float64(i / n))) tmp = 0.0 if (i <= -1e-13) tmp = t_0; elseif (i <= 9.2e-14) tmp = Float64(n * Float64(100.0 + Float64(100.0 * Float64(i * Float64(0.5 - Float64(0.5 / n)))))); elseif (i <= 1.26e+150) tmp = t_0; else tmp = Float64(0.0 / Float64(i / n)); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(100.0 * N[(N[(Exp[i] - 1), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[i, -1e-13], t$95$0, If[LessEqual[i, 9.2e-14], N[(n * N[(100.0 + N[(100.0 * N[(i * N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[i, 1.26e+150], t$95$0, N[(0.0 / N[(i / n), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 100 \cdot \frac{\mathsf{expm1}\left(i\right)}{\frac{i}{n}}\\
\mathbf{if}\;i \leq -1 \cdot 10^{-13}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;i \leq 9.2 \cdot 10^{-14}:\\
\;\;\;\;n \cdot \left(100 + 100 \cdot \left(i \cdot \left(0.5 - \frac{0.5}{n}\right)\right)\right)\\
\mathbf{elif}\;i \leq 1.26 \cdot 10^{+150}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\frac{0}{\frac{i}{n}}\\
\end{array}
\end{array}
if i < -1e-13 or 9.19999999999999993e-14 < i < 1.26e150Initial program 36.9%
Taylor expanded in n around inf 70.2%
expm1-define70.6%
Simplified70.6%
if -1e-13 < i < 9.19999999999999993e-14Initial program 5.6%
associate-/r/6.1%
associate-*r*6.1%
*-commutative6.1%
associate-*r/6.1%
sub-neg6.1%
distribute-lft-in6.1%
metadata-eval6.1%
metadata-eval6.1%
metadata-eval6.1%
fma-define6.1%
metadata-eval6.1%
Simplified6.1%
Taylor expanded in i around 0 87.9%
*-commutative87.9%
associate-*r/87.9%
metadata-eval87.9%
Simplified87.9%
if 1.26e150 < i Initial program 59.2%
associate-*r/59.3%
sub-neg59.3%
distribute-rgt-in59.3%
metadata-eval59.3%
metadata-eval59.3%
Simplified59.3%
Taylor expanded in i around 0 49.8%
Final simplification77.5%
(FPCore (i n) :precision binary64 (if (<= i 1.15e+148) (* n (* (/ 100.0 i) (expm1 i))) (* 100.0 (+ (pow (/ i n) (+ n -1.0)) (* n (/ -1.0 i))))))
double code(double i, double n) {
double tmp;
if (i <= 1.15e+148) {
tmp = n * ((100.0 / i) * expm1(i));
} else {
tmp = 100.0 * (pow((i / n), (n + -1.0)) + (n * (-1.0 / i)));
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if (i <= 1.15e+148) {
tmp = n * ((100.0 / i) * Math.expm1(i));
} else {
tmp = 100.0 * (Math.pow((i / n), (n + -1.0)) + (n * (-1.0 / i)));
}
return tmp;
}
def code(i, n): tmp = 0 if i <= 1.15e+148: tmp = n * ((100.0 / i) * math.expm1(i)) else: tmp = 100.0 * (math.pow((i / n), (n + -1.0)) + (n * (-1.0 / i))) return tmp
function code(i, n) tmp = 0.0 if (i <= 1.15e+148) tmp = Float64(n * Float64(Float64(100.0 / i) * expm1(i))); else tmp = Float64(100.0 * Float64((Float64(i / n) ^ Float64(n + -1.0)) + Float64(n * Float64(-1.0 / i)))); end return tmp end
code[i_, n_] := If[LessEqual[i, 1.15e+148], N[(n * N[(N[(100.0 / i), $MachinePrecision] * N[(Exp[i] - 1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(100.0 * N[(N[Power[N[(i / n), $MachinePrecision], N[(n + -1.0), $MachinePrecision]], $MachinePrecision] + N[(n * N[(-1.0 / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq 1.15 \cdot 10^{+148}:\\
\;\;\;\;n \cdot \left(\frac{100}{i} \cdot \mathsf{expm1}\left(i\right)\right)\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \left({\left(\frac{i}{n}\right)}^{\left(n + -1\right)} + n \cdot \frac{-1}{i}\right)\\
\end{array}
\end{array}
if i < 1.15e148Initial program 18.6%
associate-/r/18.7%
associate-*r*18.7%
*-commutative18.7%
associate-*r/18.7%
sub-neg18.7%
distribute-lft-in18.6%
metadata-eval18.6%
metadata-eval18.6%
metadata-eval18.6%
fma-define18.7%
metadata-eval18.7%
Simplified18.7%
fma-undefine18.6%
metadata-eval18.6%
metadata-eval18.6%
distribute-lft-in18.7%
sub-neg18.7%
*-commutative18.7%
add-exp-log18.7%
expm1-define18.7%
log-pow30.2%
log1p-define82.7%
Applied egg-rr82.7%
associate-/l*82.8%
Applied egg-rr82.8%
Taylor expanded in n around inf 32.9%
expm1-define80.3%
Simplified80.3%
if 1.15e148 < i Initial program 59.2%
Taylor expanded in i around inf 59.2%
div-sub59.2%
pow159.2%
pow-div96.2%
Applied egg-rr96.2%
sub-neg96.2%
metadata-eval96.2%
associate-/r/91.9%
Simplified91.9%
Final simplification81.5%
(FPCore (i n)
:precision binary64
(if (<= n -2.9e-35)
(* n (+ 100.0 (* 100.0 (* i (- 0.5 (/ 0.5 n))))))
(if (<= n 1.05e+32)
(* 100.0 (/ i (/ i n)))
(* n (/ (* 100.0 (* i (+ 1.0 (* i (+ 0.5 (* 0.5 (/ -1.0 n))))))) i)))))
double code(double i, double n) {
double tmp;
if (n <= -2.9e-35) {
tmp = n * (100.0 + (100.0 * (i * (0.5 - (0.5 / n)))));
} else if (n <= 1.05e+32) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = n * ((100.0 * (i * (1.0 + (i * (0.5 + (0.5 * (-1.0 / n))))))) / i);
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if (n <= (-2.9d-35)) then
tmp = n * (100.0d0 + (100.0d0 * (i * (0.5d0 - (0.5d0 / n)))))
else if (n <= 1.05d+32) then
tmp = 100.0d0 * (i / (i / n))
else
tmp = n * ((100.0d0 * (i * (1.0d0 + (i * (0.5d0 + (0.5d0 * ((-1.0d0) / n))))))) / i)
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (n <= -2.9e-35) {
tmp = n * (100.0 + (100.0 * (i * (0.5 - (0.5 / n)))));
} else if (n <= 1.05e+32) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = n * ((100.0 * (i * (1.0 + (i * (0.5 + (0.5 * (-1.0 / n))))))) / i);
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -2.9e-35: tmp = n * (100.0 + (100.0 * (i * (0.5 - (0.5 / n))))) elif n <= 1.05e+32: tmp = 100.0 * (i / (i / n)) else: tmp = n * ((100.0 * (i * (1.0 + (i * (0.5 + (0.5 * (-1.0 / n))))))) / i) return tmp
function code(i, n) tmp = 0.0 if (n <= -2.9e-35) tmp = Float64(n * Float64(100.0 + Float64(100.0 * Float64(i * Float64(0.5 - Float64(0.5 / n)))))); elseif (n <= 1.05e+32) tmp = Float64(100.0 * Float64(i / Float64(i / n))); else tmp = Float64(n * Float64(Float64(100.0 * Float64(i * Float64(1.0 + Float64(i * Float64(0.5 + Float64(0.5 * Float64(-1.0 / n))))))) / i)); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (n <= -2.9e-35) tmp = n * (100.0 + (100.0 * (i * (0.5 - (0.5 / n))))); elseif (n <= 1.05e+32) tmp = 100.0 * (i / (i / n)); else tmp = n * ((100.0 * (i * (1.0 + (i * (0.5 + (0.5 * (-1.0 / n))))))) / i); end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[n, -2.9e-35], N[(n * N[(100.0 + N[(100.0 * N[(i * N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.05e+32], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(n * N[(N[(100.0 * N[(i * N[(1.0 + N[(i * N[(0.5 + N[(0.5 * N[(-1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -2.9 \cdot 10^{-35}:\\
\;\;\;\;n \cdot \left(100 + 100 \cdot \left(i \cdot \left(0.5 - \frac{0.5}{n}\right)\right)\right)\\
\mathbf{elif}\;n \leq 1.05 \cdot 10^{+32}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;n \cdot \frac{100 \cdot \left(i \cdot \left(1 + i \cdot \left(0.5 + 0.5 \cdot \frac{-1}{n}\right)\right)\right)}{i}\\
\end{array}
\end{array}
if n < -2.9000000000000002e-35Initial program 22.7%
associate-/r/23.0%
associate-*r*23.0%
*-commutative23.0%
associate-*r/23.1%
sub-neg23.1%
distribute-lft-in23.0%
metadata-eval23.0%
metadata-eval23.0%
metadata-eval23.0%
fma-define23.1%
metadata-eval23.1%
Simplified23.1%
Taylor expanded in i around 0 60.3%
*-commutative60.3%
associate-*r/60.3%
metadata-eval60.3%
Simplified60.3%
if -2.9000000000000002e-35 < n < 1.05e32Initial program 26.7%
Taylor expanded in i around 0 61.3%
if 1.05e32 < n Initial program 17.1%
associate-/r/17.6%
associate-*r*17.6%
*-commutative17.6%
associate-*r/17.6%
sub-neg17.6%
distribute-lft-in17.6%
metadata-eval17.6%
metadata-eval17.6%
metadata-eval17.6%
fma-define17.6%
metadata-eval17.6%
Simplified17.6%
fma-undefine17.6%
metadata-eval17.6%
metadata-eval17.6%
distribute-lft-in17.6%
sub-neg17.6%
*-commutative17.6%
add-exp-log17.6%
expm1-define17.6%
log-pow16.0%
log1p-define74.0%
Applied egg-rr74.0%
Taylor expanded in i around 0 68.0%
Final simplification62.6%
(FPCore (i n)
:precision binary64
(let* ((t_0 (+ 100.0 (* 100.0 (* i (- 0.5 (/ 0.5 n)))))))
(if (<= n -3.2e-35)
(* n t_0)
(if (<= n 1.05e+32) (* 100.0 (/ i (/ i n))) (* n (/ (* i t_0) i))))))
double code(double i, double n) {
double t_0 = 100.0 + (100.0 * (i * (0.5 - (0.5 / n))));
double tmp;
if (n <= -3.2e-35) {
tmp = n * t_0;
} else if (n <= 1.05e+32) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = n * ((i * t_0) / i);
}
return tmp;
}
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 = 100.0d0 + (100.0d0 * (i * (0.5d0 - (0.5d0 / n))))
if (n <= (-3.2d-35)) then
tmp = n * t_0
else if (n <= 1.05d+32) then
tmp = 100.0d0 * (i / (i / n))
else
tmp = n * ((i * t_0) / i)
end if
code = tmp
end function
public static double code(double i, double n) {
double t_0 = 100.0 + (100.0 * (i * (0.5 - (0.5 / n))));
double tmp;
if (n <= -3.2e-35) {
tmp = n * t_0;
} else if (n <= 1.05e+32) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = n * ((i * t_0) / i);
}
return tmp;
}
def code(i, n): t_0 = 100.0 + (100.0 * (i * (0.5 - (0.5 / n)))) tmp = 0 if n <= -3.2e-35: tmp = n * t_0 elif n <= 1.05e+32: tmp = 100.0 * (i / (i / n)) else: tmp = n * ((i * t_0) / i) return tmp
function code(i, n) t_0 = Float64(100.0 + Float64(100.0 * Float64(i * Float64(0.5 - Float64(0.5 / n))))) tmp = 0.0 if (n <= -3.2e-35) tmp = Float64(n * t_0); elseif (n <= 1.05e+32) tmp = Float64(100.0 * Float64(i / Float64(i / n))); else tmp = Float64(n * Float64(Float64(i * t_0) / i)); end return tmp end
function tmp_2 = code(i, n) t_0 = 100.0 + (100.0 * (i * (0.5 - (0.5 / n)))); tmp = 0.0; if (n <= -3.2e-35) tmp = n * t_0; elseif (n <= 1.05e+32) tmp = 100.0 * (i / (i / n)); else tmp = n * ((i * t_0) / i); end tmp_2 = tmp; end
code[i_, n_] := Block[{t$95$0 = N[(100.0 + N[(100.0 * N[(i * N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -3.2e-35], N[(n * t$95$0), $MachinePrecision], If[LessEqual[n, 1.05e+32], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(n * N[(N[(i * t$95$0), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 100 + 100 \cdot \left(i \cdot \left(0.5 - \frac{0.5}{n}\right)\right)\\
\mathbf{if}\;n \leq -3.2 \cdot 10^{-35}:\\
\;\;\;\;n \cdot t\_0\\
\mathbf{elif}\;n \leq 1.05 \cdot 10^{+32}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;n \cdot \frac{i \cdot t\_0}{i}\\
\end{array}
\end{array}
if n < -3.1999999999999998e-35Initial program 22.7%
associate-/r/23.0%
associate-*r*23.0%
*-commutative23.0%
associate-*r/23.1%
sub-neg23.1%
distribute-lft-in23.0%
metadata-eval23.0%
metadata-eval23.0%
metadata-eval23.0%
fma-define23.1%
metadata-eval23.1%
Simplified23.1%
Taylor expanded in i around 0 60.3%
*-commutative60.3%
associate-*r/60.3%
metadata-eval60.3%
Simplified60.3%
if -3.1999999999999998e-35 < n < 1.05e32Initial program 26.7%
Taylor expanded in i around 0 61.3%
if 1.05e32 < n Initial program 17.1%
associate-/r/17.6%
associate-*r*17.6%
*-commutative17.6%
associate-*r/17.6%
sub-neg17.6%
distribute-lft-in17.6%
metadata-eval17.6%
metadata-eval17.6%
metadata-eval17.6%
fma-define17.6%
metadata-eval17.6%
Simplified17.6%
Taylor expanded in i around 0 68.0%
*-commutative68.0%
associate-*r/68.0%
metadata-eval68.0%
Simplified68.0%
Final simplification62.6%
(FPCore (i n) :precision binary64 (if (or (<= n -2.3e+28) (not (<= n 0.00048))) (* (* n 100.0) (+ 1.0 (* i 0.5))) (* 100.0 (/ i (/ i n)))))
double code(double i, double n) {
double tmp;
if ((n <= -2.3e+28) || !(n <= 0.00048)) {
tmp = (n * 100.0) * (1.0 + (i * 0.5));
} else {
tmp = 100.0 * (i / (i / n));
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if ((n <= (-2.3d+28)) .or. (.not. (n <= 0.00048d0))) then
tmp = (n * 100.0d0) * (1.0d0 + (i * 0.5d0))
else
tmp = 100.0d0 * (i / (i / n))
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if ((n <= -2.3e+28) || !(n <= 0.00048)) {
tmp = (n * 100.0) * (1.0 + (i * 0.5));
} else {
tmp = 100.0 * (i / (i / n));
}
return tmp;
}
def code(i, n): tmp = 0 if (n <= -2.3e+28) or not (n <= 0.00048): tmp = (n * 100.0) * (1.0 + (i * 0.5)) else: tmp = 100.0 * (i / (i / n)) return tmp
function code(i, n) tmp = 0.0 if ((n <= -2.3e+28) || !(n <= 0.00048)) tmp = Float64(Float64(n * 100.0) * Float64(1.0 + Float64(i * 0.5))); else tmp = Float64(100.0 * Float64(i / Float64(i / n))); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if ((n <= -2.3e+28) || ~((n <= 0.00048))) tmp = (n * 100.0) * (1.0 + (i * 0.5)); else tmp = 100.0 * (i / (i / n)); end tmp_2 = tmp; end
code[i_, n_] := If[Or[LessEqual[n, -2.3e+28], N[Not[LessEqual[n, 0.00048]], $MachinePrecision]], N[(N[(n * 100.0), $MachinePrecision] * N[(1.0 + N[(i * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -2.3 \cdot 10^{+28} \lor \neg \left(n \leq 0.00048\right):\\
\;\;\;\;\left(n \cdot 100\right) \cdot \left(1 + i \cdot 0.5\right)\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\end{array}
\end{array}
if n < -2.29999999999999984e28 or 4.80000000000000012e-4 < n Initial program 18.7%
Taylor expanded in i around 0 14.1%
associate-*r/14.1%
metadata-eval14.1%
Simplified14.1%
Taylor expanded in n around inf 64.1%
associate-*r*64.1%
*-commutative64.1%
*-commutative64.1%
Simplified64.1%
if -2.29999999999999984e28 < n < 4.80000000000000012e-4Initial program 28.3%
Taylor expanded in i around 0 59.7%
Final simplification62.2%
(FPCore (i n) :precision binary64 (if (or (<= n -9.5e+29) (not (<= n 1.5))) (* 100.0 (* n (+ 1.0 (* i 0.5)))) (* 100.0 (/ i (/ i n)))))
double code(double i, double n) {
double tmp;
if ((n <= -9.5e+29) || !(n <= 1.5)) {
tmp = 100.0 * (n * (1.0 + (i * 0.5)));
} else {
tmp = 100.0 * (i / (i / n));
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if ((n <= (-9.5d+29)) .or. (.not. (n <= 1.5d0))) then
tmp = 100.0d0 * (n * (1.0d0 + (i * 0.5d0)))
else
tmp = 100.0d0 * (i / (i / n))
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if ((n <= -9.5e+29) || !(n <= 1.5)) {
tmp = 100.0 * (n * (1.0 + (i * 0.5)));
} else {
tmp = 100.0 * (i / (i / n));
}
return tmp;
}
def code(i, n): tmp = 0 if (n <= -9.5e+29) or not (n <= 1.5): tmp = 100.0 * (n * (1.0 + (i * 0.5))) else: tmp = 100.0 * (i / (i / n)) return tmp
function code(i, n) tmp = 0.0 if ((n <= -9.5e+29) || !(n <= 1.5)) tmp = Float64(100.0 * Float64(n * Float64(1.0 + Float64(i * 0.5)))); else tmp = Float64(100.0 * Float64(i / Float64(i / n))); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if ((n <= -9.5e+29) || ~((n <= 1.5))) tmp = 100.0 * (n * (1.0 + (i * 0.5))); else tmp = 100.0 * (i / (i / n)); end tmp_2 = tmp; end
code[i_, n_] := If[Or[LessEqual[n, -9.5e+29], N[Not[LessEqual[n, 1.5]], $MachinePrecision]], N[(100.0 * N[(n * N[(1.0 + N[(i * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -9.5 \cdot 10^{+29} \lor \neg \left(n \leq 1.5\right):\\
\;\;\;\;100 \cdot \left(n \cdot \left(1 + i \cdot 0.5\right)\right)\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\end{array}
\end{array}
if n < -9.5000000000000003e29 or 1.5 < n Initial program 18.7%
Taylor expanded in i around 0 14.1%
associate-*r/14.1%
metadata-eval14.1%
Simplified14.1%
Taylor expanded in n around inf 64.1%
*-commutative64.1%
Simplified64.1%
if -9.5000000000000003e29 < n < 1.5Initial program 28.3%
Taylor expanded in i around 0 59.7%
Final simplification62.2%
(FPCore (i n) :precision binary64 (if (<= n -9.2e-35) (* n (+ 100.0 (* 100.0 (* i (- 0.5 (/ 0.5 n)))))) (if (<= n 1.5) (* 100.0 (/ i (/ i n))) (* (* n 100.0) (+ 1.0 (* i 0.5))))))
double code(double i, double n) {
double tmp;
if (n <= -9.2e-35) {
tmp = n * (100.0 + (100.0 * (i * (0.5 - (0.5 / n)))));
} else if (n <= 1.5) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = (n * 100.0) * (1.0 + (i * 0.5));
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if (n <= (-9.2d-35)) then
tmp = n * (100.0d0 + (100.0d0 * (i * (0.5d0 - (0.5d0 / n)))))
else if (n <= 1.5d0) then
tmp = 100.0d0 * (i / (i / n))
else
tmp = (n * 100.0d0) * (1.0d0 + (i * 0.5d0))
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (n <= -9.2e-35) {
tmp = n * (100.0 + (100.0 * (i * (0.5 - (0.5 / n)))));
} else if (n <= 1.5) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = (n * 100.0) * (1.0 + (i * 0.5));
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -9.2e-35: tmp = n * (100.0 + (100.0 * (i * (0.5 - (0.5 / n))))) elif n <= 1.5: tmp = 100.0 * (i / (i / n)) else: tmp = (n * 100.0) * (1.0 + (i * 0.5)) return tmp
function code(i, n) tmp = 0.0 if (n <= -9.2e-35) tmp = Float64(n * Float64(100.0 + Float64(100.0 * Float64(i * Float64(0.5 - Float64(0.5 / n)))))); elseif (n <= 1.5) tmp = Float64(100.0 * Float64(i / Float64(i / n))); else tmp = Float64(Float64(n * 100.0) * Float64(1.0 + Float64(i * 0.5))); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (n <= -9.2e-35) tmp = n * (100.0 + (100.0 * (i * (0.5 - (0.5 / n))))); elseif (n <= 1.5) tmp = 100.0 * (i / (i / n)); else tmp = (n * 100.0) * (1.0 + (i * 0.5)); end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[n, -9.2e-35], N[(n * N[(100.0 + N[(100.0 * N[(i * N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.5], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(n * 100.0), $MachinePrecision] * N[(1.0 + N[(i * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -9.2 \cdot 10^{-35}:\\
\;\;\;\;n \cdot \left(100 + 100 \cdot \left(i \cdot \left(0.5 - \frac{0.5}{n}\right)\right)\right)\\
\mathbf{elif}\;n \leq 1.5:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;\left(n \cdot 100\right) \cdot \left(1 + i \cdot 0.5\right)\\
\end{array}
\end{array}
if n < -9.1999999999999996e-35Initial program 22.7%
associate-/r/23.0%
associate-*r*23.0%
*-commutative23.0%
associate-*r/23.1%
sub-neg23.1%
distribute-lft-in23.0%
metadata-eval23.0%
metadata-eval23.0%
metadata-eval23.0%
fma-define23.1%
metadata-eval23.1%
Simplified23.1%
Taylor expanded in i around 0 60.3%
*-commutative60.3%
associate-*r/60.3%
metadata-eval60.3%
Simplified60.3%
if -9.1999999999999996e-35 < n < 1.5Initial program 28.0%
Taylor expanded in i around 0 59.2%
if 1.5 < n Initial program 16.1%
Taylor expanded in i around 0 18.6%
associate-*r/18.6%
metadata-eval18.6%
Simplified18.6%
Taylor expanded in n around inf 69.1%
associate-*r*69.1%
*-commutative69.1%
*-commutative69.1%
Simplified69.1%
Final simplification62.3%
(FPCore (i n) :precision binary64 (if (<= n -1.35e-34) (* 100.0 (+ n (* (- 0.5 (/ 0.5 n)) (* i n)))) (if (<= n 1.5) (* 100.0 (/ i (/ i n))) (* (* n 100.0) (+ 1.0 (* i 0.5))))))
double code(double i, double n) {
double tmp;
if (n <= -1.35e-34) {
tmp = 100.0 * (n + ((0.5 - (0.5 / n)) * (i * n)));
} else if (n <= 1.5) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = (n * 100.0) * (1.0 + (i * 0.5));
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if (n <= (-1.35d-34)) then
tmp = 100.0d0 * (n + ((0.5d0 - (0.5d0 / n)) * (i * n)))
else if (n <= 1.5d0) then
tmp = 100.0d0 * (i / (i / n))
else
tmp = (n * 100.0d0) * (1.0d0 + (i * 0.5d0))
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (n <= -1.35e-34) {
tmp = 100.0 * (n + ((0.5 - (0.5 / n)) * (i * n)));
} else if (n <= 1.5) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = (n * 100.0) * (1.0 + (i * 0.5));
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -1.35e-34: tmp = 100.0 * (n + ((0.5 - (0.5 / n)) * (i * n))) elif n <= 1.5: tmp = 100.0 * (i / (i / n)) else: tmp = (n * 100.0) * (1.0 + (i * 0.5)) return tmp
function code(i, n) tmp = 0.0 if (n <= -1.35e-34) tmp = Float64(100.0 * Float64(n + Float64(Float64(0.5 - Float64(0.5 / n)) * Float64(i * n)))); elseif (n <= 1.5) tmp = Float64(100.0 * Float64(i / Float64(i / n))); else tmp = Float64(Float64(n * 100.0) * Float64(1.0 + Float64(i * 0.5))); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (n <= -1.35e-34) tmp = 100.0 * (n + ((0.5 - (0.5 / n)) * (i * n))); elseif (n <= 1.5) tmp = 100.0 * (i / (i / n)); else tmp = (n * 100.0) * (1.0 + (i * 0.5)); end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[n, -1.35e-34], N[(100.0 * N[(n + N[(N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * N[(i * n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.5], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(n * 100.0), $MachinePrecision] * N[(1.0 + N[(i * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.35 \cdot 10^{-34}:\\
\;\;\;\;100 \cdot \left(n + \left(0.5 - \frac{0.5}{n}\right) \cdot \left(i \cdot n\right)\right)\\
\mathbf{elif}\;n \leq 1.5:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;\left(n \cdot 100\right) \cdot \left(1 + i \cdot 0.5\right)\\
\end{array}
\end{array}
if n < -1.35000000000000008e-34Initial program 22.7%
Taylor expanded in i around 0 60.3%
associate-*r*60.3%
*-commutative60.3%
associate-*r/60.3%
metadata-eval60.3%
Simplified60.3%
if -1.35000000000000008e-34 < n < 1.5Initial program 28.0%
Taylor expanded in i around 0 59.2%
if 1.5 < n Initial program 16.1%
Taylor expanded in i around 0 18.6%
associate-*r/18.6%
metadata-eval18.6%
Simplified18.6%
Taylor expanded in n around inf 69.1%
associate-*r*69.1%
*-commutative69.1%
*-commutative69.1%
Simplified69.1%
Final simplification62.3%
(FPCore (i n) :precision binary64 (if (or (<= n -1.3e+74) (not (<= n 1.1e+32))) (/ (* n (* i 100.0)) i) (* 100.0 (/ i (/ i n)))))
double code(double i, double n) {
double tmp;
if ((n <= -1.3e+74) || !(n <= 1.1e+32)) {
tmp = (n * (i * 100.0)) / i;
} else {
tmp = 100.0 * (i / (i / n));
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if ((n <= (-1.3d+74)) .or. (.not. (n <= 1.1d+32))) then
tmp = (n * (i * 100.0d0)) / i
else
tmp = 100.0d0 * (i / (i / n))
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if ((n <= -1.3e+74) || !(n <= 1.1e+32)) {
tmp = (n * (i * 100.0)) / i;
} else {
tmp = 100.0 * (i / (i / n));
}
return tmp;
}
def code(i, n): tmp = 0 if (n <= -1.3e+74) or not (n <= 1.1e+32): tmp = (n * (i * 100.0)) / i else: tmp = 100.0 * (i / (i / n)) return tmp
function code(i, n) tmp = 0.0 if ((n <= -1.3e+74) || !(n <= 1.1e+32)) tmp = Float64(Float64(n * Float64(i * 100.0)) / i); else tmp = Float64(100.0 * Float64(i / Float64(i / n))); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if ((n <= -1.3e+74) || ~((n <= 1.1e+32))) tmp = (n * (i * 100.0)) / i; else tmp = 100.0 * (i / (i / n)); end tmp_2 = tmp; end
code[i_, n_] := If[Or[LessEqual[n, -1.3e+74], N[Not[LessEqual[n, 1.1e+32]], $MachinePrecision]], N[(N[(n * N[(i * 100.0), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.3 \cdot 10^{+74} \lor \neg \left(n \leq 1.1 \cdot 10^{+32}\right):\\
\;\;\;\;\frac{n \cdot \left(i \cdot 100\right)}{i}\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\end{array}
\end{array}
if n < -1.3e74 or 1.1e32 < n Initial program 19.1%
associate-*r/19.1%
sub-neg19.1%
distribute-rgt-in19.1%
metadata-eval19.1%
metadata-eval19.1%
Simplified19.1%
Taylor expanded in i around 0 4.8%
*-commutative4.8%
Simplified4.8%
clear-num4.7%
associate-/r/4.8%
Applied egg-rr4.8%
*-un-lft-identity4.8%
times-frac5.3%
clear-num5.3%
/-rgt-identity5.3%
+-commutative5.3%
associate-+l+53.2%
metadata-eval53.2%
Applied egg-rr53.2%
associate-*r/62.5%
+-rgt-identity62.5%
*-commutative62.5%
Simplified62.5%
if -1.3e74 < n < 1.1e32Initial program 26.6%
Taylor expanded in i around 0 61.2%
Final simplification61.8%
(FPCore (i n) :precision binary64 (if (or (<= n -5e-35) (not (<= n 1.65e-62))) (* n 100.0) (* 100.0 (/ i (/ i n)))))
double code(double i, double n) {
double tmp;
if ((n <= -5e-35) || !(n <= 1.65e-62)) {
tmp = n * 100.0;
} else {
tmp = 100.0 * (i / (i / n));
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if ((n <= (-5d-35)) .or. (.not. (n <= 1.65d-62))) then
tmp = n * 100.0d0
else
tmp = 100.0d0 * (i / (i / n))
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if ((n <= -5e-35) || !(n <= 1.65e-62)) {
tmp = n * 100.0;
} else {
tmp = 100.0 * (i / (i / n));
}
return tmp;
}
def code(i, n): tmp = 0 if (n <= -5e-35) or not (n <= 1.65e-62): tmp = n * 100.0 else: tmp = 100.0 * (i / (i / n)) return tmp
function code(i, n) tmp = 0.0 if ((n <= -5e-35) || !(n <= 1.65e-62)) tmp = Float64(n * 100.0); else tmp = Float64(100.0 * Float64(i / Float64(i / n))); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if ((n <= -5e-35) || ~((n <= 1.65e-62))) tmp = n * 100.0; else tmp = 100.0 * (i / (i / n)); end tmp_2 = tmp; end
code[i_, n_] := If[Or[LessEqual[n, -5e-35], N[Not[LessEqual[n, 1.65e-62]], $MachinePrecision]], N[(n * 100.0), $MachinePrecision], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -5 \cdot 10^{-35} \lor \neg \left(n \leq 1.65 \cdot 10^{-62}\right):\\
\;\;\;\;n \cdot 100\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\end{array}
\end{array}
if n < -4.99999999999999964e-35 or 1.65000000000000002e-62 < n Initial program 20.9%
associate-/r/21.3%
associate-*r*21.3%
*-commutative21.3%
associate-*r/21.3%
sub-neg21.3%
distribute-lft-in21.3%
metadata-eval21.3%
metadata-eval21.3%
metadata-eval21.3%
fma-define21.3%
metadata-eval21.3%
Simplified21.3%
Taylor expanded in i around 0 54.5%
*-commutative54.5%
Simplified54.5%
if -4.99999999999999964e-35 < n < 1.65000000000000002e-62Initial program 27.2%
Taylor expanded in i around 0 63.4%
Final simplification57.2%
(FPCore (i n) :precision binary64 (* n 100.0))
double code(double i, double n) {
return n * 100.0;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
code = n * 100.0d0
end function
public static double code(double i, double n) {
return n * 100.0;
}
def code(i, n): return n * 100.0
function code(i, n) return Float64(n * 100.0) end
function tmp = code(i, n) tmp = n * 100.0; end
code[i_, n_] := N[(n * 100.0), $MachinePrecision]
\begin{array}{l}
\\
n \cdot 100
\end{array}
Initial program 22.9%
associate-/r/23.0%
associate-*r*23.0%
*-commutative23.0%
associate-*r/23.0%
sub-neg23.0%
distribute-lft-in23.0%
metadata-eval23.0%
metadata-eval23.0%
metadata-eval23.0%
fma-define23.0%
metadata-eval23.0%
Simplified23.0%
Taylor expanded in i around 0 48.8%
*-commutative48.8%
Simplified48.8%
(FPCore (i n) :precision binary64 (* i -50.0))
double code(double i, double n) {
return i * -50.0;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
code = i * (-50.0d0)
end function
public static double code(double i, double n) {
return i * -50.0;
}
def code(i, n): return i * -50.0
function code(i, n) return Float64(i * -50.0) end
function tmp = code(i, n) tmp = i * -50.0; end
code[i_, n_] := N[(i * -50.0), $MachinePrecision]
\begin{array}{l}
\\
i \cdot -50
\end{array}
Initial program 22.9%
Taylor expanded in i around 0 9.9%
associate-*r/9.9%
metadata-eval9.9%
Simplified9.9%
Taylor expanded in n around 0 2.9%
*-commutative2.9%
Simplified2.9%
(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 2024150
(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))))