
(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 22 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 -1e-61)
(* (/ -100.0 (/ i n)) (- 1.0 t_0))
(if (<= t_1 0.0)
(* 100.0 (/ (expm1 (* n (log1p (/ i n)))) (/ i n)))
(if (<= t_1 INFINITY)
(* 100.0 (/ (- (/ t_0 (/ 1.0 n)) n) i))
(* n 100.0))))))
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 <= -1e-61) {
tmp = (-100.0 / (i / n)) * (1.0 - t_0);
} else if (t_1 <= 0.0) {
tmp = 100.0 * (expm1((n * log1p((i / n)))) / (i / n));
} else if (t_1 <= ((double) INFINITY)) {
tmp = 100.0 * (((t_0 / (1.0 / n)) - n) / i);
} else {
tmp = n * 100.0;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = Math.pow((1.0 + (i / n)), n);
double t_1 = (t_0 + -1.0) / (i / n);
double tmp;
if (t_1 <= -1e-61) {
tmp = (-100.0 / (i / n)) * (1.0 - t_0);
} else if (t_1 <= 0.0) {
tmp = 100.0 * (Math.expm1((n * Math.log1p((i / n)))) / (i / n));
} else if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = 100.0 * (((t_0 / (1.0 / n)) - n) / i);
} else {
tmp = n * 100.0;
}
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 <= -1e-61: tmp = (-100.0 / (i / n)) * (1.0 - t_0) elif t_1 <= 0.0: tmp = 100.0 * (math.expm1((n * math.log1p((i / n)))) / (i / n)) elif t_1 <= math.inf: tmp = 100.0 * (((t_0 / (1.0 / n)) - n) / i) else: tmp = n * 100.0 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 <= -1e-61) tmp = Float64(Float64(-100.0 / Float64(i / n)) * Float64(1.0 - t_0)); elseif (t_1 <= 0.0) tmp = Float64(100.0 * Float64(expm1(Float64(n * log1p(Float64(i / n)))) / Float64(i / n))); elseif (t_1 <= Inf) tmp = Float64(100.0 * Float64(Float64(Float64(t_0 / Float64(1.0 / n)) - n) / i)); else tmp = Float64(n * 100.0); 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, -1e-61], N[(N[(-100.0 / N[(i / n), $MachinePrecision]), $MachinePrecision] * N[(1.0 - t$95$0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(100.0 * N[(N[(Exp[N[(n * N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]] - 1), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(100.0 * N[(N[(N[(t$95$0 / N[(1.0 / n), $MachinePrecision]), $MachinePrecision] - n), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], N[(n * 100.0), $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 -1 \cdot 10^{-61}:\\
\;\;\;\;\frac{-100}{\frac{i}{n}} \cdot \left(1 - t\_0\right)\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;100 \cdot \frac{\mathsf{expm1}\left(n \cdot \mathsf{log1p}\left(\frac{i}{n}\right)\right)}{\frac{i}{n}}\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;100 \cdot \frac{\frac{t\_0}{\frac{1}{n}} - n}{i}\\
\mathbf{else}:\\
\;\;\;\;n \cdot 100\\
\end{array}
\end{array}
if (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < -1e-61Initial program 99.9%
*-commutativeN/A
associate-*l/N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-lft-neg-outN/A
distribute-neg-fracN/A
distribute-neg-frac2N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
distribute-neg-frac2N/A
*-commutativeN/A
*-lowering-*.f64N/A
Simplified100.0%
if -1e-61 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < 0.0Initial program 21.9%
pow-to-expN/A
expm1-defineN/A
expm1-lowering-expm1.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
/-lowering-/.f6499.7%
Applied egg-rr99.7%
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 91.7%
pow-to-expN/A
expm1-defineN/A
expm1-lowering-expm1.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
/-lowering-/.f6443.3%
Applied egg-rr43.3%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f6443.3%
Applied egg-rr43.3%
div-subN/A
associate-/r/N/A
associate-/r*N/A
remove-double-divN/A
sub-divN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
exp-to-powN/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6492.2%
Applied egg-rr92.2%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
Taylor expanded in i around 0
*-lowering-*.f6476.6%
Simplified76.6%
Final simplification94.7%
(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 -2e-125)
(* t_1 100.0)
(if (<= t_1 0.0)
(* 100.0 (/ (expm1 i) (/ i n)))
(if (<= t_1 INFINITY)
(* 100.0 (/ (- (/ t_0 (/ 1.0 n)) n) i))
(* n 100.0))))))
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 <= -2e-125) {
tmp = t_1 * 100.0;
} else if (t_1 <= 0.0) {
tmp = 100.0 * (expm1(i) / (i / n));
} else if (t_1 <= ((double) INFINITY)) {
tmp = 100.0 * (((t_0 / (1.0 / n)) - n) / i);
} else {
tmp = n * 100.0;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = Math.pow((1.0 + (i / n)), n);
double t_1 = (t_0 + -1.0) / (i / n);
double tmp;
if (t_1 <= -2e-125) {
tmp = t_1 * 100.0;
} else if (t_1 <= 0.0) {
tmp = 100.0 * (Math.expm1(i) / (i / n));
} else if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = 100.0 * (((t_0 / (1.0 / n)) - n) / i);
} else {
tmp = n * 100.0;
}
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 <= -2e-125: tmp = t_1 * 100.0 elif t_1 <= 0.0: tmp = 100.0 * (math.expm1(i) / (i / n)) elif t_1 <= math.inf: tmp = 100.0 * (((t_0 / (1.0 / n)) - n) / i) else: tmp = n * 100.0 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 <= -2e-125) tmp = Float64(t_1 * 100.0); elseif (t_1 <= 0.0) tmp = Float64(100.0 * Float64(expm1(i) / Float64(i / n))); elseif (t_1 <= Inf) tmp = Float64(100.0 * Float64(Float64(Float64(t_0 / Float64(1.0 / n)) - n) / i)); else tmp = Float64(n * 100.0); 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, -2e-125], N[(t$95$1 * 100.0), $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(100.0 * N[(N[(Exp[i] - 1), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(100.0 * N[(N[(N[(t$95$0 / N[(1.0 / n), $MachinePrecision]), $MachinePrecision] - n), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], N[(n * 100.0), $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 -2 \cdot 10^{-125}:\\
\;\;\;\;t\_1 \cdot 100\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;100 \cdot \frac{\mathsf{expm1}\left(i\right)}{\frac{i}{n}}\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;100 \cdot \frac{\frac{t\_0}{\frac{1}{n}} - n}{i}\\
\mathbf{else}:\\
\;\;\;\;n \cdot 100\\
\end{array}
\end{array}
if (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < -2.00000000000000002e-125Initial program 97.4%
if -2.00000000000000002e-125 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < 0.0Initial program 20.7%
Taylor expanded in n around inf
expm1-defineN/A
expm1-lowering-expm1.f6475.4%
Simplified75.4%
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 91.7%
pow-to-expN/A
expm1-defineN/A
expm1-lowering-expm1.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
/-lowering-/.f6443.3%
Applied egg-rr43.3%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f6443.3%
Applied egg-rr43.3%
div-subN/A
associate-/r/N/A
associate-/r*N/A
remove-double-divN/A
sub-divN/A
/-lowering-/.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-commutativeN/A
exp-to-powN/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6492.2%
Applied egg-rr92.2%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
Taylor expanded in i around 0
*-lowering-*.f6476.6%
Simplified76.6%
Final simplification78.7%
(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 -2e-125)
(* t_1 100.0)
(if (<= t_1 0.0)
(* 100.0 (/ (expm1 i) (/ i n)))
(if (<= t_1 INFINITY)
(* 100.0 (- (/ t_0 (/ i n)) (/ n i)))
(* n 100.0))))))
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 <= -2e-125) {
tmp = t_1 * 100.0;
} else if (t_1 <= 0.0) {
tmp = 100.0 * (expm1(i) / (i / n));
} else if (t_1 <= ((double) INFINITY)) {
tmp = 100.0 * ((t_0 / (i / n)) - (n / i));
} else {
tmp = n * 100.0;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = Math.pow((1.0 + (i / n)), n);
double t_1 = (t_0 + -1.0) / (i / n);
double tmp;
if (t_1 <= -2e-125) {
tmp = t_1 * 100.0;
} else if (t_1 <= 0.0) {
tmp = 100.0 * (Math.expm1(i) / (i / n));
} else if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = 100.0 * ((t_0 / (i / n)) - (n / i));
} else {
tmp = n * 100.0;
}
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 <= -2e-125: tmp = t_1 * 100.0 elif t_1 <= 0.0: tmp = 100.0 * (math.expm1(i) / (i / n)) elif t_1 <= math.inf: tmp = 100.0 * ((t_0 / (i / n)) - (n / i)) else: tmp = n * 100.0 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 <= -2e-125) tmp = Float64(t_1 * 100.0); elseif (t_1 <= 0.0) tmp = Float64(100.0 * Float64(expm1(i) / Float64(i / n))); elseif (t_1 <= Inf) tmp = Float64(100.0 * Float64(Float64(t_0 / Float64(i / n)) - Float64(n / i))); else tmp = Float64(n * 100.0); 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, -2e-125], N[(t$95$1 * 100.0), $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(100.0 * N[(N[(Exp[i] - 1), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(100.0 * N[(N[(t$95$0 / N[(i / n), $MachinePrecision]), $MachinePrecision] - N[(n / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(n * 100.0), $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 -2 \cdot 10^{-125}:\\
\;\;\;\;t\_1 \cdot 100\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;100 \cdot \frac{\mathsf{expm1}\left(i\right)}{\frac{i}{n}}\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;100 \cdot \left(\frac{t\_0}{\frac{i}{n}} - \frac{n}{i}\right)\\
\mathbf{else}:\\
\;\;\;\;n \cdot 100\\
\end{array}
\end{array}
if (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < -2.00000000000000002e-125Initial program 97.4%
if -2.00000000000000002e-125 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < 0.0Initial program 20.7%
Taylor expanded in n around inf
expm1-defineN/A
expm1-lowering-expm1.f6475.4%
Simplified75.4%
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 91.7%
div-subN/A
clear-numN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6492.2%
Applied egg-rr92.2%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
Taylor expanded in i around 0
*-lowering-*.f6476.6%
Simplified76.6%
Final simplification78.7%
(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 -2e-125)
(* t_1 100.0)
(if (<= t_1 0.0)
(* 100.0 (/ (expm1 i) (/ i n)))
(if (<= t_1 INFINITY)
(* (/ -100.0 (/ i n)) (- 1.0 t_0))
(* n 100.0))))))
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 <= -2e-125) {
tmp = t_1 * 100.0;
} else if (t_1 <= 0.0) {
tmp = 100.0 * (expm1(i) / (i / n));
} else if (t_1 <= ((double) INFINITY)) {
tmp = (-100.0 / (i / n)) * (1.0 - t_0);
} else {
tmp = n * 100.0;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = Math.pow((1.0 + (i / n)), n);
double t_1 = (t_0 + -1.0) / (i / n);
double tmp;
if (t_1 <= -2e-125) {
tmp = t_1 * 100.0;
} else if (t_1 <= 0.0) {
tmp = 100.0 * (Math.expm1(i) / (i / n));
} else if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = (-100.0 / (i / n)) * (1.0 - t_0);
} else {
tmp = n * 100.0;
}
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 <= -2e-125: tmp = t_1 * 100.0 elif t_1 <= 0.0: tmp = 100.0 * (math.expm1(i) / (i / n)) elif t_1 <= math.inf: tmp = (-100.0 / (i / n)) * (1.0 - t_0) else: tmp = n * 100.0 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 <= -2e-125) tmp = Float64(t_1 * 100.0); elseif (t_1 <= 0.0) tmp = Float64(100.0 * Float64(expm1(i) / Float64(i / n))); elseif (t_1 <= Inf) tmp = Float64(Float64(-100.0 / Float64(i / n)) * Float64(1.0 - t_0)); else tmp = Float64(n * 100.0); 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, -2e-125], N[(t$95$1 * 100.0), $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(100.0 * N[(N[(Exp[i] - 1), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(N[(-100.0 / N[(i / n), $MachinePrecision]), $MachinePrecision] * N[(1.0 - t$95$0), $MachinePrecision]), $MachinePrecision], N[(n * 100.0), $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 -2 \cdot 10^{-125}:\\
\;\;\;\;t\_1 \cdot 100\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;100 \cdot \frac{\mathsf{expm1}\left(i\right)}{\frac{i}{n}}\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;\frac{-100}{\frac{i}{n}} \cdot \left(1 - t\_0\right)\\
\mathbf{else}:\\
\;\;\;\;n \cdot 100\\
\end{array}
\end{array}
if (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < -2.00000000000000002e-125Initial program 97.4%
if -2.00000000000000002e-125 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < 0.0Initial program 20.7%
Taylor expanded in n around inf
expm1-defineN/A
expm1-lowering-expm1.f6475.4%
Simplified75.4%
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 91.7%
*-commutativeN/A
associate-*l/N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-lft-neg-outN/A
distribute-neg-fracN/A
distribute-neg-frac2N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
distribute-neg-frac2N/A
*-commutativeN/A
*-lowering-*.f64N/A
Simplified91.9%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
Taylor expanded in i around 0
*-lowering-*.f6476.6%
Simplified76.6%
Final simplification78.7%
(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 -2e-125)
(* t_1 100.0)
(if (<= t_1 0.0)
(* 100.0 (/ (expm1 i) (/ i n)))
(if (<= t_1 INFINITY) (* n (* 100.0 (/ t_0 i))) (* n 100.0))))))
double code(double i, double n) {
double t_0 = pow((1.0 + (i / n)), n) + -1.0;
double t_1 = t_0 / (i / n);
double tmp;
if (t_1 <= -2e-125) {
tmp = t_1 * 100.0;
} else if (t_1 <= 0.0) {
tmp = 100.0 * (expm1(i) / (i / n));
} else if (t_1 <= ((double) INFINITY)) {
tmp = n * (100.0 * (t_0 / i));
} else {
tmp = n * 100.0;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = Math.pow((1.0 + (i / n)), n) + -1.0;
double t_1 = t_0 / (i / n);
double tmp;
if (t_1 <= -2e-125) {
tmp = t_1 * 100.0;
} else if (t_1 <= 0.0) {
tmp = 100.0 * (Math.expm1(i) / (i / n));
} else if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = n * (100.0 * (t_0 / i));
} else {
tmp = n * 100.0;
}
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 <= -2e-125: tmp = t_1 * 100.0 elif t_1 <= 0.0: tmp = 100.0 * (math.expm1(i) / (i / n)) elif t_1 <= math.inf: tmp = n * (100.0 * (t_0 / i)) else: tmp = n * 100.0 return tmp
function code(i, n) t_0 = Float64((Float64(1.0 + Float64(i / n)) ^ n) + -1.0) t_1 = Float64(t_0 / Float64(i / n)) tmp = 0.0 if (t_1 <= -2e-125) tmp = Float64(t_1 * 100.0); elseif (t_1 <= 0.0) tmp = Float64(100.0 * Float64(expm1(i) / Float64(i / n))); elseif (t_1 <= Inf) tmp = Float64(n * Float64(100.0 * Float64(t_0 / i))); else tmp = Float64(n * 100.0); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[Power[N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision], n], $MachinePrecision] + -1.0), $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 / N[(i / n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -2e-125], N[(t$95$1 * 100.0), $MachinePrecision], If[LessEqual[t$95$1, 0.0], N[(100.0 * N[(N[(Exp[i] - 1), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(n * N[(100.0 * N[(t$95$0 / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(n * 100.0), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(1 + \frac{i}{n}\right)}^{n} + -1\\
t_1 := \frac{t\_0}{\frac{i}{n}}\\
\mathbf{if}\;t\_1 \leq -2 \cdot 10^{-125}:\\
\;\;\;\;t\_1 \cdot 100\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;100 \cdot \frac{\mathsf{expm1}\left(i\right)}{\frac{i}{n}}\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;n \cdot \left(100 \cdot \frac{t\_0}{i}\right)\\
\mathbf{else}:\\
\;\;\;\;n \cdot 100\\
\end{array}
\end{array}
if (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < -2.00000000000000002e-125Initial program 97.4%
if -2.00000000000000002e-125 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < 0.0Initial program 20.7%
Taylor expanded in n around inf
expm1-defineN/A
expm1-lowering-expm1.f6475.4%
Simplified75.4%
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 91.7%
*-commutativeN/A
associate-/r/N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sub-negN/A
+-lowering-+.f64N/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
metadata-eval91.9%
Applied egg-rr91.9%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
Taylor expanded in i around 0
*-lowering-*.f6476.6%
Simplified76.6%
Final simplification78.7%
(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 -2e-125)
t_1
(if (<= t_0 0.0)
(* 100.0 (/ (expm1 i) (/ i n)))
(if (<= t_0 INFINITY) t_1 (* n 100.0))))))
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 <= -2e-125) {
tmp = t_1;
} else if (t_0 <= 0.0) {
tmp = 100.0 * (expm1(i) / (i / n));
} else if (t_0 <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = n * 100.0;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = (Math.pow((1.0 + (i / n)), n) + -1.0) / (i / n);
double t_1 = t_0 * 100.0;
double tmp;
if (t_0 <= -2e-125) {
tmp = t_1;
} else if (t_0 <= 0.0) {
tmp = 100.0 * (Math.expm1(i) / (i / n));
} else if (t_0 <= Double.POSITIVE_INFINITY) {
tmp = t_1;
} else {
tmp = n * 100.0;
}
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 <= -2e-125: tmp = t_1 elif t_0 <= 0.0: tmp = 100.0 * (math.expm1(i) / (i / n)) elif t_0 <= math.inf: tmp = t_1 else: tmp = n * 100.0 return tmp
function code(i, n) t_0 = Float64(Float64((Float64(1.0 + Float64(i / n)) ^ n) + -1.0) / Float64(i / n)) t_1 = Float64(t_0 * 100.0) tmp = 0.0 if (t_0 <= -2e-125) tmp = t_1; elseif (t_0 <= 0.0) tmp = Float64(100.0 * Float64(expm1(i) / Float64(i / n))); elseif (t_0 <= Inf) tmp = t_1; else tmp = Float64(n * 100.0); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(N[Power[N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision], n], $MachinePrecision] + -1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 * 100.0), $MachinePrecision]}, If[LessEqual[t$95$0, -2e-125], t$95$1, If[LessEqual[t$95$0, 0.0], N[(100.0 * N[(N[(Exp[i] - 1), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, Infinity], t$95$1, N[(n * 100.0), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{{\left(1 + \frac{i}{n}\right)}^{n} + -1}{\frac{i}{n}}\\
t_1 := t\_0 \cdot 100\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{-125}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 0:\\
\;\;\;\;100 \cdot \frac{\mathsf{expm1}\left(i\right)}{\frac{i}{n}}\\
\mathbf{elif}\;t\_0 \leq \infty:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;n \cdot 100\\
\end{array}
\end{array}
if (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < -2.00000000000000002e-125 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 94.2%
if -2.00000000000000002e-125 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < 0.0Initial program 20.7%
Taylor expanded in n around inf
expm1-defineN/A
expm1-lowering-expm1.f6475.4%
Simplified75.4%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
Taylor expanded in i around 0
*-lowering-*.f6476.6%
Simplified76.6%
Final simplification78.7%
(FPCore (i n)
:precision binary64
(let* ((t_0 (/ (* 100.0 (* n (expm1 i))) i)))
(if (<= n -7.6e+28)
t_0
(if (<= n -3.5e-194)
(* 100.0 (/ (expm1 i) (/ i n)))
(if (<= n 9e-217)
(* n (/ (+ -100.0 100.0) i))
(if (<= n 2.15) (* 100.0 (/ i (/ i n))) t_0))))))
double code(double i, double n) {
double t_0 = (100.0 * (n * expm1(i))) / i;
double tmp;
if (n <= -7.6e+28) {
tmp = t_0;
} else if (n <= -3.5e-194) {
tmp = 100.0 * (expm1(i) / (i / n));
} else if (n <= 9e-217) {
tmp = n * ((-100.0 + 100.0) / i);
} else if (n <= 2.15) {
tmp = 100.0 * (i / (i / n));
} 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 <= -7.6e+28) {
tmp = t_0;
} else if (n <= -3.5e-194) {
tmp = 100.0 * (Math.expm1(i) / (i / n));
} else if (n <= 9e-217) {
tmp = n * ((-100.0 + 100.0) / i);
} else if (n <= 2.15) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = (100.0 * (n * math.expm1(i))) / i tmp = 0 if n <= -7.6e+28: tmp = t_0 elif n <= -3.5e-194: tmp = 100.0 * (math.expm1(i) / (i / n)) elif n <= 9e-217: tmp = n * ((-100.0 + 100.0) / i) elif n <= 2.15: tmp = 100.0 * (i / (i / n)) else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(Float64(100.0 * Float64(n * expm1(i))) / i) tmp = 0.0 if (n <= -7.6e+28) tmp = t_0; elseif (n <= -3.5e-194) tmp = Float64(100.0 * Float64(expm1(i) / Float64(i / n))); elseif (n <= 9e-217) tmp = Float64(n * Float64(Float64(-100.0 + 100.0) / i)); elseif (n <= 2.15) 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[(100.0 * N[(n * N[(Exp[i] - 1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]}, If[LessEqual[n, -7.6e+28], t$95$0, If[LessEqual[n, -3.5e-194], N[(100.0 * N[(N[(Exp[i] - 1), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 9e-217], N[(n * N[(N[(-100.0 + 100.0), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 2.15], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{100 \cdot \left(n \cdot \mathsf{expm1}\left(i\right)\right)}{i}\\
\mathbf{if}\;n \leq -7.6 \cdot 10^{+28}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq -3.5 \cdot 10^{-194}:\\
\;\;\;\;100 \cdot \frac{\mathsf{expm1}\left(i\right)}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 9 \cdot 10^{-217}:\\
\;\;\;\;n \cdot \frac{-100 + 100}{i}\\
\mathbf{elif}\;n \leq 2.15:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -7.5999999999999998e28 or 2.14999999999999991 < n Initial program 25.6%
Taylor expanded in n around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6485.6%
Simplified85.6%
if -7.5999999999999998e28 < n < -3.5000000000000003e-194Initial program 21.1%
Taylor expanded in n around inf
expm1-defineN/A
expm1-lowering-expm1.f6473.0%
Simplified73.0%
if -3.5000000000000003e-194 < n < 8.9999999999999997e-217Initial program 63.4%
pow-to-expN/A
expm1-defineN/A
expm1-lowering-expm1.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
/-lowering-/.f6484.5%
Applied egg-rr84.5%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr63.7%
Taylor expanded in i around 0
Simplified78.8%
if 8.9999999999999997e-217 < n < 2.14999999999999991Initial program 13.2%
Taylor expanded in i around 0
Simplified69.8%
Final simplification80.0%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* (* n (expm1 i)) (/ 100.0 i))))
(if (<= n -156000.0)
t_0
(if (<= n -1.02e-199)
(* 100.0 (/ (expm1 i) (/ i n)))
(if (<= n 3e-214)
(* n (/ (+ -100.0 100.0) i))
(if (<= n 2.7) (* 100.0 (/ i (/ i n))) t_0))))))
double code(double i, double n) {
double t_0 = (n * expm1(i)) * (100.0 / i);
double tmp;
if (n <= -156000.0) {
tmp = t_0;
} else if (n <= -1.02e-199) {
tmp = 100.0 * (expm1(i) / (i / n));
} else if (n <= 3e-214) {
tmp = n * ((-100.0 + 100.0) / i);
} else if (n <= 2.7) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = t_0;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = (n * Math.expm1(i)) * (100.0 / i);
double tmp;
if (n <= -156000.0) {
tmp = t_0;
} else if (n <= -1.02e-199) {
tmp = 100.0 * (Math.expm1(i) / (i / n));
} else if (n <= 3e-214) {
tmp = n * ((-100.0 + 100.0) / i);
} else if (n <= 2.7) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = (n * math.expm1(i)) * (100.0 / i) tmp = 0 if n <= -156000.0: tmp = t_0 elif n <= -1.02e-199: tmp = 100.0 * (math.expm1(i) / (i / n)) elif n <= 3e-214: tmp = n * ((-100.0 + 100.0) / i) elif n <= 2.7: tmp = 100.0 * (i / (i / n)) else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(Float64(n * expm1(i)) * Float64(100.0 / i)) tmp = 0.0 if (n <= -156000.0) tmp = t_0; elseif (n <= -1.02e-199) tmp = Float64(100.0 * Float64(expm1(i) / Float64(i / n))); elseif (n <= 3e-214) tmp = Float64(n * Float64(Float64(-100.0 + 100.0) / i)); elseif (n <= 2.7) 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[(Exp[i] - 1), $MachinePrecision]), $MachinePrecision] * N[(100.0 / i), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -156000.0], t$95$0, If[LessEqual[n, -1.02e-199], N[(100.0 * N[(N[(Exp[i] - 1), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 3e-214], N[(n * N[(N[(-100.0 + 100.0), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 2.7], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(n \cdot \mathsf{expm1}\left(i\right)\right) \cdot \frac{100}{i}\\
\mathbf{if}\;n \leq -156000:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq -1.02 \cdot 10^{-199}:\\
\;\;\;\;100 \cdot \frac{\mathsf{expm1}\left(i\right)}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 3 \cdot 10^{-214}:\\
\;\;\;\;n \cdot \frac{-100 + 100}{i}\\
\mathbf{elif}\;n \leq 2.7:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -156000 or 2.7000000000000002 < n Initial program 26.7%
Taylor expanded in n around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6484.8%
Simplified84.8%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f64N/A
/-lowering-/.f6484.6%
Applied egg-rr84.6%
if -156000 < n < -1.02e-199Initial program 15.9%
Taylor expanded in n around inf
expm1-defineN/A
expm1-lowering-expm1.f6474.1%
Simplified74.1%
if -1.02e-199 < n < 2.99999999999999994e-214Initial program 63.4%
pow-to-expN/A
expm1-defineN/A
expm1-lowering-expm1.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
/-lowering-/.f6484.5%
Applied egg-rr84.5%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr63.7%
Taylor expanded in i around 0
Simplified78.8%
if 2.99999999999999994e-214 < n < 2.7000000000000002Initial program 13.2%
Taylor expanded in i around 0
Simplified69.8%
Final simplification79.8%
(FPCore (i n)
:precision binary64
(let* ((t_0
(+
50.0
(*
(* i -100.0)
(+ (* i -0.041666666666666664) -0.16666666666666666)))))
(if (<= n -8e+95)
(* n (+ 100.0 (* i t_0)))
(if (<= n -4.2e-195)
(* 100.0 (/ (expm1 i) (/ i n)))
(if (<= n 1.06e-178)
(* n (/ (+ -100.0 100.0) i))
(+ (* n 100.0) (* t_0 (* i n))))))))
double code(double i, double n) {
double t_0 = 50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666));
double tmp;
if (n <= -8e+95) {
tmp = n * (100.0 + (i * t_0));
} else if (n <= -4.2e-195) {
tmp = 100.0 * (expm1(i) / (i / n));
} else if (n <= 1.06e-178) {
tmp = n * ((-100.0 + 100.0) / i);
} else {
tmp = (n * 100.0) + (t_0 * (i * n));
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = 50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666));
double tmp;
if (n <= -8e+95) {
tmp = n * (100.0 + (i * t_0));
} else if (n <= -4.2e-195) {
tmp = 100.0 * (Math.expm1(i) / (i / n));
} else if (n <= 1.06e-178) {
tmp = n * ((-100.0 + 100.0) / i);
} else {
tmp = (n * 100.0) + (t_0 * (i * n));
}
return tmp;
}
def code(i, n): t_0 = 50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666)) tmp = 0 if n <= -8e+95: tmp = n * (100.0 + (i * t_0)) elif n <= -4.2e-195: tmp = 100.0 * (math.expm1(i) / (i / n)) elif n <= 1.06e-178: tmp = n * ((-100.0 + 100.0) / i) else: tmp = (n * 100.0) + (t_0 * (i * n)) return tmp
function code(i, n) t_0 = Float64(50.0 + Float64(Float64(i * -100.0) * Float64(Float64(i * -0.041666666666666664) + -0.16666666666666666))) tmp = 0.0 if (n <= -8e+95) tmp = Float64(n * Float64(100.0 + Float64(i * t_0))); elseif (n <= -4.2e-195) tmp = Float64(100.0 * Float64(expm1(i) / Float64(i / n))); elseif (n <= 1.06e-178) tmp = Float64(n * Float64(Float64(-100.0 + 100.0) / i)); else tmp = Float64(Float64(n * 100.0) + Float64(t_0 * Float64(i * n))); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(50.0 + N[(N[(i * -100.0), $MachinePrecision] * N[(N[(i * -0.041666666666666664), $MachinePrecision] + -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -8e+95], N[(n * N[(100.0 + N[(i * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, -4.2e-195], N[(100.0 * N[(N[(Exp[i] - 1), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.06e-178], N[(n * N[(N[(-100.0 + 100.0), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], N[(N[(n * 100.0), $MachinePrecision] + N[(t$95$0 * N[(i * n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 50 + \left(i \cdot -100\right) \cdot \left(i \cdot -0.041666666666666664 + -0.16666666666666666\right)\\
\mathbf{if}\;n \leq -8 \cdot 10^{+95}:\\
\;\;\;\;n \cdot \left(100 + i \cdot t\_0\right)\\
\mathbf{elif}\;n \leq -4.2 \cdot 10^{-195}:\\
\;\;\;\;100 \cdot \frac{\mathsf{expm1}\left(i\right)}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 1.06 \cdot 10^{-178}:\\
\;\;\;\;n \cdot \frac{-100 + 100}{i}\\
\mathbf{else}:\\
\;\;\;\;n \cdot 100 + t\_0 \cdot \left(i \cdot n\right)\\
\end{array}
\end{array}
if n < -8.00000000000000016e95Initial program 26.2%
*-commutativeN/A
associate-*l/N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-lft-neg-outN/A
distribute-neg-fracN/A
distribute-neg-frac2N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
distribute-neg-frac2N/A
*-commutativeN/A
*-lowering-*.f64N/A
Simplified26.3%
Taylor expanded in i around 0
Simplified63.2%
Taylor expanded in n around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6463.2%
Simplified63.2%
if -8.00000000000000016e95 < n < -4.2e-195Initial program 26.6%
Taylor expanded in n around inf
expm1-defineN/A
expm1-lowering-expm1.f6473.1%
Simplified73.1%
if -4.2e-195 < n < 1.05999999999999999e-178Initial program 56.3%
pow-to-expN/A
expm1-defineN/A
expm1-lowering-expm1.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
/-lowering-/.f6480.3%
Applied egg-rr80.3%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr56.6%
Taylor expanded in i around 0
Simplified75.9%
if 1.05999999999999999e-178 < n Initial program 18.6%
*-commutativeN/A
associate-*l/N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-lft-neg-outN/A
distribute-neg-fracN/A
distribute-neg-frac2N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
distribute-neg-frac2N/A
*-commutativeN/A
*-lowering-*.f64N/A
Simplified18.5%
Taylor expanded in i around 0
Simplified68.3%
Taylor expanded in n around inf
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6475.9%
Simplified75.9%
Final simplification72.8%
(FPCore (i n)
:precision binary64
(if (<= n -5.2e-192)
(+
(* n 100.0)
(*
i
(*
n
(-
(-
(/
(+
(* (* i -100.0) (+ 0.5 (* i 0.25)))
(+
(/
(* 100.0 (* i (+ 0.3333333333333333 (* i 0.4583333333333333))))
n)
-50.0))
n)
(* (* i -100.0) (+ 0.16666666666666666 (* i 0.041666666666666664))))
-50.0))))
(if (<= n 8.2e-179)
(* n (/ (+ -100.0 100.0) i))
(+
(* n 100.0)
(*
(+
50.0
(* (* i -100.0) (+ (* i -0.041666666666666664) -0.16666666666666666)))
(* i n))))))
double code(double i, double n) {
double tmp;
if (n <= -5.2e-192) {
tmp = (n * 100.0) + (i * (n * ((((((i * -100.0) * (0.5 + (i * 0.25))) + (((100.0 * (i * (0.3333333333333333 + (i * 0.4583333333333333)))) / n) + -50.0)) / n) - ((i * -100.0) * (0.16666666666666666 + (i * 0.041666666666666664)))) - -50.0)));
} else if (n <= 8.2e-179) {
tmp = n * ((-100.0 + 100.0) / i);
} else {
tmp = (n * 100.0) + ((50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666))) * (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 <= (-5.2d-192)) then
tmp = (n * 100.0d0) + (i * (n * ((((((i * (-100.0d0)) * (0.5d0 + (i * 0.25d0))) + (((100.0d0 * (i * (0.3333333333333333d0 + (i * 0.4583333333333333d0)))) / n) + (-50.0d0))) / n) - ((i * (-100.0d0)) * (0.16666666666666666d0 + (i * 0.041666666666666664d0)))) - (-50.0d0))))
else if (n <= 8.2d-179) then
tmp = n * (((-100.0d0) + 100.0d0) / i)
else
tmp = (n * 100.0d0) + ((50.0d0 + ((i * (-100.0d0)) * ((i * (-0.041666666666666664d0)) + (-0.16666666666666666d0)))) * (i * n))
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (n <= -5.2e-192) {
tmp = (n * 100.0) + (i * (n * ((((((i * -100.0) * (0.5 + (i * 0.25))) + (((100.0 * (i * (0.3333333333333333 + (i * 0.4583333333333333)))) / n) + -50.0)) / n) - ((i * -100.0) * (0.16666666666666666 + (i * 0.041666666666666664)))) - -50.0)));
} else if (n <= 8.2e-179) {
tmp = n * ((-100.0 + 100.0) / i);
} else {
tmp = (n * 100.0) + ((50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666))) * (i * n));
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -5.2e-192: tmp = (n * 100.0) + (i * (n * ((((((i * -100.0) * (0.5 + (i * 0.25))) + (((100.0 * (i * (0.3333333333333333 + (i * 0.4583333333333333)))) / n) + -50.0)) / n) - ((i * -100.0) * (0.16666666666666666 + (i * 0.041666666666666664)))) - -50.0))) elif n <= 8.2e-179: tmp = n * ((-100.0 + 100.0) / i) else: tmp = (n * 100.0) + ((50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666))) * (i * n)) return tmp
function code(i, n) tmp = 0.0 if (n <= -5.2e-192) tmp = Float64(Float64(n * 100.0) + Float64(i * Float64(n * Float64(Float64(Float64(Float64(Float64(Float64(i * -100.0) * Float64(0.5 + Float64(i * 0.25))) + Float64(Float64(Float64(100.0 * Float64(i * Float64(0.3333333333333333 + Float64(i * 0.4583333333333333)))) / n) + -50.0)) / n) - Float64(Float64(i * -100.0) * Float64(0.16666666666666666 + Float64(i * 0.041666666666666664)))) - -50.0)))); elseif (n <= 8.2e-179) tmp = Float64(n * Float64(Float64(-100.0 + 100.0) / i)); else tmp = Float64(Float64(n * 100.0) + Float64(Float64(50.0 + Float64(Float64(i * -100.0) * Float64(Float64(i * -0.041666666666666664) + -0.16666666666666666))) * Float64(i * n))); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (n <= -5.2e-192) tmp = (n * 100.0) + (i * (n * ((((((i * -100.0) * (0.5 + (i * 0.25))) + (((100.0 * (i * (0.3333333333333333 + (i * 0.4583333333333333)))) / n) + -50.0)) / n) - ((i * -100.0) * (0.16666666666666666 + (i * 0.041666666666666664)))) - -50.0))); elseif (n <= 8.2e-179) tmp = n * ((-100.0 + 100.0) / i); else tmp = (n * 100.0) + ((50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666))) * (i * n)); end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[n, -5.2e-192], N[(N[(n * 100.0), $MachinePrecision] + N[(i * N[(n * N[(N[(N[(N[(N[(N[(i * -100.0), $MachinePrecision] * N[(0.5 + N[(i * 0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(100.0 * N[(i * N[(0.3333333333333333 + N[(i * 0.4583333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision] + -50.0), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision] - N[(N[(i * -100.0), $MachinePrecision] * N[(0.16666666666666666 + N[(i * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - -50.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 8.2e-179], N[(n * N[(N[(-100.0 + 100.0), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], N[(N[(n * 100.0), $MachinePrecision] + N[(N[(50.0 + N[(N[(i * -100.0), $MachinePrecision] * N[(N[(i * -0.041666666666666664), $MachinePrecision] + -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(i * n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -5.2 \cdot 10^{-192}:\\
\;\;\;\;n \cdot 100 + i \cdot \left(n \cdot \left(\left(\frac{\left(i \cdot -100\right) \cdot \left(0.5 + i \cdot 0.25\right) + \left(\frac{100 \cdot \left(i \cdot \left(0.3333333333333333 + i \cdot 0.4583333333333333\right)\right)}{n} + -50\right)}{n} - \left(i \cdot -100\right) \cdot \left(0.16666666666666666 + i \cdot 0.041666666666666664\right)\right) - -50\right)\right)\\
\mathbf{elif}\;n \leq 8.2 \cdot 10^{-179}:\\
\;\;\;\;n \cdot \frac{-100 + 100}{i}\\
\mathbf{else}:\\
\;\;\;\;n \cdot 100 + \left(50 + \left(i \cdot -100\right) \cdot \left(i \cdot -0.041666666666666664 + -0.16666666666666666\right)\right) \cdot \left(i \cdot n\right)\\
\end{array}
\end{array}
if n < -5.2000000000000003e-192Initial program 26.4%
*-commutativeN/A
associate-*l/N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-lft-neg-outN/A
distribute-neg-fracN/A
distribute-neg-frac2N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
distribute-neg-frac2N/A
*-commutativeN/A
*-lowering-*.f64N/A
Simplified26.4%
Taylor expanded in i around 0
Simplified53.4%
Taylor expanded in n around -inf
Simplified62.1%
if -5.2000000000000003e-192 < n < 8.2e-179Initial program 56.3%
pow-to-expN/A
expm1-defineN/A
expm1-lowering-expm1.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
/-lowering-/.f6480.3%
Applied egg-rr80.3%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr56.6%
Taylor expanded in i around 0
Simplified75.9%
if 8.2e-179 < n Initial program 18.6%
*-commutativeN/A
associate-*l/N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-lft-neg-outN/A
distribute-neg-fracN/A
distribute-neg-frac2N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
distribute-neg-frac2N/A
*-commutativeN/A
*-lowering-*.f64N/A
Simplified18.5%
Taylor expanded in i around 0
Simplified68.3%
Taylor expanded in n around inf
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6475.9%
Simplified75.9%
Final simplification70.2%
(FPCore (i n)
:precision binary64
(if (<= n -1.65e-199)
(+
(* n 100.0)
(*
i
(+
(* (+ (/ 0.5 n) -0.5) (* n -100.0))
(*
(* i 100.0)
(* n (+ 0.16666666666666666 (* i 0.041666666666666664)))))))
(if (<= n 1.15e-178)
(* n (/ (+ -100.0 100.0) i))
(+
(* n 100.0)
(*
(+
50.0
(* (* i -100.0) (+ (* i -0.041666666666666664) -0.16666666666666666)))
(* i n))))))
double code(double i, double n) {
double tmp;
if (n <= -1.65e-199) {
tmp = (n * 100.0) + (i * ((((0.5 / n) + -0.5) * (n * -100.0)) + ((i * 100.0) * (n * (0.16666666666666666 + (i * 0.041666666666666664))))));
} else if (n <= 1.15e-178) {
tmp = n * ((-100.0 + 100.0) / i);
} else {
tmp = (n * 100.0) + ((50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666))) * (i * n));
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if (n <= (-1.65d-199)) then
tmp = (n * 100.0d0) + (i * ((((0.5d0 / n) + (-0.5d0)) * (n * (-100.0d0))) + ((i * 100.0d0) * (n * (0.16666666666666666d0 + (i * 0.041666666666666664d0))))))
else if (n <= 1.15d-178) then
tmp = n * (((-100.0d0) + 100.0d0) / i)
else
tmp = (n * 100.0d0) + ((50.0d0 + ((i * (-100.0d0)) * ((i * (-0.041666666666666664d0)) + (-0.16666666666666666d0)))) * (i * n))
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (n <= -1.65e-199) {
tmp = (n * 100.0) + (i * ((((0.5 / n) + -0.5) * (n * -100.0)) + ((i * 100.0) * (n * (0.16666666666666666 + (i * 0.041666666666666664))))));
} else if (n <= 1.15e-178) {
tmp = n * ((-100.0 + 100.0) / i);
} else {
tmp = (n * 100.0) + ((50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666))) * (i * n));
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -1.65e-199: tmp = (n * 100.0) + (i * ((((0.5 / n) + -0.5) * (n * -100.0)) + ((i * 100.0) * (n * (0.16666666666666666 + (i * 0.041666666666666664)))))) elif n <= 1.15e-178: tmp = n * ((-100.0 + 100.0) / i) else: tmp = (n * 100.0) + ((50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666))) * (i * n)) return tmp
function code(i, n) tmp = 0.0 if (n <= -1.65e-199) tmp = Float64(Float64(n * 100.0) + Float64(i * Float64(Float64(Float64(Float64(0.5 / n) + -0.5) * Float64(n * -100.0)) + Float64(Float64(i * 100.0) * Float64(n * Float64(0.16666666666666666 + Float64(i * 0.041666666666666664))))))); elseif (n <= 1.15e-178) tmp = Float64(n * Float64(Float64(-100.0 + 100.0) / i)); else tmp = Float64(Float64(n * 100.0) + Float64(Float64(50.0 + Float64(Float64(i * -100.0) * Float64(Float64(i * -0.041666666666666664) + -0.16666666666666666))) * Float64(i * n))); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (n <= -1.65e-199) tmp = (n * 100.0) + (i * ((((0.5 / n) + -0.5) * (n * -100.0)) + ((i * 100.0) * (n * (0.16666666666666666 + (i * 0.041666666666666664)))))); elseif (n <= 1.15e-178) tmp = n * ((-100.0 + 100.0) / i); else tmp = (n * 100.0) + ((50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666))) * (i * n)); end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[n, -1.65e-199], N[(N[(n * 100.0), $MachinePrecision] + N[(i * N[(N[(N[(N[(0.5 / n), $MachinePrecision] + -0.5), $MachinePrecision] * N[(n * -100.0), $MachinePrecision]), $MachinePrecision] + N[(N[(i * 100.0), $MachinePrecision] * N[(n * N[(0.16666666666666666 + N[(i * 0.041666666666666664), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.15e-178], N[(n * N[(N[(-100.0 + 100.0), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], N[(N[(n * 100.0), $MachinePrecision] + N[(N[(50.0 + N[(N[(i * -100.0), $MachinePrecision] * N[(N[(i * -0.041666666666666664), $MachinePrecision] + -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(i * n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.65 \cdot 10^{-199}:\\
\;\;\;\;n \cdot 100 + i \cdot \left(\left(\frac{0.5}{n} + -0.5\right) \cdot \left(n \cdot -100\right) + \left(i \cdot 100\right) \cdot \left(n \cdot \left(0.16666666666666666 + i \cdot 0.041666666666666664\right)\right)\right)\\
\mathbf{elif}\;n \leq 1.15 \cdot 10^{-178}:\\
\;\;\;\;n \cdot \frac{-100 + 100}{i}\\
\mathbf{else}:\\
\;\;\;\;n \cdot 100 + \left(50 + \left(i \cdot -100\right) \cdot \left(i \cdot -0.041666666666666664 + -0.16666666666666666\right)\right) \cdot \left(i \cdot n\right)\\
\end{array}
\end{array}
if n < -1.6500000000000001e-199Initial program 26.4%
*-commutativeN/A
associate-*l/N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-lft-neg-outN/A
distribute-neg-fracN/A
distribute-neg-frac2N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
distribute-neg-frac2N/A
*-commutativeN/A
*-lowering-*.f64N/A
Simplified26.4%
Taylor expanded in i around 0
Simplified53.4%
Taylor expanded in n around -inf
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6461.7%
Simplified61.7%
if -1.6500000000000001e-199 < n < 1.14999999999999997e-178Initial program 56.3%
pow-to-expN/A
expm1-defineN/A
expm1-lowering-expm1.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
/-lowering-/.f6480.3%
Applied egg-rr80.3%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr56.6%
Taylor expanded in i around 0
Simplified75.9%
if 1.14999999999999997e-178 < n Initial program 18.6%
*-commutativeN/A
associate-*l/N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-lft-neg-outN/A
distribute-neg-fracN/A
distribute-neg-frac2N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
distribute-neg-frac2N/A
*-commutativeN/A
*-lowering-*.f64N/A
Simplified18.5%
Taylor expanded in i around 0
Simplified68.3%
Taylor expanded in n around inf
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6475.9%
Simplified75.9%
Final simplification70.0%
(FPCore (i n)
:precision binary64
(if (<= n -2.7e-193)
(+
(* n 100.0)
(*
i
(+
(* n 50.0)
(* i (* n (+ (* i 4.166666666666667) 16.666666666666668))))))
(if (<= n 9.4e-180)
(* n (/ (+ -100.0 100.0) i))
(+
(* n 100.0)
(*
(+
50.0
(* (* i -100.0) (+ (* i -0.041666666666666664) -0.16666666666666666)))
(* i n))))))
double code(double i, double n) {
double tmp;
if (n <= -2.7e-193) {
tmp = (n * 100.0) + (i * ((n * 50.0) + (i * (n * ((i * 4.166666666666667) + 16.666666666666668)))));
} else if (n <= 9.4e-180) {
tmp = n * ((-100.0 + 100.0) / i);
} else {
tmp = (n * 100.0) + ((50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666))) * (i * n));
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if (n <= (-2.7d-193)) then
tmp = (n * 100.0d0) + (i * ((n * 50.0d0) + (i * (n * ((i * 4.166666666666667d0) + 16.666666666666668d0)))))
else if (n <= 9.4d-180) then
tmp = n * (((-100.0d0) + 100.0d0) / i)
else
tmp = (n * 100.0d0) + ((50.0d0 + ((i * (-100.0d0)) * ((i * (-0.041666666666666664d0)) + (-0.16666666666666666d0)))) * (i * n))
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (n <= -2.7e-193) {
tmp = (n * 100.0) + (i * ((n * 50.0) + (i * (n * ((i * 4.166666666666667) + 16.666666666666668)))));
} else if (n <= 9.4e-180) {
tmp = n * ((-100.0 + 100.0) / i);
} else {
tmp = (n * 100.0) + ((50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666))) * (i * n));
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -2.7e-193: tmp = (n * 100.0) + (i * ((n * 50.0) + (i * (n * ((i * 4.166666666666667) + 16.666666666666668))))) elif n <= 9.4e-180: tmp = n * ((-100.0 + 100.0) / i) else: tmp = (n * 100.0) + ((50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666))) * (i * n)) return tmp
function code(i, n) tmp = 0.0 if (n <= -2.7e-193) tmp = Float64(Float64(n * 100.0) + Float64(i * Float64(Float64(n * 50.0) + Float64(i * Float64(n * Float64(Float64(i * 4.166666666666667) + 16.666666666666668)))))); elseif (n <= 9.4e-180) tmp = Float64(n * Float64(Float64(-100.0 + 100.0) / i)); else tmp = Float64(Float64(n * 100.0) + Float64(Float64(50.0 + Float64(Float64(i * -100.0) * Float64(Float64(i * -0.041666666666666664) + -0.16666666666666666))) * Float64(i * n))); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (n <= -2.7e-193) tmp = (n * 100.0) + (i * ((n * 50.0) + (i * (n * ((i * 4.166666666666667) + 16.666666666666668))))); elseif (n <= 9.4e-180) tmp = n * ((-100.0 + 100.0) / i); else tmp = (n * 100.0) + ((50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666))) * (i * n)); end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[n, -2.7e-193], N[(N[(n * 100.0), $MachinePrecision] + N[(i * N[(N[(n * 50.0), $MachinePrecision] + N[(i * N[(n * N[(N[(i * 4.166666666666667), $MachinePrecision] + 16.666666666666668), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 9.4e-180], N[(n * N[(N[(-100.0 + 100.0), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], N[(N[(n * 100.0), $MachinePrecision] + N[(N[(50.0 + N[(N[(i * -100.0), $MachinePrecision] * N[(N[(i * -0.041666666666666664), $MachinePrecision] + -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(i * n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -2.7 \cdot 10^{-193}:\\
\;\;\;\;n \cdot 100 + i \cdot \left(n \cdot 50 + i \cdot \left(n \cdot \left(i \cdot 4.166666666666667 + 16.666666666666668\right)\right)\right)\\
\mathbf{elif}\;n \leq 9.4 \cdot 10^{-180}:\\
\;\;\;\;n \cdot \frac{-100 + 100}{i}\\
\mathbf{else}:\\
\;\;\;\;n \cdot 100 + \left(50 + \left(i \cdot -100\right) \cdot \left(i \cdot -0.041666666666666664 + -0.16666666666666666\right)\right) \cdot \left(i \cdot n\right)\\
\end{array}
\end{array}
if n < -2.6999999999999999e-193Initial program 26.4%
Taylor expanded in n around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6471.0%
Simplified71.0%
Taylor expanded in i around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6461.7%
Simplified61.7%
if -2.6999999999999999e-193 < n < 9.39999999999999951e-180Initial program 56.3%
pow-to-expN/A
expm1-defineN/A
expm1-lowering-expm1.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
/-lowering-/.f6480.3%
Applied egg-rr80.3%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr56.6%
Taylor expanded in i around 0
Simplified75.9%
if 9.39999999999999951e-180 < n Initial program 18.6%
*-commutativeN/A
associate-*l/N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-lft-neg-outN/A
distribute-neg-fracN/A
distribute-neg-frac2N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
distribute-neg-frac2N/A
*-commutativeN/A
*-lowering-*.f64N/A
Simplified18.5%
Taylor expanded in i around 0
Simplified68.3%
Taylor expanded in n around inf
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6475.9%
Simplified75.9%
Final simplification70.0%
(FPCore (i n)
:precision binary64
(if (<= n -2.1e-196)
(+
(* n 100.0)
(*
i
(+
(* n 50.0)
(* i (* n (+ (* i 4.166666666666667) 16.666666666666668))))))
(if (<= n 1.12e-179)
(* n (/ (+ -100.0 100.0) i))
(*
n
(+
100.0
(*
i
(+
50.0
(*
(* i -100.0)
(+ (* i -0.041666666666666664) -0.16666666666666666)))))))))
double code(double i, double n) {
double tmp;
if (n <= -2.1e-196) {
tmp = (n * 100.0) + (i * ((n * 50.0) + (i * (n * ((i * 4.166666666666667) + 16.666666666666668)))));
} else if (n <= 1.12e-179) {
tmp = n * ((-100.0 + 100.0) / i);
} else {
tmp = n * (100.0 + (i * (50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666)))));
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if (n <= (-2.1d-196)) then
tmp = (n * 100.0d0) + (i * ((n * 50.0d0) + (i * (n * ((i * 4.166666666666667d0) + 16.666666666666668d0)))))
else if (n <= 1.12d-179) then
tmp = n * (((-100.0d0) + 100.0d0) / i)
else
tmp = n * (100.0d0 + (i * (50.0d0 + ((i * (-100.0d0)) * ((i * (-0.041666666666666664d0)) + (-0.16666666666666666d0))))))
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (n <= -2.1e-196) {
tmp = (n * 100.0) + (i * ((n * 50.0) + (i * (n * ((i * 4.166666666666667) + 16.666666666666668)))));
} else if (n <= 1.12e-179) {
tmp = n * ((-100.0 + 100.0) / i);
} else {
tmp = n * (100.0 + (i * (50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666)))));
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -2.1e-196: tmp = (n * 100.0) + (i * ((n * 50.0) + (i * (n * ((i * 4.166666666666667) + 16.666666666666668))))) elif n <= 1.12e-179: tmp = n * ((-100.0 + 100.0) / i) else: tmp = n * (100.0 + (i * (50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666))))) return tmp
function code(i, n) tmp = 0.0 if (n <= -2.1e-196) tmp = Float64(Float64(n * 100.0) + Float64(i * Float64(Float64(n * 50.0) + Float64(i * Float64(n * Float64(Float64(i * 4.166666666666667) + 16.666666666666668)))))); elseif (n <= 1.12e-179) tmp = Float64(n * Float64(Float64(-100.0 + 100.0) / i)); else tmp = Float64(n * Float64(100.0 + Float64(i * Float64(50.0 + Float64(Float64(i * -100.0) * Float64(Float64(i * -0.041666666666666664) + -0.16666666666666666)))))); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (n <= -2.1e-196) tmp = (n * 100.0) + (i * ((n * 50.0) + (i * (n * ((i * 4.166666666666667) + 16.666666666666668))))); elseif (n <= 1.12e-179) tmp = n * ((-100.0 + 100.0) / i); else tmp = n * (100.0 + (i * (50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666))))); end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[n, -2.1e-196], N[(N[(n * 100.0), $MachinePrecision] + N[(i * N[(N[(n * 50.0), $MachinePrecision] + N[(i * N[(n * N[(N[(i * 4.166666666666667), $MachinePrecision] + 16.666666666666668), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.12e-179], N[(n * N[(N[(-100.0 + 100.0), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], N[(n * N[(100.0 + N[(i * N[(50.0 + N[(N[(i * -100.0), $MachinePrecision] * N[(N[(i * -0.041666666666666664), $MachinePrecision] + -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -2.1 \cdot 10^{-196}:\\
\;\;\;\;n \cdot 100 + i \cdot \left(n \cdot 50 + i \cdot \left(n \cdot \left(i \cdot 4.166666666666667 + 16.666666666666668\right)\right)\right)\\
\mathbf{elif}\;n \leq 1.12 \cdot 10^{-179}:\\
\;\;\;\;n \cdot \frac{-100 + 100}{i}\\
\mathbf{else}:\\
\;\;\;\;n \cdot \left(100 + i \cdot \left(50 + \left(i \cdot -100\right) \cdot \left(i \cdot -0.041666666666666664 + -0.16666666666666666\right)\right)\right)\\
\end{array}
\end{array}
if n < -2.09999999999999988e-196Initial program 26.4%
Taylor expanded in n around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6471.0%
Simplified71.0%
Taylor expanded in i around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6461.7%
Simplified61.7%
if -2.09999999999999988e-196 < n < 1.11999999999999999e-179Initial program 56.3%
pow-to-expN/A
expm1-defineN/A
expm1-lowering-expm1.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
/-lowering-/.f6480.3%
Applied egg-rr80.3%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr56.6%
Taylor expanded in i around 0
Simplified75.9%
if 1.11999999999999999e-179 < n Initial program 18.6%
*-commutativeN/A
associate-*l/N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-lft-neg-outN/A
distribute-neg-fracN/A
distribute-neg-frac2N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
distribute-neg-frac2N/A
*-commutativeN/A
*-lowering-*.f64N/A
Simplified18.5%
Taylor expanded in i around 0
Simplified68.3%
Taylor expanded in n around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6475.9%
Simplified75.9%
Final simplification70.0%
(FPCore (i n)
:precision binary64
(let* ((t_0
(*
n
(+
100.0
(*
i
(+
50.0
(*
(* i -100.0)
(+ (* i -0.041666666666666664) -0.16666666666666666))))))))
(if (<= n -5.2e-197)
t_0
(if (<= n 2.75e-179) (* n (/ (+ -100.0 100.0) i)) t_0))))
double code(double i, double n) {
double t_0 = n * (100.0 + (i * (50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666)))));
double tmp;
if (n <= -5.2e-197) {
tmp = t_0;
} else if (n <= 2.75e-179) {
tmp = n * ((-100.0 + 100.0) / i);
} 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 * (-100.0d0)) * ((i * (-0.041666666666666664d0)) + (-0.16666666666666666d0))))))
if (n <= (-5.2d-197)) then
tmp = t_0
else if (n <= 2.75d-179) then
tmp = n * (((-100.0d0) + 100.0d0) / i)
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 * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666)))));
double tmp;
if (n <= -5.2e-197) {
tmp = t_0;
} else if (n <= 2.75e-179) {
tmp = n * ((-100.0 + 100.0) / i);
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = n * (100.0 + (i * (50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666))))) tmp = 0 if n <= -5.2e-197: tmp = t_0 elif n <= 2.75e-179: tmp = n * ((-100.0 + 100.0) / i) else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(n * Float64(100.0 + Float64(i * Float64(50.0 + Float64(Float64(i * -100.0) * Float64(Float64(i * -0.041666666666666664) + -0.16666666666666666)))))) tmp = 0.0 if (n <= -5.2e-197) tmp = t_0; elseif (n <= 2.75e-179) tmp = Float64(n * Float64(Float64(-100.0 + 100.0) / i)); else tmp = t_0; end return tmp end
function tmp_2 = code(i, n) t_0 = n * (100.0 + (i * (50.0 + ((i * -100.0) * ((i * -0.041666666666666664) + -0.16666666666666666))))); tmp = 0.0; if (n <= -5.2e-197) tmp = t_0; elseif (n <= 2.75e-179) tmp = n * ((-100.0 + 100.0) / i); 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[(N[(i * -100.0), $MachinePrecision] * N[(N[(i * -0.041666666666666664), $MachinePrecision] + -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -5.2e-197], t$95$0, If[LessEqual[n, 2.75e-179], N[(n * N[(N[(-100.0 + 100.0), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := n \cdot \left(100 + i \cdot \left(50 + \left(i \cdot -100\right) \cdot \left(i \cdot -0.041666666666666664 + -0.16666666666666666\right)\right)\right)\\
\mathbf{if}\;n \leq -5.2 \cdot 10^{-197}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 2.75 \cdot 10^{-179}:\\
\;\;\;\;n \cdot \frac{-100 + 100}{i}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -5.2000000000000003e-197 or 2.7500000000000001e-179 < n Initial program 22.5%
*-commutativeN/A
associate-*l/N/A
sub-negN/A
remove-double-negN/A
distribute-neg-inN/A
+-commutativeN/A
sub-negN/A
distribute-lft-neg-outN/A
distribute-neg-fracN/A
distribute-neg-frac2N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
distribute-neg-frac2N/A
*-commutativeN/A
*-lowering-*.f64N/A
Simplified22.5%
Taylor expanded in i around 0
Simplified60.8%
Taylor expanded in n around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6468.7%
Simplified68.7%
if -5.2000000000000003e-197 < n < 2.7500000000000001e-179Initial program 56.3%
pow-to-expN/A
expm1-defineN/A
expm1-lowering-expm1.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
/-lowering-/.f6480.3%
Applied egg-rr80.3%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr56.6%
Taylor expanded in i around 0
Simplified75.9%
Final simplification70.0%
(FPCore (i n)
:precision binary64
(if (<= n -4.3e-198)
(+ (* n 100.0) (* i (* n (+ 50.0 (* i 16.666666666666668)))))
(if (<= n 1.4e-219)
(* n (/ (+ -100.0 100.0) i))
(if (<= n 1.55)
(* 100.0 (/ i (/ i n)))
(/ (* i (* n (+ 100.0 (* i 50.0)))) i)))))
double code(double i, double n) {
double tmp;
if (n <= -4.3e-198) {
tmp = (n * 100.0) + (i * (n * (50.0 + (i * 16.666666666666668))));
} else if (n <= 1.4e-219) {
tmp = n * ((-100.0 + 100.0) / i);
} else if (n <= 1.55) {
tmp = 100.0 * (i / (i / n));
} 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) :: tmp
if (n <= (-4.3d-198)) then
tmp = (n * 100.0d0) + (i * (n * (50.0d0 + (i * 16.666666666666668d0))))
else if (n <= 1.4d-219) then
tmp = n * (((-100.0d0) + 100.0d0) / i)
else if (n <= 1.55d0) then
tmp = 100.0d0 * (i / (i / n))
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 tmp;
if (n <= -4.3e-198) {
tmp = (n * 100.0) + (i * (n * (50.0 + (i * 16.666666666666668))));
} else if (n <= 1.4e-219) {
tmp = n * ((-100.0 + 100.0) / i);
} else if (n <= 1.55) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = (i * (n * (100.0 + (i * 50.0)))) / i;
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -4.3e-198: tmp = (n * 100.0) + (i * (n * (50.0 + (i * 16.666666666666668)))) elif n <= 1.4e-219: tmp = n * ((-100.0 + 100.0) / i) elif n <= 1.55: tmp = 100.0 * (i / (i / n)) else: tmp = (i * (n * (100.0 + (i * 50.0)))) / i return tmp
function code(i, n) tmp = 0.0 if (n <= -4.3e-198) tmp = Float64(Float64(n * 100.0) + Float64(i * Float64(n * Float64(50.0 + Float64(i * 16.666666666666668))))); elseif (n <= 1.4e-219) tmp = Float64(n * Float64(Float64(-100.0 + 100.0) / i)); elseif (n <= 1.55) tmp = Float64(100.0 * Float64(i / Float64(i / n))); 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) tmp = 0.0; if (n <= -4.3e-198) tmp = (n * 100.0) + (i * (n * (50.0 + (i * 16.666666666666668)))); elseif (n <= 1.4e-219) tmp = n * ((-100.0 + 100.0) / i); elseif (n <= 1.55) tmp = 100.0 * (i / (i / n)); else tmp = (i * (n * (100.0 + (i * 50.0)))) / i; end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[n, -4.3e-198], N[(N[(n * 100.0), $MachinePrecision] + N[(i * N[(n * N[(50.0 + N[(i * 16.666666666666668), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.4e-219], N[(n * N[(N[(-100.0 + 100.0), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.55], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(i * N[(n * N[(100.0 + N[(i * 50.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -4.3 \cdot 10^{-198}:\\
\;\;\;\;n \cdot 100 + i \cdot \left(n \cdot \left(50 + i \cdot 16.666666666666668\right)\right)\\
\mathbf{elif}\;n \leq 1.4 \cdot 10^{-219}:\\
\;\;\;\;n \cdot \frac{-100 + 100}{i}\\
\mathbf{elif}\;n \leq 1.55:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;\frac{i \cdot \left(n \cdot \left(100 + i \cdot 50\right)\right)}{i}\\
\end{array}
\end{array}
if n < -4.3000000000000003e-198Initial program 26.4%
Taylor expanded in n around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6471.0%
Simplified71.0%
Taylor expanded in i around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6461.3%
Simplified61.3%
if -4.3000000000000003e-198 < n < 1.39999999999999995e-219Initial program 63.4%
pow-to-expN/A
expm1-defineN/A
expm1-lowering-expm1.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
/-lowering-/.f6484.5%
Applied egg-rr84.5%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr63.7%
Taylor expanded in i around 0
Simplified78.8%
if 1.39999999999999995e-219 < n < 1.55000000000000004Initial program 13.2%
Taylor expanded in i around 0
Simplified69.8%
if 1.55000000000000004 < n Initial program 21.8%
Taylor expanded in n around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6488.0%
Simplified88.0%
Taylor expanded in i around 0
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6479.9%
Simplified79.9%
Final simplification70.4%
(FPCore (i n)
:precision binary64
(if (<= n -7.2e-194)
(* 100.0 (* n (+ 1.0 (* 0.16666666666666666 (* i i)))))
(if (<= n 5.6e-219)
(* n (/ (+ -100.0 100.0) i))
(if (<= n 1.5)
(* 100.0 (/ i (/ i n)))
(/ (* i (* n (+ 100.0 (* i 50.0)))) i)))))
double code(double i, double n) {
double tmp;
if (n <= -7.2e-194) {
tmp = 100.0 * (n * (1.0 + (0.16666666666666666 * (i * i))));
} else if (n <= 5.6e-219) {
tmp = n * ((-100.0 + 100.0) / i);
} else if (n <= 1.5) {
tmp = 100.0 * (i / (i / n));
} 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) :: tmp
if (n <= (-7.2d-194)) then
tmp = 100.0d0 * (n * (1.0d0 + (0.16666666666666666d0 * (i * i))))
else if (n <= 5.6d-219) then
tmp = n * (((-100.0d0) + 100.0d0) / i)
else if (n <= 1.5d0) then
tmp = 100.0d0 * (i / (i / n))
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 tmp;
if (n <= -7.2e-194) {
tmp = 100.0 * (n * (1.0 + (0.16666666666666666 * (i * i))));
} else if (n <= 5.6e-219) {
tmp = n * ((-100.0 + 100.0) / i);
} else if (n <= 1.5) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = (i * (n * (100.0 + (i * 50.0)))) / i;
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -7.2e-194: tmp = 100.0 * (n * (1.0 + (0.16666666666666666 * (i * i)))) elif n <= 5.6e-219: tmp = n * ((-100.0 + 100.0) / i) elif n <= 1.5: tmp = 100.0 * (i / (i / n)) else: tmp = (i * (n * (100.0 + (i * 50.0)))) / i return tmp
function code(i, n) tmp = 0.0 if (n <= -7.2e-194) tmp = Float64(100.0 * Float64(n * Float64(1.0 + Float64(0.16666666666666666 * Float64(i * i))))); elseif (n <= 5.6e-219) tmp = Float64(n * Float64(Float64(-100.0 + 100.0) / i)); elseif (n <= 1.5) tmp = Float64(100.0 * Float64(i / Float64(i / n))); 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) tmp = 0.0; if (n <= -7.2e-194) tmp = 100.0 * (n * (1.0 + (0.16666666666666666 * (i * i)))); elseif (n <= 5.6e-219) tmp = n * ((-100.0 + 100.0) / i); elseif (n <= 1.5) tmp = 100.0 * (i / (i / n)); else tmp = (i * (n * (100.0 + (i * 50.0)))) / i; end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[n, -7.2e-194], N[(100.0 * N[(n * N[(1.0 + N[(0.16666666666666666 * N[(i * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 5.6e-219], N[(n * N[(N[(-100.0 + 100.0), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.5], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(i * N[(n * N[(100.0 + N[(i * 50.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -7.2 \cdot 10^{-194}:\\
\;\;\;\;100 \cdot \left(n \cdot \left(1 + 0.16666666666666666 \cdot \left(i \cdot i\right)\right)\right)\\
\mathbf{elif}\;n \leq 5.6 \cdot 10^{-219}:\\
\;\;\;\;n \cdot \frac{-100 + 100}{i}\\
\mathbf{elif}\;n \leq 1.5:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;\frac{i \cdot \left(n \cdot \left(100 + i \cdot 50\right)\right)}{i}\\
\end{array}
\end{array}
if n < -7.2e-194Initial program 26.4%
Taylor expanded in i around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
Simplified57.9%
Taylor expanded in n around 0
Simplified57.9%
Taylor expanded in n around inf
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6461.3%
Simplified61.3%
if -7.2e-194 < n < 5.5999999999999998e-219Initial program 63.4%
pow-to-expN/A
expm1-defineN/A
expm1-lowering-expm1.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
/-lowering-/.f6484.5%
Applied egg-rr84.5%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr63.7%
Taylor expanded in i around 0
Simplified78.8%
if 5.5999999999999998e-219 < n < 1.5Initial program 13.2%
Taylor expanded in i around 0
Simplified69.8%
if 1.5 < n Initial program 21.8%
Taylor expanded in n around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6488.0%
Simplified88.0%
Taylor expanded in i around 0
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6479.9%
Simplified79.9%
Final simplification70.4%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* 100.0 (* n (+ 1.0 (* 0.16666666666666666 (* i i)))))))
(if (<= n -1.26e-191)
t_0
(if (<= n 2.9e-179) (* n (/ (+ -100.0 100.0) i)) t_0))))
double code(double i, double n) {
double t_0 = 100.0 * (n * (1.0 + (0.16666666666666666 * (i * i))));
double tmp;
if (n <= -1.26e-191) {
tmp = t_0;
} else if (n <= 2.9e-179) {
tmp = n * ((-100.0 + 100.0) / i);
} 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 + (0.16666666666666666d0 * (i * i))))
if (n <= (-1.26d-191)) then
tmp = t_0
else if (n <= 2.9d-179) then
tmp = n * (((-100.0d0) + 100.0d0) / i)
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 + (0.16666666666666666 * (i * i))));
double tmp;
if (n <= -1.26e-191) {
tmp = t_0;
} else if (n <= 2.9e-179) {
tmp = n * ((-100.0 + 100.0) / i);
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = 100.0 * (n * (1.0 + (0.16666666666666666 * (i * i)))) tmp = 0 if n <= -1.26e-191: tmp = t_0 elif n <= 2.9e-179: tmp = n * ((-100.0 + 100.0) / i) else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(100.0 * Float64(n * Float64(1.0 + Float64(0.16666666666666666 * Float64(i * i))))) tmp = 0.0 if (n <= -1.26e-191) tmp = t_0; elseif (n <= 2.9e-179) tmp = Float64(n * Float64(Float64(-100.0 + 100.0) / i)); else tmp = t_0; end return tmp end
function tmp_2 = code(i, n) t_0 = 100.0 * (n * (1.0 + (0.16666666666666666 * (i * i)))); tmp = 0.0; if (n <= -1.26e-191) tmp = t_0; elseif (n <= 2.9e-179) tmp = n * ((-100.0 + 100.0) / i); 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[(0.16666666666666666 * N[(i * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -1.26e-191], t$95$0, If[LessEqual[n, 2.9e-179], N[(n * N[(N[(-100.0 + 100.0), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 100 \cdot \left(n \cdot \left(1 + 0.16666666666666666 \cdot \left(i \cdot i\right)\right)\right)\\
\mathbf{if}\;n \leq -1.26 \cdot 10^{-191}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 2.9 \cdot 10^{-179}:\\
\;\;\;\;n \cdot \frac{-100 + 100}{i}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -1.26e-191 or 2.8999999999999999e-179 < n Initial program 22.5%
Taylor expanded in i around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
Simplified61.4%
Taylor expanded in n around 0
Simplified61.4%
Taylor expanded in n around inf
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6464.9%
Simplified64.9%
if -1.26e-191 < n < 2.8999999999999999e-179Initial program 56.3%
pow-to-expN/A
expm1-defineN/A
expm1-lowering-expm1.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
/-lowering-/.f6480.3%
Applied egg-rr80.3%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr56.6%
Taylor expanded in i around 0
Simplified75.9%
Final simplification66.9%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* n (+ 100.0 (* i 50.0)))))
(if (<= n -1.2e-199)
t_0
(if (<= n 8e-180) (* n (/ (+ -100.0 100.0) i)) t_0))))
double code(double i, double n) {
double t_0 = n * (100.0 + (i * 50.0));
double tmp;
if (n <= -1.2e-199) {
tmp = t_0;
} else if (n <= 8e-180) {
tmp = n * ((-100.0 + 100.0) / i);
} 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))
if (n <= (-1.2d-199)) then
tmp = t_0
else if (n <= 8d-180) then
tmp = n * (((-100.0d0) + 100.0d0) / i)
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 tmp;
if (n <= -1.2e-199) {
tmp = t_0;
} else if (n <= 8e-180) {
tmp = n * ((-100.0 + 100.0) / i);
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = n * (100.0 + (i * 50.0)) tmp = 0 if n <= -1.2e-199: tmp = t_0 elif n <= 8e-180: tmp = n * ((-100.0 + 100.0) / i) else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(n * Float64(100.0 + Float64(i * 50.0))) tmp = 0.0 if (n <= -1.2e-199) tmp = t_0; elseif (n <= 8e-180) tmp = Float64(n * Float64(Float64(-100.0 + 100.0) / i)); else tmp = t_0; end return tmp end
function tmp_2 = code(i, n) t_0 = n * (100.0 + (i * 50.0)); tmp = 0.0; if (n <= -1.2e-199) tmp = t_0; elseif (n <= 8e-180) tmp = n * ((-100.0 + 100.0) / i); 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]}, If[LessEqual[n, -1.2e-199], t$95$0, If[LessEqual[n, 8e-180], N[(n * N[(N[(-100.0 + 100.0), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := n \cdot \left(100 + i \cdot 50\right)\\
\mathbf{if}\;n \leq -1.2 \cdot 10^{-199}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 8 \cdot 10^{-180}:\\
\;\;\;\;n \cdot \frac{-100 + 100}{i}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -1.19999999999999998e-199 or 8.0000000000000002e-180 < n Initial program 22.5%
Taylor expanded in n around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6471.0%
Simplified71.0%
Taylor expanded in i around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6462.2%
Simplified62.2%
if -1.19999999999999998e-199 < n < 8.0000000000000002e-180Initial program 56.3%
pow-to-expN/A
expm1-defineN/A
expm1-lowering-expm1.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
log1p-defineN/A
log1p-lowering-log1p.f64N/A
/-lowering-/.f6480.3%
Applied egg-rr80.3%
associate-*r/N/A
associate-/r/N/A
*-lowering-*.f64N/A
Applied egg-rr56.6%
Taylor expanded in i around 0
Simplified75.9%
Final simplification64.6%
(FPCore (i n) :precision binary64 (let* ((t_0 (* n (+ 100.0 (* i 50.0))))) (if (<= n -2.7e+66) t_0 (if (<= n 1.42) (* 100.0 (/ i (/ i n))) t_0))))
double code(double i, double n) {
double t_0 = n * (100.0 + (i * 50.0));
double tmp;
if (n <= -2.7e+66) {
tmp = t_0;
} else if (n <= 1.42) {
tmp = 100.0 * (i / (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))
if (n <= (-2.7d+66)) then
tmp = t_0
else if (n <= 1.42d0) then
tmp = 100.0d0 * (i / (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));
double tmp;
if (n <= -2.7e+66) {
tmp = t_0;
} else if (n <= 1.42) {
tmp = 100.0 * (i / (i / n));
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = n * (100.0 + (i * 50.0)) tmp = 0 if n <= -2.7e+66: tmp = t_0 elif n <= 1.42: tmp = 100.0 * (i / (i / n)) else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(n * Float64(100.0 + Float64(i * 50.0))) tmp = 0.0 if (n <= -2.7e+66) tmp = t_0; elseif (n <= 1.42) tmp = Float64(100.0 * Float64(i / 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)); tmp = 0.0; if (n <= -2.7e+66) tmp = t_0; elseif (n <= 1.42) tmp = 100.0 * (i / (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 * 50.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -2.7e+66], t$95$0, If[LessEqual[n, 1.42], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := n \cdot \left(100 + i \cdot 50\right)\\
\mathbf{if}\;n \leq -2.7 \cdot 10^{+66}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 1.42:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -2.7e66 or 1.4199999999999999 < n Initial program 24.6%
Taylor expanded in n around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6485.6%
Simplified85.6%
Taylor expanded in i around 0
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6462.0%
Simplified62.0%
if -2.7e66 < n < 1.4199999999999999Initial program 32.5%
Taylor expanded in i around 0
Simplified60.7%
Final simplification61.3%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* 100.0 (/ i (/ i n)))))
(if (<= i -1.75e-26)
t_0
(if (<= i 1.22e-20) (* 100.0 (+ n (* i -0.5))) t_0))))
double code(double i, double n) {
double t_0 = 100.0 * (i / (i / n));
double tmp;
if (i <= -1.75e-26) {
tmp = t_0;
} else if (i <= 1.22e-20) {
tmp = 100.0 * (n + (i * -0.5));
} 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))
if (i <= (-1.75d-26)) then
tmp = t_0
else if (i <= 1.22d-20) then
tmp = 100.0d0 * (n + (i * (-0.5d0)))
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));
double tmp;
if (i <= -1.75e-26) {
tmp = t_0;
} else if (i <= 1.22e-20) {
tmp = 100.0 * (n + (i * -0.5));
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = 100.0 * (i / (i / n)) tmp = 0 if i <= -1.75e-26: tmp = t_0 elif i <= 1.22e-20: tmp = 100.0 * (n + (i * -0.5)) else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(100.0 * Float64(i / Float64(i / n))) tmp = 0.0 if (i <= -1.75e-26) tmp = t_0; elseif (i <= 1.22e-20) tmp = Float64(100.0 * Float64(n + Float64(i * -0.5))); else tmp = t_0; end return tmp end
function tmp_2 = code(i, n) t_0 = 100.0 * (i / (i / n)); tmp = 0.0; if (i <= -1.75e-26) tmp = t_0; elseif (i <= 1.22e-20) tmp = 100.0 * (n + (i * -0.5)); else tmp = t_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[i, -1.75e-26], t$95$0, If[LessEqual[i, 1.22e-20], N[(100.0 * N[(n + N[(i * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 100 \cdot \frac{i}{\frac{i}{n}}\\
\mathbf{if}\;i \leq -1.75 \cdot 10^{-26}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;i \leq 1.22 \cdot 10^{-20}:\\
\;\;\;\;100 \cdot \left(n + i \cdot -0.5\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if i < -1.74999999999999992e-26 or 1.22000000000000003e-20 < i Initial program 50.1%
Taylor expanded in i around 0
Simplified20.4%
if -1.74999999999999992e-26 < i < 1.22000000000000003e-20Initial program 9.9%
Taylor expanded in i around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
Simplified75.5%
Taylor expanded in n around 0
Simplified75.5%
Taylor expanded in i around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6487.0%
Simplified87.0%
(FPCore (i n) :precision binary64 (if (<= i 1.95e-16) (* n 100.0) (* 50.0 (* i n))))
double code(double i, double n) {
double tmp;
if (i <= 1.95e-16) {
tmp = n * 100.0;
} else {
tmp = 50.0 * (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 <= 1.95d-16) then
tmp = n * 100.0d0
else
tmp = 50.0d0 * (i * n)
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (i <= 1.95e-16) {
tmp = n * 100.0;
} else {
tmp = 50.0 * (i * n);
}
return tmp;
}
def code(i, n): tmp = 0 if i <= 1.95e-16: tmp = n * 100.0 else: tmp = 50.0 * (i * n) return tmp
function code(i, n) tmp = 0.0 if (i <= 1.95e-16) tmp = Float64(n * 100.0); else tmp = Float64(50.0 * Float64(i * n)); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (i <= 1.95e-16) tmp = n * 100.0; else tmp = 50.0 * (i * n); end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[i, 1.95e-16], N[(n * 100.0), $MachinePrecision], N[(50.0 * N[(i * n), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq 1.95 \cdot 10^{-16}:\\
\;\;\;\;n \cdot 100\\
\mathbf{else}:\\
\;\;\;\;50 \cdot \left(i \cdot n\right)\\
\end{array}
\end{array}
if i < 1.94999999999999989e-16Initial program 21.9%
Taylor expanded in i around 0
*-lowering-*.f6465.9%
Simplified65.9%
if 1.94999999999999989e-16 < i Initial program 47.1%
Taylor expanded in n around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
expm1-defineN/A
expm1-lowering-expm1.f6442.5%
Simplified42.5%
Taylor expanded in i around 0
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6433.6%
Simplified33.6%
Taylor expanded in i around inf
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6421.4%
Simplified21.4%
Final simplification54.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 28.6%
Taylor expanded in i around 0
*-lowering-*.f6449.8%
Simplified49.8%
Final simplification49.8%
(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 2024156
(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))))