
(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 16 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 (+ (/ i n) 1.0) n)) (t_1 (/ (- t_0 1.0) (/ i n))))
(if (<= t_1 0.0)
(/ (* (expm1 (* (log1p (/ i n)) n)) 100.0) (/ i n))
(if (<= t_1 INFINITY)
(* (* (+ (/ -1.0 i) (/ t_0 i)) 100.0) n)
(/ 1.0 (fma (* 0.01 i) (- (/ 0.5 (* n n)) (/ 0.5 n)) (/ 0.01 n)))))))
double code(double i, double n) {
double t_0 = pow(((i / n) + 1.0), n);
double t_1 = (t_0 - 1.0) / (i / n);
double tmp;
if (t_1 <= 0.0) {
tmp = (expm1((log1p((i / n)) * n)) * 100.0) / (i / n);
} else if (t_1 <= ((double) INFINITY)) {
tmp = (((-1.0 / i) + (t_0 / i)) * 100.0) * n;
} else {
tmp = 1.0 / fma((0.01 * i), ((0.5 / (n * n)) - (0.5 / n)), (0.01 / n));
}
return tmp;
}
function code(i, n) t_0 = Float64(Float64(i / n) + 1.0) ^ n t_1 = Float64(Float64(t_0 - 1.0) / Float64(i / n)) tmp = 0.0 if (t_1 <= 0.0) tmp = Float64(Float64(expm1(Float64(log1p(Float64(i / n)) * n)) * 100.0) / Float64(i / n)); elseif (t_1 <= Inf) tmp = Float64(Float64(Float64(Float64(-1.0 / i) + Float64(t_0 / i)) * 100.0) * n); else tmp = Float64(1.0 / fma(Float64(0.01 * i), Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), Float64(0.01 / n))); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[Power[N[(N[(i / n), $MachinePrecision] + 1.0), $MachinePrecision], n], $MachinePrecision]}, Block[{t$95$1 = N[(N[(t$95$0 - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, 0.0], N[(N[(N[(Exp[N[(N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision] * n), $MachinePrecision]] - 1), $MachinePrecision] * 100.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(N[(N[(N[(-1.0 / i), $MachinePrecision] + N[(t$95$0 / i), $MachinePrecision]), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], N[(1.0 / N[(N[(0.01 * i), $MachinePrecision] * N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(\frac{i}{n} + 1\right)}^{n}\\
t_1 := \frac{t\_0 - 1}{\frac{i}{n}}\\
\mathbf{if}\;t\_1 \leq 0:\\
\;\;\;\;\frac{\mathsf{expm1}\left(\mathsf{log1p}\left(\frac{i}{n}\right) \cdot n\right) \cdot 100}{\frac{i}{n}}\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;\left(\left(\frac{-1}{i} + \frac{t\_0}{i}\right) \cdot 100\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(0.01 \cdot i, \frac{0.5}{n \cdot n} - \frac{0.5}{n}, \frac{0.01}{n}\right)}\\
\end{array}
\end{array}
if (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < -0.0Initial program 28.5%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6428.5
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6498.7
Applied rewrites98.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 96.7%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6496.8
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6445.8
Applied rewrites45.8%
lift-/.f64N/A
lift-*.f64N/A
lift-expm1.f64N/A
lift-*.f64N/A
lift-log1p.f64N/A
pow-to-expN/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
*-commutativeN/A
associate-/r/N/A
associate-*l*N/A
*-commutativeN/A
lift-*.f64N/A
lower-*.f64N/A
Applied rewrites45.8%
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6445.8
Applied rewrites45.8%
lift-/.f64N/A
lift-expm1.f64N/A
div-subN/A
lift-*.f64N/A
lift-log1p.f64N/A
+-commutativeN/A
lift-+.f64N/A
pow-to-expN/A
lift-pow.f64N/A
lift-/.f64N/A
sub-negN/A
mul-1-negN/A
div-invN/A
lift-/.f64N/A
lower-+.f6497.5
lift-+.f64N/A
+-commutativeN/A
lower-+.f6497.5
Applied rewrites97.5%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f640.0
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f640.0
Applied rewrites0.0%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f6424.3
Applied rewrites24.3%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6424.3
Applied rewrites24.3%
Taylor expanded in i around 0
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6499.6
Applied rewrites99.6%
Final simplification98.7%
(FPCore (i n)
:precision binary64
(let* ((t_0 (pow (+ (/ i n) 1.0) n)) (t_1 (/ (- t_0 1.0) (/ i n))))
(if (<= t_1 0.0)
(* (* (/ 100.0 i) (expm1 (* (log1p (/ i n)) n))) n)
(if (<= t_1 INFINITY)
(* (* (+ (/ -1.0 i) (/ t_0 i)) 100.0) n)
(/ 1.0 (fma (* 0.01 i) (- (/ 0.5 (* n n)) (/ 0.5 n)) (/ 0.01 n)))))))
double code(double i, double n) {
double t_0 = pow(((i / n) + 1.0), n);
double t_1 = (t_0 - 1.0) / (i / n);
double tmp;
if (t_1 <= 0.0) {
tmp = ((100.0 / i) * expm1((log1p((i / n)) * n))) * n;
} else if (t_1 <= ((double) INFINITY)) {
tmp = (((-1.0 / i) + (t_0 / i)) * 100.0) * n;
} else {
tmp = 1.0 / fma((0.01 * i), ((0.5 / (n * n)) - (0.5 / n)), (0.01 / n));
}
return tmp;
}
function code(i, n) t_0 = Float64(Float64(i / n) + 1.0) ^ n t_1 = Float64(Float64(t_0 - 1.0) / Float64(i / n)) tmp = 0.0 if (t_1 <= 0.0) tmp = Float64(Float64(Float64(100.0 / i) * expm1(Float64(log1p(Float64(i / n)) * n))) * n); elseif (t_1 <= Inf) tmp = Float64(Float64(Float64(Float64(-1.0 / i) + Float64(t_0 / i)) * 100.0) * n); else tmp = Float64(1.0 / fma(Float64(0.01 * i), Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), Float64(0.01 / n))); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[Power[N[(N[(i / n), $MachinePrecision] + 1.0), $MachinePrecision], n], $MachinePrecision]}, Block[{t$95$1 = N[(N[(t$95$0 - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, 0.0], N[(N[(N[(100.0 / i), $MachinePrecision] * N[(Exp[N[(N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision] * n), $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(N[(N[(N[(-1.0 / i), $MachinePrecision] + N[(t$95$0 / i), $MachinePrecision]), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], N[(1.0 / N[(N[(0.01 * i), $MachinePrecision] * N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(\frac{i}{n} + 1\right)}^{n}\\
t_1 := \frac{t\_0 - 1}{\frac{i}{n}}\\
\mathbf{if}\;t\_1 \leq 0:\\
\;\;\;\;\left(\frac{100}{i} \cdot \mathsf{expm1}\left(\mathsf{log1p}\left(\frac{i}{n}\right) \cdot n\right)\right) \cdot n\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;\left(\left(\frac{-1}{i} + \frac{t\_0}{i}\right) \cdot 100\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(0.01 \cdot i, \frac{0.5}{n \cdot n} - \frac{0.5}{n}, \frac{0.01}{n}\right)}\\
\end{array}
\end{array}
if (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < -0.0Initial program 28.5%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f64N/A
lower-/.f6498.3
Applied rewrites98.3%
if -0.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < +inf.0Initial program 96.7%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6496.8
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6445.8
Applied rewrites45.8%
lift-/.f64N/A
lift-*.f64N/A
lift-expm1.f64N/A
lift-*.f64N/A
lift-log1p.f64N/A
pow-to-expN/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
*-commutativeN/A
associate-/r/N/A
associate-*l*N/A
*-commutativeN/A
lift-*.f64N/A
lower-*.f64N/A
Applied rewrites45.8%
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6445.8
Applied rewrites45.8%
lift-/.f64N/A
lift-expm1.f64N/A
div-subN/A
lift-*.f64N/A
lift-log1p.f64N/A
+-commutativeN/A
lift-+.f64N/A
pow-to-expN/A
lift-pow.f64N/A
lift-/.f64N/A
sub-negN/A
mul-1-negN/A
div-invN/A
lift-/.f64N/A
lower-+.f6497.5
lift-+.f64N/A
+-commutativeN/A
lower-+.f6497.5
Applied rewrites97.5%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f640.0
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f640.0
Applied rewrites0.0%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f6424.3
Applied rewrites24.3%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6424.3
Applied rewrites24.3%
Taylor expanded in i around 0
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6499.6
Applied rewrites99.6%
Final simplification98.4%
(FPCore (i n)
:precision binary64
(let* ((t_0 (pow (+ (/ i n) 1.0) n)) (t_1 (/ (- t_0 1.0) (/ i n))))
(if (<= t_1 0.0)
(* (* (/ (expm1 i) i) 100.0) n)
(if (<= t_1 INFINITY)
(* (* (+ (/ -1.0 i) (/ t_0 i)) 100.0) n)
(/ 1.0 (fma (* 0.01 i) (- (/ 0.5 (* n n)) (/ 0.5 n)) (/ 0.01 n)))))))
double code(double i, double n) {
double t_0 = pow(((i / n) + 1.0), n);
double t_1 = (t_0 - 1.0) / (i / n);
double tmp;
if (t_1 <= 0.0) {
tmp = ((expm1(i) / i) * 100.0) * n;
} else if (t_1 <= ((double) INFINITY)) {
tmp = (((-1.0 / i) + (t_0 / i)) * 100.0) * n;
} else {
tmp = 1.0 / fma((0.01 * i), ((0.5 / (n * n)) - (0.5 / n)), (0.01 / n));
}
return tmp;
}
function code(i, n) t_0 = Float64(Float64(i / n) + 1.0) ^ n t_1 = Float64(Float64(t_0 - 1.0) / Float64(i / n)) tmp = 0.0 if (t_1 <= 0.0) tmp = Float64(Float64(Float64(expm1(i) / i) * 100.0) * n); elseif (t_1 <= Inf) tmp = Float64(Float64(Float64(Float64(-1.0 / i) + Float64(t_0 / i)) * 100.0) * n); else tmp = Float64(1.0 / fma(Float64(0.01 * i), Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), Float64(0.01 / n))); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[Power[N[(N[(i / n), $MachinePrecision] + 1.0), $MachinePrecision], n], $MachinePrecision]}, Block[{t$95$1 = N[(N[(t$95$0 - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, 0.0], N[(N[(N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(N[(N[(N[(-1.0 / i), $MachinePrecision] + N[(t$95$0 / i), $MachinePrecision]), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], N[(1.0 / N[(N[(0.01 * i), $MachinePrecision] * N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(\frac{i}{n} + 1\right)}^{n}\\
t_1 := \frac{t\_0 - 1}{\frac{i}{n}}\\
\mathbf{if}\;t\_1 \leq 0:\\
\;\;\;\;\left(\frac{\mathsf{expm1}\left(i\right)}{i} \cdot 100\right) \cdot n\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;\left(\left(\frac{-1}{i} + \frac{t\_0}{i}\right) \cdot 100\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(0.01 \cdot i, \frac{0.5}{n \cdot n} - \frac{0.5}{n}, \frac{0.01}{n}\right)}\\
\end{array}
\end{array}
if (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < -0.0Initial program 28.5%
Taylor expanded in n around inf
associate-/l*N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6479.2
Applied rewrites79.2%
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 96.7%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6496.8
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6445.8
Applied rewrites45.8%
lift-/.f64N/A
lift-*.f64N/A
lift-expm1.f64N/A
lift-*.f64N/A
lift-log1p.f64N/A
pow-to-expN/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
*-commutativeN/A
associate-/r/N/A
associate-*l*N/A
*-commutativeN/A
lift-*.f64N/A
lower-*.f64N/A
Applied rewrites45.8%
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6445.8
Applied rewrites45.8%
lift-/.f64N/A
lift-expm1.f64N/A
div-subN/A
lift-*.f64N/A
lift-log1p.f64N/A
+-commutativeN/A
lift-+.f64N/A
pow-to-expN/A
lift-pow.f64N/A
lift-/.f64N/A
sub-negN/A
mul-1-negN/A
div-invN/A
lift-/.f64N/A
lower-+.f6497.5
lift-+.f64N/A
+-commutativeN/A
lower-+.f6497.5
Applied rewrites97.5%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f640.0
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f640.0
Applied rewrites0.0%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f6424.3
Applied rewrites24.3%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6424.3
Applied rewrites24.3%
Taylor expanded in i around 0
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6499.6
Applied rewrites99.6%
Final simplification83.9%
(FPCore (i n)
:precision binary64
(let* ((t_0 (- (pow (+ (/ i n) 1.0) n) 1.0)) (t_1 (/ t_0 (/ i n))))
(if (<= t_1 0.0)
(* (* (/ (expm1 i) i) 100.0) n)
(if (<= t_1 INFINITY)
(* (* (/ t_0 i) 100.0) n)
(/ 1.0 (fma (* 0.01 i) (- (/ 0.5 (* n n)) (/ 0.5 n)) (/ 0.01 n)))))))
double code(double i, double n) {
double t_0 = pow(((i / n) + 1.0), n) - 1.0;
double t_1 = t_0 / (i / n);
double tmp;
if (t_1 <= 0.0) {
tmp = ((expm1(i) / i) * 100.0) * n;
} else if (t_1 <= ((double) INFINITY)) {
tmp = ((t_0 / i) * 100.0) * n;
} else {
tmp = 1.0 / fma((0.01 * i), ((0.5 / (n * n)) - (0.5 / n)), (0.01 / n));
}
return tmp;
}
function code(i, n) t_0 = Float64((Float64(Float64(i / n) + 1.0) ^ n) - 1.0) t_1 = Float64(t_0 / Float64(i / n)) tmp = 0.0 if (t_1 <= 0.0) tmp = Float64(Float64(Float64(expm1(i) / i) * 100.0) * n); elseif (t_1 <= Inf) tmp = Float64(Float64(Float64(t_0 / i) * 100.0) * n); else tmp = Float64(1.0 / fma(Float64(0.01 * i), Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), Float64(0.01 / n))); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[Power[N[(N[(i / n), $MachinePrecision] + 1.0), $MachinePrecision], n], $MachinePrecision] - 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 / N[(i / n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, 0.0], N[(N[(N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(N[(N[(t$95$0 / i), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], N[(1.0 / N[(N[(0.01 * i), $MachinePrecision] * N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := {\left(\frac{i}{n} + 1\right)}^{n} - 1\\
t_1 := \frac{t\_0}{\frac{i}{n}}\\
\mathbf{if}\;t\_1 \leq 0:\\
\;\;\;\;\left(\frac{\mathsf{expm1}\left(i\right)}{i} \cdot 100\right) \cdot n\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;\left(\frac{t\_0}{i} \cdot 100\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(0.01 \cdot i, \frac{0.5}{n \cdot n} - \frac{0.5}{n}, \frac{0.01}{n}\right)}\\
\end{array}
\end{array}
if (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < -0.0Initial program 28.5%
Taylor expanded in n around inf
associate-/l*N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6479.2
Applied rewrites79.2%
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 96.7%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6496.8
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6445.8
Applied rewrites45.8%
lift-/.f64N/A
lift-*.f64N/A
lift-expm1.f64N/A
lift-*.f64N/A
lift-log1p.f64N/A
pow-to-expN/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
*-commutativeN/A
associate-/r/N/A
associate-*l*N/A
*-commutativeN/A
lift-*.f64N/A
lower-*.f64N/A
Applied rewrites45.8%
lift-*.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6445.8
Applied rewrites45.8%
lift-expm1.f64N/A
lower--.f64N/A
lift-*.f64N/A
lift-log1p.f64N/A
lift-/.f64N/A
pow-to-expN/A
lower-pow.f64N/A
lift-/.f64N/A
lower-+.f6496.9
Applied rewrites96.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%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f640.0
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f640.0
Applied rewrites0.0%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f6424.3
Applied rewrites24.3%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6424.3
Applied rewrites24.3%
Taylor expanded in i around 0
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6499.6
Applied rewrites99.6%
Final simplification83.8%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* (* (/ (expm1 i) i) 100.0) n)))
(if (<= n -8e-251)
t_0
(if (<= n 0.39)
(/ 1.0 (fma (* 0.01 i) (- (/ 0.5 (* n n)) (/ 0.5 n)) (/ 0.01 n)))
t_0))))
double code(double i, double n) {
double t_0 = ((expm1(i) / i) * 100.0) * n;
double tmp;
if (n <= -8e-251) {
tmp = t_0;
} else if (n <= 0.39) {
tmp = 1.0 / fma((0.01 * i), ((0.5 / (n * n)) - (0.5 / n)), (0.01 / n));
} else {
tmp = t_0;
}
return tmp;
}
function code(i, n) t_0 = Float64(Float64(Float64(expm1(i) / i) * 100.0) * n) tmp = 0.0 if (n <= -8e-251) tmp = t_0; elseif (n <= 0.39) tmp = Float64(1.0 / fma(Float64(0.01 * i), Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), Float64(0.01 / n))); else tmp = t_0; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision]}, If[LessEqual[n, -8e-251], t$95$0, If[LessEqual[n, 0.39], N[(1.0 / N[(N[(0.01 * i), $MachinePrecision] * N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\frac{\mathsf{expm1}\left(i\right)}{i} \cdot 100\right) \cdot n\\
\mathbf{if}\;n \leq -8 \cdot 10^{-251}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 0.39:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(0.01 \cdot i, \frac{0.5}{n \cdot n} - \frac{0.5}{n}, \frac{0.01}{n}\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -8.00000000000000012e-251 or 0.39000000000000001 < n Initial program 29.0%
Taylor expanded in n around inf
associate-/l*N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower-expm1.f6487.1
Applied rewrites87.1%
if -8.00000000000000012e-251 < n < 0.39000000000000001Initial program 36.6%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6436.6
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6486.0
Applied rewrites86.0%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f6449.5
Applied rewrites49.5%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6449.4
Applied rewrites49.4%
Taylor expanded in i around 0
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6470.0
Applied rewrites70.0%
(FPCore (i n)
:precision binary64
(let* ((t_0
(*
(/
(*
(*
(fma
(fma
(+ (/ (- (/ 0.3333333333333333 n) 0.5) n) 0.16666666666666666)
i
(- 0.5 (/ 0.5 n)))
i
1.0)
i)
n)
i)
100.0)))
(if (<= n -1.66e+167)
t_0
(if (<= n 4.2e-11)
(/ 1.0 (fma (* 0.01 i) (- (/ 0.5 (* n n)) (/ 0.5 n)) (/ 0.01 n)))
t_0))))
double code(double i, double n) {
double t_0 = (((fma(fma(((((0.3333333333333333 / n) - 0.5) / n) + 0.16666666666666666), i, (0.5 - (0.5 / n))), i, 1.0) * i) * n) / i) * 100.0;
double tmp;
if (n <= -1.66e+167) {
tmp = t_0;
} else if (n <= 4.2e-11) {
tmp = 1.0 / fma((0.01 * i), ((0.5 / (n * n)) - (0.5 / n)), (0.01 / n));
} else {
tmp = t_0;
}
return tmp;
}
function code(i, n) t_0 = Float64(Float64(Float64(Float64(fma(fma(Float64(Float64(Float64(Float64(0.3333333333333333 / n) - 0.5) / n) + 0.16666666666666666), i, Float64(0.5 - Float64(0.5 / n))), i, 1.0) * i) * n) / i) * 100.0) tmp = 0.0 if (n <= -1.66e+167) tmp = t_0; elseif (n <= 4.2e-11) tmp = Float64(1.0 / fma(Float64(0.01 * i), Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), Float64(0.01 / n))); else tmp = t_0; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(N[(N[(N[(N[(N[(N[(N[(N[(0.3333333333333333 / n), $MachinePrecision] - 0.5), $MachinePrecision] / n), $MachinePrecision] + 0.16666666666666666), $MachinePrecision] * i + N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * i + 1.0), $MachinePrecision] * i), $MachinePrecision] * n), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision]}, If[LessEqual[n, -1.66e+167], t$95$0, If[LessEqual[n, 4.2e-11], N[(1.0 / N[(N[(0.01 * i), $MachinePrecision] * N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{\left(\mathsf{fma}\left(\mathsf{fma}\left(\frac{\frac{0.3333333333333333}{n} - 0.5}{n} + 0.16666666666666666, i, 0.5 - \frac{0.5}{n}\right), i, 1\right) \cdot i\right) \cdot n}{i} \cdot 100\\
\mathbf{if}\;n \leq -1.66 \cdot 10^{+167}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 4.2 \cdot 10^{-11}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(0.01 \cdot i, \frac{0.5}{n \cdot n} - \frac{0.5}{n}, \frac{0.01}{n}\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -1.6600000000000001e167 or 4.1999999999999997e-11 < n Initial program 15.3%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f64N/A
Applied rewrites52.8%
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
associate-*l/N/A
lower-/.f64N/A
Applied rewrites75.9%
if -1.6600000000000001e167 < n < 4.1999999999999997e-11Initial program 43.8%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6443.8
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6483.2
Applied rewrites83.2%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f6449.4
Applied rewrites49.4%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6449.3
Applied rewrites49.3%
Taylor expanded in i around 0
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6463.3
Applied rewrites63.3%
Final simplification69.1%
(FPCore (i n)
:precision binary64
(if (<= n -1.66e+167)
(* (fma (* n i) (fma 0.16666666666666666 i 0.5) n) 100.0)
(if (<= n 4.2e-11)
(/ 1.0 (fma (* 0.01 i) (- (/ 0.5 (* n n)) (/ 0.5 n)) (/ 0.01 n)))
(*
(* 100.0 n)
(fma
(fma
(- (+ (/ 0.3333333333333333 (* n n)) 0.16666666666666666) (/ 0.5 n))
i
(- 0.5 (/ 0.5 n)))
i
1.0)))))
double code(double i, double n) {
double tmp;
if (n <= -1.66e+167) {
tmp = fma((n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0;
} else if (n <= 4.2e-11) {
tmp = 1.0 / fma((0.01 * i), ((0.5 / (n * n)) - (0.5 / n)), (0.01 / n));
} else {
tmp = (100.0 * n) * fma(fma((((0.3333333333333333 / (n * n)) + 0.16666666666666666) - (0.5 / n)), i, (0.5 - (0.5 / n))), i, 1.0);
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -1.66e+167) tmp = Float64(fma(Float64(n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0); elseif (n <= 4.2e-11) tmp = Float64(1.0 / fma(Float64(0.01 * i), Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), Float64(0.01 / n))); else tmp = Float64(Float64(100.0 * n) * fma(fma(Float64(Float64(Float64(0.3333333333333333 / Float64(n * n)) + 0.16666666666666666) - Float64(0.5 / n)), i, Float64(0.5 - Float64(0.5 / n))), i, 1.0)); end return tmp end
code[i_, n_] := If[LessEqual[n, -1.66e+167], N[(N[(N[(n * i), $MachinePrecision] * N[(0.16666666666666666 * i + 0.5), $MachinePrecision] + n), $MachinePrecision] * 100.0), $MachinePrecision], If[LessEqual[n, 4.2e-11], N[(1.0 / N[(N[(0.01 * i), $MachinePrecision] * N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(100.0 * n), $MachinePrecision] * N[(N[(N[(N[(N[(0.3333333333333333 / N[(n * n), $MachinePrecision]), $MachinePrecision] + 0.16666666666666666), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * i + N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * i + 1.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.66 \cdot 10^{+167}:\\
\;\;\;\;\mathsf{fma}\left(n \cdot i, \mathsf{fma}\left(0.16666666666666666, i, 0.5\right), n\right) \cdot 100\\
\mathbf{elif}\;n \leq 4.2 \cdot 10^{-11}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(0.01 \cdot i, \frac{0.5}{n \cdot n} - \frac{0.5}{n}, \frac{0.01}{n}\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(100 \cdot n\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(\left(\frac{0.3333333333333333}{n \cdot n} + 0.16666666666666666\right) - \frac{0.5}{n}, i, 0.5 - \frac{0.5}{n}\right), i, 1\right)\\
\end{array}
\end{array}
if n < -1.6600000000000001e167Initial program 11.6%
Taylor expanded in i around 0
+-commutativeN/A
Applied rewrites63.6%
Taylor expanded in n around inf
Applied rewrites63.6%
if -1.6600000000000001e167 < n < 4.1999999999999997e-11Initial program 43.8%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6443.8
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6483.2
Applied rewrites83.2%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f6449.4
Applied rewrites49.4%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6449.3
Applied rewrites49.3%
Taylor expanded in i around 0
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6463.3
Applied rewrites63.3%
if 4.1999999999999997e-11 < n Initial program 17.1%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6417.1
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6472.9
Applied rewrites72.9%
lift-/.f64N/A
lift-*.f64N/A
lift-expm1.f64N/A
lift-*.f64N/A
lift-log1p.f64N/A
pow-to-expN/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
*-commutativeN/A
associate-/r/N/A
associate-*l*N/A
*-commutativeN/A
lift-*.f64N/A
lower-*.f64N/A
Applied rewrites73.1%
Taylor expanded in i around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites78.7%
Final simplification68.1%
(FPCore (i n)
:precision binary64
(if (<= n -1.66e+167)
(* (fma (* n i) (fma 0.16666666666666666 i 0.5) n) 100.0)
(if (<= n 4.2e-11)
(/ 1.0 (fma (* 0.01 i) (- (/ 0.5 (* n n)) (/ 0.5 n)) (/ 0.01 n)))
(*
(fma
(* n i)
(fma
(- (+ (/ 0.3333333333333333 (* n n)) 0.16666666666666666) (/ 0.5 n))
i
(- 0.5 (/ 0.5 n)))
n)
100.0))))
double code(double i, double n) {
double tmp;
if (n <= -1.66e+167) {
tmp = fma((n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0;
} else if (n <= 4.2e-11) {
tmp = 1.0 / fma((0.01 * i), ((0.5 / (n * n)) - (0.5 / n)), (0.01 / n));
} else {
tmp = fma((n * i), fma((((0.3333333333333333 / (n * n)) + 0.16666666666666666) - (0.5 / n)), i, (0.5 - (0.5 / n))), n) * 100.0;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -1.66e+167) tmp = Float64(fma(Float64(n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0); elseif (n <= 4.2e-11) tmp = Float64(1.0 / fma(Float64(0.01 * i), Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), Float64(0.01 / n))); else tmp = Float64(fma(Float64(n * i), fma(Float64(Float64(Float64(0.3333333333333333 / Float64(n * n)) + 0.16666666666666666) - Float64(0.5 / n)), i, Float64(0.5 - Float64(0.5 / n))), n) * 100.0); end return tmp end
code[i_, n_] := If[LessEqual[n, -1.66e+167], N[(N[(N[(n * i), $MachinePrecision] * N[(0.16666666666666666 * i + 0.5), $MachinePrecision] + n), $MachinePrecision] * 100.0), $MachinePrecision], If[LessEqual[n, 4.2e-11], N[(1.0 / N[(N[(0.01 * i), $MachinePrecision] * N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(n * i), $MachinePrecision] * N[(N[(N[(N[(0.3333333333333333 / N[(n * n), $MachinePrecision]), $MachinePrecision] + 0.16666666666666666), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * i + N[(0.5 - N[(0.5 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + n), $MachinePrecision] * 100.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.66 \cdot 10^{+167}:\\
\;\;\;\;\mathsf{fma}\left(n \cdot i, \mathsf{fma}\left(0.16666666666666666, i, 0.5\right), n\right) \cdot 100\\
\mathbf{elif}\;n \leq 4.2 \cdot 10^{-11}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(0.01 \cdot i, \frac{0.5}{n \cdot n} - \frac{0.5}{n}, \frac{0.01}{n}\right)}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(n \cdot i, \mathsf{fma}\left(\left(\frac{0.3333333333333333}{n \cdot n} + 0.16666666666666666\right) - \frac{0.5}{n}, i, 0.5 - \frac{0.5}{n}\right), n\right) \cdot 100\\
\end{array}
\end{array}
if n < -1.6600000000000001e167Initial program 11.6%
Taylor expanded in i around 0
+-commutativeN/A
Applied rewrites63.6%
Taylor expanded in n around inf
Applied rewrites63.6%
if -1.6600000000000001e167 < n < 4.1999999999999997e-11Initial program 43.8%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6443.8
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6483.2
Applied rewrites83.2%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f6449.4
Applied rewrites49.4%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6449.3
Applied rewrites49.3%
Taylor expanded in i around 0
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6463.3
Applied rewrites63.3%
if 4.1999999999999997e-11 < n Initial program 17.1%
Taylor expanded in i around 0
+-commutativeN/A
Applied rewrites78.7%
Final simplification68.1%
(FPCore (i n)
:precision binary64
(if (<= n -1.66e+167)
(* (fma (* n i) (fma 0.16666666666666666 i 0.5) n) 100.0)
(if (<= n 1.3e+54)
(/ 1.0 (fma (* 0.01 i) (- (/ 0.5 (* n n)) (/ 0.5 n)) (/ 0.01 n)))
(* (* (fma (fma 0.16666666666666666 i 0.5) i 1.0) n) 100.0))))
double code(double i, double n) {
double tmp;
if (n <= -1.66e+167) {
tmp = fma((n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0;
} else if (n <= 1.3e+54) {
tmp = 1.0 / fma((0.01 * i), ((0.5 / (n * n)) - (0.5 / n)), (0.01 / n));
} else {
tmp = (fma(fma(0.16666666666666666, i, 0.5), i, 1.0) * n) * 100.0;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -1.66e+167) tmp = Float64(fma(Float64(n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0); elseif (n <= 1.3e+54) tmp = Float64(1.0 / fma(Float64(0.01 * i), Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), Float64(0.01 / n))); else tmp = Float64(Float64(fma(fma(0.16666666666666666, i, 0.5), i, 1.0) * n) * 100.0); end return tmp end
code[i_, n_] := If[LessEqual[n, -1.66e+167], N[(N[(N[(n * i), $MachinePrecision] * N[(0.16666666666666666 * i + 0.5), $MachinePrecision] + n), $MachinePrecision] * 100.0), $MachinePrecision], If[LessEqual[n, 1.3e+54], N[(1.0 / N[(N[(0.01 * i), $MachinePrecision] * N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] + N[(0.01 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(0.16666666666666666 * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * n), $MachinePrecision] * 100.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.66 \cdot 10^{+167}:\\
\;\;\;\;\mathsf{fma}\left(n \cdot i, \mathsf{fma}\left(0.16666666666666666, i, 0.5\right), n\right) \cdot 100\\
\mathbf{elif}\;n \leq 1.3 \cdot 10^{+54}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(0.01 \cdot i, \frac{0.5}{n \cdot n} - \frac{0.5}{n}, \frac{0.01}{n}\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, i, 0.5\right), i, 1\right) \cdot n\right) \cdot 100\\
\end{array}
\end{array}
if n < -1.6600000000000001e167Initial program 11.6%
Taylor expanded in i around 0
+-commutativeN/A
Applied rewrites63.6%
Taylor expanded in n around inf
Applied rewrites63.6%
if -1.6600000000000001e167 < n < 1.30000000000000003e54Initial program 40.6%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6440.7
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6484.8
Applied rewrites84.8%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f6452.0
Applied rewrites52.0%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6451.9
Applied rewrites51.9%
Taylor expanded in i around 0
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
lower--.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6465.0
Applied rewrites65.0%
if 1.30000000000000003e54 < n Initial program 18.6%
Taylor expanded in i around 0
+-commutativeN/A
Applied rewrites78.2%
Taylor expanded in n around inf
Applied rewrites78.2%
Final simplification68.0%
(FPCore (i n)
:precision binary64
(if (<= i -7.2e+31)
(/ 0.0 i)
(if (<= i 9200000000000.0)
(* (fma (* n i) (fma 0.16666666666666666 i 0.5) n) 100.0)
(*
(/
(fma
(fma
(fma -0.5 i -0.5)
i
(* (fma (fma 0.16666666666666666 i 0.5) i 1.0) n))
n
(* (* i i) 0.3333333333333333))
n)
100.0))))
double code(double i, double n) {
double tmp;
if (i <= -7.2e+31) {
tmp = 0.0 / i;
} else if (i <= 9200000000000.0) {
tmp = fma((n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0;
} else {
tmp = (fma(fma(fma(-0.5, i, -0.5), i, (fma(fma(0.16666666666666666, i, 0.5), i, 1.0) * n)), n, ((i * i) * 0.3333333333333333)) / n) * 100.0;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (i <= -7.2e+31) tmp = Float64(0.0 / i); elseif (i <= 9200000000000.0) tmp = Float64(fma(Float64(n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0); else tmp = Float64(Float64(fma(fma(fma(-0.5, i, -0.5), i, Float64(fma(fma(0.16666666666666666, i, 0.5), i, 1.0) * n)), n, Float64(Float64(i * i) * 0.3333333333333333)) / n) * 100.0); end return tmp end
code[i_, n_] := If[LessEqual[i, -7.2e+31], N[(0.0 / i), $MachinePrecision], If[LessEqual[i, 9200000000000.0], N[(N[(N[(n * i), $MachinePrecision] * N[(0.16666666666666666 * i + 0.5), $MachinePrecision] + n), $MachinePrecision] * 100.0), $MachinePrecision], N[(N[(N[(N[(N[(-0.5 * i + -0.5), $MachinePrecision] * i + N[(N[(N[(0.16666666666666666 * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * n), $MachinePrecision]), $MachinePrecision] * n + N[(N[(i * i), $MachinePrecision] * 0.3333333333333333), $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision] * 100.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -7.2 \cdot 10^{+31}:\\
\;\;\;\;\frac{0}{i}\\
\mathbf{elif}\;i \leq 9200000000000:\\
\;\;\;\;\mathsf{fma}\left(n \cdot i, \mathsf{fma}\left(0.16666666666666666, i, 0.5\right), n\right) \cdot 100\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-0.5, i, -0.5\right), i, \mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, i, 0.5\right), i, 1\right) \cdot n\right), n, \left(i \cdot i\right) \cdot 0.3333333333333333\right)}{n} \cdot 100\\
\end{array}
\end{array}
if i < -7.19999999999999992e31Initial program 74.9%
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
+-commutativeN/A
clear-numN/A
associate-/r/N/A
distribute-lft-neg-inN/A
distribute-frac-neg2N/A
lower-fma.f64N/A
frac-2negN/A
remove-double-negN/A
lower-/.f64N/A
metadata-evalN/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
Applied rewrites66.4%
Taylor expanded in i around 0
associate-*r/N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
lower-/.f6431.5
Applied rewrites31.5%
if -7.19999999999999992e31 < i < 9.2e12Initial program 12.0%
Taylor expanded in i around 0
+-commutativeN/A
Applied rewrites68.7%
Taylor expanded in n around inf
Applied rewrites78.0%
if 9.2e12 < i Initial program 44.7%
Taylor expanded in i around 0
+-commutativeN/A
Applied rewrites39.9%
Taylor expanded in n around 0
Applied rewrites49.2%
Final simplification63.1%
(FPCore (i n)
:precision binary64
(if (<= n -1.8e+56)
(* (fma (* n i) (fma 0.16666666666666666 i 0.5) n) 100.0)
(if (<= n -8e-251)
(/ (* 100.0 i) (/ i n))
(if (<= n 5.2e-250)
(/ 0.0 i)
(if (<= n 1.45e+21)
(* (* 100.0 i) (/ n i))
(* (* (fma (fma 0.16666666666666666 i 0.5) i 1.0) n) 100.0))))))
double code(double i, double n) {
double tmp;
if (n <= -1.8e+56) {
tmp = fma((n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0;
} else if (n <= -8e-251) {
tmp = (100.0 * i) / (i / n);
} else if (n <= 5.2e-250) {
tmp = 0.0 / i;
} else if (n <= 1.45e+21) {
tmp = (100.0 * i) * (n / i);
} else {
tmp = (fma(fma(0.16666666666666666, i, 0.5), i, 1.0) * n) * 100.0;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -1.8e+56) tmp = Float64(fma(Float64(n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0); elseif (n <= -8e-251) tmp = Float64(Float64(100.0 * i) / Float64(i / n)); elseif (n <= 5.2e-250) tmp = Float64(0.0 / i); elseif (n <= 1.45e+21) tmp = Float64(Float64(100.0 * i) * Float64(n / i)); else tmp = Float64(Float64(fma(fma(0.16666666666666666, i, 0.5), i, 1.0) * n) * 100.0); end return tmp end
code[i_, n_] := If[LessEqual[n, -1.8e+56], N[(N[(N[(n * i), $MachinePrecision] * N[(0.16666666666666666 * i + 0.5), $MachinePrecision] + n), $MachinePrecision] * 100.0), $MachinePrecision], If[LessEqual[n, -8e-251], N[(N[(100.0 * i), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 5.2e-250], N[(0.0 / i), $MachinePrecision], If[LessEqual[n, 1.45e+21], N[(N[(100.0 * i), $MachinePrecision] * N[(n / i), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(0.16666666666666666 * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * n), $MachinePrecision] * 100.0), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.8 \cdot 10^{+56}:\\
\;\;\;\;\mathsf{fma}\left(n \cdot i, \mathsf{fma}\left(0.16666666666666666, i, 0.5\right), n\right) \cdot 100\\
\mathbf{elif}\;n \leq -8 \cdot 10^{-251}:\\
\;\;\;\;\frac{100 \cdot i}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 5.2 \cdot 10^{-250}:\\
\;\;\;\;\frac{0}{i}\\
\mathbf{elif}\;n \leq 1.45 \cdot 10^{+21}:\\
\;\;\;\;\left(100 \cdot i\right) \cdot \frac{n}{i}\\
\mathbf{else}:\\
\;\;\;\;\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, i, 0.5\right), i, 1\right) \cdot n\right) \cdot 100\\
\end{array}
\end{array}
if n < -1.79999999999999999e56Initial program 26.5%
Taylor expanded in i around 0
+-commutativeN/A
Applied rewrites58.4%
Taylor expanded in n around inf
Applied rewrites58.4%
if -1.79999999999999999e56 < n < -8.00000000000000012e-251Initial program 47.7%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6447.8
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6487.2
Applied rewrites87.2%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f6456.9
Applied rewrites56.9%
if -8.00000000000000012e-251 < n < 5.20000000000000016e-250Initial program 84.2%
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
+-commutativeN/A
clear-numN/A
associate-/r/N/A
distribute-lft-neg-inN/A
distribute-frac-neg2N/A
lower-fma.f64N/A
frac-2negN/A
remove-double-negN/A
lower-/.f64N/A
metadata-evalN/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
Applied rewrites16.3%
Taylor expanded in i around 0
associate-*r/N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
lower-/.f6484.2
Applied rewrites84.2%
if 5.20000000000000016e-250 < n < 1.45e21Initial program 16.0%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6416.0
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6482.0
Applied rewrites82.0%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f6459.9
Applied rewrites59.9%
lift-/.f64N/A
clear-numN/A
associate-/r/N/A
lift-/.f64N/A
clear-numN/A
lower-*.f64N/A
lower-/.f6460.0
Applied rewrites60.0%
if 1.45e21 < n Initial program 17.0%
Taylor expanded in i around 0
+-commutativeN/A
Applied rewrites79.0%
Taylor expanded in n around inf
Applied rewrites79.1%
Final simplification65.8%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* (* 100.0 i) (/ n i))))
(if (<= n -1.4e+56)
(* (fma (* n i) (fma 0.16666666666666666 i 0.5) n) 100.0)
(if (<= n -8.2e-251)
t_0
(if (<= n 5.2e-250)
(/ 0.0 i)
(if (<= n 1.45e+21)
t_0
(* (* (fma (fma 0.16666666666666666 i 0.5) i 1.0) n) 100.0)))))))
double code(double i, double n) {
double t_0 = (100.0 * i) * (n / i);
double tmp;
if (n <= -1.4e+56) {
tmp = fma((n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0;
} else if (n <= -8.2e-251) {
tmp = t_0;
} else if (n <= 5.2e-250) {
tmp = 0.0 / i;
} else if (n <= 1.45e+21) {
tmp = t_0;
} else {
tmp = (fma(fma(0.16666666666666666, i, 0.5), i, 1.0) * n) * 100.0;
}
return tmp;
}
function code(i, n) t_0 = Float64(Float64(100.0 * i) * Float64(n / i)) tmp = 0.0 if (n <= -1.4e+56) tmp = Float64(fma(Float64(n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0); elseif (n <= -8.2e-251) tmp = t_0; elseif (n <= 5.2e-250) tmp = Float64(0.0 / i); elseif (n <= 1.45e+21) tmp = t_0; else tmp = Float64(Float64(fma(fma(0.16666666666666666, i, 0.5), i, 1.0) * n) * 100.0); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(100.0 * i), $MachinePrecision] * N[(n / i), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -1.4e+56], N[(N[(N[(n * i), $MachinePrecision] * N[(0.16666666666666666 * i + 0.5), $MachinePrecision] + n), $MachinePrecision] * 100.0), $MachinePrecision], If[LessEqual[n, -8.2e-251], t$95$0, If[LessEqual[n, 5.2e-250], N[(0.0 / i), $MachinePrecision], If[LessEqual[n, 1.45e+21], t$95$0, N[(N[(N[(N[(0.16666666666666666 * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * n), $MachinePrecision] * 100.0), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(100 \cdot i\right) \cdot \frac{n}{i}\\
\mathbf{if}\;n \leq -1.4 \cdot 10^{+56}:\\
\;\;\;\;\mathsf{fma}\left(n \cdot i, \mathsf{fma}\left(0.16666666666666666, i, 0.5\right), n\right) \cdot 100\\
\mathbf{elif}\;n \leq -8.2 \cdot 10^{-251}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 5.2 \cdot 10^{-250}:\\
\;\;\;\;\frac{0}{i}\\
\mathbf{elif}\;n \leq 1.45 \cdot 10^{+21}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, i, 0.5\right), i, 1\right) \cdot n\right) \cdot 100\\
\end{array}
\end{array}
if n < -1.40000000000000004e56Initial program 26.5%
Taylor expanded in i around 0
+-commutativeN/A
Applied rewrites58.4%
Taylor expanded in n around inf
Applied rewrites58.4%
if -1.40000000000000004e56 < n < -8.1999999999999997e-251 or 5.20000000000000016e-250 < n < 1.45e21Initial program 33.6%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6433.7
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6484.9
Applied rewrites84.9%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f6458.3
Applied rewrites58.3%
lift-/.f64N/A
clear-numN/A
associate-/r/N/A
lift-/.f64N/A
clear-numN/A
lower-*.f64N/A
lower-/.f6458.2
Applied rewrites58.2%
if -8.1999999999999997e-251 < n < 5.20000000000000016e-250Initial program 84.2%
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
+-commutativeN/A
clear-numN/A
associate-/r/N/A
distribute-lft-neg-inN/A
distribute-frac-neg2N/A
lower-fma.f64N/A
frac-2negN/A
remove-double-negN/A
lower-/.f64N/A
metadata-evalN/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
Applied rewrites16.3%
Taylor expanded in i around 0
associate-*r/N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
lower-/.f6484.2
Applied rewrites84.2%
if 1.45e21 < n Initial program 17.0%
Taylor expanded in i around 0
+-commutativeN/A
Applied rewrites79.0%
Taylor expanded in n around inf
Applied rewrites79.1%
Final simplification65.8%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* (* 100.0 i) (/ n i)))
(t_1 (* (* (fma (fma 0.16666666666666666 i 0.5) i 1.0) n) 100.0)))
(if (<= n -1.4e+56)
t_1
(if (<= n -8.2e-251)
t_0
(if (<= n 5.2e-250) (/ 0.0 i) (if (<= n 1.45e+21) t_0 t_1))))))
double code(double i, double n) {
double t_0 = (100.0 * i) * (n / i);
double t_1 = (fma(fma(0.16666666666666666, i, 0.5), i, 1.0) * n) * 100.0;
double tmp;
if (n <= -1.4e+56) {
tmp = t_1;
} else if (n <= -8.2e-251) {
tmp = t_0;
} else if (n <= 5.2e-250) {
tmp = 0.0 / i;
} else if (n <= 1.45e+21) {
tmp = t_0;
} else {
tmp = t_1;
}
return tmp;
}
function code(i, n) t_0 = Float64(Float64(100.0 * i) * Float64(n / i)) t_1 = Float64(Float64(fma(fma(0.16666666666666666, i, 0.5), i, 1.0) * n) * 100.0) tmp = 0.0 if (n <= -1.4e+56) tmp = t_1; elseif (n <= -8.2e-251) tmp = t_0; elseif (n <= 5.2e-250) tmp = Float64(0.0 / i); elseif (n <= 1.45e+21) tmp = t_0; else tmp = t_1; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(100.0 * i), $MachinePrecision] * N[(n / i), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(N[(0.16666666666666666 * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * n), $MachinePrecision] * 100.0), $MachinePrecision]}, If[LessEqual[n, -1.4e+56], t$95$1, If[LessEqual[n, -8.2e-251], t$95$0, If[LessEqual[n, 5.2e-250], N[(0.0 / i), $MachinePrecision], If[LessEqual[n, 1.45e+21], t$95$0, t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(100 \cdot i\right) \cdot \frac{n}{i}\\
t_1 := \left(\mathsf{fma}\left(\mathsf{fma}\left(0.16666666666666666, i, 0.5\right), i, 1\right) \cdot n\right) \cdot 100\\
\mathbf{if}\;n \leq -1.4 \cdot 10^{+56}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;n \leq -8.2 \cdot 10^{-251}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 5.2 \cdot 10^{-250}:\\
\;\;\;\;\frac{0}{i}\\
\mathbf{elif}\;n \leq 1.45 \cdot 10^{+21}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if n < -1.40000000000000004e56 or 1.45e21 < n Initial program 21.6%
Taylor expanded in i around 0
+-commutativeN/A
Applied rewrites68.9%
Taylor expanded in n around inf
Applied rewrites68.9%
if -1.40000000000000004e56 < n < -8.1999999999999997e-251 or 5.20000000000000016e-250 < n < 1.45e21Initial program 33.6%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6433.7
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6484.9
Applied rewrites84.9%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f6458.3
Applied rewrites58.3%
lift-/.f64N/A
clear-numN/A
associate-/r/N/A
lift-/.f64N/A
clear-numN/A
lower-*.f64N/A
lower-/.f6458.2
Applied rewrites58.2%
if -8.1999999999999997e-251 < n < 5.20000000000000016e-250Initial program 84.2%
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
+-commutativeN/A
clear-numN/A
associate-/r/N/A
distribute-lft-neg-inN/A
distribute-frac-neg2N/A
lower-fma.f64N/A
frac-2negN/A
remove-double-negN/A
lower-/.f64N/A
metadata-evalN/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
Applied rewrites16.3%
Taylor expanded in i around 0
associate-*r/N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
lower-/.f6484.2
Applied rewrites84.2%
Final simplification65.8%
(FPCore (i n)
:precision binary64
(if (<= i -7.2e+31)
(/ 0.0 i)
(if (<= i 3.1e-201)
(* 100.0 n)
(if (<= i 7.4e+171) (/ (* (* 100.0 i) n) i) (/ 0.0 i)))))
double code(double i, double n) {
double tmp;
if (i <= -7.2e+31) {
tmp = 0.0 / i;
} else if (i <= 3.1e-201) {
tmp = 100.0 * n;
} else if (i <= 7.4e+171) {
tmp = ((100.0 * i) * n) / i;
} else {
tmp = 0.0 / i;
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if (i <= (-7.2d+31)) then
tmp = 0.0d0 / i
else if (i <= 3.1d-201) then
tmp = 100.0d0 * n
else if (i <= 7.4d+171) then
tmp = ((100.0d0 * i) * n) / i
else
tmp = 0.0d0 / i
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (i <= -7.2e+31) {
tmp = 0.0 / i;
} else if (i <= 3.1e-201) {
tmp = 100.0 * n;
} else if (i <= 7.4e+171) {
tmp = ((100.0 * i) * n) / i;
} else {
tmp = 0.0 / i;
}
return tmp;
}
def code(i, n): tmp = 0 if i <= -7.2e+31: tmp = 0.0 / i elif i <= 3.1e-201: tmp = 100.0 * n elif i <= 7.4e+171: tmp = ((100.0 * i) * n) / i else: tmp = 0.0 / i return tmp
function code(i, n) tmp = 0.0 if (i <= -7.2e+31) tmp = Float64(0.0 / i); elseif (i <= 3.1e-201) tmp = Float64(100.0 * n); elseif (i <= 7.4e+171) tmp = Float64(Float64(Float64(100.0 * i) * n) / i); else tmp = Float64(0.0 / i); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (i <= -7.2e+31) tmp = 0.0 / i; elseif (i <= 3.1e-201) tmp = 100.0 * n; elseif (i <= 7.4e+171) tmp = ((100.0 * i) * n) / i; else tmp = 0.0 / i; end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[i, -7.2e+31], N[(0.0 / i), $MachinePrecision], If[LessEqual[i, 3.1e-201], N[(100.0 * n), $MachinePrecision], If[LessEqual[i, 7.4e+171], N[(N[(N[(100.0 * i), $MachinePrecision] * n), $MachinePrecision] / i), $MachinePrecision], N[(0.0 / i), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -7.2 \cdot 10^{+31}:\\
\;\;\;\;\frac{0}{i}\\
\mathbf{elif}\;i \leq 3.1 \cdot 10^{-201}:\\
\;\;\;\;100 \cdot n\\
\mathbf{elif}\;i \leq 7.4 \cdot 10^{+171}:\\
\;\;\;\;\frac{\left(100 \cdot i\right) \cdot n}{i}\\
\mathbf{else}:\\
\;\;\;\;\frac{0}{i}\\
\end{array}
\end{array}
if i < -7.19999999999999992e31 or 7.39999999999999996e171 < i Initial program 68.2%
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
+-commutativeN/A
clear-numN/A
associate-/r/N/A
distribute-lft-neg-inN/A
distribute-frac-neg2N/A
lower-fma.f64N/A
frac-2negN/A
remove-double-negN/A
lower-/.f64N/A
metadata-evalN/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
Applied rewrites59.9%
Taylor expanded in i around 0
associate-*r/N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
lower-/.f6436.6
Applied rewrites36.6%
if -7.19999999999999992e31 < i < 3.0999999999999999e-201Initial program 8.8%
Taylor expanded in i around 0
lower-*.f6481.9
Applied rewrites81.9%
if 3.0999999999999999e-201 < i < 7.39999999999999996e171Initial program 26.2%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6426.2
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6481.9
Applied rewrites81.9%
Taylor expanded in i around 0
*-commutativeN/A
lower-*.f6431.9
Applied rewrites31.9%
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
associate-*l/N/A
lower-/.f64N/A
lower-*.f6454.7
Applied rewrites54.7%
(FPCore (i n) :precision binary64 (if (<= i -7.2e+31) (/ 0.0 i) (if (<= i 2e-16) (* 100.0 n) (/ 0.0 i))))
double code(double i, double n) {
double tmp;
if (i <= -7.2e+31) {
tmp = 0.0 / i;
} else if (i <= 2e-16) {
tmp = 100.0 * n;
} else {
tmp = 0.0 / i;
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if (i <= (-7.2d+31)) then
tmp = 0.0d0 / i
else if (i <= 2d-16) then
tmp = 100.0d0 * n
else
tmp = 0.0d0 / i
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (i <= -7.2e+31) {
tmp = 0.0 / i;
} else if (i <= 2e-16) {
tmp = 100.0 * n;
} else {
tmp = 0.0 / i;
}
return tmp;
}
def code(i, n): tmp = 0 if i <= -7.2e+31: tmp = 0.0 / i elif i <= 2e-16: tmp = 100.0 * n else: tmp = 0.0 / i return tmp
function code(i, n) tmp = 0.0 if (i <= -7.2e+31) tmp = Float64(0.0 / i); elseif (i <= 2e-16) tmp = Float64(100.0 * n); else tmp = Float64(0.0 / i); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (i <= -7.2e+31) tmp = 0.0 / i; elseif (i <= 2e-16) tmp = 100.0 * n; else tmp = 0.0 / i; end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[i, -7.2e+31], N[(0.0 / i), $MachinePrecision], If[LessEqual[i, 2e-16], N[(100.0 * n), $MachinePrecision], N[(0.0 / i), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -7.2 \cdot 10^{+31}:\\
\;\;\;\;\frac{0}{i}\\
\mathbf{elif}\;i \leq 2 \cdot 10^{-16}:\\
\;\;\;\;100 \cdot n\\
\mathbf{else}:\\
\;\;\;\;\frac{0}{i}\\
\end{array}
\end{array}
if i < -7.19999999999999992e31 or 2e-16 < i Initial program 56.8%
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
+-commutativeN/A
clear-numN/A
associate-/r/N/A
distribute-lft-neg-inN/A
distribute-frac-neg2N/A
lower-fma.f64N/A
frac-2negN/A
remove-double-negN/A
lower-/.f64N/A
metadata-evalN/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
Applied rewrites48.4%
Taylor expanded in i around 0
associate-*r/N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
lower-/.f6427.5
Applied rewrites27.5%
if -7.19999999999999992e31 < i < 2e-16Initial program 10.5%
Taylor expanded in i around 0
lower-*.f6479.8
Applied rewrites79.8%
(FPCore (i n) :precision binary64 (* 100.0 n))
double code(double i, double n) {
return 100.0 * n;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
code = 100.0d0 * n
end function
public static double code(double i, double n) {
return 100.0 * n;
}
def code(i, n): return 100.0 * n
function code(i, n) return Float64(100.0 * n) end
function tmp = code(i, n) tmp = 100.0 * n; end
code[i_, n_] := N[(100.0 * n), $MachinePrecision]
\begin{array}{l}
\\
100 \cdot n
\end{array}
Initial program 30.7%
Taylor expanded in i around 0
lower-*.f6447.3
Applied rewrites47.3%
(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 2024268
(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))))