
(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 18 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)) (t_1 (/ (+ t_0 -1.0) (/ i n))))
(if (<= t_1 0.0)
(* (/ (expm1 (* n (log1p (/ i n)))) i) (* n 100.0))
(if (<= t_1 INFINITY)
(* 100.0 (- (* t_0 (/ n i)) (/ n i)))
(/ -1.0 (- (* 0.01 (* i (/ (+ 0.5 (/ -0.5 n)) n))) (/ 0.01 n)))))))
double code(double i, double n) {
double t_0 = pow((1.0 + (i / n)), n);
double t_1 = (t_0 + -1.0) / (i / n);
double tmp;
if (t_1 <= 0.0) {
tmp = (expm1((n * log1p((i / n)))) / i) * (n * 100.0);
} else if (t_1 <= ((double) INFINITY)) {
tmp = 100.0 * ((t_0 * (n / i)) - (n / i));
} else {
tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n));
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = Math.pow((1.0 + (i / n)), n);
double t_1 = (t_0 + -1.0) / (i / n);
double tmp;
if (t_1 <= 0.0) {
tmp = (Math.expm1((n * Math.log1p((i / n)))) / i) * (n * 100.0);
} else if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = 100.0 * ((t_0 * (n / i)) - (n / i));
} else {
tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n));
}
return tmp;
}
def code(i, n): t_0 = math.pow((1.0 + (i / n)), n) t_1 = (t_0 + -1.0) / (i / n) tmp = 0 if t_1 <= 0.0: tmp = (math.expm1((n * math.log1p((i / n)))) / i) * (n * 100.0) elif t_1 <= math.inf: tmp = 100.0 * ((t_0 * (n / i)) - (n / i)) else: tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n)) return tmp
function code(i, n) t_0 = Float64(1.0 + Float64(i / n)) ^ n t_1 = Float64(Float64(t_0 + -1.0) / Float64(i / n)) tmp = 0.0 if (t_1 <= 0.0) tmp = Float64(Float64(expm1(Float64(n * log1p(Float64(i / n)))) / i) * Float64(n * 100.0)); elseif (t_1 <= Inf) tmp = Float64(100.0 * Float64(Float64(t_0 * Float64(n / i)) - Float64(n / i))); else tmp = Float64(-1.0 / Float64(Float64(0.01 * Float64(i * Float64(Float64(0.5 + Float64(-0.5 / n)) / n))) - Float64(0.01 / n))); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[Power[N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision], n], $MachinePrecision]}, Block[{t$95$1 = N[(N[(t$95$0 + -1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, 0.0], N[(N[(N[(Exp[N[(n * N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]] - 1), $MachinePrecision] / i), $MachinePrecision] * N[(n * 100.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(100.0 * N[(N[(t$95$0 * N[(n / i), $MachinePrecision]), $MachinePrecision] - N[(n / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(-1.0 / N[(N[(0.01 * N[(i * N[(N[(0.5 + N[(-0.5 / n), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(0.01 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(1 + \frac{i}{n}\right)}^{n}\\
t_1 := \frac{t\_0 + -1}{\frac{i}{n}}\\
\mathbf{if}\;t\_1 \leq 0:\\
\;\;\;\;\frac{\mathsf{expm1}\left(n \cdot \mathsf{log1p}\left(\frac{i}{n}\right)\right)}{i} \cdot \left(n \cdot 100\right)\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;100 \cdot \left(t\_0 \cdot \frac{n}{i} - \frac{n}{i}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{-1}{0.01 \cdot \left(i \cdot \frac{0.5 + \frac{-0.5}{n}}{n}\right) - \frac{0.01}{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)) < 0.0Initial program 24.8%
associate-*r/24.8%
sub-neg24.8%
distribute-rgt-in24.8%
metadata-eval24.8%
metadata-eval24.8%
Simplified24.8%
metadata-eval24.8%
metadata-eval24.8%
distribute-rgt-in24.8%
sub-neg24.8%
*-commutative24.8%
associate-*l/24.8%
associate-/r/24.8%
associate-*l*24.7%
add-exp-log24.7%
expm1-define24.7%
log-pow32.1%
log1p-define97.9%
Applied egg-rr97.9%
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 98.5%
div-sub98.8%
clear-num98.4%
sub-neg98.4%
div-inv98.4%
clear-num98.7%
Applied egg-rr98.7%
sub-neg98.7%
Simplified98.7%
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%
*-commutative1.9%
fma-undefine1.9%
*-commutative1.9%
associate-/r/0.0%
clear-num0.0%
frac-2neg0.0%
metadata-eval0.0%
fma-define0.0%
Applied egg-rr0.0%
distribute-neg-frac20.0%
fma-define0.0%
+-commutative0.0%
distribute-neg-in0.0%
metadata-eval0.0%
*-commutative0.0%
distribute-lft-neg-in0.0%
metadata-eval0.0%
Simplified0.0%
Taylor expanded in i around 0 99.6%
associate-/l*99.6%
cancel-sign-sub-inv99.6%
metadata-eval99.6%
associate-*r/99.6%
metadata-eval99.6%
associate-*r/99.7%
metadata-eval99.7%
Simplified99.7%
Final simplification98.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 -5e-188)
(* t_0 (* 100.0 (/ n i)))
(if (<= t_1 0.0)
(/
-1.0
(+
(* 0.01 (/ (* i (+ 0.5 (* 0.5 (/ -1.0 n)))) n))
(* 0.01 (/ -1.0 n))))
(if (<= t_1 INFINITY)
(* t_1 100.0)
(/ -1.0 (- (* 0.01 (* i (/ (+ 0.5 (/ -0.5 n)) n))) (/ 0.01 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 tmp;
if (t_1 <= -5e-188) {
tmp = t_0 * (100.0 * (n / i));
} else if (t_1 <= 0.0) {
tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n)));
} else if (t_1 <= ((double) INFINITY)) {
tmp = t_1 * 100.0;
} else {
tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / 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 tmp;
if (t_1 <= -5e-188) {
tmp = t_0 * (100.0 * (n / i));
} else if (t_1 <= 0.0) {
tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n)));
} else if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = t_1 * 100.0;
} else {
tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n));
}
return tmp;
}
def code(i, n): t_0 = math.pow((1.0 + (i / n)), n) + -1.0 t_1 = t_0 / (i / n) tmp = 0 if t_1 <= -5e-188: tmp = t_0 * (100.0 * (n / i)) elif t_1 <= 0.0: tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n))) elif t_1 <= math.inf: tmp = t_1 * 100.0 else: tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / 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)) tmp = 0.0 if (t_1 <= -5e-188) tmp = Float64(t_0 * Float64(100.0 * Float64(n / i))); elseif (t_1 <= 0.0) tmp = Float64(-1.0 / Float64(Float64(0.01 * Float64(Float64(i * Float64(0.5 + Float64(0.5 * Float64(-1.0 / n)))) / n)) + Float64(0.01 * Float64(-1.0 / n)))); elseif (t_1 <= Inf) tmp = Float64(t_1 * 100.0); else tmp = Float64(-1.0 / Float64(Float64(0.01 * Float64(i * Float64(Float64(0.5 + Float64(-0.5 / n)) / n))) - Float64(0.01 / n))); end return tmp end
function tmp_2 = code(i, n) t_0 = ((1.0 + (i / n)) ^ n) + -1.0; t_1 = t_0 / (i / n); tmp = 0.0; if (t_1 <= -5e-188) tmp = t_0 * (100.0 * (n / i)); elseif (t_1 <= 0.0) tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n))); elseif (t_1 <= Inf) tmp = t_1 * 100.0; else tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n)); end tmp_2 = 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, -5e-188], N[(t$95$0 * N[(100.0 * N[(n / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(-1.0 / N[(N[(0.01 * N[(N[(i * N[(0.5 + N[(0.5 * N[(-1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 * N[(-1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(t$95$1 * 100.0), $MachinePrecision], N[(-1.0 / N[(N[(0.01 * N[(i * N[(N[(0.5 + N[(-0.5 / n), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(0.01 / n), $MachinePrecision]), $MachinePrecision]), $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 -5 \cdot 10^{-188}:\\
\;\;\;\;t\_0 \cdot \left(100 \cdot \frac{n}{i}\right)\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;\frac{-1}{0.01 \cdot \frac{i \cdot \left(0.5 + 0.5 \cdot \frac{-1}{n}\right)}{n} + 0.01 \cdot \frac{-1}{n}}\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;t\_1 \cdot 100\\
\mathbf{else}:\\
\;\;\;\;\frac{-1}{0.01 \cdot \left(i \cdot \frac{0.5 + \frac{-0.5}{n}}{n}\right) - \frac{0.01}{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)) < -5.0000000000000001e-188Initial program 99.3%
associate-*r/99.3%
sub-neg99.3%
distribute-rgt-in99.3%
metadata-eval99.3%
metadata-eval99.3%
Simplified99.3%
metadata-eval99.3%
metadata-eval99.3%
distribute-rgt-in99.3%
sub-neg99.3%
associate-*r/99.3%
*-commutative99.3%
div-inv99.3%
clear-num99.6%
associate-*l*99.6%
add-exp-log99.6%
expm1-define99.6%
log-pow90.5%
log1p-define90.5%
Applied egg-rr90.5%
expm1-undefine90.5%
*-commutative90.5%
log1p-undefine90.5%
exp-to-pow99.6%
Applied egg-rr99.6%
if -5.0000000000000001e-188 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < 0.0Initial program 19.8%
associate-/r/19.8%
associate-*r*19.8%
*-commutative19.8%
associate-*r/19.8%
sub-neg19.8%
distribute-lft-in19.8%
metadata-eval19.8%
metadata-eval19.8%
metadata-eval19.8%
fma-define19.8%
metadata-eval19.8%
Simplified19.8%
*-commutative19.8%
fma-undefine19.8%
*-commutative19.8%
associate-/r/19.8%
clear-num19.8%
frac-2neg19.8%
metadata-eval19.8%
fma-define19.8%
Applied egg-rr19.8%
distribute-neg-frac219.8%
fma-define19.8%
+-commutative19.8%
distribute-neg-in19.8%
metadata-eval19.8%
*-commutative19.8%
distribute-lft-neg-in19.8%
metadata-eval19.8%
Simplified19.8%
Taylor expanded in i around 0 70.5%
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 98.5%
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%
*-commutative1.9%
fma-undefine1.9%
*-commutative1.9%
associate-/r/0.0%
clear-num0.0%
frac-2neg0.0%
metadata-eval0.0%
fma-define0.0%
Applied egg-rr0.0%
distribute-neg-frac20.0%
fma-define0.0%
+-commutative0.0%
distribute-neg-in0.0%
metadata-eval0.0%
*-commutative0.0%
distribute-lft-neg-in0.0%
metadata-eval0.0%
Simplified0.0%
Taylor expanded in i around 0 99.6%
associate-/l*99.6%
cancel-sign-sub-inv99.6%
metadata-eval99.6%
associate-*r/99.6%
metadata-eval99.6%
associate-*r/99.7%
metadata-eval99.7%
Simplified99.7%
Final simplification80.8%
(FPCore (i n)
:precision binary64
(let* ((t_0 (pow (+ 1.0 (/ i n)) n)) (t_1 (/ (+ t_0 -1.0) (/ i n))))
(if (<= t_1 -5e-188)
(* n (/ (+ (* t_0 100.0) -100.0) i))
(if (<= t_1 0.0)
(/
-1.0
(+
(* 0.01 (/ (* i (+ 0.5 (* 0.5 (/ -1.0 n)))) n))
(* 0.01 (/ -1.0 n))))
(if (<= t_1 INFINITY)
(* t_1 100.0)
(/ -1.0 (- (* 0.01 (* i (/ (+ 0.5 (/ -0.5 n)) n))) (/ 0.01 n))))))))
double code(double i, double n) {
double t_0 = pow((1.0 + (i / n)), n);
double t_1 = (t_0 + -1.0) / (i / n);
double tmp;
if (t_1 <= -5e-188) {
tmp = n * (((t_0 * 100.0) + -100.0) / i);
} else if (t_1 <= 0.0) {
tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n)));
} else if (t_1 <= ((double) INFINITY)) {
tmp = t_1 * 100.0;
} else {
tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n));
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = Math.pow((1.0 + (i / n)), n);
double t_1 = (t_0 + -1.0) / (i / n);
double tmp;
if (t_1 <= -5e-188) {
tmp = n * (((t_0 * 100.0) + -100.0) / i);
} else if (t_1 <= 0.0) {
tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n)));
} else if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = t_1 * 100.0;
} else {
tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n));
}
return tmp;
}
def code(i, n): t_0 = math.pow((1.0 + (i / n)), n) t_1 = (t_0 + -1.0) / (i / n) tmp = 0 if t_1 <= -5e-188: tmp = n * (((t_0 * 100.0) + -100.0) / i) elif t_1 <= 0.0: tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n))) elif t_1 <= math.inf: tmp = t_1 * 100.0 else: tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n)) return tmp
function code(i, n) t_0 = Float64(1.0 + Float64(i / n)) ^ n t_1 = Float64(Float64(t_0 + -1.0) / Float64(i / n)) tmp = 0.0 if (t_1 <= -5e-188) tmp = Float64(n * Float64(Float64(Float64(t_0 * 100.0) + -100.0) / i)); elseif (t_1 <= 0.0) tmp = Float64(-1.0 / Float64(Float64(0.01 * Float64(Float64(i * Float64(0.5 + Float64(0.5 * Float64(-1.0 / n)))) / n)) + Float64(0.01 * Float64(-1.0 / n)))); elseif (t_1 <= Inf) tmp = Float64(t_1 * 100.0); else tmp = Float64(-1.0 / Float64(Float64(0.01 * Float64(i * Float64(Float64(0.5 + Float64(-0.5 / n)) / n))) - Float64(0.01 / n))); end return tmp end
function tmp_2 = code(i, n) t_0 = (1.0 + (i / n)) ^ n; t_1 = (t_0 + -1.0) / (i / n); tmp = 0.0; if (t_1 <= -5e-188) tmp = n * (((t_0 * 100.0) + -100.0) / i); elseif (t_1 <= 0.0) tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n))); elseif (t_1 <= Inf) tmp = t_1 * 100.0; else tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n)); end tmp_2 = tmp; end
code[i_, n_] := Block[{t$95$0 = N[Power[N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision], n], $MachinePrecision]}, Block[{t$95$1 = N[(N[(t$95$0 + -1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -5e-188], N[(n * N[(N[(N[(t$95$0 * 100.0), $MachinePrecision] + -100.0), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(-1.0 / N[(N[(0.01 * N[(N[(i * N[(0.5 + N[(0.5 * N[(-1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 * N[(-1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(t$95$1 * 100.0), $MachinePrecision], N[(-1.0 / N[(N[(0.01 * N[(i * N[(N[(0.5 + N[(-0.5 / n), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(0.01 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(1 + \frac{i}{n}\right)}^{n}\\
t_1 := \frac{t\_0 + -1}{\frac{i}{n}}\\
\mathbf{if}\;t\_1 \leq -5 \cdot 10^{-188}:\\
\;\;\;\;n \cdot \frac{t\_0 \cdot 100 + -100}{i}\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;\frac{-1}{0.01 \cdot \frac{i \cdot \left(0.5 + 0.5 \cdot \frac{-1}{n}\right)}{n} + 0.01 \cdot \frac{-1}{n}}\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;t\_1 \cdot 100\\
\mathbf{else}:\\
\;\;\;\;\frac{-1}{0.01 \cdot \left(i \cdot \frac{0.5 + \frac{-0.5}{n}}{n}\right) - \frac{0.01}{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)) < -5.0000000000000001e-188Initial program 99.3%
associate-/r/99.3%
associate-*r*99.1%
*-commutative99.1%
associate-*r/99.4%
sub-neg99.4%
distribute-lft-in99.4%
metadata-eval99.4%
metadata-eval99.4%
metadata-eval99.4%
fma-define99.4%
metadata-eval99.4%
Simplified99.4%
fma-undefine99.4%
*-commutative99.4%
Applied egg-rr99.4%
if -5.0000000000000001e-188 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < 0.0Initial program 19.8%
associate-/r/19.8%
associate-*r*19.8%
*-commutative19.8%
associate-*r/19.8%
sub-neg19.8%
distribute-lft-in19.8%
metadata-eval19.8%
metadata-eval19.8%
metadata-eval19.8%
fma-define19.8%
metadata-eval19.8%
Simplified19.8%
*-commutative19.8%
fma-undefine19.8%
*-commutative19.8%
associate-/r/19.8%
clear-num19.8%
frac-2neg19.8%
metadata-eval19.8%
fma-define19.8%
Applied egg-rr19.8%
distribute-neg-frac219.8%
fma-define19.8%
+-commutative19.8%
distribute-neg-in19.8%
metadata-eval19.8%
*-commutative19.8%
distribute-lft-neg-in19.8%
metadata-eval19.8%
Simplified19.8%
Taylor expanded in i around 0 70.5%
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 98.5%
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%
*-commutative1.9%
fma-undefine1.9%
*-commutative1.9%
associate-/r/0.0%
clear-num0.0%
frac-2neg0.0%
metadata-eval0.0%
fma-define0.0%
Applied egg-rr0.0%
distribute-neg-frac20.0%
fma-define0.0%
+-commutative0.0%
distribute-neg-in0.0%
metadata-eval0.0%
*-commutative0.0%
distribute-lft-neg-in0.0%
metadata-eval0.0%
Simplified0.0%
Taylor expanded in i around 0 99.6%
associate-/l*99.6%
cancel-sign-sub-inv99.6%
metadata-eval99.6%
associate-*r/99.6%
metadata-eval99.6%
associate-*r/99.7%
metadata-eval99.7%
Simplified99.7%
Final simplification80.8%
(FPCore (i n)
:precision binary64
(let* ((t_0 (/ (+ (pow (+ 1.0 (/ i n)) n) -1.0) (/ i n))) (t_1 (* t_0 100.0)))
(if (<= t_0 -5e-188)
t_1
(if (<= t_0 0.0)
(/
-1.0
(+
(* 0.01 (/ (* i (+ 0.5 (* 0.5 (/ -1.0 n)))) n))
(* 0.01 (/ -1.0 n))))
(if (<= t_0 INFINITY)
t_1
(/ -1.0 (- (* 0.01 (* i (/ (+ 0.5 (/ -0.5 n)) n))) (/ 0.01 n))))))))
double code(double i, double n) {
double t_0 = (pow((1.0 + (i / n)), n) + -1.0) / (i / n);
double t_1 = t_0 * 100.0;
double tmp;
if (t_0 <= -5e-188) {
tmp = t_1;
} else if (t_0 <= 0.0) {
tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n)));
} else if (t_0 <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n));
}
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 t_1 = t_0 * 100.0;
double tmp;
if (t_0 <= -5e-188) {
tmp = t_1;
} else if (t_0 <= 0.0) {
tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n)));
} else if (t_0 <= Double.POSITIVE_INFINITY) {
tmp = t_1;
} else {
tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n));
}
return tmp;
}
def code(i, n): t_0 = (math.pow((1.0 + (i / n)), n) + -1.0) / (i / n) t_1 = t_0 * 100.0 tmp = 0 if t_0 <= -5e-188: tmp = t_1 elif t_0 <= 0.0: tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n))) elif t_0 <= math.inf: tmp = t_1 else: tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n)) return tmp
function code(i, n) t_0 = Float64(Float64((Float64(1.0 + Float64(i / n)) ^ n) + -1.0) / Float64(i / n)) t_1 = Float64(t_0 * 100.0) tmp = 0.0 if (t_0 <= -5e-188) tmp = t_1; elseif (t_0 <= 0.0) tmp = Float64(-1.0 / Float64(Float64(0.01 * Float64(Float64(i * Float64(0.5 + Float64(0.5 * Float64(-1.0 / n)))) / n)) + Float64(0.01 * Float64(-1.0 / n)))); elseif (t_0 <= Inf) tmp = t_1; else tmp = Float64(-1.0 / Float64(Float64(0.01 * Float64(i * Float64(Float64(0.5 + Float64(-0.5 / n)) / n))) - Float64(0.01 / n))); end return tmp end
function tmp_2 = code(i, n) t_0 = (((1.0 + (i / n)) ^ n) + -1.0) / (i / n); t_1 = t_0 * 100.0; tmp = 0.0; if (t_0 <= -5e-188) tmp = t_1; elseif (t_0 <= 0.0) tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n))); elseif (t_0 <= Inf) tmp = t_1; else tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n)); end tmp_2 = 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]}, Block[{t$95$1 = N[(t$95$0 * 100.0), $MachinePrecision]}, If[LessEqual[t$95$0, -5e-188], t$95$1, If[LessEqual[t$95$0, 0.0], N[(-1.0 / N[(N[(0.01 * N[(N[(i * N[(0.5 + N[(0.5 * N[(-1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 * N[(-1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, Infinity], t$95$1, N[(-1.0 / N[(N[(0.01 * N[(i * N[(N[(0.5 + N[(-0.5 / n), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(0.01 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{{\left(1 + \frac{i}{n}\right)}^{n} + -1}{\frac{i}{n}}\\
t_1 := t\_0 \cdot 100\\
\mathbf{if}\;t\_0 \leq -5 \cdot 10^{-188}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;\frac{-1}{0.01 \cdot \frac{i \cdot \left(0.5 + 0.5 \cdot \frac{-1}{n}\right)}{n} + 0.01 \cdot \frac{-1}{n}}\\
\mathbf{elif}\;t\_0 \leq \infty:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{-1}{0.01 \cdot \left(i \cdot \frac{0.5 + \frac{-0.5}{n}}{n}\right) - \frac{0.01}{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)) < -5.0000000000000001e-188 or 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 98.8%
if -5.0000000000000001e-188 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < 0.0Initial program 19.8%
associate-/r/19.8%
associate-*r*19.8%
*-commutative19.8%
associate-*r/19.8%
sub-neg19.8%
distribute-lft-in19.8%
metadata-eval19.8%
metadata-eval19.8%
metadata-eval19.8%
fma-define19.8%
metadata-eval19.8%
Simplified19.8%
*-commutative19.8%
fma-undefine19.8%
*-commutative19.8%
associate-/r/19.8%
clear-num19.8%
frac-2neg19.8%
metadata-eval19.8%
fma-define19.8%
Applied egg-rr19.8%
distribute-neg-frac219.8%
fma-define19.8%
+-commutative19.8%
distribute-neg-in19.8%
metadata-eval19.8%
*-commutative19.8%
distribute-lft-neg-in19.8%
metadata-eval19.8%
Simplified19.8%
Taylor expanded in i around 0 70.5%
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%
*-commutative1.9%
fma-undefine1.9%
*-commutative1.9%
associate-/r/0.0%
clear-num0.0%
frac-2neg0.0%
metadata-eval0.0%
fma-define0.0%
Applied egg-rr0.0%
distribute-neg-frac20.0%
fma-define0.0%
+-commutative0.0%
distribute-neg-in0.0%
metadata-eval0.0%
*-commutative0.0%
distribute-lft-neg-in0.0%
metadata-eval0.0%
Simplified0.0%
Taylor expanded in i around 0 99.6%
associate-/l*99.6%
cancel-sign-sub-inv99.6%
metadata-eval99.6%
associate-*r/99.6%
metadata-eval99.6%
associate-*r/99.7%
metadata-eval99.7%
Simplified99.7%
Final simplification80.8%
(FPCore (i n)
:precision binary64
(let* ((t_0 (pow (+ 1.0 (/ i n)) n)) (t_1 (/ (+ t_0 -1.0) (/ i n))))
(if (<= t_1 0.0)
(* n (* (expm1 (* n (log1p (/ i n)))) (/ 100.0 i)))
(if (<= t_1 INFINITY)
(* 100.0 (- (* t_0 (/ n i)) (/ n i)))
(/ -1.0 (- (* 0.01 (* i (/ (+ 0.5 (/ -0.5 n)) n))) (/ 0.01 n)))))))
double code(double i, double n) {
double t_0 = pow((1.0 + (i / n)), n);
double t_1 = (t_0 + -1.0) / (i / n);
double tmp;
if (t_1 <= 0.0) {
tmp = n * (expm1((n * log1p((i / n)))) * (100.0 / i));
} else if (t_1 <= ((double) INFINITY)) {
tmp = 100.0 * ((t_0 * (n / i)) - (n / i));
} else {
tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n));
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = Math.pow((1.0 + (i / n)), n);
double t_1 = (t_0 + -1.0) / (i / n);
double tmp;
if (t_1 <= 0.0) {
tmp = n * (Math.expm1((n * Math.log1p((i / n)))) * (100.0 / i));
} else if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = 100.0 * ((t_0 * (n / i)) - (n / i));
} else {
tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n));
}
return tmp;
}
def code(i, n): t_0 = math.pow((1.0 + (i / n)), n) t_1 = (t_0 + -1.0) / (i / n) tmp = 0 if t_1 <= 0.0: tmp = n * (math.expm1((n * math.log1p((i / n)))) * (100.0 / i)) elif t_1 <= math.inf: tmp = 100.0 * ((t_0 * (n / i)) - (n / i)) else: tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n)) return tmp
function code(i, n) t_0 = Float64(1.0 + Float64(i / n)) ^ n t_1 = Float64(Float64(t_0 + -1.0) / Float64(i / n)) tmp = 0.0 if (t_1 <= 0.0) tmp = Float64(n * Float64(expm1(Float64(n * log1p(Float64(i / n)))) * Float64(100.0 / i))); elseif (t_1 <= Inf) tmp = Float64(100.0 * Float64(Float64(t_0 * Float64(n / i)) - Float64(n / i))); else tmp = Float64(-1.0 / Float64(Float64(0.01 * Float64(i * Float64(Float64(0.5 + Float64(-0.5 / n)) / n))) - Float64(0.01 / n))); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[Power[N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision], n], $MachinePrecision]}, Block[{t$95$1 = N[(N[(t$95$0 + -1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, 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$1, Infinity], N[(100.0 * N[(N[(t$95$0 * N[(n / i), $MachinePrecision]), $MachinePrecision] - N[(n / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(-1.0 / N[(N[(0.01 * N[(i * N[(N[(0.5 + N[(-0.5 / n), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(0.01 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(1 + \frac{i}{n}\right)}^{n}\\
t_1 := \frac{t\_0 + -1}{\frac{i}{n}}\\
\mathbf{if}\;t\_1 \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\_1 \leq \infty:\\
\;\;\;\;100 \cdot \left(t\_0 \cdot \frac{n}{i} - \frac{n}{i}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{-1}{0.01 \cdot \left(i \cdot \frac{0.5 + \frac{-0.5}{n}}{n}\right) - \frac{0.01}{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)) < 0.0Initial program 24.8%
associate-/r/24.8%
associate-*r*24.7%
*-commutative24.7%
associate-*r/24.7%
sub-neg24.7%
distribute-lft-in24.7%
metadata-eval24.7%
metadata-eval24.7%
metadata-eval24.7%
fma-define24.7%
metadata-eval24.7%
Simplified24.7%
fma-undefine24.7%
metadata-eval24.7%
metadata-eval24.7%
distribute-lft-in24.7%
sub-neg24.7%
*-commutative24.7%
add-exp-log24.7%
expm1-define24.7%
log-pow32.1%
log1p-define97.8%
Applied egg-rr97.8%
associate-/l*97.8%
Applied egg-rr97.8%
metadata-eval97.8%
associate-*r/97.7%
*-commutative97.7%
associate-*r/97.8%
metadata-eval97.8%
Simplified97.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 98.5%
div-sub98.8%
clear-num98.4%
sub-neg98.4%
div-inv98.4%
clear-num98.7%
Applied egg-rr98.7%
sub-neg98.7%
Simplified98.7%
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%
*-commutative1.9%
fma-undefine1.9%
*-commutative1.9%
associate-/r/0.0%
clear-num0.0%
frac-2neg0.0%
metadata-eval0.0%
fma-define0.0%
Applied egg-rr0.0%
distribute-neg-frac20.0%
fma-define0.0%
+-commutative0.0%
distribute-neg-in0.0%
metadata-eval0.0%
*-commutative0.0%
distribute-lft-neg-in0.0%
metadata-eval0.0%
Simplified0.0%
Taylor expanded in i around 0 99.6%
associate-/l*99.6%
cancel-sign-sub-inv99.6%
metadata-eval99.6%
associate-*r/99.6%
metadata-eval99.6%
associate-*r/99.7%
metadata-eval99.7%
Simplified99.7%
Final simplification98.3%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* (* n 100.0) (/ (expm1 (* i (+ 1.0 (/ (* i -0.5) n)))) i))))
(if (<= n -6.5e+131)
t_0
(if (<= n 2.55e-9)
(/
-1.0
(+
(* 0.01 (/ (* i (+ 0.5 (* 0.5 (/ -1.0 n)))) n))
(* 0.01 (/ -1.0 n))))
(if (<= n 1.3e+103)
(*
n
(/
(*
i
(+
100.0
(*
i
(*
100.0
(+
(*
i
(+
0.16666666666666666
(- (/ 0.3333333333333333 (pow n 2.0)) (/ 0.5 n))))
(- 0.5 (/ 0.5 n)))))))
i))
t_0)))))
double code(double i, double n) {
double t_0 = (n * 100.0) * (expm1((i * (1.0 + ((i * -0.5) / n)))) / i);
double tmp;
if (n <= -6.5e+131) {
tmp = t_0;
} else if (n <= 2.55e-9) {
tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n)));
} else if (n <= 1.3e+103) {
tmp = n * ((i * (100.0 + (i * (100.0 * ((i * (0.16666666666666666 + ((0.3333333333333333 / pow(n, 2.0)) - (0.5 / n)))) + (0.5 - (0.5 / n))))))) / i);
} else {
tmp = t_0;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = (n * 100.0) * (Math.expm1((i * (1.0 + ((i * -0.5) / n)))) / i);
double tmp;
if (n <= -6.5e+131) {
tmp = t_0;
} else if (n <= 2.55e-9) {
tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n)));
} else if (n <= 1.3e+103) {
tmp = n * ((i * (100.0 + (i * (100.0 * ((i * (0.16666666666666666 + ((0.3333333333333333 / Math.pow(n, 2.0)) - (0.5 / n)))) + (0.5 - (0.5 / n))))))) / i);
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = (n * 100.0) * (math.expm1((i * (1.0 + ((i * -0.5) / n)))) / i) tmp = 0 if n <= -6.5e+131: tmp = t_0 elif n <= 2.55e-9: tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n))) elif n <= 1.3e+103: tmp = n * ((i * (100.0 + (i * (100.0 * ((i * (0.16666666666666666 + ((0.3333333333333333 / math.pow(n, 2.0)) - (0.5 / n)))) + (0.5 - (0.5 / n))))))) / i) else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(Float64(n * 100.0) * Float64(expm1(Float64(i * Float64(1.0 + Float64(Float64(i * -0.5) / n)))) / i)) tmp = 0.0 if (n <= -6.5e+131) tmp = t_0; elseif (n <= 2.55e-9) tmp = Float64(-1.0 / Float64(Float64(0.01 * Float64(Float64(i * Float64(0.5 + Float64(0.5 * Float64(-1.0 / n)))) / n)) + Float64(0.01 * Float64(-1.0 / n)))); elseif (n <= 1.3e+103) tmp = Float64(n * Float64(Float64(i * Float64(100.0 + Float64(i * Float64(100.0 * Float64(Float64(i * Float64(0.16666666666666666 + Float64(Float64(0.3333333333333333 / (n ^ 2.0)) - Float64(0.5 / n)))) + Float64(0.5 - Float64(0.5 / n))))))) / i)); else tmp = t_0; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(n * 100.0), $MachinePrecision] * N[(N[(Exp[N[(i * N[(1.0 + N[(N[(i * -0.5), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]] - 1), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -6.5e+131], t$95$0, If[LessEqual[n, 2.55e-9], N[(-1.0 / N[(N[(0.01 * N[(N[(i * N[(0.5 + N[(0.5 * N[(-1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 * N[(-1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.3e+103], N[(n * N[(N[(i * N[(100.0 + N[(i * N[(100.0 * N[(N[(i * N[(0.16666666666666666 + N[(N[(0.3333333333333333 / N[Power[n, 2.0], $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(n \cdot 100\right) \cdot \frac{\mathsf{expm1}\left(i \cdot \left(1 + \frac{i \cdot -0.5}{n}\right)\right)}{i}\\
\mathbf{if}\;n \leq -6.5 \cdot 10^{+131}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 2.55 \cdot 10^{-9}:\\
\;\;\;\;\frac{-1}{0.01 \cdot \frac{i \cdot \left(0.5 + 0.5 \cdot \frac{-1}{n}\right)}{n} + 0.01 \cdot \frac{-1}{n}}\\
\mathbf{elif}\;n \leq 1.3 \cdot 10^{+103}:\\
\;\;\;\;n \cdot \frac{i \cdot \left(100 + i \cdot \left(100 \cdot \left(i \cdot \left(0.16666666666666666 + \left(\frac{0.3333333333333333}{{n}^{2}} - \frac{0.5}{n}\right)\right) + \left(0.5 - \frac{0.5}{n}\right)\right)\right)\right)}{i}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -6.5e131 or 1.3000000000000001e103 < n Initial program 14.2%
associate-*r/14.2%
sub-neg14.2%
distribute-rgt-in14.2%
metadata-eval14.2%
metadata-eval14.2%
Simplified14.2%
metadata-eval14.2%
metadata-eval14.2%
distribute-rgt-in14.2%
sub-neg14.2%
*-commutative14.2%
associate-*l/14.2%
associate-/r/14.9%
associate-*l*14.9%
add-exp-log14.9%
expm1-define14.9%
log-pow12.1%
log1p-define62.2%
Applied egg-rr62.2%
Taylor expanded in i around 0 88.7%
associate-*r/88.7%
*-commutative88.7%
Simplified88.7%
if -6.5e131 < n < 2.55000000000000009e-9Initial program 35.0%
associate-/r/35.2%
associate-*r*35.2%
*-commutative35.2%
associate-*r/35.2%
sub-neg35.2%
distribute-lft-in35.2%
metadata-eval35.2%
metadata-eval35.2%
metadata-eval35.2%
fma-define35.2%
metadata-eval35.2%
Simplified35.2%
*-commutative35.2%
fma-undefine35.2%
*-commutative35.2%
associate-/r/35.0%
clear-num35.0%
frac-2neg35.0%
metadata-eval35.0%
fma-define35.0%
Applied egg-rr35.0%
distribute-neg-frac235.0%
fma-define35.0%
+-commutative35.0%
distribute-neg-in35.0%
metadata-eval35.0%
*-commutative35.0%
distribute-lft-neg-in35.0%
metadata-eval35.0%
Simplified35.0%
Taylor expanded in i around 0 80.2%
if 2.55000000000000009e-9 < n < 1.3000000000000001e103Initial program 42.3%
associate-/r/42.5%
associate-*r*42.5%
*-commutative42.5%
associate-*r/42.5%
sub-neg42.5%
distribute-lft-in42.5%
metadata-eval42.5%
metadata-eval42.5%
metadata-eval42.5%
fma-define42.5%
metadata-eval42.5%
Simplified42.5%
Taylor expanded in i around 0 91.1%
distribute-lft-out91.1%
associate--l+91.1%
associate-*r/91.1%
metadata-eval91.1%
associate-*r/91.1%
metadata-eval91.1%
associate-*r/91.1%
metadata-eval91.1%
Simplified91.1%
Final simplification84.6%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* (* n 100.0) (/ (expm1 (* i (+ 1.0 (/ (* i -0.5) n)))) i))))
(if (<= n -6.5e+131)
t_0
(if (<= n 2.55e-9)
(/
-1.0
(+
(* 0.01 (/ (* i (+ 0.5 (* 0.5 (/ -1.0 n)))) n))
(* 0.01 (/ -1.0 n))))
(if (<= n 2.7e+101)
(*
n
(+
100.0
(*
i
(*
100.0
(+
(*
i
(+
0.16666666666666666
(- (/ 0.3333333333333333 (pow n 2.0)) (/ 0.5 n))))
(- 0.5 (/ 0.5 n)))))))
t_0)))))
double code(double i, double n) {
double t_0 = (n * 100.0) * (expm1((i * (1.0 + ((i * -0.5) / n)))) / i);
double tmp;
if (n <= -6.5e+131) {
tmp = t_0;
} else if (n <= 2.55e-9) {
tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n)));
} else if (n <= 2.7e+101) {
tmp = n * (100.0 + (i * (100.0 * ((i * (0.16666666666666666 + ((0.3333333333333333 / pow(n, 2.0)) - (0.5 / n)))) + (0.5 - (0.5 / n))))));
} else {
tmp = t_0;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = (n * 100.0) * (Math.expm1((i * (1.0 + ((i * -0.5) / n)))) / i);
double tmp;
if (n <= -6.5e+131) {
tmp = t_0;
} else if (n <= 2.55e-9) {
tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n)));
} else if (n <= 2.7e+101) {
tmp = n * (100.0 + (i * (100.0 * ((i * (0.16666666666666666 + ((0.3333333333333333 / Math.pow(n, 2.0)) - (0.5 / n)))) + (0.5 - (0.5 / n))))));
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = (n * 100.0) * (math.expm1((i * (1.0 + ((i * -0.5) / n)))) / i) tmp = 0 if n <= -6.5e+131: tmp = t_0 elif n <= 2.55e-9: tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n))) elif n <= 2.7e+101: tmp = n * (100.0 + (i * (100.0 * ((i * (0.16666666666666666 + ((0.3333333333333333 / math.pow(n, 2.0)) - (0.5 / n)))) + (0.5 - (0.5 / n)))))) else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(Float64(n * 100.0) * Float64(expm1(Float64(i * Float64(1.0 + Float64(Float64(i * -0.5) / n)))) / i)) tmp = 0.0 if (n <= -6.5e+131) tmp = t_0; elseif (n <= 2.55e-9) tmp = Float64(-1.0 / Float64(Float64(0.01 * Float64(Float64(i * Float64(0.5 + Float64(0.5 * Float64(-1.0 / n)))) / n)) + Float64(0.01 * Float64(-1.0 / n)))); elseif (n <= 2.7e+101) tmp = Float64(n * Float64(100.0 + Float64(i * Float64(100.0 * Float64(Float64(i * Float64(0.16666666666666666 + Float64(Float64(0.3333333333333333 / (n ^ 2.0)) - Float64(0.5 / n)))) + Float64(0.5 - Float64(0.5 / n))))))); else tmp = t_0; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(n * 100.0), $MachinePrecision] * N[(N[(Exp[N[(i * N[(1.0 + N[(N[(i * -0.5), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]] - 1), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -6.5e+131], t$95$0, If[LessEqual[n, 2.55e-9], N[(-1.0 / N[(N[(0.01 * N[(N[(i * N[(0.5 + N[(0.5 * N[(-1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 * N[(-1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 2.7e+101], N[(n * N[(100.0 + N[(i * N[(100.0 * N[(N[(i * N[(0.16666666666666666 + N[(N[(0.3333333333333333 / N[Power[n, 2.0], $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(n \cdot 100\right) \cdot \frac{\mathsf{expm1}\left(i \cdot \left(1 + \frac{i \cdot -0.5}{n}\right)\right)}{i}\\
\mathbf{if}\;n \leq -6.5 \cdot 10^{+131}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 2.55 \cdot 10^{-9}:\\
\;\;\;\;\frac{-1}{0.01 \cdot \frac{i \cdot \left(0.5 + 0.5 \cdot \frac{-1}{n}\right)}{n} + 0.01 \cdot \frac{-1}{n}}\\
\mathbf{elif}\;n \leq 2.7 \cdot 10^{+101}:\\
\;\;\;\;n \cdot \left(100 + i \cdot \left(100 \cdot \left(i \cdot \left(0.16666666666666666 + \left(\frac{0.3333333333333333}{{n}^{2}} - \frac{0.5}{n}\right)\right) + \left(0.5 - \frac{0.5}{n}\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -6.5e131 or 2.70000000000000006e101 < n Initial program 14.2%
associate-*r/14.2%
sub-neg14.2%
distribute-rgt-in14.2%
metadata-eval14.2%
metadata-eval14.2%
Simplified14.2%
metadata-eval14.2%
metadata-eval14.2%
distribute-rgt-in14.2%
sub-neg14.2%
*-commutative14.2%
associate-*l/14.2%
associate-/r/14.9%
associate-*l*14.9%
add-exp-log14.9%
expm1-define14.9%
log-pow12.1%
log1p-define62.2%
Applied egg-rr62.2%
Taylor expanded in i around 0 88.7%
associate-*r/88.7%
*-commutative88.7%
Simplified88.7%
if -6.5e131 < n < 2.55000000000000009e-9Initial program 35.0%
associate-/r/35.2%
associate-*r*35.2%
*-commutative35.2%
associate-*r/35.2%
sub-neg35.2%
distribute-lft-in35.2%
metadata-eval35.2%
metadata-eval35.2%
metadata-eval35.2%
fma-define35.2%
metadata-eval35.2%
Simplified35.2%
*-commutative35.2%
fma-undefine35.2%
*-commutative35.2%
associate-/r/35.0%
clear-num35.0%
frac-2neg35.0%
metadata-eval35.0%
fma-define35.0%
Applied egg-rr35.0%
distribute-neg-frac235.0%
fma-define35.0%
+-commutative35.0%
distribute-neg-in35.0%
metadata-eval35.0%
*-commutative35.0%
distribute-lft-neg-in35.0%
metadata-eval35.0%
Simplified35.0%
Taylor expanded in i around 0 80.2%
if 2.55000000000000009e-9 < n < 2.70000000000000006e101Initial program 42.3%
associate-/r/42.5%
associate-*r*42.5%
*-commutative42.5%
associate-*r/42.5%
sub-neg42.5%
distribute-lft-in42.5%
metadata-eval42.5%
metadata-eval42.5%
metadata-eval42.5%
fma-define42.5%
metadata-eval42.5%
Simplified42.5%
Taylor expanded in i around 0 87.2%
distribute-lft-out87.2%
associate--l+87.2%
associate-*r/87.2%
metadata-eval87.2%
associate-*r/87.2%
metadata-eval87.2%
associate-*r/87.2%
metadata-eval87.2%
Simplified87.2%
Final simplification84.3%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* (* n 100.0) (/ (expm1 (* i (+ 1.0 (/ (* i -0.5) n)))) i))))
(if (<= n -6.5e+131)
t_0
(if (<= n 2.55e-9)
(/
-1.0
(+
(* 0.01 (/ (* i (+ 0.5 (* 0.5 (/ -1.0 n)))) n))
(* 0.01 (/ -1.0 n))))
(if (<= n 5.6e+103)
(* n (/ (* i (+ 100.0 (* 100.0 (* i (- 0.5 (/ 0.5 n)))))) i))
t_0)))))
double code(double i, double n) {
double t_0 = (n * 100.0) * (expm1((i * (1.0 + ((i * -0.5) / n)))) / i);
double tmp;
if (n <= -6.5e+131) {
tmp = t_0;
} else if (n <= 2.55e-9) {
tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n)));
} else if (n <= 5.6e+103) {
tmp = n * ((i * (100.0 + (100.0 * (i * (0.5 - (0.5 / n)))))) / i);
} else {
tmp = t_0;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = (n * 100.0) * (Math.expm1((i * (1.0 + ((i * -0.5) / n)))) / i);
double tmp;
if (n <= -6.5e+131) {
tmp = t_0;
} else if (n <= 2.55e-9) {
tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n)));
} else if (n <= 5.6e+103) {
tmp = n * ((i * (100.0 + (100.0 * (i * (0.5 - (0.5 / n)))))) / i);
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = (n * 100.0) * (math.expm1((i * (1.0 + ((i * -0.5) / n)))) / i) tmp = 0 if n <= -6.5e+131: tmp = t_0 elif n <= 2.55e-9: tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n))) elif n <= 5.6e+103: tmp = n * ((i * (100.0 + (100.0 * (i * (0.5 - (0.5 / n)))))) / i) else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(Float64(n * 100.0) * Float64(expm1(Float64(i * Float64(1.0 + Float64(Float64(i * -0.5) / n)))) / i)) tmp = 0.0 if (n <= -6.5e+131) tmp = t_0; elseif (n <= 2.55e-9) tmp = Float64(-1.0 / Float64(Float64(0.01 * Float64(Float64(i * Float64(0.5 + Float64(0.5 * Float64(-1.0 / n)))) / n)) + Float64(0.01 * Float64(-1.0 / n)))); elseif (n <= 5.6e+103) tmp = Float64(n * Float64(Float64(i * Float64(100.0 + Float64(100.0 * Float64(i * Float64(0.5 - Float64(0.5 / n)))))) / i)); else tmp = t_0; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(n * 100.0), $MachinePrecision] * N[(N[(Exp[N[(i * N[(1.0 + N[(N[(i * -0.5), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]] - 1), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -6.5e+131], t$95$0, If[LessEqual[n, 2.55e-9], N[(-1.0 / N[(N[(0.01 * N[(N[(i * N[(0.5 + N[(0.5 * N[(-1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 * N[(-1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 5.6e+103], N[(n * N[(N[(i * N[(100.0 + N[(100.0 * N[(i * N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(n \cdot 100\right) \cdot \frac{\mathsf{expm1}\left(i \cdot \left(1 + \frac{i \cdot -0.5}{n}\right)\right)}{i}\\
\mathbf{if}\;n \leq -6.5 \cdot 10^{+131}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 2.55 \cdot 10^{-9}:\\
\;\;\;\;\frac{-1}{0.01 \cdot \frac{i \cdot \left(0.5 + 0.5 \cdot \frac{-1}{n}\right)}{n} + 0.01 \cdot \frac{-1}{n}}\\
\mathbf{elif}\;n \leq 5.6 \cdot 10^{+103}:\\
\;\;\;\;n \cdot \frac{i \cdot \left(100 + 100 \cdot \left(i \cdot \left(0.5 - \frac{0.5}{n}\right)\right)\right)}{i}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -6.5e131 or 5.60000000000000017e103 < n Initial program 14.2%
associate-*r/14.2%
sub-neg14.2%
distribute-rgt-in14.2%
metadata-eval14.2%
metadata-eval14.2%
Simplified14.2%
metadata-eval14.2%
metadata-eval14.2%
distribute-rgt-in14.2%
sub-neg14.2%
*-commutative14.2%
associate-*l/14.2%
associate-/r/14.9%
associate-*l*14.9%
add-exp-log14.9%
expm1-define14.9%
log-pow12.1%
log1p-define62.2%
Applied egg-rr62.2%
Taylor expanded in i around 0 88.7%
associate-*r/88.7%
*-commutative88.7%
Simplified88.7%
if -6.5e131 < n < 2.55000000000000009e-9Initial program 35.0%
associate-/r/35.2%
associate-*r*35.2%
*-commutative35.2%
associate-*r/35.2%
sub-neg35.2%
distribute-lft-in35.2%
metadata-eval35.2%
metadata-eval35.2%
metadata-eval35.2%
fma-define35.2%
metadata-eval35.2%
Simplified35.2%
*-commutative35.2%
fma-undefine35.2%
*-commutative35.2%
associate-/r/35.0%
clear-num35.0%
frac-2neg35.0%
metadata-eval35.0%
fma-define35.0%
Applied egg-rr35.0%
distribute-neg-frac235.0%
fma-define35.0%
+-commutative35.0%
distribute-neg-in35.0%
metadata-eval35.0%
*-commutative35.0%
distribute-lft-neg-in35.0%
metadata-eval35.0%
Simplified35.0%
Taylor expanded in i around 0 80.2%
if 2.55000000000000009e-9 < n < 5.60000000000000017e103Initial program 42.3%
associate-/r/42.5%
associate-*r*42.5%
*-commutative42.5%
associate-*r/42.5%
sub-neg42.5%
distribute-lft-in42.5%
metadata-eval42.5%
metadata-eval42.5%
metadata-eval42.5%
fma-define42.5%
metadata-eval42.5%
Simplified42.5%
Taylor expanded in i around 0 86.8%
*-commutative86.8%
associate-*r/86.8%
metadata-eval86.8%
Simplified86.8%
Final simplification84.3%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* n (/ (* 100.0 (expm1 (* i (+ 1.0 (/ (* i -0.5) n))))) i))))
(if (<= n -6.5e+131)
t_0
(if (<= n 2.55e-9)
(/
-1.0
(+
(* 0.01 (/ (* i (+ 0.5 (* 0.5 (/ -1.0 n)))) n))
(* 0.01 (/ -1.0 n))))
(if (<= n 1.15e+99)
(* n (/ (* i (+ 100.0 (* 100.0 (* i (- 0.5 (/ 0.5 n)))))) i))
t_0)))))
double code(double i, double n) {
double t_0 = n * ((100.0 * expm1((i * (1.0 + ((i * -0.5) / n))))) / i);
double tmp;
if (n <= -6.5e+131) {
tmp = t_0;
} else if (n <= 2.55e-9) {
tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n)));
} else if (n <= 1.15e+99) {
tmp = n * ((i * (100.0 + (100.0 * (i * (0.5 - (0.5 / n)))))) / i);
} else {
tmp = t_0;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = n * ((100.0 * Math.expm1((i * (1.0 + ((i * -0.5) / n))))) / i);
double tmp;
if (n <= -6.5e+131) {
tmp = t_0;
} else if (n <= 2.55e-9) {
tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n)));
} else if (n <= 1.15e+99) {
tmp = n * ((i * (100.0 + (100.0 * (i * (0.5 - (0.5 / n)))))) / i);
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = n * ((100.0 * math.expm1((i * (1.0 + ((i * -0.5) / n))))) / i) tmp = 0 if n <= -6.5e+131: tmp = t_0 elif n <= 2.55e-9: tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n))) elif n <= 1.15e+99: tmp = n * ((i * (100.0 + (100.0 * (i * (0.5 - (0.5 / n)))))) / i) else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(n * Float64(Float64(100.0 * expm1(Float64(i * Float64(1.0 + Float64(Float64(i * -0.5) / n))))) / i)) tmp = 0.0 if (n <= -6.5e+131) tmp = t_0; elseif (n <= 2.55e-9) tmp = Float64(-1.0 / Float64(Float64(0.01 * Float64(Float64(i * Float64(0.5 + Float64(0.5 * Float64(-1.0 / n)))) / n)) + Float64(0.01 * Float64(-1.0 / n)))); elseif (n <= 1.15e+99) tmp = Float64(n * Float64(Float64(i * Float64(100.0 + Float64(100.0 * Float64(i * Float64(0.5 - Float64(0.5 / n)))))) / i)); else tmp = t_0; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(n * N[(N[(100.0 * N[(Exp[N[(i * N[(1.0 + N[(N[(i * -0.5), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -6.5e+131], t$95$0, If[LessEqual[n, 2.55e-9], N[(-1.0 / N[(N[(0.01 * N[(N[(i * N[(0.5 + N[(0.5 * N[(-1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 * N[(-1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.15e+99], N[(n * N[(N[(i * N[(100.0 + N[(100.0 * N[(i * N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := n \cdot \frac{100 \cdot \mathsf{expm1}\left(i \cdot \left(1 + \frac{i \cdot -0.5}{n}\right)\right)}{i}\\
\mathbf{if}\;n \leq -6.5 \cdot 10^{+131}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 2.55 \cdot 10^{-9}:\\
\;\;\;\;\frac{-1}{0.01 \cdot \frac{i \cdot \left(0.5 + 0.5 \cdot \frac{-1}{n}\right)}{n} + 0.01 \cdot \frac{-1}{n}}\\
\mathbf{elif}\;n \leq 1.15 \cdot 10^{+99}:\\
\;\;\;\;n \cdot \frac{i \cdot \left(100 + 100 \cdot \left(i \cdot \left(0.5 - \frac{0.5}{n}\right)\right)\right)}{i}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -6.5e131 or 1.1500000000000001e99 < n Initial program 14.2%
associate-/r/14.9%
associate-*r*14.9%
*-commutative14.9%
associate-*r/14.9%
sub-neg14.9%
distribute-lft-in14.9%
metadata-eval14.9%
metadata-eval14.9%
metadata-eval14.9%
fma-define14.9%
metadata-eval14.9%
Simplified14.9%
fma-undefine14.9%
metadata-eval14.9%
metadata-eval14.9%
distribute-lft-in14.9%
sub-neg14.9%
*-commutative14.9%
add-exp-log14.9%
expm1-define14.9%
log-pow12.1%
log1p-define62.1%
Applied egg-rr62.1%
Taylor expanded in i around 0 88.5%
associate-*r/88.7%
*-commutative88.7%
Simplified88.5%
if -6.5e131 < n < 2.55000000000000009e-9Initial program 35.0%
associate-/r/35.2%
associate-*r*35.2%
*-commutative35.2%
associate-*r/35.2%
sub-neg35.2%
distribute-lft-in35.2%
metadata-eval35.2%
metadata-eval35.2%
metadata-eval35.2%
fma-define35.2%
metadata-eval35.2%
Simplified35.2%
*-commutative35.2%
fma-undefine35.2%
*-commutative35.2%
associate-/r/35.0%
clear-num35.0%
frac-2neg35.0%
metadata-eval35.0%
fma-define35.0%
Applied egg-rr35.0%
distribute-neg-frac235.0%
fma-define35.0%
+-commutative35.0%
distribute-neg-in35.0%
metadata-eval35.0%
*-commutative35.0%
distribute-lft-neg-in35.0%
metadata-eval35.0%
Simplified35.0%
Taylor expanded in i around 0 80.2%
if 2.55000000000000009e-9 < n < 1.1500000000000001e99Initial program 42.3%
associate-/r/42.5%
associate-*r*42.5%
*-commutative42.5%
associate-*r/42.5%
sub-neg42.5%
distribute-lft-in42.5%
metadata-eval42.5%
metadata-eval42.5%
metadata-eval42.5%
fma-define42.5%
metadata-eval42.5%
Simplified42.5%
Taylor expanded in i around 0 86.8%
*-commutative86.8%
associate-*r/86.8%
metadata-eval86.8%
Simplified86.8%
Final simplification84.2%
(FPCore (i n)
:precision binary64
(let* ((t_0 (- 0.5 (/ 0.5 n))))
(if (<= n -7.8e-66)
(* 100.0 (+ n (* t_0 (* i n))))
(if (<= n 2.5e-9)
(* 100.0 (/ i (/ i n)))
(* n (/ (* i (+ 100.0 (* 100.0 (* i t_0)))) i))))))
double code(double i, double n) {
double t_0 = 0.5 - (0.5 / n);
double tmp;
if (n <= -7.8e-66) {
tmp = 100.0 * (n + (t_0 * (i * n)));
} else if (n <= 2.5e-9) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = n * ((i * (100.0 + (100.0 * (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 = 0.5d0 - (0.5d0 / n)
if (n <= (-7.8d-66)) then
tmp = 100.0d0 * (n + (t_0 * (i * n)))
else if (n <= 2.5d-9) then
tmp = 100.0d0 * (i / (i / n))
else
tmp = n * ((i * (100.0d0 + (100.0d0 * (i * t_0)))) / i)
end if
code = tmp
end function
public static double code(double i, double n) {
double t_0 = 0.5 - (0.5 / n);
double tmp;
if (n <= -7.8e-66) {
tmp = 100.0 * (n + (t_0 * (i * n)));
} else if (n <= 2.5e-9) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = n * ((i * (100.0 + (100.0 * (i * t_0)))) / i);
}
return tmp;
}
def code(i, n): t_0 = 0.5 - (0.5 / n) tmp = 0 if n <= -7.8e-66: tmp = 100.0 * (n + (t_0 * (i * n))) elif n <= 2.5e-9: tmp = 100.0 * (i / (i / n)) else: tmp = n * ((i * (100.0 + (100.0 * (i * t_0)))) / i) return tmp
function code(i, n) t_0 = Float64(0.5 - Float64(0.5 / n)) tmp = 0.0 if (n <= -7.8e-66) tmp = Float64(100.0 * Float64(n + Float64(t_0 * Float64(i * n)))); elseif (n <= 2.5e-9) tmp = Float64(100.0 * Float64(i / Float64(i / n))); else tmp = Float64(n * Float64(Float64(i * Float64(100.0 + Float64(100.0 * Float64(i * t_0)))) / i)); end return tmp end
function tmp_2 = code(i, n) t_0 = 0.5 - (0.5 / n); tmp = 0.0; if (n <= -7.8e-66) tmp = 100.0 * (n + (t_0 * (i * n))); elseif (n <= 2.5e-9) tmp = 100.0 * (i / (i / n)); else tmp = n * ((i * (100.0 + (100.0 * (i * t_0)))) / i); end tmp_2 = tmp; end
code[i_, n_] := Block[{t$95$0 = N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -7.8e-66], N[(100.0 * N[(n + N[(t$95$0 * N[(i * n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 2.5e-9], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(n * N[(N[(i * N[(100.0 + N[(100.0 * N[(i * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.5 - \frac{0.5}{n}\\
\mathbf{if}\;n \leq -7.8 \cdot 10^{-66}:\\
\;\;\;\;100 \cdot \left(n + t\_0 \cdot \left(i \cdot n\right)\right)\\
\mathbf{elif}\;n \leq 2.5 \cdot 10^{-9}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;n \cdot \frac{i \cdot \left(100 + 100 \cdot \left(i \cdot t\_0\right)\right)}{i}\\
\end{array}
\end{array}
if n < -7.79999999999999965e-66Initial program 22.3%
Taylor expanded in i around 0 64.7%
associate-*r*64.7%
*-commutative64.7%
associate-*r/64.7%
metadata-eval64.7%
Simplified64.7%
if -7.79999999999999965e-66 < n < 2.5000000000000001e-9Initial program 39.2%
Taylor expanded in i around 0 40.2%
+-commutative40.2%
Simplified40.2%
Taylor expanded in i around 0 66.7%
if 2.5000000000000001e-9 < n Initial program 21.4%
associate-/r/22.0%
associate-*r*22.0%
*-commutative22.0%
associate-*r/22.0%
sub-neg22.0%
distribute-lft-in22.0%
metadata-eval22.0%
metadata-eval22.0%
metadata-eval22.0%
fma-define22.0%
metadata-eval22.0%
Simplified22.0%
Taylor expanded in i around 0 71.5%
*-commutative71.5%
associate-*r/71.5%
metadata-eval71.5%
Simplified71.5%
Final simplification67.4%
(FPCore (i n)
:precision binary64
(if (<= n 2.55e-9)
(/
-1.0
(+ (* 0.01 (/ (* i (+ 0.5 (* 0.5 (/ -1.0 n)))) n)) (* 0.01 (/ -1.0 n))))
(* n (/ (* i (+ 100.0 (* 100.0 (* i (- 0.5 (/ 0.5 n)))))) i))))
double code(double i, double n) {
double tmp;
if (n <= 2.55e-9) {
tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n)));
} else {
tmp = n * ((i * (100.0 + (100.0 * (i * (0.5 - (0.5 / 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.55d-9) then
tmp = (-1.0d0) / ((0.01d0 * ((i * (0.5d0 + (0.5d0 * ((-1.0d0) / n)))) / n)) + (0.01d0 * ((-1.0d0) / n)))
else
tmp = n * ((i * (100.0d0 + (100.0d0 * (i * (0.5d0 - (0.5d0 / n)))))) / i)
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (n <= 2.55e-9) {
tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n)));
} else {
tmp = n * ((i * (100.0 + (100.0 * (i * (0.5 - (0.5 / n)))))) / i);
}
return tmp;
}
def code(i, n): tmp = 0 if n <= 2.55e-9: tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n))) else: tmp = n * ((i * (100.0 + (100.0 * (i * (0.5 - (0.5 / n)))))) / i) return tmp
function code(i, n) tmp = 0.0 if (n <= 2.55e-9) tmp = Float64(-1.0 / Float64(Float64(0.01 * Float64(Float64(i * Float64(0.5 + Float64(0.5 * Float64(-1.0 / n)))) / n)) + Float64(0.01 * Float64(-1.0 / n)))); else tmp = Float64(n * Float64(Float64(i * Float64(100.0 + Float64(100.0 * Float64(i * Float64(0.5 - Float64(0.5 / n)))))) / i)); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (n <= 2.55e-9) tmp = -1.0 / ((0.01 * ((i * (0.5 + (0.5 * (-1.0 / n)))) / n)) + (0.01 * (-1.0 / n))); else tmp = n * ((i * (100.0 + (100.0 * (i * (0.5 - (0.5 / n)))))) / i); end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[n, 2.55e-9], N[(-1.0 / N[(N[(0.01 * N[(N[(i * N[(0.5 + N[(0.5 * N[(-1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 * N[(-1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(n * N[(N[(i * N[(100.0 + N[(100.0 * N[(i * N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq 2.55 \cdot 10^{-9}:\\
\;\;\;\;\frac{-1}{0.01 \cdot \frac{i \cdot \left(0.5 + 0.5 \cdot \frac{-1}{n}\right)}{n} + 0.01 \cdot \frac{-1}{n}}\\
\mathbf{else}:\\
\;\;\;\;n \cdot \frac{i \cdot \left(100 + 100 \cdot \left(i \cdot \left(0.5 - \frac{0.5}{n}\right)\right)\right)}{i}\\
\end{array}
\end{array}
if n < 2.55000000000000009e-9Initial program 29.5%
associate-/r/29.9%
associate-*r*29.9%
*-commutative29.9%
associate-*r/29.9%
sub-neg29.9%
distribute-lft-in29.9%
metadata-eval29.9%
metadata-eval29.9%
metadata-eval29.9%
fma-define29.9%
metadata-eval29.9%
Simplified29.9%
*-commutative29.9%
fma-undefine29.9%
*-commutative29.9%
associate-/r/29.5%
clear-num29.5%
frac-2neg29.5%
metadata-eval29.5%
fma-define29.5%
Applied egg-rr29.5%
distribute-neg-frac229.5%
fma-define29.5%
+-commutative29.5%
distribute-neg-in29.5%
metadata-eval29.5%
*-commutative29.5%
distribute-lft-neg-in29.5%
metadata-eval29.5%
Simplified29.5%
Taylor expanded in i around 0 75.5%
if 2.55000000000000009e-9 < n Initial program 21.4%
associate-/r/22.0%
associate-*r*22.0%
*-commutative22.0%
associate-*r/22.0%
sub-neg22.0%
distribute-lft-in22.0%
metadata-eval22.0%
metadata-eval22.0%
metadata-eval22.0%
fma-define22.0%
metadata-eval22.0%
Simplified22.0%
Taylor expanded in i around 0 71.5%
*-commutative71.5%
associate-*r/71.5%
metadata-eval71.5%
Simplified71.5%
Final simplification74.3%
(FPCore (i n) :precision binary64 (if (<= n 2.3e-9) (/ -1.0 (- (* 0.01 (* i (/ (+ 0.5 (/ -0.5 n)) n))) (/ 0.01 n))) (* n (/ (* i (+ 100.0 (* 100.0 (* i (- 0.5 (/ 0.5 n)))))) i))))
double code(double i, double n) {
double tmp;
if (n <= 2.3e-9) {
tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n));
} else {
tmp = n * ((i * (100.0 + (100.0 * (i * (0.5 - (0.5 / 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.3d-9) then
tmp = (-1.0d0) / ((0.01d0 * (i * ((0.5d0 + ((-0.5d0) / n)) / n))) - (0.01d0 / n))
else
tmp = n * ((i * (100.0d0 + (100.0d0 * (i * (0.5d0 - (0.5d0 / n)))))) / i)
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (n <= 2.3e-9) {
tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n));
} else {
tmp = n * ((i * (100.0 + (100.0 * (i * (0.5 - (0.5 / n)))))) / i);
}
return tmp;
}
def code(i, n): tmp = 0 if n <= 2.3e-9: tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n)) else: tmp = n * ((i * (100.0 + (100.0 * (i * (0.5 - (0.5 / n)))))) / i) return tmp
function code(i, n) tmp = 0.0 if (n <= 2.3e-9) tmp = Float64(-1.0 / Float64(Float64(0.01 * Float64(i * Float64(Float64(0.5 + Float64(-0.5 / n)) / n))) - Float64(0.01 / n))); else tmp = Float64(n * Float64(Float64(i * Float64(100.0 + Float64(100.0 * Float64(i * Float64(0.5 - Float64(0.5 / n)))))) / i)); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (n <= 2.3e-9) tmp = -1.0 / ((0.01 * (i * ((0.5 + (-0.5 / n)) / n))) - (0.01 / n)); else tmp = n * ((i * (100.0 + (100.0 * (i * (0.5 - (0.5 / n)))))) / i); end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[n, 2.3e-9], N[(-1.0 / N[(N[(0.01 * N[(i * N[(N[(0.5 + N[(-0.5 / n), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(0.01 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(n * N[(N[(i * N[(100.0 + N[(100.0 * N[(i * N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq 2.3 \cdot 10^{-9}:\\
\;\;\;\;\frac{-1}{0.01 \cdot \left(i \cdot \frac{0.5 + \frac{-0.5}{n}}{n}\right) - \frac{0.01}{n}}\\
\mathbf{else}:\\
\;\;\;\;n \cdot \frac{i \cdot \left(100 + 100 \cdot \left(i \cdot \left(0.5 - \frac{0.5}{n}\right)\right)\right)}{i}\\
\end{array}
\end{array}
if n < 2.2999999999999999e-9Initial program 29.5%
associate-/r/29.9%
associate-*r*29.9%
*-commutative29.9%
associate-*r/29.9%
sub-neg29.9%
distribute-lft-in29.9%
metadata-eval29.9%
metadata-eval29.9%
metadata-eval29.9%
fma-define29.9%
metadata-eval29.9%
Simplified29.9%
*-commutative29.9%
fma-undefine29.9%
*-commutative29.9%
associate-/r/29.5%
clear-num29.5%
frac-2neg29.5%
metadata-eval29.5%
fma-define29.5%
Applied egg-rr29.5%
distribute-neg-frac229.5%
fma-define29.5%
+-commutative29.5%
distribute-neg-in29.5%
metadata-eval29.5%
*-commutative29.5%
distribute-lft-neg-in29.5%
metadata-eval29.5%
Simplified29.5%
Taylor expanded in i around 0 75.5%
associate-/l*71.6%
cancel-sign-sub-inv71.6%
metadata-eval71.6%
associate-*r/71.6%
metadata-eval71.6%
associate-*r/71.7%
metadata-eval71.7%
Simplified71.7%
if 2.2999999999999999e-9 < n Initial program 21.4%
associate-/r/22.0%
associate-*r*22.0%
*-commutative22.0%
associate-*r/22.0%
sub-neg22.0%
distribute-lft-in22.0%
metadata-eval22.0%
metadata-eval22.0%
metadata-eval22.0%
fma-define22.0%
metadata-eval22.0%
Simplified22.0%
Taylor expanded in i around 0 71.5%
*-commutative71.5%
associate-*r/71.5%
metadata-eval71.5%
Simplified71.5%
Final simplification71.6%
(FPCore (i n) :precision binary64 (if (<= n -6.6e-76) (* 100.0 (+ n (* (- 0.5 (/ 0.5 n)) (* i n)))) (if (<= n 8.2e-38) (* 100.0 (/ i (/ i n))) (* 100.0 (/ (* i n) i)))))
double code(double i, double n) {
double tmp;
if (n <= -6.6e-76) {
tmp = 100.0 * (n + ((0.5 - (0.5 / n)) * (i * n)));
} else if (n <= 8.2e-38) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = 100.0 * ((i * 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 <= (-6.6d-76)) then
tmp = 100.0d0 * (n + ((0.5d0 - (0.5d0 / n)) * (i * n)))
else if (n <= 8.2d-38) then
tmp = 100.0d0 * (i / (i / n))
else
tmp = 100.0d0 * ((i * n) / i)
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (n <= -6.6e-76) {
tmp = 100.0 * (n + ((0.5 - (0.5 / n)) * (i * n)));
} else if (n <= 8.2e-38) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = 100.0 * ((i * n) / i);
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -6.6e-76: tmp = 100.0 * (n + ((0.5 - (0.5 / n)) * (i * n))) elif n <= 8.2e-38: tmp = 100.0 * (i / (i / n)) else: tmp = 100.0 * ((i * n) / i) return tmp
function code(i, n) tmp = 0.0 if (n <= -6.6e-76) tmp = Float64(100.0 * Float64(n + Float64(Float64(0.5 - Float64(0.5 / n)) * Float64(i * n)))); elseif (n <= 8.2e-38) tmp = Float64(100.0 * Float64(i / Float64(i / n))); else tmp = Float64(100.0 * Float64(Float64(i * n) / i)); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (n <= -6.6e-76) tmp = 100.0 * (n + ((0.5 - (0.5 / n)) * (i * n))); elseif (n <= 8.2e-38) tmp = 100.0 * (i / (i / n)); else tmp = 100.0 * ((i * n) / i); end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[n, -6.6e-76], 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, 8.2e-38], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(100.0 * N[(N[(i * n), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -6.6 \cdot 10^{-76}:\\
\;\;\;\;100 \cdot \left(n + \left(0.5 - \frac{0.5}{n}\right) \cdot \left(i \cdot n\right)\right)\\
\mathbf{elif}\;n \leq 8.2 \cdot 10^{-38}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \frac{i \cdot n}{i}\\
\end{array}
\end{array}
if n < -6.59999999999999967e-76Initial program 22.3%
Taylor expanded in i around 0 64.7%
associate-*r*64.7%
*-commutative64.7%
associate-*r/64.7%
metadata-eval64.7%
Simplified64.7%
if -6.59999999999999967e-76 < n < 8.1999999999999996e-38Initial program 41.9%
Taylor expanded in i around 0 43.9%
+-commutative43.9%
Simplified43.9%
Taylor expanded in i around 0 67.4%
if 8.1999999999999996e-38 < n Initial program 20.7%
Taylor expanded in i around 0 3.6%
+-commutative3.6%
Simplified3.6%
clear-num3.6%
inv-pow3.6%
Applied egg-rr3.6%
unpow-13.6%
Simplified3.6%
clear-num3.6%
div-inv3.6%
add-exp-log2.3%
expm1-define2.3%
+-commutative2.3%
log1p-define27.4%
expm1-log1p-u28.6%
clear-num28.7%
Applied egg-rr28.7%
associate-*r/67.0%
Applied egg-rr67.0%
Final simplification66.2%
(FPCore (i n) :precision binary64 (if (or (<= n -1e+21) (not (<= n 8.5e-41))) (* 100.0 (/ (* i n) i)) (* 100.0 (/ i (/ i n)))))
double code(double i, double n) {
double tmp;
if ((n <= -1e+21) || !(n <= 8.5e-41)) {
tmp = 100.0 * ((i * n) / 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 <= (-1d+21)) .or. (.not. (n <= 8.5d-41))) then
tmp = 100.0d0 * ((i * n) / 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 <= -1e+21) || !(n <= 8.5e-41)) {
tmp = 100.0 * ((i * n) / i);
} else {
tmp = 100.0 * (i / (i / n));
}
return tmp;
}
def code(i, n): tmp = 0 if (n <= -1e+21) or not (n <= 8.5e-41): tmp = 100.0 * ((i * n) / i) else: tmp = 100.0 * (i / (i / n)) return tmp
function code(i, n) tmp = 0.0 if ((n <= -1e+21) || !(n <= 8.5e-41)) tmp = Float64(100.0 * Float64(Float64(i * n) / 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 <= -1e+21) || ~((n <= 8.5e-41))) tmp = 100.0 * ((i * n) / i); else tmp = 100.0 * (i / (i / n)); end tmp_2 = tmp; end
code[i_, n_] := If[Or[LessEqual[n, -1e+21], N[Not[LessEqual[n, 8.5e-41]], $MachinePrecision]], N[(100.0 * N[(N[(i * n), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1 \cdot 10^{+21} \lor \neg \left(n \leq 8.5 \cdot 10^{-41}\right):\\
\;\;\;\;100 \cdot \frac{i \cdot n}{i}\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\end{array}
\end{array}
if n < -1e21 or 8.4999999999999996e-41 < n Initial program 22.1%
Taylor expanded in i around 0 3.7%
+-commutative3.7%
Simplified3.7%
clear-num3.6%
inv-pow3.6%
Applied egg-rr3.6%
unpow-13.6%
Simplified3.6%
clear-num3.7%
div-inv3.6%
add-exp-log2.4%
expm1-define2.4%
+-commutative2.4%
log1p-define28.1%
expm1-log1p-u29.3%
clear-num29.4%
Applied egg-rr29.4%
associate-*r/66.0%
Applied egg-rr66.0%
if -1e21 < n < 8.4999999999999996e-41Initial program 37.0%
Taylor expanded in i around 0 36.3%
+-commutative36.3%
Simplified36.3%
Taylor expanded in i around 0 66.6%
Final simplification66.2%
(FPCore (i n) :precision binary64 (if (or (<= i -9e+151) (not (<= i 5e-50))) (* 100.0 (* i (/ n i))) (* n 100.0)))
double code(double i, double n) {
double tmp;
if ((i <= -9e+151) || !(i <= 5e-50)) {
tmp = 100.0 * (i * (n / i));
} else {
tmp = n * 100.0;
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if ((i <= (-9d+151)) .or. (.not. (i <= 5d-50))) then
tmp = 100.0d0 * (i * (n / i))
else
tmp = n * 100.0d0
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if ((i <= -9e+151) || !(i <= 5e-50)) {
tmp = 100.0 * (i * (n / i));
} else {
tmp = n * 100.0;
}
return tmp;
}
def code(i, n): tmp = 0 if (i <= -9e+151) or not (i <= 5e-50): tmp = 100.0 * (i * (n / i)) else: tmp = n * 100.0 return tmp
function code(i, n) tmp = 0.0 if ((i <= -9e+151) || !(i <= 5e-50)) tmp = Float64(100.0 * Float64(i * Float64(n / i))); else tmp = Float64(n * 100.0); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if ((i <= -9e+151) || ~((i <= 5e-50))) tmp = 100.0 * (i * (n / i)); else tmp = n * 100.0; end tmp_2 = tmp; end
code[i_, n_] := If[Or[LessEqual[i, -9e+151], N[Not[LessEqual[i, 5e-50]], $MachinePrecision]], N[(100.0 * N[(i * N[(n / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(n * 100.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -9 \cdot 10^{+151} \lor \neg \left(i \leq 5 \cdot 10^{-50}\right):\\
\;\;\;\;100 \cdot \left(i \cdot \frac{n}{i}\right)\\
\mathbf{else}:\\
\;\;\;\;n \cdot 100\\
\end{array}
\end{array}
if i < -8.9999999999999997e151 or 4.99999999999999968e-50 < i Initial program 56.5%
Taylor expanded in i around 0 26.5%
+-commutative26.5%
Simplified26.5%
clear-num26.5%
inv-pow26.5%
Applied egg-rr26.5%
unpow-126.5%
Simplified26.5%
clear-num26.5%
div-inv26.5%
add-exp-log13.6%
expm1-define13.6%
+-commutative13.6%
log1p-define20.8%
expm1-log1p-u33.7%
clear-num32.7%
Applied egg-rr32.7%
if -8.9999999999999997e151 < i < 4.99999999999999968e-50Initial program 8.7%
Taylor expanded in i around 0 7.1%
+-commutative7.1%
Simplified7.1%
Taylor expanded in i around 0 75.9%
*-commutative75.9%
Simplified75.9%
Final simplification59.4%
(FPCore (i n) :precision binary64 (if (<= n -1.7e+21) (* 100.0 (/ 1.0 (* i (/ (/ 1.0 n) i)))) (if (<= n 5e-47) (* 100.0 (/ i (/ i n))) (* 100.0 (/ (* i n) i)))))
double code(double i, double n) {
double tmp;
if (n <= -1.7e+21) {
tmp = 100.0 * (1.0 / (i * ((1.0 / n) / i)));
} else if (n <= 5e-47) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = 100.0 * ((i * 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 <= (-1.7d+21)) then
tmp = 100.0d0 * (1.0d0 / (i * ((1.0d0 / n) / i)))
else if (n <= 5d-47) then
tmp = 100.0d0 * (i / (i / n))
else
tmp = 100.0d0 * ((i * n) / i)
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (n <= -1.7e+21) {
tmp = 100.0 * (1.0 / (i * ((1.0 / n) / i)));
} else if (n <= 5e-47) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = 100.0 * ((i * n) / i);
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -1.7e+21: tmp = 100.0 * (1.0 / (i * ((1.0 / n) / i))) elif n <= 5e-47: tmp = 100.0 * (i / (i / n)) else: tmp = 100.0 * ((i * n) / i) return tmp
function code(i, n) tmp = 0.0 if (n <= -1.7e+21) tmp = Float64(100.0 * Float64(1.0 / Float64(i * Float64(Float64(1.0 / n) / i)))); elseif (n <= 5e-47) tmp = Float64(100.0 * Float64(i / Float64(i / n))); else tmp = Float64(100.0 * Float64(Float64(i * n) / i)); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (n <= -1.7e+21) tmp = 100.0 * (1.0 / (i * ((1.0 / n) / i))); elseif (n <= 5e-47) tmp = 100.0 * (i / (i / n)); else tmp = 100.0 * ((i * n) / i); end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[n, -1.7e+21], N[(100.0 * N[(1.0 / N[(i * N[(N[(1.0 / n), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 5e-47], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(100.0 * N[(N[(i * n), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.7 \cdot 10^{+21}:\\
\;\;\;\;100 \cdot \frac{1}{i \cdot \frac{\frac{1}{n}}{i}}\\
\mathbf{elif}\;n \leq 5 \cdot 10^{-47}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \frac{i \cdot n}{i}\\
\end{array}
\end{array}
if n < -1.7e21Initial program 23.4%
Taylor expanded in i around 0 3.8%
+-commutative3.8%
Simplified3.8%
clear-num3.6%
inv-pow3.6%
Applied egg-rr3.6%
unpow-13.6%
Simplified3.6%
clear-num3.8%
div-inv3.6%
add-exp-log2.6%
expm1-define2.6%
+-commutative2.6%
log1p-define28.8%
expm1-log1p-u29.9%
clear-num30.0%
Applied egg-rr30.0%
associate-*r/64.9%
pow164.9%
metadata-eval64.9%
pow-flip64.9%
inv-pow64.9%
div-inv64.7%
*-un-lft-identity64.7%
associate-*l/64.7%
*-commutative64.7%
clear-num64.6%
frac-times65.0%
metadata-eval65.0%
Applied egg-rr65.0%
if -1.7e21 < n < 5.00000000000000011e-47Initial program 37.0%
Taylor expanded in i around 0 36.3%
+-commutative36.3%
Simplified36.3%
Taylor expanded in i around 0 66.6%
if 5.00000000000000011e-47 < n Initial program 20.7%
Taylor expanded in i around 0 3.6%
+-commutative3.6%
Simplified3.6%
clear-num3.6%
inv-pow3.6%
Applied egg-rr3.6%
unpow-13.6%
Simplified3.6%
clear-num3.6%
div-inv3.6%
add-exp-log2.3%
expm1-define2.3%
+-commutative2.3%
log1p-define27.4%
expm1-log1p-u28.6%
clear-num28.7%
Applied egg-rr28.7%
associate-*r/67.0%
Applied egg-rr67.0%
Final simplification66.2%
(FPCore (i n) :precision binary64 (if (<= i -2e+154) (* 100.0 (* i (/ n i))) (if (<= i 1e-49) (* n 100.0) (* 100.0 (/ i (/ i n))))))
double code(double i, double n) {
double tmp;
if (i <= -2e+154) {
tmp = 100.0 * (i * (n / i));
} else if (i <= 1e-49) {
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 (i <= (-2d+154)) then
tmp = 100.0d0 * (i * (n / i))
else if (i <= 1d-49) 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 (i <= -2e+154) {
tmp = 100.0 * (i * (n / i));
} else if (i <= 1e-49) {
tmp = n * 100.0;
} else {
tmp = 100.0 * (i / (i / n));
}
return tmp;
}
def code(i, n): tmp = 0 if i <= -2e+154: tmp = 100.0 * (i * (n / i)) elif i <= 1e-49: tmp = n * 100.0 else: tmp = 100.0 * (i / (i / n)) return tmp
function code(i, n) tmp = 0.0 if (i <= -2e+154) tmp = Float64(100.0 * Float64(i * Float64(n / i))); elseif (i <= 1e-49) 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 (i <= -2e+154) tmp = 100.0 * (i * (n / i)); elseif (i <= 1e-49) tmp = n * 100.0; else tmp = 100.0 * (i / (i / n)); end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[i, -2e+154], N[(100.0 * N[(i * N[(n / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[i, 1e-49], N[(n * 100.0), $MachinePrecision], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -2 \cdot 10^{+154}:\\
\;\;\;\;100 \cdot \left(i \cdot \frac{n}{i}\right)\\
\mathbf{elif}\;i \leq 10^{-49}:\\
\;\;\;\;n \cdot 100\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\end{array}
\end{array}
if i < -2.00000000000000007e154Initial program 79.3%
Taylor expanded in i around 0 45.1%
+-commutative45.1%
Simplified45.1%
clear-num45.1%
inv-pow45.1%
Applied egg-rr45.1%
unpow-145.1%
Simplified45.1%
clear-num45.1%
div-inv45.1%
add-exp-log0.0%
expm1-define0.0%
+-commutative0.0%
log1p-define0.0%
expm1-log1p-u45.1%
clear-num45.1%
Applied egg-rr45.1%
if -2.00000000000000007e154 < i < 9.99999999999999936e-50Initial program 8.7%
Taylor expanded in i around 0 7.0%
+-commutative7.0%
Simplified7.0%
Taylor expanded in i around 0 76.0%
*-commutative76.0%
Simplified76.0%
if 9.99999999999999936e-50 < i Initial program 48.1%
Taylor expanded in i around 0 19.3%
+-commutative19.3%
Simplified19.3%
Taylor expanded in i around 0 28.1%
(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 27.0%
Taylor expanded in i around 0 14.5%
+-commutative14.5%
Simplified14.5%
Taylor expanded in i around 0 51.7%
*-commutative51.7%
Simplified51.7%
(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 2024179
(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))))