
(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 13 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 2e-170)
(/ (* (expm1 (* n (log1p (/ i n)))) 100.0) (/ i n))
(if (<= t_1 INFINITY)
(* 100.0 (+ (* t_0 (/ n i)) (* n (/ -1.0 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-170) {
tmp = (expm1((n * log1p((i / n)))) * 100.0) / (i / n);
} else if (t_1 <= ((double) INFINITY)) {
tmp = 100.0 * ((t_0 * (n / i)) + (n * (-1.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);
double t_1 = (t_0 + -1.0) / (i / n);
double tmp;
if (t_1 <= 2e-170) {
tmp = (Math.expm1((n * Math.log1p((i / n)))) * 100.0) / (i / n);
} else if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = 100.0 * ((t_0 * (n / i)) + (n * (-1.0 / 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-170: tmp = (math.expm1((n * math.log1p((i / n)))) * 100.0) / (i / n) elif t_1 <= math.inf: tmp = 100.0 * ((t_0 * (n / i)) + (n * (-1.0 / 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-170) tmp = Float64(Float64(expm1(Float64(n * log1p(Float64(i / n)))) * 100.0) / Float64(i / n)); elseif (t_1 <= Inf) tmp = Float64(100.0 * Float64(Float64(t_0 * Float64(n / i)) + Float64(n * Float64(-1.0 / 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-170], N[(N[(N[(Exp[N[(n * N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]] - 1), $MachinePrecision] * 100.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(100.0 * N[(N[(t$95$0 * N[(n / i), $MachinePrecision]), $MachinePrecision] + N[(n * N[(-1.0 / i), $MachinePrecision]), $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^{-170}:\\
\;\;\;\;\frac{\mathsf{expm1}\left(n \cdot \mathsf{log1p}\left(\frac{i}{n}\right)\right) \cdot 100}{\frac{i}{n}}\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;100 \cdot \left(t\_0 \cdot \frac{n}{i} + n \cdot \frac{-1}{i}\right)\\
\mathbf{else}:\\
\;\;\;\;n \cdot 100\\
\end{array}
\end{array}
if (/.f64 (-.f64 (pow.f64 (+.f64 1 (/.f64 i n)) n) 1) (/.f64 i n)) < 1.99999999999999997e-170Initial program 28.4%
associate-*r/28.4%
*-commutative28.4%
add-exp-log28.4%
expm1-define28.4%
log-pow39.4%
log1p-define97.5%
Applied egg-rr97.5%
if 1.99999999999999997e-170 < (/.f64 (-.f64 (pow.f64 (+.f64 1 (/.f64 i n)) n) 1) (/.f64 i n)) < +inf.0Initial program 99.7%
div-sub99.7%
clear-num99.7%
div-inv99.8%
cancel-sign-sub-inv99.8%
div-inv99.8%
clear-num99.8%
Applied egg-rr99.8%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 1 (/.f64 i n)) n) 1) (/.f64 i n)) Initial program 0.0%
Taylor expanded in i around 0 81.0%
*-commutative81.0%
Simplified81.0%
Final simplification94.5%
(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-170)
(* (expm1 (* n (log1p (/ i n)))) (* 100.0 (/ n i)))
(if (<= t_1 INFINITY)
(* 100.0 (+ (* t_0 (/ n i)) (* n (/ -1.0 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-170) {
tmp = expm1((n * log1p((i / n)))) * (100.0 * (n / i));
} else if (t_1 <= ((double) INFINITY)) {
tmp = 100.0 * ((t_0 * (n / i)) + (n * (-1.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);
double t_1 = (t_0 + -1.0) / (i / n);
double tmp;
if (t_1 <= 2e-170) {
tmp = Math.expm1((n * Math.log1p((i / n)))) * (100.0 * (n / i));
} else if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = 100.0 * ((t_0 * (n / i)) + (n * (-1.0 / 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-170: tmp = math.expm1((n * math.log1p((i / n)))) * (100.0 * (n / i)) elif t_1 <= math.inf: tmp = 100.0 * ((t_0 * (n / i)) + (n * (-1.0 / 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-170) tmp = Float64(expm1(Float64(n * log1p(Float64(i / n)))) * Float64(100.0 * Float64(n / i))); elseif (t_1 <= Inf) tmp = Float64(100.0 * Float64(Float64(t_0 * Float64(n / i)) + Float64(n * Float64(-1.0 / 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-170], N[(N[(Exp[N[(n * N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]] - 1), $MachinePrecision] * N[(100.0 * N[(n / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(100.0 * N[(N[(t$95$0 * N[(n / i), $MachinePrecision]), $MachinePrecision] + N[(n * N[(-1.0 / i), $MachinePrecision]), $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^{-170}:\\
\;\;\;\;\mathsf{expm1}\left(n \cdot \mathsf{log1p}\left(\frac{i}{n}\right)\right) \cdot \left(100 \cdot \frac{n}{i}\right)\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;100 \cdot \left(t\_0 \cdot \frac{n}{i} + n \cdot \frac{-1}{i}\right)\\
\mathbf{else}:\\
\;\;\;\;n \cdot 100\\
\end{array}
\end{array}
if (/.f64 (-.f64 (pow.f64 (+.f64 1 (/.f64 i n)) n) 1) (/.f64 i n)) < 1.99999999999999997e-170Initial program 28.4%
associate-*r/28.4%
sub-neg28.4%
distribute-lft-in28.4%
metadata-eval28.4%
metadata-eval28.4%
Simplified28.4%
metadata-eval28.4%
metadata-eval28.4%
distribute-lft-in28.4%
sub-neg28.4%
associate-*r/28.4%
*-commutative28.4%
div-inv28.4%
clear-num28.4%
associate-*l*28.3%
add-exp-log28.3%
expm1-define28.3%
log-pow39.5%
log1p-define96.5%
Applied egg-rr96.5%
if 1.99999999999999997e-170 < (/.f64 (-.f64 (pow.f64 (+.f64 1 (/.f64 i n)) n) 1) (/.f64 i n)) < +inf.0Initial program 99.7%
div-sub99.7%
clear-num99.7%
div-inv99.8%
cancel-sign-sub-inv99.8%
div-inv99.8%
clear-num99.8%
Applied egg-rr99.8%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 1 (/.f64 i n)) n) 1) (/.f64 i n)) Initial program 0.0%
Taylor expanded in i around 0 81.0%
*-commutative81.0%
Simplified81.0%
Final simplification93.8%
(FPCore (i n)
:precision binary64
(let* ((t_0 (log (/ i n))) (t_1 (* 100.0 (* n (/ (expm1 i) i)))))
(if (<= n -3.9e-83)
t_1
(if (<= n -9.5e-189)
(* 100.0 (/ (* n t_0) (/ i n)))
(if (<= n 1.08e-120)
(/ 0.0 (/ i n))
(if (<= n 1.38e-27) (* 100.0 (* (pow n 2.0) (/ t_0 i))) t_1))))))
double code(double i, double n) {
double t_0 = log((i / n));
double t_1 = 100.0 * (n * (expm1(i) / i));
double tmp;
if (n <= -3.9e-83) {
tmp = t_1;
} else if (n <= -9.5e-189) {
tmp = 100.0 * ((n * t_0) / (i / n));
} else if (n <= 1.08e-120) {
tmp = 0.0 / (i / n);
} else if (n <= 1.38e-27) {
tmp = 100.0 * (pow(n, 2.0) * (t_0 / i));
} else {
tmp = t_1;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = Math.log((i / n));
double t_1 = 100.0 * (n * (Math.expm1(i) / i));
double tmp;
if (n <= -3.9e-83) {
tmp = t_1;
} else if (n <= -9.5e-189) {
tmp = 100.0 * ((n * t_0) / (i / n));
} else if (n <= 1.08e-120) {
tmp = 0.0 / (i / n);
} else if (n <= 1.38e-27) {
tmp = 100.0 * (Math.pow(n, 2.0) * (t_0 / i));
} else {
tmp = t_1;
}
return tmp;
}
def code(i, n): t_0 = math.log((i / n)) t_1 = 100.0 * (n * (math.expm1(i) / i)) tmp = 0 if n <= -3.9e-83: tmp = t_1 elif n <= -9.5e-189: tmp = 100.0 * ((n * t_0) / (i / n)) elif n <= 1.08e-120: tmp = 0.0 / (i / n) elif n <= 1.38e-27: tmp = 100.0 * (math.pow(n, 2.0) * (t_0 / i)) else: tmp = t_1 return tmp
function code(i, n) t_0 = log(Float64(i / n)) t_1 = Float64(100.0 * Float64(n * Float64(expm1(i) / i))) tmp = 0.0 if (n <= -3.9e-83) tmp = t_1; elseif (n <= -9.5e-189) tmp = Float64(100.0 * Float64(Float64(n * t_0) / Float64(i / n))); elseif (n <= 1.08e-120) tmp = Float64(0.0 / Float64(i / n)); elseif (n <= 1.38e-27) tmp = Float64(100.0 * Float64((n ^ 2.0) * Float64(t_0 / i))); else tmp = t_1; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[Log[N[(i / n), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(100.0 * N[(n * N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -3.9e-83], t$95$1, If[LessEqual[n, -9.5e-189], N[(100.0 * N[(N[(n * t$95$0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.08e-120], N[(0.0 / N[(i / n), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.38e-27], N[(100.0 * N[(N[Power[n, 2.0], $MachinePrecision] * N[(t$95$0 / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \log \left(\frac{i}{n}\right)\\
t_1 := 100 \cdot \left(n \cdot \frac{\mathsf{expm1}\left(i\right)}{i}\right)\\
\mathbf{if}\;n \leq -3.9 \cdot 10^{-83}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;n \leq -9.5 \cdot 10^{-189}:\\
\;\;\;\;100 \cdot \frac{n \cdot t\_0}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 1.08 \cdot 10^{-120}:\\
\;\;\;\;\frac{0}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 1.38 \cdot 10^{-27}:\\
\;\;\;\;100 \cdot \left({n}^{2} \cdot \frac{t\_0}{i}\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if n < -3.9e-83 or 1.38e-27 < n Initial program 21.7%
Taylor expanded in n around inf 38.6%
*-commutative38.6%
associate-/l*38.5%
expm1-define88.7%
Simplified88.7%
if -3.9e-83 < n < -9.499999999999999e-189Initial program 36.8%
Taylor expanded in i around inf 36.5%
Taylor expanded in n around 0 0.0%
neg-mul-10.0%
sub-neg0.0%
log-div84.7%
Simplified84.7%
if -9.499999999999999e-189 < n < 1.0800000000000001e-120Initial program 59.7%
associate-*r/59.7%
sub-neg59.7%
distribute-lft-in59.7%
metadata-eval59.7%
metadata-eval59.7%
Simplified59.7%
Taylor expanded in i around 0 77.2%
if 1.0800000000000001e-120 < n < 1.38e-27Initial program 14.5%
associate-*r/14.5%
sub-neg14.5%
distribute-lft-in14.5%
metadata-eval14.5%
metadata-eval14.5%
Simplified14.5%
metadata-eval14.5%
metadata-eval14.5%
distribute-lft-in14.5%
sub-neg14.5%
associate-*r/14.5%
*-commutative14.5%
div-inv14.5%
clear-num14.6%
associate-*l*14.6%
add-exp-log14.6%
expm1-define14.6%
log-pow63.6%
log1p-define89.1%
Applied egg-rr89.1%
*-commutative89.1%
associate-*r*89.1%
associate-*r/84.2%
Applied egg-rr84.2%
Taylor expanded in n around 0 68.1%
associate-/l*68.0%
neg-mul-168.0%
sub-neg68.0%
log-div57.8%
Simplified57.8%
Final simplification84.4%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* 100.0 (* n (/ (expm1 i) i))))
(t_1 (* 100.0 (/ (* n (log (/ i n))) (/ i n)))))
(if (<= n -4.3e-83)
t_0
(if (<= n -1.7e-190)
t_1
(if (<= n 3.8e-120) (/ 0.0 (/ i n)) (if (<= n 1.3e-27) t_1 t_0))))))
double code(double i, double n) {
double t_0 = 100.0 * (n * (expm1(i) / i));
double t_1 = 100.0 * ((n * log((i / n))) / (i / n));
double tmp;
if (n <= -4.3e-83) {
tmp = t_0;
} else if (n <= -1.7e-190) {
tmp = t_1;
} else if (n <= 3.8e-120) {
tmp = 0.0 / (i / n);
} else if (n <= 1.3e-27) {
tmp = t_1;
} 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 t_1 = 100.0 * ((n * Math.log((i / n))) / (i / n));
double tmp;
if (n <= -4.3e-83) {
tmp = t_0;
} else if (n <= -1.7e-190) {
tmp = t_1;
} else if (n <= 3.8e-120) {
tmp = 0.0 / (i / n);
} else if (n <= 1.3e-27) {
tmp = t_1;
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = 100.0 * (n * (math.expm1(i) / i)) t_1 = 100.0 * ((n * math.log((i / n))) / (i / n)) tmp = 0 if n <= -4.3e-83: tmp = t_0 elif n <= -1.7e-190: tmp = t_1 elif n <= 3.8e-120: tmp = 0.0 / (i / n) elif n <= 1.3e-27: tmp = t_1 else: tmp = t_0 return tmp
function code(i, n) t_0 = Float64(100.0 * Float64(n * Float64(expm1(i) / i))) t_1 = Float64(100.0 * Float64(Float64(n * log(Float64(i / n))) / Float64(i / n))) tmp = 0.0 if (n <= -4.3e-83) tmp = t_0; elseif (n <= -1.7e-190) tmp = t_1; elseif (n <= 3.8e-120) tmp = Float64(0.0 / Float64(i / n)); elseif (n <= 1.3e-27) tmp = t_1; else tmp = t_0; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(100.0 * N[(n * N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(100.0 * N[(N[(n * N[Log[N[(i / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -4.3e-83], t$95$0, If[LessEqual[n, -1.7e-190], t$95$1, If[LessEqual[n, 3.8e-120], N[(0.0 / N[(i / n), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 1.3e-27], t$95$1, t$95$0]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 100 \cdot \left(n \cdot \frac{\mathsf{expm1}\left(i\right)}{i}\right)\\
t_1 := 100 \cdot \frac{n \cdot \log \left(\frac{i}{n}\right)}{\frac{i}{n}}\\
\mathbf{if}\;n \leq -4.3 \cdot 10^{-83}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq -1.7 \cdot 10^{-190}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;n \leq 3.8 \cdot 10^{-120}:\\
\;\;\;\;\frac{0}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 1.3 \cdot 10^{-27}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -4.30000000000000033e-83 or 1.30000000000000009e-27 < n Initial program 21.7%
Taylor expanded in n around inf 38.6%
*-commutative38.6%
associate-/l*38.5%
expm1-define88.7%
Simplified88.7%
if -4.30000000000000033e-83 < n < -1.69999999999999991e-190 or 3.7999999999999997e-120 < n < 1.30000000000000009e-27Initial program 24.0%
Taylor expanded in i around inf 23.4%
Taylor expanded in n around 0 39.2%
neg-mul-139.2%
sub-neg39.2%
log-div69.1%
Simplified69.1%
if -1.69999999999999991e-190 < n < 3.7999999999999997e-120Initial program 59.7%
associate-*r/59.7%
sub-neg59.7%
distribute-lft-in59.7%
metadata-eval59.7%
metadata-eval59.7%
Simplified59.7%
Taylor expanded in i around 0 77.2%
Final simplification84.4%
(FPCore (i n) :precision binary64 (if (or (<= i -2.1e-8) (not (<= i 2e-20))) (* 100.0 (/ (expm1 i) (/ i n))) (* n (+ 100.0 (* i 50.0)))))
double code(double i, double n) {
double tmp;
if ((i <= -2.1e-8) || !(i <= 2e-20)) {
tmp = 100.0 * (expm1(i) / (i / n));
} else {
tmp = n * (100.0 + (i * 50.0));
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if ((i <= -2.1e-8) || !(i <= 2e-20)) {
tmp = 100.0 * (Math.expm1(i) / (i / n));
} else {
tmp = n * (100.0 + (i * 50.0));
}
return tmp;
}
def code(i, n): tmp = 0 if (i <= -2.1e-8) or not (i <= 2e-20): tmp = 100.0 * (math.expm1(i) / (i / n)) else: tmp = n * (100.0 + (i * 50.0)) return tmp
function code(i, n) tmp = 0.0 if ((i <= -2.1e-8) || !(i <= 2e-20)) tmp = Float64(100.0 * Float64(expm1(i) / Float64(i / n))); else tmp = Float64(n * Float64(100.0 + Float64(i * 50.0))); end return tmp end
code[i_, n_] := If[Or[LessEqual[i, -2.1e-8], N[Not[LessEqual[i, 2e-20]], $MachinePrecision]], N[(100.0 * N[(N[(Exp[i] - 1), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(n * N[(100.0 + N[(i * 50.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -2.1 \cdot 10^{-8} \lor \neg \left(i \leq 2 \cdot 10^{-20}\right):\\
\;\;\;\;100 \cdot \frac{\mathsf{expm1}\left(i\right)}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;n \cdot \left(100 + i \cdot 50\right)\\
\end{array}
\end{array}
if i < -2.09999999999999994e-8 or 1.99999999999999989e-20 < i Initial program 45.8%
Taylor expanded in n around inf 59.2%
expm1-define59.5%
Simplified59.5%
if -2.09999999999999994e-8 < i < 1.99999999999999989e-20Initial program 11.4%
Taylor expanded in i around 0 80.7%
Taylor expanded in n around inf 81.4%
*-commutative81.4%
Simplified81.4%
Final simplification70.9%
(FPCore (i n) :precision binary64 (if (or (<= n -6.3e-156) (not (<= n 2.5e-87))) (* n (/ (* 100.0 (expm1 i)) i)) (/ 0.0 (/ i n))))
double code(double i, double n) {
double tmp;
if ((n <= -6.3e-156) || !(n <= 2.5e-87)) {
tmp = n * ((100.0 * expm1(i)) / i);
} else {
tmp = 0.0 / (i / n);
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if ((n <= -6.3e-156) || !(n <= 2.5e-87)) {
tmp = n * ((100.0 * Math.expm1(i)) / i);
} else {
tmp = 0.0 / (i / n);
}
return tmp;
}
def code(i, n): tmp = 0 if (n <= -6.3e-156) or not (n <= 2.5e-87): tmp = n * ((100.0 * math.expm1(i)) / i) else: tmp = 0.0 / (i / n) return tmp
function code(i, n) tmp = 0.0 if ((n <= -6.3e-156) || !(n <= 2.5e-87)) tmp = Float64(n * Float64(Float64(100.0 * expm1(i)) / i)); else tmp = Float64(0.0 / Float64(i / n)); end return tmp end
code[i_, n_] := If[Or[LessEqual[n, -6.3e-156], N[Not[LessEqual[n, 2.5e-87]], $MachinePrecision]], N[(n * N[(N[(100.0 * N[(Exp[i] - 1), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], N[(0.0 / N[(i / n), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -6.3 \cdot 10^{-156} \lor \neg \left(n \leq 2.5 \cdot 10^{-87}\right):\\
\;\;\;\;n \cdot \frac{100 \cdot \mathsf{expm1}\left(i\right)}{i}\\
\mathbf{else}:\\
\;\;\;\;\frac{0}{\frac{i}{n}}\\
\end{array}
\end{array}
if n < -6.3000000000000001e-156 or 2.50000000000000021e-87 < n Initial program 21.2%
associate-*r/21.2%
sub-neg21.2%
distribute-lft-in21.2%
metadata-eval21.2%
metadata-eval21.2%
Simplified21.2%
Taylor expanded in n around inf 35.1%
*-commutative35.1%
Simplified35.1%
Taylor expanded in i around inf 35.5%
fma-neg35.4%
metadata-eval35.4%
associate-*r/35.5%
fma-undefine35.5%
metadata-eval35.5%
distribute-lft-in35.5%
metadata-eval35.5%
sub-neg35.5%
expm1-define83.0%
Simplified83.0%
if -6.3000000000000001e-156 < n < 2.50000000000000021e-87Initial program 54.4%
associate-*r/54.4%
sub-neg54.4%
distribute-lft-in54.4%
metadata-eval54.4%
metadata-eval54.4%
Simplified54.4%
Taylor expanded in i around 0 69.8%
Final simplification80.3%
(FPCore (i n) :precision binary64 (if (or (<= n -1e-158) (not (<= n 2.8e-87))) (* 100.0 (* n (/ (expm1 i) i))) (/ 0.0 (/ i n))))
double code(double i, double n) {
double tmp;
if ((n <= -1e-158) || !(n <= 2.8e-87)) {
tmp = 100.0 * (n * (expm1(i) / i));
} else {
tmp = 0.0 / (i / n);
}
return tmp;
}
public static double code(double i, double n) {
double tmp;
if ((n <= -1e-158) || !(n <= 2.8e-87)) {
tmp = 100.0 * (n * (Math.expm1(i) / i));
} else {
tmp = 0.0 / (i / n);
}
return tmp;
}
def code(i, n): tmp = 0 if (n <= -1e-158) or not (n <= 2.8e-87): tmp = 100.0 * (n * (math.expm1(i) / i)) else: tmp = 0.0 / (i / n) return tmp
function code(i, n) tmp = 0.0 if ((n <= -1e-158) || !(n <= 2.8e-87)) tmp = Float64(100.0 * Float64(n * Float64(expm1(i) / i))); else tmp = Float64(0.0 / Float64(i / n)); end return tmp end
code[i_, n_] := If[Or[LessEqual[n, -1e-158], N[Not[LessEqual[n, 2.8e-87]], $MachinePrecision]], N[(100.0 * N[(n * N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.0 / N[(i / n), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1 \cdot 10^{-158} \lor \neg \left(n \leq 2.8 \cdot 10^{-87}\right):\\
\;\;\;\;100 \cdot \left(n \cdot \frac{\mathsf{expm1}\left(i\right)}{i}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{0}{\frac{i}{n}}\\
\end{array}
\end{array}
if n < -1.00000000000000006e-158 or 2.8000000000000001e-87 < n Initial program 21.2%
Taylor expanded in n around inf 35.4%
*-commutative35.4%
associate-/l*35.4%
expm1-define83.0%
Simplified83.0%
if -1.00000000000000006e-158 < n < 2.8000000000000001e-87Initial program 54.4%
associate-*r/54.4%
sub-neg54.4%
distribute-lft-in54.4%
metadata-eval54.4%
metadata-eval54.4%
Simplified54.4%
Taylor expanded in i around 0 69.8%
Final simplification80.3%
(FPCore (i n) :precision binary64 (if (<= n -1.7e-73) (* n (+ 100.0 (* i 50.0))) (if (<= n 2.8e-87) (/ 0.0 (/ i n)) (+ (* n 100.0) (* 50.0 (* i n))))))
double code(double i, double n) {
double tmp;
if (n <= -1.7e-73) {
tmp = n * (100.0 + (i * 50.0));
} else if (n <= 2.8e-87) {
tmp = 0.0 / (i / n);
} else {
tmp = (n * 100.0) + (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 (n <= (-1.7d-73)) then
tmp = n * (100.0d0 + (i * 50.0d0))
else if (n <= 2.8d-87) then
tmp = 0.0d0 / (i / n)
else
tmp = (n * 100.0d0) + (50.0d0 * (i * n))
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (n <= -1.7e-73) {
tmp = n * (100.0 + (i * 50.0));
} else if (n <= 2.8e-87) {
tmp = 0.0 / (i / n);
} else {
tmp = (n * 100.0) + (50.0 * (i * n));
}
return tmp;
}
def code(i, n): tmp = 0 if n <= -1.7e-73: tmp = n * (100.0 + (i * 50.0)) elif n <= 2.8e-87: tmp = 0.0 / (i / n) else: tmp = (n * 100.0) + (50.0 * (i * n)) return tmp
function code(i, n) tmp = 0.0 if (n <= -1.7e-73) tmp = Float64(n * Float64(100.0 + Float64(i * 50.0))); elseif (n <= 2.8e-87) tmp = Float64(0.0 / Float64(i / n)); else tmp = Float64(Float64(n * 100.0) + Float64(50.0 * Float64(i * n))); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (n <= -1.7e-73) tmp = n * (100.0 + (i * 50.0)); elseif (n <= 2.8e-87) tmp = 0.0 / (i / n); else tmp = (n * 100.0) + (50.0 * (i * n)); end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[n, -1.7e-73], N[(n * N[(100.0 + N[(i * 50.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 2.8e-87], N[(0.0 / N[(i / n), $MachinePrecision]), $MachinePrecision], N[(N[(n * 100.0), $MachinePrecision] + N[(50.0 * N[(i * n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.7 \cdot 10^{-73}:\\
\;\;\;\;n \cdot \left(100 + i \cdot 50\right)\\
\mathbf{elif}\;n \leq 2.8 \cdot 10^{-87}:\\
\;\;\;\;\frac{0}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;n \cdot 100 + 50 \cdot \left(i \cdot n\right)\\
\end{array}
\end{array}
if n < -1.7000000000000001e-73Initial program 24.2%
Taylor expanded in i around 0 58.6%
Taylor expanded in n around inf 58.6%
*-commutative58.6%
Simplified58.6%
if -1.7000000000000001e-73 < n < 2.8000000000000001e-87Initial program 51.2%
associate-*r/51.2%
sub-neg51.2%
distribute-lft-in51.2%
metadata-eval51.2%
metadata-eval51.2%
Simplified51.2%
Taylor expanded in i around 0 64.3%
if 2.8000000000000001e-87 < n Initial program 17.7%
Taylor expanded in i around 0 64.0%
Taylor expanded in n around inf 64.0%
*-commutative64.0%
Simplified64.0%
Final simplification62.2%
(FPCore (i n) :precision binary64 (if (or (<= i -6.2e+234) (not (<= i 5e-14))) (* 100.0 (* i (/ n i))) (* n 100.0)))
double code(double i, double n) {
double tmp;
if ((i <= -6.2e+234) || !(i <= 5e-14)) {
tmp = 100.0 * (i * (n / i));
} else {
tmp = n * 100.0;
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if ((i <= (-6.2d+234)) .or. (.not. (i <= 5d-14))) then
tmp = 100.0d0 * (i * (n / i))
else
tmp = n * 100.0d0
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if ((i <= -6.2e+234) || !(i <= 5e-14)) {
tmp = 100.0 * (i * (n / i));
} else {
tmp = n * 100.0;
}
return tmp;
}
def code(i, n): tmp = 0 if (i <= -6.2e+234) or not (i <= 5e-14): tmp = 100.0 * (i * (n / i)) else: tmp = n * 100.0 return tmp
function code(i, n) tmp = 0.0 if ((i <= -6.2e+234) || !(i <= 5e-14)) tmp = Float64(100.0 * Float64(i * Float64(n / i))); else tmp = Float64(n * 100.0); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if ((i <= -6.2e+234) || ~((i <= 5e-14))) tmp = 100.0 * (i * (n / i)); else tmp = n * 100.0; end tmp_2 = tmp; end
code[i_, n_] := If[Or[LessEqual[i, -6.2e+234], N[Not[LessEqual[i, 5e-14]], $MachinePrecision]], N[(100.0 * N[(i * N[(n / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(n * 100.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -6.2 \cdot 10^{+234} \lor \neg \left(i \leq 5 \cdot 10^{-14}\right):\\
\;\;\;\;100 \cdot \left(i \cdot \frac{n}{i}\right)\\
\mathbf{else}:\\
\;\;\;\;n \cdot 100\\
\end{array}
\end{array}
if i < -6.19999999999999979e234 or 5.0000000000000002e-14 < i Initial program 47.2%
Taylor expanded in i around 0 22.0%
clear-num22.0%
associate-/r/22.0%
clear-num18.0%
Applied egg-rr18.0%
if -6.19999999999999979e234 < i < 5.0000000000000002e-14Initial program 18.2%
Taylor expanded in i around 0 64.8%
*-commutative64.8%
Simplified64.8%
Final simplification49.1%
(FPCore (i n) :precision binary64 (if (or (<= n -1.45e+33) (not (<= n 0.23))) (* n (+ 100.0 (* i 50.0))) (* 100.0 (/ i (/ i n)))))
double code(double i, double n) {
double tmp;
if ((n <= -1.45e+33) || !(n <= 0.23)) {
tmp = n * (100.0 + (i * 50.0));
} else {
tmp = 100.0 * (i / (i / n));
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if ((n <= (-1.45d+33)) .or. (.not. (n <= 0.23d0))) then
tmp = n * (100.0d0 + (i * 50.0d0))
else
tmp = 100.0d0 * (i / (i / n))
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if ((n <= -1.45e+33) || !(n <= 0.23)) {
tmp = n * (100.0 + (i * 50.0));
} else {
tmp = 100.0 * (i / (i / n));
}
return tmp;
}
def code(i, n): tmp = 0 if (n <= -1.45e+33) or not (n <= 0.23): tmp = n * (100.0 + (i * 50.0)) else: tmp = 100.0 * (i / (i / n)) return tmp
function code(i, n) tmp = 0.0 if ((n <= -1.45e+33) || !(n <= 0.23)) tmp = Float64(n * Float64(100.0 + Float64(i * 50.0))); else tmp = Float64(100.0 * Float64(i / Float64(i / n))); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if ((n <= -1.45e+33) || ~((n <= 0.23))) tmp = n * (100.0 + (i * 50.0)); else tmp = 100.0 * (i / (i / n)); end tmp_2 = tmp; end
code[i_, n_] := If[Or[LessEqual[n, -1.45e+33], N[Not[LessEqual[n, 0.23]], $MachinePrecision]], N[(n * N[(100.0 + N[(i * 50.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.45 \cdot 10^{+33} \lor \neg \left(n \leq 0.23\right):\\
\;\;\;\;n \cdot \left(100 + i \cdot 50\right)\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\end{array}
\end{array}
if n < -1.45000000000000012e33 or 0.23000000000000001 < n Initial program 19.8%
Taylor expanded in i around 0 65.5%
Taylor expanded in n around inf 65.5%
*-commutative65.5%
Simplified65.5%
if -1.45000000000000012e33 < n < 0.23000000000000001Initial program 40.7%
Taylor expanded in i around 0 43.5%
Final simplification56.9%
(FPCore (i n) :precision binary64 (if (or (<= n -1.7e-73) (not (<= n 2.5e-87))) (* n (+ 100.0 (* i 50.0))) (/ 0.0 (/ i n))))
double code(double i, double n) {
double tmp;
if ((n <= -1.7e-73) || !(n <= 2.5e-87)) {
tmp = n * (100.0 + (i * 50.0));
} else {
tmp = 0.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 ((n <= (-1.7d-73)) .or. (.not. (n <= 2.5d-87))) then
tmp = n * (100.0d0 + (i * 50.0d0))
else
tmp = 0.0d0 / (i / n)
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if ((n <= -1.7e-73) || !(n <= 2.5e-87)) {
tmp = n * (100.0 + (i * 50.0));
} else {
tmp = 0.0 / (i / n);
}
return tmp;
}
def code(i, n): tmp = 0 if (n <= -1.7e-73) or not (n <= 2.5e-87): tmp = n * (100.0 + (i * 50.0)) else: tmp = 0.0 / (i / n) return tmp
function code(i, n) tmp = 0.0 if ((n <= -1.7e-73) || !(n <= 2.5e-87)) tmp = Float64(n * Float64(100.0 + Float64(i * 50.0))); else tmp = Float64(0.0 / Float64(i / n)); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if ((n <= -1.7e-73) || ~((n <= 2.5e-87))) tmp = n * (100.0 + (i * 50.0)); else tmp = 0.0 / (i / n); end tmp_2 = tmp; end
code[i_, n_] := If[Or[LessEqual[n, -1.7e-73], N[Not[LessEqual[n, 2.5e-87]], $MachinePrecision]], N[(n * N[(100.0 + N[(i * 50.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.0 / N[(i / n), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.7 \cdot 10^{-73} \lor \neg \left(n \leq 2.5 \cdot 10^{-87}\right):\\
\;\;\;\;n \cdot \left(100 + i \cdot 50\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{0}{\frac{i}{n}}\\
\end{array}
\end{array}
if n < -1.7000000000000001e-73 or 2.50000000000000021e-87 < n Initial program 20.7%
Taylor expanded in i around 0 61.5%
Taylor expanded in n around inf 61.6%
*-commutative61.6%
Simplified61.6%
if -1.7000000000000001e-73 < n < 2.50000000000000021e-87Initial program 51.2%
associate-*r/51.2%
sub-neg51.2%
distribute-lft-in51.2%
metadata-eval51.2%
metadata-eval51.2%
Simplified51.2%
Taylor expanded in i around 0 64.3%
Final simplification62.2%
(FPCore (i n) :precision binary64 (if (<= i -5e+69) (* 100.0 (* i (/ n i))) (if (<= i 1e-6) (* n 100.0) (* 100.0 (/ i (/ i n))))))
double code(double i, double n) {
double tmp;
if (i <= -5e+69) {
tmp = 100.0 * (i * (n / i));
} else if (i <= 1e-6) {
tmp = n * 100.0;
} else {
tmp = 100.0 * (i / (i / n));
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if (i <= (-5d+69)) then
tmp = 100.0d0 * (i * (n / i))
else if (i <= 1d-6) then
tmp = n * 100.0d0
else
tmp = 100.0d0 * (i / (i / n))
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (i <= -5e+69) {
tmp = 100.0 * (i * (n / i));
} else if (i <= 1e-6) {
tmp = n * 100.0;
} else {
tmp = 100.0 * (i / (i / n));
}
return tmp;
}
def code(i, n): tmp = 0 if i <= -5e+69: tmp = 100.0 * (i * (n / i)) elif i <= 1e-6: tmp = n * 100.0 else: tmp = 100.0 * (i / (i / n)) return tmp
function code(i, n) tmp = 0.0 if (i <= -5e+69) tmp = Float64(100.0 * Float64(i * Float64(n / i))); elseif (i <= 1e-6) tmp = Float64(n * 100.0); else tmp = Float64(100.0 * Float64(i / Float64(i / n))); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (i <= -5e+69) tmp = 100.0 * (i * (n / i)); elseif (i <= 1e-6) tmp = n * 100.0; else tmp = 100.0 * (i / (i / n)); end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[i, -5e+69], N[(100.0 * N[(i * N[(n / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[i, 1e-6], N[(n * 100.0), $MachinePrecision], N[(100.0 * N[(i / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -5 \cdot 10^{+69}:\\
\;\;\;\;100 \cdot \left(i \cdot \frac{n}{i}\right)\\
\mathbf{elif}\;i \leq 10^{-6}:\\
\;\;\;\;n \cdot 100\\
\mathbf{else}:\\
\;\;\;\;100 \cdot \frac{i}{\frac{i}{n}}\\
\end{array}
\end{array}
if i < -5.00000000000000036e69Initial program 56.8%
Taylor expanded in i around 0 17.7%
clear-num17.7%
associate-/r/17.7%
clear-num17.7%
Applied egg-rr17.7%
if -5.00000000000000036e69 < i < 9.99999999999999955e-7Initial program 13.8%
Taylor expanded in i around 0 74.3%
*-commutative74.3%
Simplified74.3%
if 9.99999999999999955e-7 < i Initial program 42.5%
Taylor expanded in i around 0 18.5%
Final simplification50.4%
(FPCore (i n) :precision binary64 (* n 100.0))
double code(double i, double n) {
return n * 100.0;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
code = n * 100.0d0
end function
public static double code(double i, double n) {
return n * 100.0;
}
def code(i, n): return n * 100.0
function code(i, n) return Float64(n * 100.0) end
function tmp = code(i, n) tmp = n * 100.0; end
code[i_, n_] := N[(n * 100.0), $MachinePrecision]
\begin{array}{l}
\\
n \cdot 100
\end{array}
Initial program 27.9%
Taylor expanded in i around 0 45.1%
*-commutative45.1%
Simplified45.1%
Final simplification45.1%
(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 2024039
(FPCore (i n)
:name "Compound Interest"
:precision binary64
:herbie-target
(* 100.0 (/ (- (exp (* n (if (== (+ 1.0 (/ i n)) 1.0) (/ i n) (/ (* (/ i n) (log (+ 1.0 (/ i n)))) (- (+ (/ i n) 1.0) 1.0))))) 1.0) (/ i n)))
(* 100.0 (/ (- (pow (+ 1.0 (/ i n)) n) 1.0) (/ i n))))