
(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 (if (<= (/ (+ (pow (+ 1.0 (/ i n)) n) -1.0) (/ i n)) 0.0) (* n (* 100.0 (/ (expm1 (* n (log1p (/ i n)))) i))) (* n (/ (* 100.0 (expm1 i)) i))))
double code(double i, double n) {
double tmp;
if (((pow((1.0 + (i / n)), n) + -1.0) / (i / n)) <= 0.0) {
tmp = n * (100.0 * (expm1((n * log1p((i / n)))) / i));
} else {
tmp = n * ((100.0 * expm1(i)) / i);
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if (((Math.pow((1.0 + (i / n)), n) + -1.0) / (i / n)) <= 0.0) {
tmp = n * (100.0 * (Math.expm1((n * Math.log1p((i / n)))) / i));
} else {
tmp = n * ((100.0 * Math.expm1(i)) / i);
}
return tmp;
}
def code(i, n): tmp = 0 if ((math.pow((1.0 + (i / n)), n) + -1.0) / (i / n)) <= 0.0: tmp = n * (100.0 * (math.expm1((n * math.log1p((i / n)))) / i)) else: tmp = n * ((100.0 * math.expm1(i)) / i) return tmp
function code(i, n) tmp = 0.0 if (Float64(Float64((Float64(1.0 + Float64(i / n)) ^ n) + -1.0) / Float64(i / n)) <= 0.0) tmp = Float64(n * Float64(100.0 * Float64(expm1(Float64(n * log1p(Float64(i / n)))) / i))); else tmp = Float64(n * Float64(Float64(100.0 * expm1(i)) / i)); end return tmp end
code[i_, n_] := If[LessEqual[N[(N[(N[Power[N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision], n], $MachinePrecision] + -1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision], 0.0], N[(n * N[(100.0 * N[(N[(Exp[N[(n * N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]] - 1), $MachinePrecision] / i), $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}\;\frac{{\left(1 + \frac{i}{n}\right)}^{n} + -1}{\frac{i}{n}} \leq 0:\\
\;\;\;\;n \cdot \left(100 \cdot \frac{\mathsf{expm1}\left(n \cdot \mathsf{log1p}\left(\frac{i}{n}\right)\right)}{i}\right)\\
\mathbf{else}:\\
\;\;\;\;n \cdot \frac{100 \cdot \mathsf{expm1}\left(i\right)}{i}\\
\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 29.7%
associate-*r/29.7%
sub-neg29.7%
distribute-rgt-in29.7%
metadata-eval29.7%
metadata-eval29.7%
Simplified29.7%
metadata-eval29.7%
metadata-eval29.7%
distribute-rgt-in29.7%
sub-neg29.7%
associate-*r/29.7%
associate-/r/29.2%
associate-*r*29.2%
add-exp-log29.2%
expm1-define29.2%
log-pow36.1%
log1p-define96.7%
Applied egg-rr96.7%
if -0.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 26.2%
associate-/r/27.5%
associate-*r*27.5%
*-commutative27.5%
associate-*r/27.5%
sub-neg27.5%
distribute-lft-in27.5%
metadata-eval27.5%
metadata-eval27.5%
metadata-eval27.5%
fma-define27.5%
metadata-eval27.5%
Simplified27.5%
Taylor expanded in n around inf 21.1%
sub-neg21.1%
metadata-eval21.1%
metadata-eval21.1%
distribute-lft-in21.1%
metadata-eval21.1%
sub-neg21.1%
expm1-define72.5%
Simplified72.5%
Final simplification90.9%
(FPCore (i n) :precision binary64 (if (<= (/ (+ (pow (+ 1.0 (/ i n)) n) -1.0) (/ i n)) 0.0) (* 100.0 (* n (/ (expm1 (* n (log1p (/ i n)))) i))) (* n (/ (* 100.0 (expm1 i)) i))))
double code(double i, double n) {
double tmp;
if (((pow((1.0 + (i / n)), n) + -1.0) / (i / n)) <= 0.0) {
tmp = 100.0 * (n * (expm1((n * log1p((i / n)))) / i));
} else {
tmp = n * ((100.0 * expm1(i)) / i);
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if (((Math.pow((1.0 + (i / n)), n) + -1.0) / (i / n)) <= 0.0) {
tmp = 100.0 * (n * (Math.expm1((n * Math.log1p((i / n)))) / i));
} else {
tmp = n * ((100.0 * Math.expm1(i)) / i);
}
return tmp;
}
def code(i, n): tmp = 0 if ((math.pow((1.0 + (i / n)), n) + -1.0) / (i / n)) <= 0.0: tmp = 100.0 * (n * (math.expm1((n * math.log1p((i / n)))) / i)) else: tmp = n * ((100.0 * math.expm1(i)) / i) return tmp
function code(i, n) tmp = 0.0 if (Float64(Float64((Float64(1.0 + Float64(i / n)) ^ n) + -1.0) / Float64(i / n)) <= 0.0) tmp = Float64(100.0 * Float64(n * Float64(expm1(Float64(n * log1p(Float64(i / n)))) / i))); else tmp = Float64(n * Float64(Float64(100.0 * expm1(i)) / i)); end return tmp end
code[i_, n_] := If[LessEqual[N[(N[(N[Power[N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision], n], $MachinePrecision] + -1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision], 0.0], N[(100.0 * N[(n * N[(N[(Exp[N[(n * N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]] - 1), $MachinePrecision] / i), $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}\;\frac{{\left(1 + \frac{i}{n}\right)}^{n} + -1}{\frac{i}{n}} \leq 0:\\
\;\;\;\;100 \cdot \left(n \cdot \frac{\mathsf{expm1}\left(n \cdot \mathsf{log1p}\left(\frac{i}{n}\right)\right)}{i}\right)\\
\mathbf{else}:\\
\;\;\;\;n \cdot \frac{100 \cdot \mathsf{expm1}\left(i\right)}{i}\\
\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 29.7%
associate-/r/29.2%
add-exp-log29.2%
expm1-define29.2%
log-pow36.0%
log1p-define96.7%
Applied egg-rr96.7%
if -0.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 26.2%
associate-/r/27.5%
associate-*r*27.5%
*-commutative27.5%
associate-*r/27.5%
sub-neg27.5%
distribute-lft-in27.5%
metadata-eval27.5%
metadata-eval27.5%
metadata-eval27.5%
fma-define27.5%
metadata-eval27.5%
Simplified27.5%
Taylor expanded in n around inf 21.1%
sub-neg21.1%
metadata-eval21.1%
metadata-eval21.1%
distribute-lft-in21.1%
metadata-eval21.1%
sub-neg21.1%
expm1-define72.5%
Simplified72.5%
Final simplification90.9%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* 100.0 (/ (* n (expm1 i)) i))))
(if (<= n -2e+28)
t_0
(if (<= n -1.95e-145)
(* 100.0 (* (expm1 i) (/ n i)))
(if (<= n 2.45e-234)
0.0
(if (<= n 7.5e-41)
(* 100.0 (/ i (/ i n)))
(if (<= n 2.8e+37)
(+
(* n 100.0)
(*
i
(*
n
(+
50.0
(* i (+ 16.666666666666668 (* i 4.166666666666667)))))))
t_0)))))))
double code(double i, double n) {
double t_0 = 100.0 * ((n * expm1(i)) / i);
double tmp;
if (n <= -2e+28) {
tmp = t_0;
} else if (n <= -1.95e-145) {
tmp = 100.0 * (expm1(i) * (n / i));
} else if (n <= 2.45e-234) {
tmp = 0.0;
} else if (n <= 7.5e-41) {
tmp = 100.0 * (i / (i / n));
} else if (n <= 2.8e+37) {
tmp = (n * 100.0) + (i * (n * (50.0 + (i * (16.666666666666668 + (i * 4.166666666666667))))));
} else {
tmp = t_0;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = 100.0 * ((n * Math.expm1(i)) / i);
double tmp;
if (n <= -2e+28) {
tmp = t_0;
} else if (n <= -1.95e-145) {
tmp = 100.0 * (Math.expm1(i) * (n / i));
} else if (n <= 2.45e-234) {
tmp = 0.0;
} else if (n <= 7.5e-41) {
tmp = 100.0 * (i / (i / n));
} else if (n <= 2.8e+37) {
tmp = (n * 100.0) + (i * (n * (50.0 + (i * (16.666666666666668 + (i * 4.166666666666667))))));
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = 100.0 * ((n * math.expm1(i)) / i) tmp = 0 if n <= -2e+28: tmp = t_0 elif n <= -1.95e-145: tmp = 100.0 * (math.expm1(i) * (n / i)) elif n <= 2.45e-234: tmp = 0.0 elif n <= 7.5e-41: tmp = 100.0 * (i / (i / n)) elif n <= 2.8e+37: tmp = (n * 100.0) + (i * (n * (50.0 + (i * (16.666666666666668 + (i * 4.166666666666667)))))) else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(100.0 * Float64(Float64(n * expm1(i)) / i)) tmp = 0.0 if (n <= -2e+28) tmp = t_0; elseif (n <= -1.95e-145) tmp = Float64(100.0 * Float64(expm1(i) * Float64(n / i))); elseif (n <= 2.45e-234) tmp = 0.0; elseif (n <= 7.5e-41) tmp = Float64(100.0 * Float64(i / Float64(i / n))); elseif (n <= 2.8e+37) tmp = Float64(Float64(n * 100.0) + Float64(i * Float64(n * Float64(50.0 + Float64(i * Float64(16.666666666666668 + Float64(i * 4.166666666666667))))))); else tmp = t_0; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(100.0 * N[(N[(n * N[(Exp[i] - 1), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -2e+28], t$95$0, If[LessEqual[n, -1.95e-145], N[(100.0 * N[(N[(Exp[i] - 1), $MachinePrecision] * N[(n / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 2.45e-234], 0.0, If[LessEqual[n, 7.5e-41], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 2.8e+37], N[(N[(n * 100.0), $MachinePrecision] + N[(i * N[(n * N[(50.0 + N[(i * N[(16.666666666666668 + N[(i * 4.166666666666667), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 100 \cdot \frac{n \cdot \mathsf{expm1}\left(i\right)}{i}\\
\mathbf{if}\;n \leq -2 \cdot 10^{+28}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq -1.95 \cdot 10^{-145}:\\
\;\;\;\;100 \cdot \left(\mathsf{expm1}\left(i\right) \cdot \frac{n}{i}\right)\\
\mathbf{elif}\;n \leq 2.45 \cdot 10^{-234}:\\
\;\;\;\;0\\
\mathbf{elif}\;n \leq 7.5 \cdot 10^{-41}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 2.8 \cdot 10^{+37}:\\
\;\;\;\;n \cdot 100 + i \cdot \left(n \cdot \left(50 + i \cdot \left(16.666666666666668 + i \cdot 4.166666666666667\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -1.99999999999999992e28 or 2.7999999999999998e37 < n Initial program 21.9%
Taylor expanded in n around inf 49.6%
expm1-define73.4%
Simplified73.4%
clear-num73.3%
un-div-inv73.4%
associate-/l/95.9%
Applied egg-rr95.9%
/-rgt-identity95.9%
associate-/r/95.9%
associate-/r/95.9%
associate-/r/96.0%
/-rgt-identity96.0%
*-commutative96.0%
Simplified96.0%
Taylor expanded in i around inf 50.1%
associate-*r/49.4%
expm1-define95.2%
associate-*r*95.4%
Simplified95.4%
Taylor expanded in n around 0 50.1%
expm1-define96.1%
Simplified96.1%
if -1.99999999999999992e28 < n < -1.95000000000000015e-145Initial program 15.4%
associate-/r/15.4%
add-exp-log15.4%
expm1-define15.4%
log-pow38.6%
log1p-define99.4%
Applied egg-rr99.4%
Taylor expanded in n around inf 8.5%
expm1-define56.0%
associate-*l/70.5%
Simplified70.5%
if -1.95000000000000015e-145 < n < 2.45000000000000004e-234Initial program 73.2%
associate-*r/73.2%
sub-neg73.2%
distribute-rgt-in73.2%
metadata-eval73.2%
metadata-eval73.2%
Simplified73.2%
Taylor expanded in i around 0 80.9%
Taylor expanded in i around 0 80.9%
if 2.45000000000000004e-234 < n < 7.50000000000000049e-41Initial program 23.2%
Taylor expanded in i around 0 66.3%
if 7.50000000000000049e-41 < n < 2.7999999999999998e37Initial program 17.6%
associate-*r/17.6%
sub-neg17.6%
distribute-rgt-in17.6%
metadata-eval17.6%
metadata-eval17.6%
Simplified17.6%
Taylor expanded in n around inf 15.6%
Taylor expanded in i around 0 95.8%
Taylor expanded in n around 0 95.8%
Final simplification87.0%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* 100.0 (/ (expm1 i) (/ i n))))
(t_1 (/ (* i (* n (+ 100.0 (* i 50.0)))) i)))
(if (<= i -5.2e-8)
t_0
(if (<= i -2.35e-177)
t_1
(if (<= i 2.9e-227)
(+ (* n 100.0) (* i -50.0))
(if (<= i 1.9e-53)
t_1
(if (<= i 1.22e+117)
t_0
(if (<= i 1.06e+158)
0.0
(if (<= i 2.2e+204)
(/ (* i (+ (* n 100.0) (* 50.0 (* i n)))) i)
(if (<= i 1.5e+273) 0.0 (* (* i n) (/ 100.0 i))))))))))))
double code(double i, double n) {
double t_0 = 100.0 * (expm1(i) / (i / n));
double t_1 = (i * (n * (100.0 + (i * 50.0)))) / i;
double tmp;
if (i <= -5.2e-8) {
tmp = t_0;
} else if (i <= -2.35e-177) {
tmp = t_1;
} else if (i <= 2.9e-227) {
tmp = (n * 100.0) + (i * -50.0);
} else if (i <= 1.9e-53) {
tmp = t_1;
} else if (i <= 1.22e+117) {
tmp = t_0;
} else if (i <= 1.06e+158) {
tmp = 0.0;
} else if (i <= 2.2e+204) {
tmp = (i * ((n * 100.0) + (50.0 * (i * n)))) / i;
} else if (i <= 1.5e+273) {
tmp = 0.0;
} else {
tmp = (i * n) * (100.0 / i);
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = 100.0 * (Math.expm1(i) / (i / n));
double t_1 = (i * (n * (100.0 + (i * 50.0)))) / i;
double tmp;
if (i <= -5.2e-8) {
tmp = t_0;
} else if (i <= -2.35e-177) {
tmp = t_1;
} else if (i <= 2.9e-227) {
tmp = (n * 100.0) + (i * -50.0);
} else if (i <= 1.9e-53) {
tmp = t_1;
} else if (i <= 1.22e+117) {
tmp = t_0;
} else if (i <= 1.06e+158) {
tmp = 0.0;
} else if (i <= 2.2e+204) {
tmp = (i * ((n * 100.0) + (50.0 * (i * n)))) / i;
} else if (i <= 1.5e+273) {
tmp = 0.0;
} else {
tmp = (i * n) * (100.0 / i);
}
return tmp;
}
def code(i, n): t_0 = 100.0 * (math.expm1(i) / (i / n)) t_1 = (i * (n * (100.0 + (i * 50.0)))) / i tmp = 0 if i <= -5.2e-8: tmp = t_0 elif i <= -2.35e-177: tmp = t_1 elif i <= 2.9e-227: tmp = (n * 100.0) + (i * -50.0) elif i <= 1.9e-53: tmp = t_1 elif i <= 1.22e+117: tmp = t_0 elif i <= 1.06e+158: tmp = 0.0 elif i <= 2.2e+204: tmp = (i * ((n * 100.0) + (50.0 * (i * n)))) / i elif i <= 1.5e+273: tmp = 0.0 else: tmp = (i * n) * (100.0 / i) return tmp
function code(i, n) t_0 = Float64(100.0 * Float64(expm1(i) / Float64(i / n))) t_1 = Float64(Float64(i * Float64(n * Float64(100.0 + Float64(i * 50.0)))) / i) tmp = 0.0 if (i <= -5.2e-8) tmp = t_0; elseif (i <= -2.35e-177) tmp = t_1; elseif (i <= 2.9e-227) tmp = Float64(Float64(n * 100.0) + Float64(i * -50.0)); elseif (i <= 1.9e-53) tmp = t_1; elseif (i <= 1.22e+117) tmp = t_0; elseif (i <= 1.06e+158) tmp = 0.0; elseif (i <= 2.2e+204) tmp = Float64(Float64(i * Float64(Float64(n * 100.0) + Float64(50.0 * Float64(i * n)))) / i); elseif (i <= 1.5e+273) tmp = 0.0; else tmp = Float64(Float64(i * n) * Float64(100.0 / i)); 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]}, Block[{t$95$1 = N[(N[(i * N[(n * N[(100.0 + N[(i * 50.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]}, If[LessEqual[i, -5.2e-8], t$95$0, If[LessEqual[i, -2.35e-177], t$95$1, If[LessEqual[i, 2.9e-227], N[(N[(n * 100.0), $MachinePrecision] + N[(i * -50.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[i, 1.9e-53], t$95$1, If[LessEqual[i, 1.22e+117], t$95$0, If[LessEqual[i, 1.06e+158], 0.0, If[LessEqual[i, 2.2e+204], N[(N[(i * N[(N[(n * 100.0), $MachinePrecision] + N[(50.0 * N[(i * n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision], If[LessEqual[i, 1.5e+273], 0.0, N[(N[(i * n), $MachinePrecision] * N[(100.0 / i), $MachinePrecision]), $MachinePrecision]]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 100 \cdot \frac{\mathsf{expm1}\left(i\right)}{\frac{i}{n}}\\
t_1 := \frac{i \cdot \left(n \cdot \left(100 + i \cdot 50\right)\right)}{i}\\
\mathbf{if}\;i \leq -5.2 \cdot 10^{-8}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;i \leq -2.35 \cdot 10^{-177}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;i \leq 2.9 \cdot 10^{-227}:\\
\;\;\;\;n \cdot 100 + i \cdot -50\\
\mathbf{elif}\;i \leq 1.9 \cdot 10^{-53}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;i \leq 1.22 \cdot 10^{+117}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;i \leq 1.06 \cdot 10^{+158}:\\
\;\;\;\;0\\
\mathbf{elif}\;i \leq 2.2 \cdot 10^{+204}:\\
\;\;\;\;\frac{i \cdot \left(n \cdot 100 + 50 \cdot \left(i \cdot n\right)\right)}{i}\\
\mathbf{elif}\;i \leq 1.5 \cdot 10^{+273}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;\left(i \cdot n\right) \cdot \frac{100}{i}\\
\end{array}
\end{array}
if i < -5.2000000000000002e-8 or 1.8999999999999999e-53 < i < 1.22000000000000004e117Initial program 43.4%
Taylor expanded in n around inf 71.3%
expm1-define76.6%
Simplified76.6%
if -5.2000000000000002e-8 < i < -2.34999999999999983e-177 or 2.90000000000000011e-227 < i < 1.8999999999999999e-53Initial program 15.3%
Taylor expanded in n around inf 16.2%
expm1-define64.6%
Simplified64.6%
clear-num64.5%
un-div-inv64.4%
associate-/l/88.9%
Applied egg-rr88.9%
/-rgt-identity88.9%
associate-/r/88.9%
associate-/r/88.8%
associate-/r/89.0%
/-rgt-identity89.0%
*-commutative89.0%
Simplified89.0%
Taylor expanded in i around inf 16.4%
associate-*r/16.4%
expm1-define88.9%
associate-*r*89.2%
Simplified89.2%
Taylor expanded in i around 0 89.2%
associate-*r*89.2%
distribute-rgt-out89.2%
*-commutative89.2%
Applied egg-rr89.2%
if -2.34999999999999983e-177 < i < 2.90000000000000011e-227Initial program 4.5%
associate-*r/4.5%
sub-neg4.5%
distribute-rgt-in4.5%
metadata-eval4.5%
metadata-eval4.5%
Simplified4.5%
metadata-eval4.5%
metadata-eval4.5%
distribute-rgt-in4.5%
sub-neg4.5%
associate-*r/4.5%
associate-/r/5.3%
associate-*r*5.3%
add-exp-log5.3%
expm1-define5.3%
log-pow11.6%
log1p-define59.9%
Applied egg-rr59.9%
Taylor expanded in i around 0 93.4%
sub-neg93.4%
associate-*r/93.4%
metadata-eval93.4%
sub-neg93.4%
Simplified93.4%
Taylor expanded in n around 0 93.4%
associate-*r/93.4%
*-commutative93.4%
associate-/l*93.4%
Simplified93.4%
Taylor expanded in i around 0 93.4%
if 1.22000000000000004e117 < i < 1.06e158 or 2.20000000000000011e204 < i < 1.5e273Initial program 38.5%
associate-*r/38.5%
sub-neg38.5%
distribute-rgt-in38.5%
metadata-eval38.5%
metadata-eval38.5%
Simplified38.5%
Taylor expanded in i around 0 67.5%
Taylor expanded in i around 0 67.5%
if 1.06e158 < i < 2.20000000000000011e204Initial program 62.9%
Taylor expanded in n around inf 75.3%
expm1-define75.3%
Simplified75.3%
clear-num75.3%
un-div-inv75.3%
associate-/l/75.3%
Applied egg-rr75.3%
/-rgt-identity75.3%
associate-/r/75.3%
associate-/r/75.3%
associate-/r/75.3%
/-rgt-identity75.3%
*-commutative75.3%
Simplified75.3%
Taylor expanded in i around inf 75.3%
associate-*r/75.3%
expm1-define75.3%
associate-*r*75.3%
Simplified75.3%
Taylor expanded in i around 0 75.4%
if 1.5e273 < i Initial program 75.0%
Taylor expanded in n around inf 50.2%
expm1-define50.2%
Simplified50.2%
clear-num50.2%
un-div-inv50.2%
associate-/l/50.6%
Applied egg-rr50.6%
/-rgt-identity50.6%
associate-/r/50.6%
associate-/r/50.6%
associate-/r/50.6%
/-rgt-identity50.6%
*-commutative50.6%
Simplified50.6%
Taylor expanded in i around 0 51.5%
*-commutative51.5%
Simplified51.5%
Final simplification82.0%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* 100.0 (* (expm1 i) (/ n i))))
(t_1 (/ (* i (+ (* n 100.0) (* 50.0 (* i n)))) i)))
(if (<= i -0.235)
t_0
(if (<= i -2.15e-178)
(*
(/ 100.0 i)
(*
i
(+
n
(*
i
(+
(* n 0.5)
(* (* i n) (+ 0.16666666666666666 (* i 0.041666666666666664))))))))
(if (<= i 8.5e-237)
(+ (* n 100.0) (* i -50.0))
(if (<= i 3e-36)
t_1
(if (<= i 1.22e+117)
t_0
(if (<= i 1.06e+158)
0.0
(if (<= i 2.25e+204)
t_1
(if (<= i 2.8e+273) 0.0 (* (* i n) (/ 100.0 i))))))))))))
double code(double i, double n) {
double t_0 = 100.0 * (expm1(i) * (n / i));
double t_1 = (i * ((n * 100.0) + (50.0 * (i * n)))) / i;
double tmp;
if (i <= -0.235) {
tmp = t_0;
} else if (i <= -2.15e-178) {
tmp = (100.0 / i) * (i * (n + (i * ((n * 0.5) + ((i * n) * (0.16666666666666666 + (i * 0.041666666666666664)))))));
} else if (i <= 8.5e-237) {
tmp = (n * 100.0) + (i * -50.0);
} else if (i <= 3e-36) {
tmp = t_1;
} else if (i <= 1.22e+117) {
tmp = t_0;
} else if (i <= 1.06e+158) {
tmp = 0.0;
} else if (i <= 2.25e+204) {
tmp = t_1;
} else if (i <= 2.8e+273) {
tmp = 0.0;
} else {
tmp = (i * n) * (100.0 / i);
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = 100.0 * (Math.expm1(i) * (n / i));
double t_1 = (i * ((n * 100.0) + (50.0 * (i * n)))) / i;
double tmp;
if (i <= -0.235) {
tmp = t_0;
} else if (i <= -2.15e-178) {
tmp = (100.0 / i) * (i * (n + (i * ((n * 0.5) + ((i * n) * (0.16666666666666666 + (i * 0.041666666666666664)))))));
} else if (i <= 8.5e-237) {
tmp = (n * 100.0) + (i * -50.0);
} else if (i <= 3e-36) {
tmp = t_1;
} else if (i <= 1.22e+117) {
tmp = t_0;
} else if (i <= 1.06e+158) {
tmp = 0.0;
} else if (i <= 2.25e+204) {
tmp = t_1;
} else if (i <= 2.8e+273) {
tmp = 0.0;
} else {
tmp = (i * n) * (100.0 / i);
}
return tmp;
}
def code(i, n): t_0 = 100.0 * (math.expm1(i) * (n / i)) t_1 = (i * ((n * 100.0) + (50.0 * (i * n)))) / i tmp = 0 if i <= -0.235: tmp = t_0 elif i <= -2.15e-178: tmp = (100.0 / i) * (i * (n + (i * ((n * 0.5) + ((i * n) * (0.16666666666666666 + (i * 0.041666666666666664))))))) elif i <= 8.5e-237: tmp = (n * 100.0) + (i * -50.0) elif i <= 3e-36: tmp = t_1 elif i <= 1.22e+117: tmp = t_0 elif i <= 1.06e+158: tmp = 0.0 elif i <= 2.25e+204: tmp = t_1 elif i <= 2.8e+273: tmp = 0.0 else: tmp = (i * n) * (100.0 / i) return tmp
function code(i, n) t_0 = Float64(100.0 * Float64(expm1(i) * Float64(n / i))) t_1 = Float64(Float64(i * Float64(Float64(n * 100.0) + Float64(50.0 * Float64(i * n)))) / i) tmp = 0.0 if (i <= -0.235) tmp = t_0; elseif (i <= -2.15e-178) tmp = Float64(Float64(100.0 / i) * Float64(i * Float64(n + Float64(i * Float64(Float64(n * 0.5) + Float64(Float64(i * n) * Float64(0.16666666666666666 + Float64(i * 0.041666666666666664)))))))); elseif (i <= 8.5e-237) tmp = Float64(Float64(n * 100.0) + Float64(i * -50.0)); elseif (i <= 3e-36) tmp = t_1; elseif (i <= 1.22e+117) tmp = t_0; elseif (i <= 1.06e+158) tmp = 0.0; elseif (i <= 2.25e+204) tmp = t_1; elseif (i <= 2.8e+273) tmp = 0.0; else tmp = Float64(Float64(i * n) * Float64(100.0 / i)); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(100.0 * N[(N[(Exp[i] - 1), $MachinePrecision] * N[(n / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(i * N[(N[(n * 100.0), $MachinePrecision] + N[(50.0 * N[(i * n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]}, If[LessEqual[i, -0.235], t$95$0, If[LessEqual[i, -2.15e-178], N[(N[(100.0 / i), $MachinePrecision] * N[(i * N[(n + N[(i * N[(N[(n * 0.5), $MachinePrecision] + N[(N[(i * n), $MachinePrecision] * N[(0.16666666666666666 + N[(i * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[i, 8.5e-237], N[(N[(n * 100.0), $MachinePrecision] + N[(i * -50.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[i, 3e-36], t$95$1, If[LessEqual[i, 1.22e+117], t$95$0, If[LessEqual[i, 1.06e+158], 0.0, If[LessEqual[i, 2.25e+204], t$95$1, If[LessEqual[i, 2.8e+273], 0.0, N[(N[(i * n), $MachinePrecision] * N[(100.0 / i), $MachinePrecision]), $MachinePrecision]]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 100 \cdot \left(\mathsf{expm1}\left(i\right) \cdot \frac{n}{i}\right)\\
t_1 := \frac{i \cdot \left(n \cdot 100 + 50 \cdot \left(i \cdot n\right)\right)}{i}\\
\mathbf{if}\;i \leq -0.235:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;i \leq -2.15 \cdot 10^{-178}:\\
\;\;\;\;\frac{100}{i} \cdot \left(i \cdot \left(n + i \cdot \left(n \cdot 0.5 + \left(i \cdot n\right) \cdot \left(0.16666666666666666 + i \cdot 0.041666666666666664\right)\right)\right)\right)\\
\mathbf{elif}\;i \leq 8.5 \cdot 10^{-237}:\\
\;\;\;\;n \cdot 100 + i \cdot -50\\
\mathbf{elif}\;i \leq 3 \cdot 10^{-36}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;i \leq 1.22 \cdot 10^{+117}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;i \leq 1.06 \cdot 10^{+158}:\\
\;\;\;\;0\\
\mathbf{elif}\;i \leq 2.25 \cdot 10^{+204}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;i \leq 2.8 \cdot 10^{+273}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;\left(i \cdot n\right) \cdot \frac{100}{i}\\
\end{array}
\end{array}
if i < -0.23499999999999999 or 3.0000000000000002e-36 < i < 1.22000000000000004e117Initial program 45.3%
associate-/r/44.1%
add-exp-log44.1%
expm1-define44.1%
log-pow46.8%
log1p-define94.3%
Applied egg-rr94.3%
Taylor expanded in n around inf 74.5%
expm1-define77.5%
associate-*l/76.5%
Simplified76.5%
if -0.23499999999999999 < i < -2.15e-178Initial program 15.5%
Taylor expanded in n around inf 19.5%
expm1-define63.7%
Simplified63.7%
clear-num63.5%
un-div-inv63.5%
associate-/l/87.5%
Applied egg-rr87.5%
/-rgt-identity87.5%
associate-/r/87.5%
associate-/r/87.4%
associate-/r/87.7%
/-rgt-identity87.7%
*-commutative87.7%
Simplified87.7%
Taylor expanded in i around 0 87.0%
Taylor expanded in n around 0 87.0%
associate-*r*87.0%
*-commutative87.0%
*-commutative87.0%
Simplified87.0%
if -2.15e-178 < i < 8.49999999999999951e-237Initial program 4.5%
associate-*r/4.5%
sub-neg4.5%
distribute-rgt-in4.5%
metadata-eval4.5%
metadata-eval4.5%
Simplified4.5%
metadata-eval4.5%
metadata-eval4.5%
distribute-rgt-in4.5%
sub-neg4.5%
associate-*r/4.5%
associate-/r/5.3%
associate-*r*5.3%
add-exp-log5.3%
expm1-define5.3%
log-pow10.2%
log1p-define58.3%
Applied egg-rr58.3%
Taylor expanded in i around 0 94.5%
sub-neg94.5%
associate-*r/94.5%
metadata-eval94.5%
sub-neg94.5%
Simplified94.5%
Taylor expanded in n around 0 94.5%
associate-*r/94.5%
*-commutative94.5%
associate-/l*94.5%
Simplified94.5%
Taylor expanded in i around 0 94.5%
if 8.49999999999999951e-237 < i < 3.0000000000000002e-36 or 1.06e158 < i < 2.25000000000000001e204Initial program 24.8%
Taylor expanded in n around inf 26.7%
expm1-define66.9%
Simplified66.9%
clear-num66.8%
un-div-inv66.7%
associate-/l/83.2%
Applied egg-rr83.2%
/-rgt-identity83.2%
associate-/r/83.2%
associate-/r/83.2%
associate-/r/83.2%
/-rgt-identity83.2%
*-commutative83.2%
Simplified83.2%
Taylor expanded in i around inf 26.9%
associate-*r/26.9%
expm1-define83.4%
associate-*r*83.5%
Simplified83.5%
Taylor expanded in i around 0 83.5%
if 1.22000000000000004e117 < i < 1.06e158 or 2.25000000000000001e204 < i < 2.80000000000000018e273Initial program 38.5%
associate-*r/38.5%
sub-neg38.5%
distribute-rgt-in38.5%
metadata-eval38.5%
metadata-eval38.5%
Simplified38.5%
Taylor expanded in i around 0 67.5%
Taylor expanded in i around 0 67.5%
if 2.80000000000000018e273 < i Initial program 75.0%
Taylor expanded in n around inf 50.2%
expm1-define50.2%
Simplified50.2%
clear-num50.2%
un-div-inv50.2%
associate-/l/50.6%
Applied egg-rr50.6%
/-rgt-identity50.6%
associate-/r/50.6%
associate-/r/50.6%
associate-/r/50.6%
/-rgt-identity50.6%
*-commutative50.6%
Simplified50.6%
Taylor expanded in i around 0 51.5%
*-commutative51.5%
Simplified51.5%
Final simplification81.6%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* n (/ (* 100.0 (expm1 i)) i))))
(if (<= n -1.9e-238)
t_0
(if (<= n 2.7e-216)
(/ 0.0 (/ i n))
(if (<= n 7.5e-41) (* 100.0 (/ i (/ i n))) t_0)))))
double code(double i, double n) {
double t_0 = n * ((100.0 * expm1(i)) / i);
double tmp;
if (n <= -1.9e-238) {
tmp = t_0;
} else if (n <= 2.7e-216) {
tmp = 0.0 / (i / n);
} else if (n <= 7.5e-41) {
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 * Math.expm1(i)) / i);
double tmp;
if (n <= -1.9e-238) {
tmp = t_0;
} else if (n <= 2.7e-216) {
tmp = 0.0 / (i / n);
} else if (n <= 7.5e-41) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = n * ((100.0 * math.expm1(i)) / i) tmp = 0 if n <= -1.9e-238: tmp = t_0 elif n <= 2.7e-216: tmp = 0.0 / (i / n) elif n <= 7.5e-41: tmp = 100.0 * (i / (i / n)) else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(n * Float64(Float64(100.0 * expm1(i)) / i)) tmp = 0.0 if (n <= -1.9e-238) tmp = t_0; elseif (n <= 2.7e-216) tmp = Float64(0.0 / Float64(i / n)); elseif (n <= 7.5e-41) 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 * N[(Exp[i] - 1), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -1.9e-238], t$95$0, If[LessEqual[n, 2.7e-216], N[(0.0 / N[(i / n), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 7.5e-41], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := n \cdot \frac{100 \cdot \mathsf{expm1}\left(i\right)}{i}\\
\mathbf{if}\;n \leq -1.9 \cdot 10^{-238}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 2.7 \cdot 10^{-216}:\\
\;\;\;\;\frac{0}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 7.5 \cdot 10^{-41}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -1.8999999999999998e-238 or 7.50000000000000049e-41 < n Initial program 22.1%
associate-/r/22.3%
associate-*r*22.3%
*-commutative22.3%
associate-*r/22.3%
sub-neg22.3%
distribute-lft-in22.3%
metadata-eval22.3%
metadata-eval22.3%
metadata-eval22.3%
fma-define22.3%
metadata-eval22.3%
Simplified22.3%
Taylor expanded in n around inf 39.4%
sub-neg39.4%
metadata-eval39.4%
metadata-eval39.4%
distribute-lft-in39.5%
metadata-eval39.5%
sub-neg39.5%
expm1-define89.5%
Simplified89.5%
if -1.8999999999999998e-238 < n < 2.6999999999999999e-216Initial program 70.8%
associate-*r/70.8%
sub-neg70.8%
distribute-rgt-in70.8%
metadata-eval70.8%
metadata-eval70.8%
Simplified70.8%
Taylor expanded in i around 0 86.0%
if 2.6999999999999999e-216 < n < 7.50000000000000049e-41Initial program 26.6%
Taylor expanded in i around 0 63.9%
Final simplification86.1%
(FPCore (i n)
:precision binary64
(let* ((t_0
(*
(/ 100.0 i)
(*
i
(+
n
(*
i
(+
(* n 0.5)
(*
(* i n)
(+ 0.16666666666666666 (* i 0.041666666666666664))))))))))
(if (<= n -6.2e+77)
t_0
(if (<= n -4.8e-218)
(* 100.0 (/ i (/ i n)))
(if (<= n 6e-118) 0.0 t_0)))))
double code(double i, double n) {
double t_0 = (100.0 / i) * (i * (n + (i * ((n * 0.5) + ((i * n) * (0.16666666666666666 + (i * 0.041666666666666664)))))));
double tmp;
if (n <= -6.2e+77) {
tmp = t_0;
} else if (n <= -4.8e-218) {
tmp = 100.0 * (i / (i / n));
} else if (n <= 6e-118) {
tmp = 0.0;
} else {
tmp = t_0;
}
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 / i) * (i * (n + (i * ((n * 0.5d0) + ((i * n) * (0.16666666666666666d0 + (i * 0.041666666666666664d0)))))))
if (n <= (-6.2d+77)) then
tmp = t_0
else if (n <= (-4.8d-218)) then
tmp = 100.0d0 * (i / (i / n))
else if (n <= 6d-118) then
tmp = 0.0d0
else
tmp = t_0
end if
code = tmp
end function
public static double code(double i, double n) {
double t_0 = (100.0 / i) * (i * (n + (i * ((n * 0.5) + ((i * n) * (0.16666666666666666 + (i * 0.041666666666666664)))))));
double tmp;
if (n <= -6.2e+77) {
tmp = t_0;
} else if (n <= -4.8e-218) {
tmp = 100.0 * (i / (i / n));
} else if (n <= 6e-118) {
tmp = 0.0;
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = (100.0 / i) * (i * (n + (i * ((n * 0.5) + ((i * n) * (0.16666666666666666 + (i * 0.041666666666666664))))))) tmp = 0 if n <= -6.2e+77: tmp = t_0 elif n <= -4.8e-218: tmp = 100.0 * (i / (i / n)) elif n <= 6e-118: tmp = 0.0 else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(Float64(100.0 / i) * Float64(i * Float64(n + Float64(i * Float64(Float64(n * 0.5) + Float64(Float64(i * n) * Float64(0.16666666666666666 + Float64(i * 0.041666666666666664)))))))) tmp = 0.0 if (n <= -6.2e+77) tmp = t_0; elseif (n <= -4.8e-218) tmp = Float64(100.0 * Float64(i / Float64(i / n))); elseif (n <= 6e-118) tmp = 0.0; else tmp = t_0; end return tmp end
function tmp_2 = code(i, n) t_0 = (100.0 / i) * (i * (n + (i * ((n * 0.5) + ((i * n) * (0.16666666666666666 + (i * 0.041666666666666664))))))); tmp = 0.0; if (n <= -6.2e+77) tmp = t_0; elseif (n <= -4.8e-218) tmp = 100.0 * (i / (i / n)); elseif (n <= 6e-118) tmp = 0.0; else tmp = t_0; end tmp_2 = tmp; end
code[i_, n_] := Block[{t$95$0 = N[(N[(100.0 / i), $MachinePrecision] * N[(i * N[(n + N[(i * N[(N[(n * 0.5), $MachinePrecision] + N[(N[(i * n), $MachinePrecision] * N[(0.16666666666666666 + N[(i * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -6.2e+77], t$95$0, If[LessEqual[n, -4.8e-218], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 6e-118], 0.0, t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{100}{i} \cdot \left(i \cdot \left(n + i \cdot \left(n \cdot 0.5 + \left(i \cdot n\right) \cdot \left(0.16666666666666666 + i \cdot 0.041666666666666664\right)\right)\right)\right)\\
\mathbf{if}\;n \leq -6.2 \cdot 10^{+77}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq -4.8 \cdot 10^{-218}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 6 \cdot 10^{-118}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -6.19999999999999997e77 or 6.00000000000000035e-118 < n Initial program 19.0%
Taylor expanded in n around inf 39.8%
expm1-define71.1%
Simplified71.1%
clear-num71.0%
un-div-inv71.0%
associate-/l/90.2%
Applied egg-rr90.2%
/-rgt-identity90.2%
associate-/r/90.2%
associate-/r/90.1%
associate-/r/90.3%
/-rgt-identity90.3%
*-commutative90.3%
Simplified90.3%
Taylor expanded in i around 0 73.1%
Taylor expanded in n around 0 73.1%
associate-*r*73.1%
*-commutative73.1%
*-commutative73.1%
Simplified73.1%
if -6.19999999999999997e77 < n < -4.8000000000000002e-218Initial program 32.2%
Taylor expanded in i around 0 55.8%
if -4.8000000000000002e-218 < n < 6.00000000000000035e-118Initial program 55.6%
associate-*r/55.6%
sub-neg55.6%
distribute-rgt-in55.6%
metadata-eval55.6%
metadata-eval55.6%
Simplified55.6%
Taylor expanded in i around 0 77.6%
Taylor expanded in i around 0 77.6%
Final simplification70.5%
(FPCore (i n)
:precision binary64
(if (<= n -1.9e-148)
(*
100.0
(+
n
(*
i
(+
(* n 0.5)
(* i (+ (* (* i n) 0.041666666666666664) (* n 0.16666666666666666)))))))
(if (<= n -8.2e-209)
0.0
(if (<= n -4.7e-225)
(* 100.0 (/ i (/ i n)))
(if (<= n 2.4e-83)
(/ 0.0 (/ i n))
(*
100.0
(*
n
(+
1.0
(*
i
(+
0.5
(*
i
(+ 0.16666666666666666 (* i 0.041666666666666664)))))))))))))
double code(double i, double n) {
double tmp;
if (n <= -1.9e-148) {
tmp = 100.0 * (n + (i * ((n * 0.5) + (i * (((i * n) * 0.041666666666666664) + (n * 0.16666666666666666))))));
} else if (n <= -8.2e-209) {
tmp = 0.0;
} else if (n <= -4.7e-225) {
tmp = 100.0 * (i / (i / n));
} else if (n <= 2.4e-83) {
tmp = 0.0 / (i / n);
} else {
tmp = 100.0 * (n * (1.0 + (i * (0.5 + (i * (0.16666666666666666 + (i * 0.041666666666666664)))))));
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if (n <= (-1.9d-148)) then
tmp = 100.0d0 * (n + (i * ((n * 0.5d0) + (i * (((i * n) * 0.041666666666666664d0) + (n * 0.16666666666666666d0))))))
else if (n <= (-8.2d-209)) then
tmp = 0.0d0
else if (n <= (-4.7d-225)) then
tmp = 100.0d0 * (i / (i / n))
else if (n <= 2.4d-83) then
tmp = 0.0d0 / (i / n)
else
tmp = 100.0d0 * (n * (1.0d0 + (i * (0.5d0 + (i * (0.16666666666666666d0 + (i * 0.041666666666666664d0)))))))
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (n <= -1.9e-148) {
tmp = 100.0 * (n + (i * ((n * 0.5) + (i * (((i * n) * 0.041666666666666664) + (n * 0.16666666666666666))))));
} else if (n <= -8.2e-209) {
tmp = 0.0;
} else if (n <= -4.7e-225) {
tmp = 100.0 * (i / (i / n));
} else if (n <= 2.4e-83) {
tmp = 0.0 / (i / n);
} else {
tmp = 100.0 * (n * (1.0 + (i * (0.5 + (i * (0.16666666666666666 + (i * 0.041666666666666664)))))));
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -1.9e-148: tmp = 100.0 * (n + (i * ((n * 0.5) + (i * (((i * n) * 0.041666666666666664) + (n * 0.16666666666666666)))))) elif n <= -8.2e-209: tmp = 0.0 elif n <= -4.7e-225: tmp = 100.0 * (i / (i / n)) elif n <= 2.4e-83: tmp = 0.0 / (i / n) else: tmp = 100.0 * (n * (1.0 + (i * (0.5 + (i * (0.16666666666666666 + (i * 0.041666666666666664))))))) return tmp
function code(i, n) tmp = 0.0 if (n <= -1.9e-148) tmp = Float64(100.0 * Float64(n + Float64(i * Float64(Float64(n * 0.5) + Float64(i * Float64(Float64(Float64(i * n) * 0.041666666666666664) + Float64(n * 0.16666666666666666))))))); elseif (n <= -8.2e-209) tmp = 0.0; elseif (n <= -4.7e-225) tmp = Float64(100.0 * Float64(i / Float64(i / n))); elseif (n <= 2.4e-83) tmp = Float64(0.0 / Float64(i / n)); else tmp = Float64(100.0 * Float64(n * Float64(1.0 + Float64(i * Float64(0.5 + Float64(i * Float64(0.16666666666666666 + Float64(i * 0.041666666666666664)))))))); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (n <= -1.9e-148) tmp = 100.0 * (n + (i * ((n * 0.5) + (i * (((i * n) * 0.041666666666666664) + (n * 0.16666666666666666)))))); elseif (n <= -8.2e-209) tmp = 0.0; elseif (n <= -4.7e-225) tmp = 100.0 * (i / (i / n)); elseif (n <= 2.4e-83) tmp = 0.0 / (i / n); else tmp = 100.0 * (n * (1.0 + (i * (0.5 + (i * (0.16666666666666666 + (i * 0.041666666666666664))))))); end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[n, -1.9e-148], N[(100.0 * N[(n + N[(i * N[(N[(n * 0.5), $MachinePrecision] + N[(i * N[(N[(N[(i * n), $MachinePrecision] * 0.041666666666666664), $MachinePrecision] + N[(n * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, -8.2e-209], 0.0, If[LessEqual[n, -4.7e-225], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 2.4e-83], N[(0.0 / N[(i / n), $MachinePrecision]), $MachinePrecision], N[(100.0 * N[(n * N[(1.0 + N[(i * N[(0.5 + N[(i * N[(0.16666666666666666 + N[(i * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.9 \cdot 10^{-148}:\\
\;\;\;\;100 \cdot \left(n + i \cdot \left(n \cdot 0.5 + i \cdot \left(\left(i \cdot n\right) \cdot 0.041666666666666664 + n \cdot 0.16666666666666666\right)\right)\right)\\
\mathbf{elif}\;n \leq -8.2 \cdot 10^{-209}:\\
\;\;\;\;0\\
\mathbf{elif}\;n \leq -4.7 \cdot 10^{-225}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 2.4 \cdot 10^{-83}:\\
\;\;\;\;\frac{0}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \left(n \cdot \left(1 + i \cdot \left(0.5 + i \cdot \left(0.16666666666666666 + i \cdot 0.041666666666666664\right)\right)\right)\right)\\
\end{array}
\end{array}
if n < -1.90000000000000007e-148Initial program 20.6%
Taylor expanded in n around inf 36.6%
expm1-define71.8%
Simplified71.8%
Taylor expanded in i around 0 59.1%
if -1.90000000000000007e-148 < n < -8.19999999999999955e-209Initial program 73.1%
associate-*r/73.1%
sub-neg73.1%
distribute-rgt-in73.1%
metadata-eval73.1%
metadata-eval73.1%
Simplified73.1%
Taylor expanded in i around 0 73.1%
Taylor expanded in i around 0 73.1%
if -8.19999999999999955e-209 < n < -4.70000000000000014e-225Initial program 6.3%
Taylor expanded in i around 0 69.0%
if -4.70000000000000014e-225 < n < 2.4000000000000001e-83Initial program 53.3%
associate-*r/53.3%
sub-neg53.3%
distribute-rgt-in53.3%
metadata-eval53.3%
metadata-eval53.3%
Simplified53.3%
Taylor expanded in i around 0 76.6%
if 2.4000000000000001e-83 < n Initial program 20.1%
Taylor expanded in n around inf 40.3%
expm1-define75.9%
Simplified75.9%
clear-num75.8%
un-div-inv75.8%
associate-/l/92.9%
Applied egg-rr92.9%
/-rgt-identity92.9%
associate-/r/92.9%
associate-/r/92.8%
associate-/r/92.9%
/-rgt-identity92.9%
*-commutative92.9%
Simplified92.9%
Taylor expanded in i around 0 82.9%
Taylor expanded in n around 0 81.0%
Final simplification70.9%
(FPCore (i n)
:precision binary64
(let* ((t_0
(*
100.0
(*
n
(+
1.0
(*
i
(+
0.5
(* i (+ 0.16666666666666666 (* i 0.041666666666666664))))))))))
(if (<= n -1.22e-148)
t_0
(if (<= n -4.8e-209)
0.0
(if (<= n -7e-224)
(* 100.0 (/ i (/ i n)))
(if (<= n 2.4e-83) (/ 0.0 (/ i n)) t_0))))))
double code(double i, double n) {
double t_0 = 100.0 * (n * (1.0 + (i * (0.5 + (i * (0.16666666666666666 + (i * 0.041666666666666664)))))));
double tmp;
if (n <= -1.22e-148) {
tmp = t_0;
} else if (n <= -4.8e-209) {
tmp = 0.0;
} else if (n <= -7e-224) {
tmp = 100.0 * (i / (i / n));
} else if (n <= 2.4e-83) {
tmp = 0.0 / (i / n);
} else {
tmp = t_0;
}
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 * (n * (1.0d0 + (i * (0.5d0 + (i * (0.16666666666666666d0 + (i * 0.041666666666666664d0)))))))
if (n <= (-1.22d-148)) then
tmp = t_0
else if (n <= (-4.8d-209)) then
tmp = 0.0d0
else if (n <= (-7d-224)) then
tmp = 100.0d0 * (i / (i / n))
else if (n <= 2.4d-83) then
tmp = 0.0d0 / (i / n)
else
tmp = t_0
end if
code = tmp
end function
public static double code(double i, double n) {
double t_0 = 100.0 * (n * (1.0 + (i * (0.5 + (i * (0.16666666666666666 + (i * 0.041666666666666664)))))));
double tmp;
if (n <= -1.22e-148) {
tmp = t_0;
} else if (n <= -4.8e-209) {
tmp = 0.0;
} else if (n <= -7e-224) {
tmp = 100.0 * (i / (i / n));
} else if (n <= 2.4e-83) {
tmp = 0.0 / (i / n);
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = 100.0 * (n * (1.0 + (i * (0.5 + (i * (0.16666666666666666 + (i * 0.041666666666666664))))))) tmp = 0 if n <= -1.22e-148: tmp = t_0 elif n <= -4.8e-209: tmp = 0.0 elif n <= -7e-224: tmp = 100.0 * (i / (i / n)) elif n <= 2.4e-83: tmp = 0.0 / (i / n) else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(100.0 * Float64(n * Float64(1.0 + Float64(i * Float64(0.5 + Float64(i * Float64(0.16666666666666666 + Float64(i * 0.041666666666666664)))))))) tmp = 0.0 if (n <= -1.22e-148) tmp = t_0; elseif (n <= -4.8e-209) tmp = 0.0; elseif (n <= -7e-224) tmp = Float64(100.0 * Float64(i / Float64(i / n))); elseif (n <= 2.4e-83) tmp = Float64(0.0 / Float64(i / n)); else tmp = t_0; end return tmp end
function tmp_2 = code(i, n) t_0 = 100.0 * (n * (1.0 + (i * (0.5 + (i * (0.16666666666666666 + (i * 0.041666666666666664))))))); tmp = 0.0; if (n <= -1.22e-148) tmp = t_0; elseif (n <= -4.8e-209) tmp = 0.0; elseif (n <= -7e-224) tmp = 100.0 * (i / (i / n)); elseif (n <= 2.4e-83) tmp = 0.0 / (i / n); else tmp = t_0; end tmp_2 = tmp; end
code[i_, n_] := Block[{t$95$0 = N[(100.0 * N[(n * N[(1.0 + N[(i * N[(0.5 + N[(i * N[(0.16666666666666666 + N[(i * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -1.22e-148], t$95$0, If[LessEqual[n, -4.8e-209], 0.0, If[LessEqual[n, -7e-224], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 2.4e-83], N[(0.0 / N[(i / n), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 100 \cdot \left(n \cdot \left(1 + i \cdot \left(0.5 + i \cdot \left(0.16666666666666666 + i \cdot 0.041666666666666664\right)\right)\right)\right)\\
\mathbf{if}\;n \leq -1.22 \cdot 10^{-148}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq -4.8 \cdot 10^{-209}:\\
\;\;\;\;0\\
\mathbf{elif}\;n \leq -7 \cdot 10^{-224}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 2.4 \cdot 10^{-83}:\\
\;\;\;\;\frac{0}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -1.21999999999999992e-148 or 2.4000000000000001e-83 < n Initial program 20.4%
Taylor expanded in n around inf 38.3%
expm1-define73.7%
Simplified73.7%
clear-num73.6%
un-div-inv73.6%
associate-/l/87.5%
Applied egg-rr87.5%
/-rgt-identity87.5%
associate-/r/87.5%
associate-/r/87.0%
associate-/r/87.2%
/-rgt-identity87.2%
*-commutative87.2%
Simplified87.2%
Taylor expanded in i around 0 67.8%
Taylor expanded in n around 0 69.1%
if -1.21999999999999992e-148 < n < -4.8000000000000002e-209Initial program 73.1%
associate-*r/73.1%
sub-neg73.1%
distribute-rgt-in73.1%
metadata-eval73.1%
metadata-eval73.1%
Simplified73.1%
Taylor expanded in i around 0 73.1%
Taylor expanded in i around 0 73.1%
if -4.8000000000000002e-209 < n < -7.00000000000000037e-224Initial program 6.3%
Taylor expanded in i around 0 69.0%
if -7.00000000000000037e-224 < n < 2.4000000000000001e-83Initial program 53.3%
associate-*r/53.3%
sub-neg53.3%
distribute-rgt-in53.3%
metadata-eval53.3%
metadata-eval53.3%
Simplified53.3%
Taylor expanded in i around 0 76.6%
Final simplification70.9%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* 100.0 (/ i (/ i n)))))
(if (<= n -2.4e-148)
(* n (+ 100.0 (* i (+ 50.0 (* i 16.666666666666668)))))
(if (<= n -9.5e-209)
0.0
(if (<= n -3.6e-221)
t_0
(if (<= n 4e-252)
(/ 0.0 (/ i n))
(if (<= n 1.45e-92)
t_0
(/ (* i (* n (+ 100.0 (* i 50.0)))) i))))))))
double code(double i, double n) {
double t_0 = 100.0 * (i / (i / n));
double tmp;
if (n <= -2.4e-148) {
tmp = n * (100.0 + (i * (50.0 + (i * 16.666666666666668))));
} else if (n <= -9.5e-209) {
tmp = 0.0;
} else if (n <= -3.6e-221) {
tmp = t_0;
} else if (n <= 4e-252) {
tmp = 0.0 / (i / n);
} else if (n <= 1.45e-92) {
tmp = t_0;
} else {
tmp = (i * (n * (100.0 + (i * 50.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 * (i / (i / n))
if (n <= (-2.4d-148)) then
tmp = n * (100.0d0 + (i * (50.0d0 + (i * 16.666666666666668d0))))
else if (n <= (-9.5d-209)) then
tmp = 0.0d0
else if (n <= (-3.6d-221)) then
tmp = t_0
else if (n <= 4d-252) then
tmp = 0.0d0 / (i / n)
else if (n <= 1.45d-92) then
tmp = t_0
else
tmp = (i * (n * (100.0d0 + (i * 50.0d0)))) / i
end if
code = tmp
end function
public static double code(double i, double n) {
double t_0 = 100.0 * (i / (i / n));
double tmp;
if (n <= -2.4e-148) {
tmp = n * (100.0 + (i * (50.0 + (i * 16.666666666666668))));
} else if (n <= -9.5e-209) {
tmp = 0.0;
} else if (n <= -3.6e-221) {
tmp = t_0;
} else if (n <= 4e-252) {
tmp = 0.0 / (i / n);
} else if (n <= 1.45e-92) {
tmp = t_0;
} else {
tmp = (i * (n * (100.0 + (i * 50.0)))) / i;
}
return tmp;
}
def code(i, n): t_0 = 100.0 * (i / (i / n)) tmp = 0 if n <= -2.4e-148: tmp = n * (100.0 + (i * (50.0 + (i * 16.666666666666668)))) elif n <= -9.5e-209: tmp = 0.0 elif n <= -3.6e-221: tmp = t_0 elif n <= 4e-252: tmp = 0.0 / (i / n) elif n <= 1.45e-92: tmp = t_0 else: tmp = (i * (n * (100.0 + (i * 50.0)))) / i return tmp
function code(i, n) t_0 = Float64(100.0 * Float64(i / Float64(i / n))) tmp = 0.0 if (n <= -2.4e-148) tmp = Float64(n * Float64(100.0 + Float64(i * Float64(50.0 + Float64(i * 16.666666666666668))))); elseif (n <= -9.5e-209) tmp = 0.0; elseif (n <= -3.6e-221) tmp = t_0; elseif (n <= 4e-252) tmp = Float64(0.0 / Float64(i / n)); elseif (n <= 1.45e-92) tmp = t_0; else tmp = Float64(Float64(i * Float64(n * Float64(100.0 + Float64(i * 50.0)))) / i); end return tmp end
function tmp_2 = code(i, n) t_0 = 100.0 * (i / (i / n)); tmp = 0.0; if (n <= -2.4e-148) tmp = n * (100.0 + (i * (50.0 + (i * 16.666666666666668)))); elseif (n <= -9.5e-209) tmp = 0.0; elseif (n <= -3.6e-221) tmp = t_0; elseif (n <= 4e-252) tmp = 0.0 / (i / n); elseif (n <= 1.45e-92) tmp = t_0; else tmp = (i * (n * (100.0 + (i * 50.0)))) / i; end tmp_2 = tmp; end
code[i_, n_] := Block[{t$95$0 = N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -2.4e-148], N[(n * N[(100.0 + N[(i * N[(50.0 + N[(i * 16.666666666666668), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, -9.5e-209], 0.0, If[LessEqual[n, -3.6e-221], t$95$0, If[LessEqual[n, 4e-252], N[(0.0 / N[(i / n), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.45e-92], t$95$0, N[(N[(i * N[(n * N[(100.0 + N[(i * 50.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{if}\;n \leq -2.4 \cdot 10^{-148}:\\
\;\;\;\;n \cdot \left(100 + i \cdot \left(50 + i \cdot 16.666666666666668\right)\right)\\
\mathbf{elif}\;n \leq -9.5 \cdot 10^{-209}:\\
\;\;\;\;0\\
\mathbf{elif}\;n \leq -3.6 \cdot 10^{-221}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 4 \cdot 10^{-252}:\\
\;\;\;\;\frac{0}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 1.45 \cdot 10^{-92}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\frac{i \cdot \left(n \cdot \left(100 + i \cdot 50\right)\right)}{i}\\
\end{array}
\end{array}
if n < -2.4000000000000001e-148Initial program 20.6%
associate-*r/20.6%
sub-neg20.6%
distribute-rgt-in20.7%
metadata-eval20.7%
metadata-eval20.7%
Simplified20.7%
Taylor expanded in n around inf 36.6%
Taylor expanded in i around 0 58.0%
Taylor expanded in n around 0 59.1%
if -2.4000000000000001e-148 < n < -9.50000000000000028e-209Initial program 73.1%
associate-*r/73.1%
sub-neg73.1%
distribute-rgt-in73.1%
metadata-eval73.1%
metadata-eval73.1%
Simplified73.1%
Taylor expanded in i around 0 73.1%
Taylor expanded in i around 0 73.1%
if -9.50000000000000028e-209 < n < -3.60000000000000011e-221 or 3.99999999999999977e-252 < n < 1.44999999999999992e-92Initial program 28.6%
Taylor expanded in i around 0 63.3%
if -3.60000000000000011e-221 < n < 3.99999999999999977e-252Initial program 80.8%
associate-*r/80.8%
sub-neg80.8%
distribute-rgt-in80.8%
metadata-eval80.8%
metadata-eval80.8%
Simplified80.8%
Taylor expanded in i around 0 92.8%
if 1.44999999999999992e-92 < n Initial program 20.5%
Taylor expanded in n around inf 39.0%
expm1-define73.4%
Simplified73.4%
clear-num73.3%
un-div-inv73.3%
associate-/l/89.8%
Applied egg-rr89.8%
/-rgt-identity89.8%
associate-/r/89.9%
associate-/r/89.7%
associate-/r/89.9%
/-rgt-identity89.9%
*-commutative89.9%
Simplified89.9%
Taylor expanded in i around inf 39.4%
associate-*r/39.3%
expm1-define89.8%
associate-*r*90.0%
Simplified90.0%
Taylor expanded in i around 0 77.1%
associate-*r*77.1%
distribute-rgt-out77.1%
*-commutative77.1%
Applied egg-rr77.1%
Final simplification69.6%
(FPCore (i n)
:precision binary64
(let* ((t_0
(*
n
(+
100.0
(*
i
(+ 50.0 (* i (+ 16.666666666666668 (* i 4.166666666666667)))))))))
(if (<= n -9.2e-149)
t_0
(if (<= n -4.9e-209)
0.0
(if (<= n -2.5e-233)
(* 100.0 (/ i (/ i n)))
(if (<= n 2.4e-83) (/ 0.0 (/ i n)) t_0))))))
double code(double i, double n) {
double t_0 = n * (100.0 + (i * (50.0 + (i * (16.666666666666668 + (i * 4.166666666666667))))));
double tmp;
if (n <= -9.2e-149) {
tmp = t_0;
} else if (n <= -4.9e-209) {
tmp = 0.0;
} else if (n <= -2.5e-233) {
tmp = 100.0 * (i / (i / n));
} else if (n <= 2.4e-83) {
tmp = 0.0 / (i / n);
} else {
tmp = t_0;
}
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 = n * (100.0d0 + (i * (50.0d0 + (i * (16.666666666666668d0 + (i * 4.166666666666667d0))))))
if (n <= (-9.2d-149)) then
tmp = t_0
else if (n <= (-4.9d-209)) then
tmp = 0.0d0
else if (n <= (-2.5d-233)) then
tmp = 100.0d0 * (i / (i / n))
else if (n <= 2.4d-83) then
tmp = 0.0d0 / (i / n)
else
tmp = t_0
end if
code = tmp
end function
public static double code(double i, double n) {
double t_0 = n * (100.0 + (i * (50.0 + (i * (16.666666666666668 + (i * 4.166666666666667))))));
double tmp;
if (n <= -9.2e-149) {
tmp = t_0;
} else if (n <= -4.9e-209) {
tmp = 0.0;
} else if (n <= -2.5e-233) {
tmp = 100.0 * (i / (i / n));
} else if (n <= 2.4e-83) {
tmp = 0.0 / (i / n);
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = n * (100.0 + (i * (50.0 + (i * (16.666666666666668 + (i * 4.166666666666667)))))) tmp = 0 if n <= -9.2e-149: tmp = t_0 elif n <= -4.9e-209: tmp = 0.0 elif n <= -2.5e-233: tmp = 100.0 * (i / (i / n)) elif n <= 2.4e-83: tmp = 0.0 / (i / n) else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(n * Float64(100.0 + Float64(i * Float64(50.0 + Float64(i * Float64(16.666666666666668 + Float64(i * 4.166666666666667))))))) tmp = 0.0 if (n <= -9.2e-149) tmp = t_0; elseif (n <= -4.9e-209) tmp = 0.0; elseif (n <= -2.5e-233) tmp = Float64(100.0 * Float64(i / Float64(i / n))); elseif (n <= 2.4e-83) tmp = Float64(0.0 / Float64(i / n)); else tmp = t_0; end return tmp end
function tmp_2 = code(i, n) t_0 = n * (100.0 + (i * (50.0 + (i * (16.666666666666668 + (i * 4.166666666666667)))))); tmp = 0.0; if (n <= -9.2e-149) tmp = t_0; elseif (n <= -4.9e-209) tmp = 0.0; elseif (n <= -2.5e-233) tmp = 100.0 * (i / (i / n)); elseif (n <= 2.4e-83) tmp = 0.0 / (i / n); else tmp = t_0; end tmp_2 = tmp; end
code[i_, n_] := Block[{t$95$0 = N[(n * N[(100.0 + N[(i * N[(50.0 + N[(i * N[(16.666666666666668 + N[(i * 4.166666666666667), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -9.2e-149], t$95$0, If[LessEqual[n, -4.9e-209], 0.0, If[LessEqual[n, -2.5e-233], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 2.4e-83], N[(0.0 / N[(i / n), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := n \cdot \left(100 + i \cdot \left(50 + i \cdot \left(16.666666666666668 + i \cdot 4.166666666666667\right)\right)\right)\\
\mathbf{if}\;n \leq -9.2 \cdot 10^{-149}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq -4.9 \cdot 10^{-209}:\\
\;\;\;\;0\\
\mathbf{elif}\;n \leq -2.5 \cdot 10^{-233}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 2.4 \cdot 10^{-83}:\\
\;\;\;\;\frac{0}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -9.1999999999999999e-149 or 2.4000000000000001e-83 < n Initial program 20.4%
associate-*r/20.4%
sub-neg20.4%
distribute-rgt-in20.4%
metadata-eval20.4%
metadata-eval20.4%
Simplified20.4%
Taylor expanded in n around inf 38.3%
Taylor expanded in i around 0 53.3%
*-commutative53.3%
Simplified53.3%
Taylor expanded in n around 0 69.1%
if -9.1999999999999999e-149 < n < -4.90000000000000035e-209Initial program 73.1%
associate-*r/73.1%
sub-neg73.1%
distribute-rgt-in73.1%
metadata-eval73.1%
metadata-eval73.1%
Simplified73.1%
Taylor expanded in i around 0 73.1%
Taylor expanded in i around 0 73.1%
if -4.90000000000000035e-209 < n < -2.50000000000000006e-233Initial program 6.3%
Taylor expanded in i around 0 69.0%
if -2.50000000000000006e-233 < n < 2.4000000000000001e-83Initial program 53.3%
associate-*r/53.3%
sub-neg53.3%
distribute-rgt-in53.3%
metadata-eval53.3%
metadata-eval53.3%
Simplified53.3%
Taylor expanded in i around 0 76.6%
Final simplification70.9%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* n (+ 100.0 (* i (+ 50.0 (* i 16.666666666666668)))))))
(if (<= n -8.5e-149)
t_0
(if (<= n -1.12e-208)
0.0
(if (<= n -2.2e-232)
(* 100.0 (/ i (/ i n)))
(if (<= n 1.32e-86) (/ 0.0 (/ i n)) t_0))))))
double code(double i, double n) {
double t_0 = n * (100.0 + (i * (50.0 + (i * 16.666666666666668))));
double tmp;
if (n <= -8.5e-149) {
tmp = t_0;
} else if (n <= -1.12e-208) {
tmp = 0.0;
} else if (n <= -2.2e-232) {
tmp = 100.0 * (i / (i / n));
} else if (n <= 1.32e-86) {
tmp = 0.0 / (i / n);
} else {
tmp = t_0;
}
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 = n * (100.0d0 + (i * (50.0d0 + (i * 16.666666666666668d0))))
if (n <= (-8.5d-149)) then
tmp = t_0
else if (n <= (-1.12d-208)) then
tmp = 0.0d0
else if (n <= (-2.2d-232)) then
tmp = 100.0d0 * (i / (i / n))
else if (n <= 1.32d-86) then
tmp = 0.0d0 / (i / n)
else
tmp = t_0
end if
code = tmp
end function
public static double code(double i, double n) {
double t_0 = n * (100.0 + (i * (50.0 + (i * 16.666666666666668))));
double tmp;
if (n <= -8.5e-149) {
tmp = t_0;
} else if (n <= -1.12e-208) {
tmp = 0.0;
} else if (n <= -2.2e-232) {
tmp = 100.0 * (i / (i / n));
} else if (n <= 1.32e-86) {
tmp = 0.0 / (i / n);
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = n * (100.0 + (i * (50.0 + (i * 16.666666666666668)))) tmp = 0 if n <= -8.5e-149: tmp = t_0 elif n <= -1.12e-208: tmp = 0.0 elif n <= -2.2e-232: tmp = 100.0 * (i / (i / n)) elif n <= 1.32e-86: tmp = 0.0 / (i / n) else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(n * Float64(100.0 + Float64(i * Float64(50.0 + Float64(i * 16.666666666666668))))) tmp = 0.0 if (n <= -8.5e-149) tmp = t_0; elseif (n <= -1.12e-208) tmp = 0.0; elseif (n <= -2.2e-232) tmp = Float64(100.0 * Float64(i / Float64(i / n))); elseif (n <= 1.32e-86) tmp = Float64(0.0 / Float64(i / n)); else tmp = t_0; end return tmp end
function tmp_2 = code(i, n) t_0 = n * (100.0 + (i * (50.0 + (i * 16.666666666666668)))); tmp = 0.0; if (n <= -8.5e-149) tmp = t_0; elseif (n <= -1.12e-208) tmp = 0.0; elseif (n <= -2.2e-232) tmp = 100.0 * (i / (i / n)); elseif (n <= 1.32e-86) tmp = 0.0 / (i / n); else tmp = t_0; end tmp_2 = tmp; end
code[i_, n_] := Block[{t$95$0 = N[(n * N[(100.0 + N[(i * N[(50.0 + N[(i * 16.666666666666668), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -8.5e-149], t$95$0, If[LessEqual[n, -1.12e-208], 0.0, If[LessEqual[n, -2.2e-232], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.32e-86], N[(0.0 / N[(i / n), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := n \cdot \left(100 + i \cdot \left(50 + i \cdot 16.666666666666668\right)\right)\\
\mathbf{if}\;n \leq -8.5 \cdot 10^{-149}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq -1.12 \cdot 10^{-208}:\\
\;\;\;\;0\\
\mathbf{elif}\;n \leq -2.2 \cdot 10^{-232}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 1.32 \cdot 10^{-86}:\\
\;\;\;\;\frac{0}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -8.5000000000000006e-149 or 1.32e-86 < n Initial program 20.3%
associate-*r/20.3%
sub-neg20.3%
distribute-rgt-in20.3%
metadata-eval20.3%
metadata-eval20.3%
Simplified20.3%
Taylor expanded in n around inf 38.1%
Taylor expanded in i around 0 66.9%
Taylor expanded in n around 0 67.4%
if -8.5000000000000006e-149 < n < -1.12000000000000005e-208Initial program 73.1%
associate-*r/73.1%
sub-neg73.1%
distribute-rgt-in73.1%
metadata-eval73.1%
metadata-eval73.1%
Simplified73.1%
Taylor expanded in i around 0 73.1%
Taylor expanded in i around 0 73.1%
if -1.12000000000000005e-208 < n < -2.20000000000000002e-232Initial program 6.3%
Taylor expanded in i around 0 69.0%
if -2.20000000000000002e-232 < n < 1.32e-86Initial program 54.3%
associate-*r/54.3%
sub-neg54.3%
distribute-rgt-in54.3%
metadata-eval54.3%
metadata-eval54.3%
Simplified54.3%
Taylor expanded in i around 0 76.1%
Final simplification69.5%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* 100.0 (/ i (/ i n)))))
(if (<= n -5e+59)
(* (* i n) (/ 100.0 i))
(if (<= n -7.8e-224)
t_0
(if (<= n 3.9e-160)
(/ 0.0 (/ i n))
(if (<= n 21000000000000.0) t_0 (* n (+ 100.0 (* i 50.0)))))))))
double code(double i, double n) {
double t_0 = 100.0 * (i / (i / n));
double tmp;
if (n <= -5e+59) {
tmp = (i * n) * (100.0 / i);
} else if (n <= -7.8e-224) {
tmp = t_0;
} else if (n <= 3.9e-160) {
tmp = 0.0 / (i / n);
} else if (n <= 21000000000000.0) {
tmp = t_0;
} else {
tmp = n * (100.0 + (i * 50.0));
}
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 * (i / (i / n))
if (n <= (-5d+59)) then
tmp = (i * n) * (100.0d0 / i)
else if (n <= (-7.8d-224)) then
tmp = t_0
else if (n <= 3.9d-160) then
tmp = 0.0d0 / (i / n)
else if (n <= 21000000000000.0d0) then
tmp = t_0
else
tmp = n * (100.0d0 + (i * 50.0d0))
end if
code = tmp
end function
public static double code(double i, double n) {
double t_0 = 100.0 * (i / (i / n));
double tmp;
if (n <= -5e+59) {
tmp = (i * n) * (100.0 / i);
} else if (n <= -7.8e-224) {
tmp = t_0;
} else if (n <= 3.9e-160) {
tmp = 0.0 / (i / n);
} else if (n <= 21000000000000.0) {
tmp = t_0;
} else {
tmp = n * (100.0 + (i * 50.0));
}
return tmp;
}
def code(i, n): t_0 = 100.0 * (i / (i / n)) tmp = 0 if n <= -5e+59: tmp = (i * n) * (100.0 / i) elif n <= -7.8e-224: tmp = t_0 elif n <= 3.9e-160: tmp = 0.0 / (i / n) elif n <= 21000000000000.0: tmp = t_0 else: tmp = n * (100.0 + (i * 50.0)) return tmp
function code(i, n) t_0 = Float64(100.0 * Float64(i / Float64(i / n))) tmp = 0.0 if (n <= -5e+59) tmp = Float64(Float64(i * n) * Float64(100.0 / i)); elseif (n <= -7.8e-224) tmp = t_0; elseif (n <= 3.9e-160) tmp = Float64(0.0 / Float64(i / n)); elseif (n <= 21000000000000.0) tmp = t_0; else tmp = Float64(n * Float64(100.0 + Float64(i * 50.0))); end return tmp end
function tmp_2 = code(i, n) t_0 = 100.0 * (i / (i / n)); tmp = 0.0; if (n <= -5e+59) tmp = (i * n) * (100.0 / i); elseif (n <= -7.8e-224) tmp = t_0; elseif (n <= 3.9e-160) tmp = 0.0 / (i / n); elseif (n <= 21000000000000.0) tmp = t_0; else tmp = n * (100.0 + (i * 50.0)); end tmp_2 = tmp; end
code[i_, n_] := Block[{t$95$0 = N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -5e+59], N[(N[(i * n), $MachinePrecision] * N[(100.0 / i), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, -7.8e-224], t$95$0, If[LessEqual[n, 3.9e-160], N[(0.0 / N[(i / n), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 21000000000000.0], t$95$0, N[(n * N[(100.0 + N[(i * 50.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{if}\;n \leq -5 \cdot 10^{+59}:\\
\;\;\;\;\left(i \cdot n\right) \cdot \frac{100}{i}\\
\mathbf{elif}\;n \leq -7.8 \cdot 10^{-224}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 3.9 \cdot 10^{-160}:\\
\;\;\;\;\frac{0}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 21000000000000:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;n \cdot \left(100 + i \cdot 50\right)\\
\end{array}
\end{array}
if n < -4.9999999999999997e59Initial program 17.1%
Taylor expanded in n around inf 44.7%
expm1-define68.7%
Simplified68.7%
clear-num68.6%
un-div-inv68.7%
associate-/l/93.3%
Applied egg-rr93.3%
/-rgt-identity93.3%
associate-/r/93.2%
associate-/r/93.2%
associate-/r/93.4%
/-rgt-identity93.4%
*-commutative93.4%
Simplified93.4%
Taylor expanded in i around 0 58.5%
*-commutative58.5%
Simplified58.5%
if -4.9999999999999997e59 < n < -7.7999999999999996e-224 or 3.89999999999999989e-160 < n < 2.1e13Initial program 26.6%
Taylor expanded in i around 0 65.3%
if -7.7999999999999996e-224 < n < 3.89999999999999989e-160Initial program 62.0%
associate-*r/62.0%
sub-neg62.0%
distribute-rgt-in62.0%
metadata-eval62.0%
metadata-eval62.0%
Simplified62.0%
Taylor expanded in i around 0 80.6%
if 2.1e13 < n Initial program 21.9%
associate-*r/21.9%
sub-neg21.9%
distribute-rgt-in21.9%
metadata-eval21.9%
metadata-eval21.9%
Simplified21.9%
Taylor expanded in n around inf 47.4%
Taylor expanded in i around 0 77.8%
Taylor expanded in i around 0 73.7%
+-commutative73.7%
associate-*r*73.7%
distribute-rgt-out73.7%
*-commutative73.7%
Simplified73.7%
Final simplification68.6%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* 100.0 (/ i (/ i n)))))
(if (<= n -3.1e+59)
(* (* i n) (/ 100.0 i))
(if (<= n -9.5e-224)
t_0
(if (<= n 2.6e-161)
0.0
(if (<= n 21000000000000.0) t_0 (* n (+ 100.0 (* i 50.0)))))))))
double code(double i, double n) {
double t_0 = 100.0 * (i / (i / n));
double tmp;
if (n <= -3.1e+59) {
tmp = (i * n) * (100.0 / i);
} else if (n <= -9.5e-224) {
tmp = t_0;
} else if (n <= 2.6e-161) {
tmp = 0.0;
} else if (n <= 21000000000000.0) {
tmp = t_0;
} else {
tmp = n * (100.0 + (i * 50.0));
}
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 * (i / (i / n))
if (n <= (-3.1d+59)) then
tmp = (i * n) * (100.0d0 / i)
else if (n <= (-9.5d-224)) then
tmp = t_0
else if (n <= 2.6d-161) then
tmp = 0.0d0
else if (n <= 21000000000000.0d0) then
tmp = t_0
else
tmp = n * (100.0d0 + (i * 50.0d0))
end if
code = tmp
end function
public static double code(double i, double n) {
double t_0 = 100.0 * (i / (i / n));
double tmp;
if (n <= -3.1e+59) {
tmp = (i * n) * (100.0 / i);
} else if (n <= -9.5e-224) {
tmp = t_0;
} else if (n <= 2.6e-161) {
tmp = 0.0;
} else if (n <= 21000000000000.0) {
tmp = t_0;
} else {
tmp = n * (100.0 + (i * 50.0));
}
return tmp;
}
def code(i, n): t_0 = 100.0 * (i / (i / n)) tmp = 0 if n <= -3.1e+59: tmp = (i * n) * (100.0 / i) elif n <= -9.5e-224: tmp = t_0 elif n <= 2.6e-161: tmp = 0.0 elif n <= 21000000000000.0: tmp = t_0 else: tmp = n * (100.0 + (i * 50.0)) return tmp
function code(i, n) t_0 = Float64(100.0 * Float64(i / Float64(i / n))) tmp = 0.0 if (n <= -3.1e+59) tmp = Float64(Float64(i * n) * Float64(100.0 / i)); elseif (n <= -9.5e-224) tmp = t_0; elseif (n <= 2.6e-161) tmp = 0.0; elseif (n <= 21000000000000.0) tmp = t_0; else tmp = Float64(n * Float64(100.0 + Float64(i * 50.0))); end return tmp end
function tmp_2 = code(i, n) t_0 = 100.0 * (i / (i / n)); tmp = 0.0; if (n <= -3.1e+59) tmp = (i * n) * (100.0 / i); elseif (n <= -9.5e-224) tmp = t_0; elseif (n <= 2.6e-161) tmp = 0.0; elseif (n <= 21000000000000.0) tmp = t_0; else tmp = n * (100.0 + (i * 50.0)); end tmp_2 = tmp; end
code[i_, n_] := Block[{t$95$0 = N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -3.1e+59], N[(N[(i * n), $MachinePrecision] * N[(100.0 / i), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, -9.5e-224], t$95$0, If[LessEqual[n, 2.6e-161], 0.0, If[LessEqual[n, 21000000000000.0], t$95$0, N[(n * N[(100.0 + N[(i * 50.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{if}\;n \leq -3.1 \cdot 10^{+59}:\\
\;\;\;\;\left(i \cdot n\right) \cdot \frac{100}{i}\\
\mathbf{elif}\;n \leq -9.5 \cdot 10^{-224}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 2.6 \cdot 10^{-161}:\\
\;\;\;\;0\\
\mathbf{elif}\;n \leq 21000000000000:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;n \cdot \left(100 + i \cdot 50\right)\\
\end{array}
\end{array}
if n < -3.10000000000000015e59Initial program 17.1%
Taylor expanded in n around inf 44.7%
expm1-define68.7%
Simplified68.7%
clear-num68.6%
un-div-inv68.7%
associate-/l/93.3%
Applied egg-rr93.3%
/-rgt-identity93.3%
associate-/r/93.2%
associate-/r/93.2%
associate-/r/93.4%
/-rgt-identity93.4%
*-commutative93.4%
Simplified93.4%
Taylor expanded in i around 0 58.5%
*-commutative58.5%
Simplified58.5%
if -3.10000000000000015e59 < n < -9.5000000000000003e-224 or 2.59999999999999995e-161 < n < 2.1e13Initial program 26.6%
Taylor expanded in i around 0 65.3%
if -9.5000000000000003e-224 < n < 2.59999999999999995e-161Initial program 62.0%
associate-*r/62.0%
sub-neg62.0%
distribute-rgt-in62.0%
metadata-eval62.0%
metadata-eval62.0%
Simplified62.0%
Taylor expanded in i around 0 80.6%
Taylor expanded in i around 0 80.6%
if 2.1e13 < n Initial program 21.9%
associate-*r/21.9%
sub-neg21.9%
distribute-rgt-in21.9%
metadata-eval21.9%
metadata-eval21.9%
Simplified21.9%
Taylor expanded in n around inf 47.4%
Taylor expanded in i around 0 77.8%
Taylor expanded in i around 0 73.7%
+-commutative73.7%
associate-*r*73.7%
distribute-rgt-out73.7%
*-commutative73.7%
Simplified73.7%
Final simplification68.6%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* n (+ 100.0 (* i 50.0)))) (t_1 (* 100.0 (/ i (/ i n)))))
(if (<= n -6.2e+77)
t_0
(if (<= n -2e-244)
t_1
(if (<= n 3.4e-237) 0.0 (if (<= n 3.2e-77) t_1 t_0))))))
double code(double i, double n) {
double t_0 = n * (100.0 + (i * 50.0));
double t_1 = 100.0 * (i / (i / n));
double tmp;
if (n <= -6.2e+77) {
tmp = t_0;
} else if (n <= -2e-244) {
tmp = t_1;
} else if (n <= 3.4e-237) {
tmp = 0.0;
} else if (n <= 3.2e-77) {
tmp = t_1;
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = n * (100.0d0 + (i * 50.0d0))
t_1 = 100.0d0 * (i / (i / n))
if (n <= (-6.2d+77)) then
tmp = t_0
else if (n <= (-2d-244)) then
tmp = t_1
else if (n <= 3.4d-237) then
tmp = 0.0d0
else if (n <= 3.2d-77) then
tmp = t_1
else
tmp = t_0
end if
code = tmp
end function
public static double code(double i, double n) {
double t_0 = n * (100.0 + (i * 50.0));
double t_1 = 100.0 * (i / (i / n));
double tmp;
if (n <= -6.2e+77) {
tmp = t_0;
} else if (n <= -2e-244) {
tmp = t_1;
} else if (n <= 3.4e-237) {
tmp = 0.0;
} else if (n <= 3.2e-77) {
tmp = t_1;
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = n * (100.0 + (i * 50.0)) t_1 = 100.0 * (i / (i / n)) tmp = 0 if n <= -6.2e+77: tmp = t_0 elif n <= -2e-244: tmp = t_1 elif n <= 3.4e-237: tmp = 0.0 elif n <= 3.2e-77: tmp = t_1 else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(n * Float64(100.0 + Float64(i * 50.0))) t_1 = Float64(100.0 * Float64(i / Float64(i / n))) tmp = 0.0 if (n <= -6.2e+77) tmp = t_0; elseif (n <= -2e-244) tmp = t_1; elseif (n <= 3.4e-237) tmp = 0.0; elseif (n <= 3.2e-77) tmp = t_1; else tmp = t_0; end return tmp end
function tmp_2 = code(i, n) t_0 = n * (100.0 + (i * 50.0)); t_1 = 100.0 * (i / (i / n)); tmp = 0.0; if (n <= -6.2e+77) tmp = t_0; elseif (n <= -2e-244) tmp = t_1; elseif (n <= 3.4e-237) tmp = 0.0; elseif (n <= 3.2e-77) tmp = t_1; else tmp = t_0; end tmp_2 = tmp; end
code[i_, n_] := Block[{t$95$0 = N[(n * N[(100.0 + N[(i * 50.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -6.2e+77], t$95$0, If[LessEqual[n, -2e-244], t$95$1, If[LessEqual[n, 3.4e-237], 0.0, If[LessEqual[n, 3.2e-77], t$95$1, t$95$0]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := n \cdot \left(100 + i \cdot 50\right)\\
t_1 := 100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{if}\;n \leq -6.2 \cdot 10^{+77}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq -2 \cdot 10^{-244}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;n \leq 3.4 \cdot 10^{-237}:\\
\;\;\;\;0\\
\mathbf{elif}\;n \leq 3.2 \cdot 10^{-77}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -6.19999999999999997e77 or 3.2e-77 < n Initial program 18.5%
associate-*r/18.5%
sub-neg18.5%
distribute-rgt-in18.5%
metadata-eval18.5%
metadata-eval18.5%
Simplified18.5%
Taylor expanded in n around inf 41.6%
Taylor expanded in i around 0 70.6%
Taylor expanded in i around 0 68.0%
+-commutative68.0%
associate-*r*68.0%
distribute-rgt-out68.0%
*-commutative68.0%
Simplified68.0%
if -6.19999999999999997e77 < n < -1.9999999999999999e-244 or 3.4000000000000002e-237 < n < 3.2e-77Initial program 29.9%
Taylor expanded in i around 0 61.4%
if -1.9999999999999999e-244 < n < 3.4000000000000002e-237Initial program 82.2%
associate-*r/82.2%
sub-neg82.2%
distribute-rgt-in82.2%
metadata-eval82.2%
metadata-eval82.2%
Simplified82.2%
Taylor expanded in i around 0 93.3%
Taylor expanded in i around 0 93.3%
(FPCore (i n)
:precision binary64
(if (<= n -6.2e+77)
(* n 100.0)
(if (<= n -5.2e-247)
(* 100.0 (/ i (/ i n)))
(if (<= n 1.32e-86) 0.0 (* n 100.0)))))
double code(double i, double n) {
double tmp;
if (n <= -6.2e+77) {
tmp = n * 100.0;
} else if (n <= -5.2e-247) {
tmp = 100.0 * (i / (i / n));
} else if (n <= 1.32e-86) {
tmp = 0.0;
} 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 (n <= (-6.2d+77)) then
tmp = n * 100.0d0
else if (n <= (-5.2d-247)) then
tmp = 100.0d0 * (i / (i / n))
else if (n <= 1.32d-86) then
tmp = 0.0d0
else
tmp = n * 100.0d0
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (n <= -6.2e+77) {
tmp = n * 100.0;
} else if (n <= -5.2e-247) {
tmp = 100.0 * (i / (i / n));
} else if (n <= 1.32e-86) {
tmp = 0.0;
} else {
tmp = n * 100.0;
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -6.2e+77: tmp = n * 100.0 elif n <= -5.2e-247: tmp = 100.0 * (i / (i / n)) elif n <= 1.32e-86: tmp = 0.0 else: tmp = n * 100.0 return tmp
function code(i, n) tmp = 0.0 if (n <= -6.2e+77) tmp = Float64(n * 100.0); elseif (n <= -5.2e-247) tmp = Float64(100.0 * Float64(i / Float64(i / n))); elseif (n <= 1.32e-86) tmp = 0.0; else tmp = Float64(n * 100.0); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (n <= -6.2e+77) tmp = n * 100.0; elseif (n <= -5.2e-247) tmp = 100.0 * (i / (i / n)); elseif (n <= 1.32e-86) tmp = 0.0; else tmp = n * 100.0; end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[n, -6.2e+77], N[(n * 100.0), $MachinePrecision], If[LessEqual[n, -5.2e-247], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.32e-86], 0.0, N[(n * 100.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -6.2 \cdot 10^{+77}:\\
\;\;\;\;n \cdot 100\\
\mathbf{elif}\;n \leq -5.2 \cdot 10^{-247}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 1.32 \cdot 10^{-86}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;n \cdot 100\\
\end{array}
\end{array}
if n < -6.19999999999999997e77 or 1.32e-86 < n Initial program 18.4%
Taylor expanded in i around 0 55.3%
*-commutative55.3%
Simplified55.3%
if -6.19999999999999997e77 < n < -5.2e-247Initial program 33.8%
Taylor expanded in i around 0 57.3%
if -5.2e-247 < n < 1.32e-86Initial program 53.4%
associate-*r/53.4%
sub-neg53.4%
distribute-rgt-in53.4%
metadata-eval53.4%
metadata-eval53.4%
Simplified53.4%
Taylor expanded in i around 0 76.5%
Taylor expanded in i around 0 76.5%
(FPCore (i n) :precision binary64 (if (or (<= n -1.45e-134) (not (<= n 2.4e-83))) (* n 100.0) 0.0))
double code(double i, double n) {
double tmp;
if ((n <= -1.45e-134) || !(n <= 2.4e-83)) {
tmp = n * 100.0;
} else {
tmp = 0.0;
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if ((n <= (-1.45d-134)) .or. (.not. (n <= 2.4d-83))) then
tmp = n * 100.0d0
else
tmp = 0.0d0
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if ((n <= -1.45e-134) || !(n <= 2.4e-83)) {
tmp = n * 100.0;
} else {
tmp = 0.0;
}
return tmp;
}
def code(i, n): tmp = 0 if (n <= -1.45e-134) or not (n <= 2.4e-83): tmp = n * 100.0 else: tmp = 0.0 return tmp
function code(i, n) tmp = 0.0 if ((n <= -1.45e-134) || !(n <= 2.4e-83)) tmp = Float64(n * 100.0); else tmp = 0.0; end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if ((n <= -1.45e-134) || ~((n <= 2.4e-83))) tmp = n * 100.0; else tmp = 0.0; end tmp_2 = tmp; end
code[i_, n_] := If[Or[LessEqual[n, -1.45e-134], N[Not[LessEqual[n, 2.4e-83]], $MachinePrecision]], N[(n * 100.0), $MachinePrecision], 0.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.45 \cdot 10^{-134} \lor \neg \left(n \leq 2.4 \cdot 10^{-83}\right):\\
\;\;\;\;n \cdot 100\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}
\end{array}
if n < -1.44999999999999997e-134 or 2.4000000000000001e-83 < n Initial program 20.6%
Taylor expanded in i around 0 55.6%
*-commutative55.6%
Simplified55.6%
if -1.44999999999999997e-134 < n < 2.4000000000000001e-83Initial program 51.2%
associate-*r/51.2%
sub-neg51.2%
distribute-rgt-in51.2%
metadata-eval51.2%
metadata-eval51.2%
Simplified51.2%
Taylor expanded in i around 0 70.0%
Taylor expanded in i around 0 70.0%
Final simplification59.5%
(FPCore (i n) :precision binary64 0.0)
double code(double i, double n) {
return 0.0;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
code = 0.0d0
end function
public static double code(double i, double n) {
return 0.0;
}
def code(i, n): return 0.0
function code(i, n) return 0.0 end
function tmp = code(i, n) tmp = 0.0; end
code[i_, n_] := 0.0
\begin{array}{l}
\\
0
\end{array}
Initial program 28.9%
associate-*r/28.9%
sub-neg28.9%
distribute-rgt-in28.9%
metadata-eval28.9%
metadata-eval28.9%
Simplified28.9%
Taylor expanded in i around 0 21.3%
Taylor expanded in i around 0 21.5%
(FPCore (i n)
:precision binary64
(let* ((t_0 (+ 1.0 (/ i n))))
(*
100.0
(/
(-
(exp
(*
n
(if (== t_0 1.0)
(/ i n)
(/ (* (/ i n) (log t_0)) (- (+ (/ i n) 1.0) 1.0)))))
1.0)
(/ i n)))))
double code(double i, double n) {
double t_0 = 1.0 + (i / n);
double tmp;
if (t_0 == 1.0) {
tmp = i / n;
} else {
tmp = ((i / n) * log(t_0)) / (((i / n) + 1.0) - 1.0);
}
return 100.0 * ((exp((n * tmp)) - 1.0) / (i / n));
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: t_0
real(8) :: tmp
t_0 = 1.0d0 + (i / n)
if (t_0 == 1.0d0) then
tmp = i / n
else
tmp = ((i / n) * log(t_0)) / (((i / n) + 1.0d0) - 1.0d0)
end if
code = 100.0d0 * ((exp((n * tmp)) - 1.0d0) / (i / n))
end function
public static double code(double i, double n) {
double t_0 = 1.0 + (i / n);
double tmp;
if (t_0 == 1.0) {
tmp = i / n;
} else {
tmp = ((i / n) * Math.log(t_0)) / (((i / n) + 1.0) - 1.0);
}
return 100.0 * ((Math.exp((n * tmp)) - 1.0) / (i / n));
}
def code(i, n): t_0 = 1.0 + (i / n) tmp = 0 if t_0 == 1.0: tmp = i / n else: tmp = ((i / n) * math.log(t_0)) / (((i / n) + 1.0) - 1.0) return 100.0 * ((math.exp((n * tmp)) - 1.0) / (i / n))
function code(i, n) t_0 = Float64(1.0 + Float64(i / n)) tmp = 0.0 if (t_0 == 1.0) tmp = Float64(i / n); else tmp = Float64(Float64(Float64(i / n) * log(t_0)) / Float64(Float64(Float64(i / n) + 1.0) - 1.0)); end return Float64(100.0 * Float64(Float64(exp(Float64(n * tmp)) - 1.0) / Float64(i / n))) end
function tmp_2 = code(i, n) t_0 = 1.0 + (i / n); tmp = 0.0; if (t_0 == 1.0) tmp = i / n; else tmp = ((i / n) * log(t_0)) / (((i / n) + 1.0) - 1.0); end tmp_2 = 100.0 * ((exp((n * tmp)) - 1.0) / (i / n)); end
code[i_, n_] := Block[{t$95$0 = N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision]}, N[(100.0 * N[(N[(N[Exp[N[(n * If[Equal[t$95$0, 1.0], N[(i / n), $MachinePrecision], N[(N[(N[(i / n), $MachinePrecision] * N[Log[t$95$0], $MachinePrecision]), $MachinePrecision] / N[(N[(N[(i / n), $MachinePrecision] + 1.0), $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]], $MachinePrecision] - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 1 + \frac{i}{n}\\
100 \cdot \frac{e^{n \cdot \begin{array}{l}
\mathbf{if}\;t\_0 = 1:\\
\;\;\;\;\frac{i}{n}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{i}{n} \cdot \log t\_0}{\left(\frac{i}{n} + 1\right) - 1}\\
\end{array}} - 1}{\frac{i}{n}}
\end{array}
\end{array}
herbie shell --seed 2024107
(FPCore (i n)
:name "Compound Interest"
:precision binary64
:alt
(* 100.0 (/ (- (exp (* n (if (== (+ 1.0 (/ i n)) 1.0) (/ i n) (/ (* (/ i n) (log (+ 1.0 (/ i n)))) (- (+ (/ i n) 1.0) 1.0))))) 1.0) (/ i n)))
(* 100.0 (/ (- (pow (+ 1.0 (/ i n)) n) 1.0) (/ i n))))