
(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 21 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 -1e-84)
(* (/ (- (* t_0 n) n) i) 100.0)
(if (<= t_1 INFINITY)
(/ (* (expm1 (* (log1p (/ i n)) n)) 100.0) (/ i n))
(/ 100.0 (fma (- (/ 0.5 (* n n)) (/ 0.5 n)) i (/ 1.0 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 <= -1e-84) {
tmp = (((t_0 * n) - n) / i) * 100.0;
} else if (t_1 <= ((double) INFINITY)) {
tmp = (expm1((log1p((i / n)) * n)) * 100.0) / (i / n);
} else {
tmp = 100.0 / fma(((0.5 / (n * n)) - (0.5 / n)), i, (1.0 / 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 <= -1e-84) tmp = Float64(Float64(Float64(Float64(t_0 * n) - n) / i) * 100.0); elseif (t_1 <= Inf) tmp = Float64(Float64(expm1(Float64(log1p(Float64(i / n)) * n)) * 100.0) / Float64(i / n)); else tmp = Float64(100.0 / fma(Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), i, Float64(1.0 / 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, -1e-84], N[(N[(N[(N[(t$95$0 * n), $MachinePrecision] - n), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision], If[LessEqual[t$95$1, Infinity], 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], N[(100.0 / N[(N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * i + N[(1.0 / 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 -1 \cdot 10^{-84}:\\
\;\;\;\;\frac{t\_0 \cdot n - n}{i} \cdot 100\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;\frac{\mathsf{expm1}\left(\mathsf{log1p}\left(\frac{i}{n}\right) \cdot n\right) \cdot 100}{\frac{i}{n}}\\
\mathbf{else}:\\
\;\;\;\;\frac{100}{\mathsf{fma}\left(\frac{0.5}{n \cdot n} - \frac{0.5}{n}, i, \frac{1}{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)) < -1e-84Initial program 99.7%
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
lower--.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f6499.9
Applied rewrites99.9%
lift--.f64N/A
lift-*.f64N/A
lift-/.f64N/A
associate-*l/N/A
lift-/.f64N/A
sub-divN/A
lower-/.f64N/A
lower--.f64N/A
lift-pow.f64N/A
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
lower-*.f64N/A
lift-/.f64N/A
+-commutativeN/A
lift-+.f64N/A
lift-pow.f6499.9
Applied rewrites99.9%
if -1e-84 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < +inf.0Initial program 24.5%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f6424.5
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6497.1
Applied rewrites97.1%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
Taylor expanded in i around 0
lower-+.f6413.3
Applied rewrites13.3%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6413.3
Applied rewrites13.3%
Taylor expanded in i around 0
*-commutativeN/A
lower-fma.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
lower-/.f6499.8
Applied rewrites99.8%
Final simplification97.7%
(FPCore (i n)
:precision binary64
(let* ((t_0 (pow (+ (/ i n) 1.0) n))
(t_1 (/ (- t_0 1.0) (/ i n)))
(t_2 (* (/ (- (* t_0 n) n) i) 100.0)))
(if (<= t_1 -1e-153)
t_2
(if (<= t_1 0.0)
(/
(*
(expm1
(fma (* (fma (/ 0.3333333333333333 n) (/ i n) (/ -0.5 n)) i) i i))
100.0)
(/ i n))
(if (<= t_1 INFINITY)
t_2
(/ 100.0 (fma (- (/ 0.5 (* n n)) (/ 0.5 n)) i (/ 1.0 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 t_2 = (((t_0 * n) - n) / i) * 100.0;
double tmp;
if (t_1 <= -1e-153) {
tmp = t_2;
} else if (t_1 <= 0.0) {
tmp = (expm1(fma((fma((0.3333333333333333 / n), (i / n), (-0.5 / n)) * i), i, i)) * 100.0) / (i / n);
} else if (t_1 <= ((double) INFINITY)) {
tmp = t_2;
} else {
tmp = 100.0 / fma(((0.5 / (n * n)) - (0.5 / n)), i, (1.0 / 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)) t_2 = Float64(Float64(Float64(Float64(t_0 * n) - n) / i) * 100.0) tmp = 0.0 if (t_1 <= -1e-153) tmp = t_2; elseif (t_1 <= 0.0) tmp = Float64(Float64(expm1(fma(Float64(fma(Float64(0.3333333333333333 / n), Float64(i / n), Float64(-0.5 / n)) * i), i, i)) * 100.0) / Float64(i / n)); elseif (t_1 <= Inf) tmp = t_2; else tmp = Float64(100.0 / fma(Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), i, Float64(1.0 / 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]}, Block[{t$95$2 = N[(N[(N[(N[(t$95$0 * n), $MachinePrecision] - n), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision]}, If[LessEqual[t$95$1, -1e-153], t$95$2, If[LessEqual[t$95$1, 0.0], N[(N[(N[(Exp[N[(N[(N[(N[(0.3333333333333333 / n), $MachinePrecision] * N[(i / n), $MachinePrecision] + N[(-0.5 / n), $MachinePrecision]), $MachinePrecision] * i), $MachinePrecision] * i + i), $MachinePrecision]] - 1), $MachinePrecision] * 100.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], t$95$2, N[(100.0 / N[(N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * i + N[(1.0 / 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}}\\
t_2 := \frac{t\_0 \cdot n - n}{i} \cdot 100\\
\mathbf{if}\;t\_1 \leq -1 \cdot 10^{-153}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;\frac{\mathsf{expm1}\left(\mathsf{fma}\left(\mathsf{fma}\left(\frac{0.3333333333333333}{n}, \frac{i}{n}, \frac{-0.5}{n}\right) \cdot i, i, i\right)\right) \cdot 100}{\frac{i}{n}}\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;\frac{100}{\mathsf{fma}\left(\frac{0.5}{n \cdot n} - \frac{0.5}{n}, i, \frac{1}{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)) < -1.00000000000000004e-153 or -0.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < +inf.0Initial program 92.8%
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
lower--.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f6492.9
Applied rewrites92.9%
lift--.f64N/A
lift-*.f64N/A
lift-/.f64N/A
associate-*l/N/A
lift-/.f64N/A
sub-divN/A
lower-/.f64N/A
lower--.f64N/A
lift-pow.f64N/A
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
lower-*.f64N/A
lift-/.f64N/A
+-commutativeN/A
lift-+.f64N/A
lift-pow.f6493.1
Applied rewrites93.1%
if -1.00000000000000004e-153 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < -0.0Initial program 15.1%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
associate-/r/N/A
lift-/.f64N/A
clear-numN/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6414.6
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6498.0
Applied rewrites98.0%
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-expm1.f64N/A
lift-*.f64N/A
lift-log1p.f64N/A
pow-to-expN/A
lift-/.f64N/A
*-commutativeN/A
lift-/.f64N/A
clear-numN/A
lift-/.f64N/A
div-invN/A
lower-/.f64N/A
Applied rewrites99.6%
Taylor expanded in i around 0
+-commutativeN/A
distribute-rgt-inN/A
*-lft-identityN/A
lower-fma.f64N/A
Applied rewrites82.9%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
Taylor expanded in i around 0
lower-+.f6413.3
Applied rewrites13.3%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6413.3
Applied rewrites13.3%
Taylor expanded in i around 0
*-commutativeN/A
lower-fma.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
lower-/.f6499.8
Applied rewrites99.8%
Final simplification87.4%
(FPCore (i n)
:precision binary64
(let* ((t_0 (pow (+ (/ i n) 1.0) n))
(t_1 (/ (- t_0 1.0) (/ i n)))
(t_2 (* (/ (- (* t_0 n) n) i) 100.0)))
(if (<= t_1 -5e-259)
t_2
(if (<= t_1 0.0)
(* (* 100.0 n) (/ (expm1 i) i))
(if (<= t_1 INFINITY)
t_2
(/ 100.0 (fma (- (/ 0.5 (* n n)) (/ 0.5 n)) i (/ 1.0 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 t_2 = (((t_0 * n) - n) / i) * 100.0;
double tmp;
if (t_1 <= -5e-259) {
tmp = t_2;
} else if (t_1 <= 0.0) {
tmp = (100.0 * n) * (expm1(i) / i);
} else if (t_1 <= ((double) INFINITY)) {
tmp = t_2;
} else {
tmp = 100.0 / fma(((0.5 / (n * n)) - (0.5 / n)), i, (1.0 / 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)) t_2 = Float64(Float64(Float64(Float64(t_0 * n) - n) / i) * 100.0) tmp = 0.0 if (t_1 <= -5e-259) tmp = t_2; elseif (t_1 <= 0.0) tmp = Float64(Float64(100.0 * n) * Float64(expm1(i) / i)); elseif (t_1 <= Inf) tmp = t_2; else tmp = Float64(100.0 / fma(Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), i, Float64(1.0 / 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]}, Block[{t$95$2 = N[(N[(N[(N[(t$95$0 * n), $MachinePrecision] - n), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision]}, If[LessEqual[t$95$1, -5e-259], t$95$2, If[LessEqual[t$95$1, 0.0], N[(N[(100.0 * n), $MachinePrecision] * N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], t$95$2, N[(100.0 / N[(N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * i + N[(1.0 / 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}}\\
t_2 := \frac{t\_0 \cdot n - n}{i} \cdot 100\\
\mathbf{if}\;t\_1 \leq -5 \cdot 10^{-259}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq 0:\\
\;\;\;\;\left(100 \cdot n\right) \cdot \frac{\mathsf{expm1}\left(i\right)}{i}\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;\frac{100}{\mathsf{fma}\left(\frac{0.5}{n \cdot n} - \frac{0.5}{n}, i, \frac{1}{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)) < -4.99999999999999977e-259 or -0.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < +inf.0Initial program 93.0%
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
lower--.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f6493.0
Applied rewrites93.0%
lift--.f64N/A
lift-*.f64N/A
lift-/.f64N/A
associate-*l/N/A
lift-/.f64N/A
sub-divN/A
lower-/.f64N/A
lower--.f64N/A
lift-pow.f64N/A
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
lower-*.f64N/A
lift-/.f64N/A
+-commutativeN/A
lift-+.f64N/A
lift-pow.f6493.3
Applied rewrites93.3%
if -4.99999999999999977e-259 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < -0.0Initial program 14.6%
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.f6482.8
Applied rewrites82.8%
Applied rewrites82.8%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
Taylor expanded in i around 0
lower-+.f6413.3
Applied rewrites13.3%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6413.3
Applied rewrites13.3%
Taylor expanded in i around 0
*-commutativeN/A
lower-fma.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
lower-/.f6499.8
Applied rewrites99.8%
Final simplification87.3%
(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 -10000.0)
(fma (* (/ t_0 i) n) 100.0 (/ -1.0 (/ i (* 100.0 n))))
(if (<= t_1 INFINITY)
(* (/ (* (expm1 (* (log1p (/ i n)) n)) 100.0) i) n)
(/ 100.0 (fma (- (/ 0.5 (* n n)) (/ 0.5 n)) i (/ 1.0 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 <= -10000.0) {
tmp = fma(((t_0 / i) * n), 100.0, (-1.0 / (i / (100.0 * n))));
} else if (t_1 <= ((double) INFINITY)) {
tmp = ((expm1((log1p((i / n)) * n)) * 100.0) / i) * n;
} else {
tmp = 100.0 / fma(((0.5 / (n * n)) - (0.5 / n)), i, (1.0 / 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 <= -10000.0) tmp = fma(Float64(Float64(t_0 / i) * n), 100.0, Float64(-1.0 / Float64(i / Float64(100.0 * n)))); elseif (t_1 <= Inf) tmp = Float64(Float64(Float64(expm1(Float64(log1p(Float64(i / n)) * n)) * 100.0) / i) * n); else tmp = Float64(100.0 / fma(Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), i, Float64(1.0 / 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, -10000.0], N[(N[(N[(t$95$0 / i), $MachinePrecision] * n), $MachinePrecision] * 100.0 + N[(-1.0 / N[(i / N[(100.0 * n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(N[(N[(N[(Exp[N[(N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision] * n), $MachinePrecision]] - 1), $MachinePrecision] * 100.0), $MachinePrecision] / i), $MachinePrecision] * n), $MachinePrecision], N[(100.0 / N[(N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * i + N[(1.0 / 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 -10000:\\
\;\;\;\;\mathsf{fma}\left(\frac{t\_0}{i} \cdot n, 100, \frac{-1}{\frac{i}{100 \cdot n}}\right)\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;\frac{\mathsf{expm1}\left(\mathsf{log1p}\left(\frac{i}{n}\right) \cdot n\right) \cdot 100}{i} \cdot n\\
\mathbf{else}:\\
\;\;\;\;\frac{100}{\mathsf{fma}\left(\frac{0.5}{n \cdot n} - \frac{0.5}{n}, i, \frac{1}{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)) < -1e4Initial program 99.8%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
distribute-rgt-inN/A
lower-fma.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-*.f64N/A
Applied rewrites99.8%
lift-*.f64N/A
lift-/.f64N/A
associate-*l/N/A
clear-numN/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f64100.0
Applied rewrites100.0%
if -1e4 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < +inf.0Initial program 26.3%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
associate-/r/N/A
lift-/.f64N/A
clear-numN/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6425.9
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6495.8
Applied rewrites95.8%
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
lift-expm1.f64N/A
lift-*.f64N/A
lift-log1p.f64N/A
pow-to-expN/A
lift-/.f64N/A
*-commutativeN/A
lift-/.f64N/A
clear-numN/A
lift-/.f64N/A
div-invN/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
Applied rewrites96.7%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
Taylor expanded in i around 0
lower-+.f6413.3
Applied rewrites13.3%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6413.3
Applied rewrites13.3%
Taylor expanded in i around 0
*-commutativeN/A
lower-fma.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
lower-/.f6499.8
Applied rewrites99.8%
Final simplification97.3%
(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 -1e-124)
(* (/ (- (* t_0 n) n) i) 100.0)
(if (<= t_1 INFINITY)
(* (* 100.0 n) (/ (expm1 (* (log1p (/ i n)) n)) i))
(/ 100.0 (fma (- (/ 0.5 (* n n)) (/ 0.5 n)) i (/ 1.0 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 <= -1e-124) {
tmp = (((t_0 * n) - n) / i) * 100.0;
} else if (t_1 <= ((double) INFINITY)) {
tmp = (100.0 * n) * (expm1((log1p((i / n)) * n)) / i);
} else {
tmp = 100.0 / fma(((0.5 / (n * n)) - (0.5 / n)), i, (1.0 / 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 <= -1e-124) tmp = Float64(Float64(Float64(Float64(t_0 * n) - n) / i) * 100.0); elseif (t_1 <= Inf) tmp = Float64(Float64(100.0 * n) * Float64(expm1(Float64(log1p(Float64(i / n)) * n)) / i)); else tmp = Float64(100.0 / fma(Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), i, Float64(1.0 / 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, -1e-124], N[(N[(N[(N[(t$95$0 * n), $MachinePrecision] - n), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(N[(100.0 * n), $MachinePrecision] * N[(N[(Exp[N[(N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision] * n), $MachinePrecision]] - 1), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision], N[(100.0 / N[(N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * i + N[(1.0 / 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 -1 \cdot 10^{-124}:\\
\;\;\;\;\frac{t\_0 \cdot n - n}{i} \cdot 100\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;\left(100 \cdot n\right) \cdot \frac{\mathsf{expm1}\left(\mathsf{log1p}\left(\frac{i}{n}\right) \cdot n\right)}{i}\\
\mathbf{else}:\\
\;\;\;\;\frac{100}{\mathsf{fma}\left(\frac{0.5}{n \cdot n} - \frac{0.5}{n}, i, \frac{1}{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)) < -9.99999999999999933e-125Initial program 99.6%
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
lower--.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-/.f6499.9
Applied rewrites99.9%
lift--.f64N/A
lift-*.f64N/A
lift-/.f64N/A
associate-*l/N/A
lift-/.f64N/A
sub-divN/A
lower-/.f64N/A
lower--.f64N/A
lift-pow.f64N/A
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
lower-*.f64N/A
lift-/.f64N/A
+-commutativeN/A
lift-+.f64N/A
lift-pow.f6499.9
Applied rewrites99.9%
if -9.99999999999999933e-125 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < +inf.0Initial program 24.1%
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
associate-*l*N/A
lower-*.f64N/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
*-commutativeN/A
lower-*.f6496.6
Applied rewrites96.6%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
Taylor expanded in i around 0
lower-+.f6413.3
Applied rewrites13.3%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6413.3
Applied rewrites13.3%
Taylor expanded in i around 0
*-commutativeN/A
lower-fma.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
lower-/.f6499.8
Applied rewrites99.8%
Final simplification97.3%
(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 -10000.0)
(fma (* (/ t_0 i) n) 100.0 (/ -1.0 (/ i (* 100.0 n))))
(if (<= t_1 INFINITY)
(* (* (/ 100.0 i) (expm1 (* (log1p (/ i n)) n))) n)
(/ 100.0 (fma (- (/ 0.5 (* n n)) (/ 0.5 n)) i (/ 1.0 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 <= -10000.0) {
tmp = fma(((t_0 / i) * n), 100.0, (-1.0 / (i / (100.0 * n))));
} else if (t_1 <= ((double) INFINITY)) {
tmp = ((100.0 / i) * expm1((log1p((i / n)) * n))) * n;
} else {
tmp = 100.0 / fma(((0.5 / (n * n)) - (0.5 / n)), i, (1.0 / 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 <= -10000.0) tmp = fma(Float64(Float64(t_0 / i) * n), 100.0, Float64(-1.0 / Float64(i / Float64(100.0 * n)))); elseif (t_1 <= Inf) tmp = Float64(Float64(Float64(100.0 / i) * expm1(Float64(log1p(Float64(i / n)) * n))) * n); else tmp = Float64(100.0 / fma(Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), i, Float64(1.0 / 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, -10000.0], N[(N[(N[(t$95$0 / i), $MachinePrecision] * n), $MachinePrecision] * 100.0 + N[(-1.0 / N[(i / N[(100.0 * n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], 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], N[(100.0 / N[(N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * i + N[(1.0 / 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 -10000:\\
\;\;\;\;\mathsf{fma}\left(\frac{t\_0}{i} \cdot n, 100, \frac{-1}{\frac{i}{100 \cdot n}}\right)\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;\left(\frac{100}{i} \cdot \mathsf{expm1}\left(\mathsf{log1p}\left(\frac{i}{n}\right) \cdot n\right)\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;\frac{100}{\mathsf{fma}\left(\frac{0.5}{n \cdot n} - \frac{0.5}{n}, i, \frac{1}{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)) < -1e4Initial program 99.8%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
distribute-rgt-inN/A
lower-fma.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-*.f64N/A
Applied rewrites99.8%
lift-*.f64N/A
lift-/.f64N/A
associate-*l/N/A
clear-numN/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f64100.0
Applied rewrites100.0%
if -1e4 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < +inf.0Initial program 26.3%
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-/.f6496.7
Applied rewrites96.7%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
Taylor expanded in i around 0
lower-+.f6413.3
Applied rewrites13.3%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6413.3
Applied rewrites13.3%
Taylor expanded in i around 0
*-commutativeN/A
lower-fma.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
lower-/.f6499.8
Applied rewrites99.8%
Final simplification97.3%
(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 -10000.0)
(fma (* (/ t_0 i) n) 100.0 (/ -1.0 (/ i (* 100.0 n))))
(if (<= t_1 INFINITY)
(* (* (/ 100.0 i) n) (expm1 (* (log1p (/ i n)) n)))
(/ 100.0 (fma (- (/ 0.5 (* n n)) (/ 0.5 n)) i (/ 1.0 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 <= -10000.0) {
tmp = fma(((t_0 / i) * n), 100.0, (-1.0 / (i / (100.0 * n))));
} else if (t_1 <= ((double) INFINITY)) {
tmp = ((100.0 / i) * n) * expm1((log1p((i / n)) * n));
} else {
tmp = 100.0 / fma(((0.5 / (n * n)) - (0.5 / n)), i, (1.0 / 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 <= -10000.0) tmp = fma(Float64(Float64(t_0 / i) * n), 100.0, Float64(-1.0 / Float64(i / Float64(100.0 * n)))); elseif (t_1 <= Inf) tmp = Float64(Float64(Float64(100.0 / i) * n) * expm1(Float64(log1p(Float64(i / n)) * n))); else tmp = Float64(100.0 / fma(Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), i, Float64(1.0 / 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, -10000.0], N[(N[(N[(t$95$0 / i), $MachinePrecision] * n), $MachinePrecision] * 100.0 + N[(-1.0 / N[(i / N[(100.0 * n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[(N[(N[(100.0 / i), $MachinePrecision] * n), $MachinePrecision] * N[(Exp[N[(N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision] * n), $MachinePrecision]] - 1), $MachinePrecision]), $MachinePrecision], N[(100.0 / N[(N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * i + N[(1.0 / 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 -10000:\\
\;\;\;\;\mathsf{fma}\left(\frac{t\_0}{i} \cdot n, 100, \frac{-1}{\frac{i}{100 \cdot n}}\right)\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;\left(\frac{100}{i} \cdot n\right) \cdot \mathsf{expm1}\left(\mathsf{log1p}\left(\frac{i}{n}\right) \cdot n\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{100}{\mathsf{fma}\left(\frac{0.5}{n \cdot n} - \frac{0.5}{n}, i, \frac{1}{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)) < -1e4Initial program 99.8%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
distribute-rgt-inN/A
lower-fma.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-*.f64N/A
Applied rewrites99.8%
lift-*.f64N/A
lift-/.f64N/A
associate-*l/N/A
clear-numN/A
lower-/.f64N/A
lower-/.f64N/A
lower-*.f64100.0
Applied rewrites100.0%
if -1e4 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < +inf.0Initial program 26.3%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
associate-/r/N/A
lift-/.f64N/A
clear-numN/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-/.f6425.9
lift--.f64N/A
lift-pow.f64N/A
pow-to-expN/A
lower-expm1.f64N/A
lower-*.f64N/A
lift-+.f64N/A
lower-log1p.f6495.8
Applied rewrites95.8%
lift-*.f64N/A
lift-/.f64N/A
associate-*l/N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6495.7
Applied rewrites95.7%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
Taylor expanded in i around 0
lower-+.f6413.3
Applied rewrites13.3%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6413.3
Applied rewrites13.3%
Taylor expanded in i around 0
*-commutativeN/A
lower-fma.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
lower-/.f6499.8
Applied rewrites99.8%
Final simplification96.6%
(FPCore (i n)
:precision binary64
(if (<= n -10200000.0)
(* (/ (* (expm1 i) n) i) 100.0)
(if (<= n 44.0)
(/ 100.0 (fma (- (/ 0.5 (* n n)) (/ 0.5 n)) i (/ 1.0 n)))
(* (* 100.0 n) (/ (expm1 i) i)))))
double code(double i, double n) {
double tmp;
if (n <= -10200000.0) {
tmp = ((expm1(i) * n) / i) * 100.0;
} else if (n <= 44.0) {
tmp = 100.0 / fma(((0.5 / (n * n)) - (0.5 / n)), i, (1.0 / n));
} else {
tmp = (100.0 * n) * (expm1(i) / i);
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -10200000.0) tmp = Float64(Float64(Float64(expm1(i) * n) / i) * 100.0); elseif (n <= 44.0) tmp = Float64(100.0 / fma(Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), i, Float64(1.0 / n))); else tmp = Float64(Float64(100.0 * n) * Float64(expm1(i) / i)); end return tmp end
code[i_, n_] := If[LessEqual[n, -10200000.0], N[(N[(N[(N[(Exp[i] - 1), $MachinePrecision] * n), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision], If[LessEqual[n, 44.0], N[(100.0 / N[(N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * i + N[(1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(100.0 * n), $MachinePrecision] * N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -10200000:\\
\;\;\;\;\frac{\mathsf{expm1}\left(i\right) \cdot n}{i} \cdot 100\\
\mathbf{elif}\;n \leq 44:\\
\;\;\;\;\frac{100}{\mathsf{fma}\left(\frac{0.5}{n \cdot n} - \frac{0.5}{n}, i, \frac{1}{n}\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(100 \cdot n\right) \cdot \frac{\mathsf{expm1}\left(i\right)}{i}\\
\end{array}
\end{array}
if n < -1.02e7Initial program 19.9%
Taylor expanded in n around inf
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-expm1.f6492.6
Applied rewrites92.6%
if -1.02e7 < n < 44Initial program 27.0%
Taylor expanded in i around 0
lower-+.f6426.1
Applied rewrites26.1%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6426.1
Applied rewrites26.1%
Taylor expanded in i around 0
*-commutativeN/A
lower-fma.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
lower-/.f6470.5
Applied rewrites70.5%
if 44 < n Initial program 23.1%
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.f6493.4
Applied rewrites93.4%
Applied rewrites93.4%
Final simplification84.7%
(FPCore (i n)
:precision binary64
(let* ((t_0 (/ (expm1 i) i)))
(if (<= n -10200000.0)
(* (* t_0 100.0) n)
(if (<= n 44.0)
(/ 100.0 (fma (- (/ 0.5 (* n n)) (/ 0.5 n)) i (/ 1.0 n)))
(* (* 100.0 n) t_0)))))
double code(double i, double n) {
double t_0 = expm1(i) / i;
double tmp;
if (n <= -10200000.0) {
tmp = (t_0 * 100.0) * n;
} else if (n <= 44.0) {
tmp = 100.0 / fma(((0.5 / (n * n)) - (0.5 / n)), i, (1.0 / n));
} else {
tmp = (100.0 * n) * t_0;
}
return tmp;
}
function code(i, n) t_0 = Float64(expm1(i) / i) tmp = 0.0 if (n <= -10200000.0) tmp = Float64(Float64(t_0 * 100.0) * n); elseif (n <= 44.0) tmp = Float64(100.0 / fma(Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), i, Float64(1.0 / n))); else tmp = Float64(Float64(100.0 * n) * t_0); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision]}, If[LessEqual[n, -10200000.0], N[(N[(t$95$0 * 100.0), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[n, 44.0], N[(100.0 / N[(N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * i + N[(1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(100.0 * n), $MachinePrecision] * t$95$0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{\mathsf{expm1}\left(i\right)}{i}\\
\mathbf{if}\;n \leq -10200000:\\
\;\;\;\;\left(t\_0 \cdot 100\right) \cdot n\\
\mathbf{elif}\;n \leq 44:\\
\;\;\;\;\frac{100}{\mathsf{fma}\left(\frac{0.5}{n \cdot n} - \frac{0.5}{n}, i, \frac{1}{n}\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(100 \cdot n\right) \cdot t\_0\\
\end{array}
\end{array}
if n < -1.02e7Initial program 19.9%
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.f6492.6
Applied rewrites92.6%
if -1.02e7 < n < 44Initial program 27.0%
Taylor expanded in i around 0
lower-+.f6426.1
Applied rewrites26.1%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6426.1
Applied rewrites26.1%
Taylor expanded in i around 0
*-commutativeN/A
lower-fma.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
lower-/.f6470.5
Applied rewrites70.5%
if 44 < n Initial program 23.1%
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.f6493.4
Applied rewrites93.4%
Applied rewrites93.4%
Final simplification84.7%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* (* (/ (expm1 i) i) 100.0) n)))
(if (<= n -10200000.0)
t_0
(if (<= n 44.0)
(/ 100.0 (fma (- (/ 0.5 (* n n)) (/ 0.5 n)) i (/ 1.0 n)))
t_0))))
double code(double i, double n) {
double t_0 = ((expm1(i) / i) * 100.0) * n;
double tmp;
if (n <= -10200000.0) {
tmp = t_0;
} else if (n <= 44.0) {
tmp = 100.0 / fma(((0.5 / (n * n)) - (0.5 / n)), i, (1.0 / 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 <= -10200000.0) tmp = t_0; elseif (n <= 44.0) tmp = Float64(100.0 / fma(Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), i, Float64(1.0 / 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, -10200000.0], t$95$0, If[LessEqual[n, 44.0], N[(100.0 / N[(N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * i + N[(1.0 / 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 -10200000:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 44:\\
\;\;\;\;\frac{100}{\mathsf{fma}\left(\frac{0.5}{n \cdot n} - \frac{0.5}{n}, i, \frac{1}{n}\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -1.02e7 or 44 < n Initial program 21.7%
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.f6493.1
Applied rewrites93.1%
if -1.02e7 < n < 44Initial program 27.0%
Taylor expanded in i around 0
lower-+.f6426.1
Applied rewrites26.1%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6426.1
Applied rewrites26.1%
Taylor expanded in i around 0
*-commutativeN/A
lower-fma.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
lower-/.f6470.5
Applied rewrites70.5%
(FPCore (i n)
:precision binary64
(if (<= n -3.3e+138)
(* (fma (fma 16.666666666666668 i 50.0) i 100.0) n)
(if (<= n 44.0)
(/ 100.0 (fma (- (/ 0.5 (* n n)) (/ 0.5 n)) i (/ 1.0 n)))
(*
(/
(*
(fma
(fma (fma 0.041666666666666664 i 0.16666666666666666) i 0.5)
i
1.0)
i)
i)
(* 100.0 n)))))
double code(double i, double n) {
double tmp;
if (n <= -3.3e+138) {
tmp = fma(fma(16.666666666666668, i, 50.0), i, 100.0) * n;
} else if (n <= 44.0) {
tmp = 100.0 / fma(((0.5 / (n * n)) - (0.5 / n)), i, (1.0 / n));
} else {
tmp = ((fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * i) / i) * (100.0 * n);
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -3.3e+138) tmp = Float64(fma(fma(16.666666666666668, i, 50.0), i, 100.0) * n); elseif (n <= 44.0) tmp = Float64(100.0 / fma(Float64(Float64(0.5 / Float64(n * n)) - Float64(0.5 / n)), i, Float64(1.0 / n))); else tmp = Float64(Float64(Float64(fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * i) / i) * Float64(100.0 * n)); end return tmp end
code[i_, n_] := If[LessEqual[n, -3.3e+138], N[(N[(N[(16.666666666666668 * i + 50.0), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[n, 44.0], N[(100.0 / N[(N[(N[(0.5 / N[(n * n), $MachinePrecision]), $MachinePrecision] - N[(0.5 / n), $MachinePrecision]), $MachinePrecision] * i + N[(1.0 / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(N[(N[(0.041666666666666664 * i + 0.16666666666666666), $MachinePrecision] * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * i), $MachinePrecision] / i), $MachinePrecision] * N[(100.0 * n), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -3.3 \cdot 10^{+138}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(16.666666666666668, i, 50\right), i, 100\right) \cdot n\\
\mathbf{elif}\;n \leq 44:\\
\;\;\;\;\frac{100}{\mathsf{fma}\left(\frac{0.5}{n \cdot n} - \frac{0.5}{n}, i, \frac{1}{n}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, i, 0.16666666666666666\right), i, 0.5\right), i, 1\right) \cdot i}{i} \cdot \left(100 \cdot n\right)\\
\end{array}
\end{array}
if n < -3.29999999999999978e138Initial program 11.8%
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.f6495.6
Applied rewrites95.6%
Taylor expanded in i around 0
Applied rewrites77.0%
if -3.29999999999999978e138 < n < 44Initial program 28.8%
Taylor expanded in i around 0
lower-+.f6422.5
Applied rewrites22.5%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6422.5
Applied rewrites22.5%
Taylor expanded in i around 0
*-commutativeN/A
lower-fma.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
lower-/.f6469.7
Applied rewrites69.7%
if 44 < n Initial program 23.1%
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.f6493.4
Applied rewrites93.4%
Applied rewrites93.4%
Taylor expanded in i around 0
Applied rewrites83.8%
Final simplification76.1%
(FPCore (i n)
:precision binary64
(if (<= n -3.7e-150)
(* (fma (* n i) (fma 0.16666666666666666 i 0.5) n) 100.0)
(if (<= n 1.5e-213)
(* (/ (- 1.0 1.0) (/ i n)) 100.0)
(*
(/
(*
(fma
(fma (fma 0.041666666666666664 i 0.16666666666666666) i 0.5)
i
1.0)
i)
i)
(* 100.0 n)))))
double code(double i, double n) {
double tmp;
if (n <= -3.7e-150) {
tmp = fma((n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0;
} else if (n <= 1.5e-213) {
tmp = ((1.0 - 1.0) / (i / n)) * 100.0;
} else {
tmp = ((fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * i) / i) * (100.0 * n);
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -3.7e-150) tmp = Float64(fma(Float64(n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0); elseif (n <= 1.5e-213) tmp = Float64(Float64(Float64(1.0 - 1.0) / Float64(i / n)) * 100.0); else tmp = Float64(Float64(Float64(fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * i) / i) * Float64(100.0 * n)); end return tmp end
code[i_, n_] := If[LessEqual[n, -3.7e-150], N[(N[(N[(n * i), $MachinePrecision] * N[(0.16666666666666666 * i + 0.5), $MachinePrecision] + n), $MachinePrecision] * 100.0), $MachinePrecision], If[LessEqual[n, 1.5e-213], N[(N[(N[(1.0 - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision] * 100.0), $MachinePrecision], N[(N[(N[(N[(N[(N[(0.041666666666666664 * i + 0.16666666666666666), $MachinePrecision] * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * i), $MachinePrecision] / i), $MachinePrecision] * N[(100.0 * n), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -3.7 \cdot 10^{-150}:\\
\;\;\;\;\mathsf{fma}\left(n \cdot i, \mathsf{fma}\left(0.16666666666666666, i, 0.5\right), n\right) \cdot 100\\
\mathbf{elif}\;n \leq 1.5 \cdot 10^{-213}:\\
\;\;\;\;\frac{1 - 1}{\frac{i}{n}} \cdot 100\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, i, 0.16666666666666666\right), i, 0.5\right), i, 1\right) \cdot i}{i} \cdot \left(100 \cdot n\right)\\
\end{array}
\end{array}
if n < -3.70000000000000001e-150Initial program 18.4%
Taylor expanded in i around 0
+-commutativeN/A
Applied rewrites72.3%
Taylor expanded in n around inf
Applied rewrites72.2%
if -3.70000000000000001e-150 < n < 1.49999999999999993e-213Initial program 59.3%
Taylor expanded in i around 0
Applied rewrites66.7%
if 1.49999999999999993e-213 < n Initial program 20.3%
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.f6483.4
Applied rewrites83.4%
Applied rewrites83.4%
Taylor expanded in i around 0
Applied rewrites76.7%
Final simplification73.9%
(FPCore (i n)
:precision binary64
(if (<= n -3.7e-150)
(* (fma (* n i) (fma 0.16666666666666666 i 0.5) n) 100.0)
(if (<= n 1.5e-213)
(* (/ (- 1.0 1.0) (/ i n)) 100.0)
(*
(fma (fma (fma 4.166666666666667 i 16.666666666666668) i 50.0) i 100.0)
n))))
double code(double i, double n) {
double tmp;
if (n <= -3.7e-150) {
tmp = fma((n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0;
} else if (n <= 1.5e-213) {
tmp = ((1.0 - 1.0) / (i / n)) * 100.0;
} else {
tmp = fma(fma(fma(4.166666666666667, i, 16.666666666666668), i, 50.0), i, 100.0) * n;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -3.7e-150) tmp = Float64(fma(Float64(n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0); elseif (n <= 1.5e-213) tmp = Float64(Float64(Float64(1.0 - 1.0) / Float64(i / n)) * 100.0); else tmp = Float64(fma(fma(fma(4.166666666666667, i, 16.666666666666668), i, 50.0), i, 100.0) * n); end return tmp end
code[i_, n_] := If[LessEqual[n, -3.7e-150], N[(N[(N[(n * i), $MachinePrecision] * N[(0.16666666666666666 * i + 0.5), $MachinePrecision] + n), $MachinePrecision] * 100.0), $MachinePrecision], If[LessEqual[n, 1.5e-213], N[(N[(N[(1.0 - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision] * 100.0), $MachinePrecision], N[(N[(N[(N[(4.166666666666667 * i + 16.666666666666668), $MachinePrecision] * i + 50.0), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -3.7 \cdot 10^{-150}:\\
\;\;\;\;\mathsf{fma}\left(n \cdot i, \mathsf{fma}\left(0.16666666666666666, i, 0.5\right), n\right) \cdot 100\\
\mathbf{elif}\;n \leq 1.5 \cdot 10^{-213}:\\
\;\;\;\;\frac{1 - 1}{\frac{i}{n}} \cdot 100\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4.166666666666667, i, 16.666666666666668\right), i, 50\right), i, 100\right) \cdot n\\
\end{array}
\end{array}
if n < -3.70000000000000001e-150Initial program 18.4%
Taylor expanded in i around 0
+-commutativeN/A
Applied rewrites72.3%
Taylor expanded in n around inf
Applied rewrites72.2%
if -3.70000000000000001e-150 < n < 1.49999999999999993e-213Initial program 59.3%
Taylor expanded in i around 0
Applied rewrites66.7%
if 1.49999999999999993e-213 < n Initial program 20.3%
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.f6483.4
Applied rewrites83.4%
Taylor expanded in i around 0
Applied rewrites76.0%
Final simplification73.6%
(FPCore (i n)
:precision binary64
(if (<= n -3.7e-150)
(* (fma (* n i) (fma 0.16666666666666666 i 0.5) n) 100.0)
(if (<= n 1.5e-213)
(/ 0.0 i)
(*
(fma (fma (fma 4.166666666666667 i 16.666666666666668) i 50.0) i 100.0)
n))))
double code(double i, double n) {
double tmp;
if (n <= -3.7e-150) {
tmp = fma((n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0;
} else if (n <= 1.5e-213) {
tmp = 0.0 / i;
} else {
tmp = fma(fma(fma(4.166666666666667, i, 16.666666666666668), i, 50.0), i, 100.0) * n;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -3.7e-150) tmp = Float64(fma(Float64(n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0); elseif (n <= 1.5e-213) tmp = Float64(0.0 / i); else tmp = Float64(fma(fma(fma(4.166666666666667, i, 16.666666666666668), i, 50.0), i, 100.0) * n); end return tmp end
code[i_, n_] := If[LessEqual[n, -3.7e-150], N[(N[(N[(n * i), $MachinePrecision] * N[(0.16666666666666666 * i + 0.5), $MachinePrecision] + n), $MachinePrecision] * 100.0), $MachinePrecision], If[LessEqual[n, 1.5e-213], N[(0.0 / i), $MachinePrecision], N[(N[(N[(N[(4.166666666666667 * i + 16.666666666666668), $MachinePrecision] * i + 50.0), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -3.7 \cdot 10^{-150}:\\
\;\;\;\;\mathsf{fma}\left(n \cdot i, \mathsf{fma}\left(0.16666666666666666, i, 0.5\right), n\right) \cdot 100\\
\mathbf{elif}\;n \leq 1.5 \cdot 10^{-213}:\\
\;\;\;\;\frac{0}{i}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4.166666666666667, i, 16.666666666666668\right), i, 50\right), i, 100\right) \cdot n\\
\end{array}
\end{array}
if n < -3.70000000000000001e-150Initial program 18.4%
Taylor expanded in i around 0
+-commutativeN/A
Applied rewrites72.3%
Taylor expanded in n around inf
Applied rewrites72.2%
if -3.70000000000000001e-150 < n < 1.49999999999999993e-213Initial program 59.3%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
distribute-rgt-inN/A
lower-fma.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-*.f64N/A
Applied rewrites18.8%
Taylor expanded in i around 0
distribute-rgt-outN/A
metadata-evalN/A
mul0-rgtN/A
lower-/.f6466.7
Applied rewrites66.7%
if 1.49999999999999993e-213 < n Initial program 20.3%
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.f6483.4
Applied rewrites83.4%
Taylor expanded in i around 0
Applied rewrites76.0%
Final simplification73.6%
(FPCore (i n) :precision binary64 (let* ((t_0 (* (fma (* n i) (fma 0.16666666666666666 i 0.5) n) 100.0))) (if (<= n -3.7e-150) t_0 (if (<= n 1.5e-213) (/ 0.0 i) t_0))))
double code(double i, double n) {
double t_0 = fma((n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0;
double tmp;
if (n <= -3.7e-150) {
tmp = t_0;
} else if (n <= 1.5e-213) {
tmp = 0.0 / i;
} else {
tmp = t_0;
}
return tmp;
}
function code(i, n) t_0 = Float64(fma(Float64(n * i), fma(0.16666666666666666, i, 0.5), n) * 100.0) tmp = 0.0 if (n <= -3.7e-150) tmp = t_0; elseif (n <= 1.5e-213) tmp = Float64(0.0 / i); else tmp = t_0; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(N[(n * i), $MachinePrecision] * N[(0.16666666666666666 * i + 0.5), $MachinePrecision] + n), $MachinePrecision] * 100.0), $MachinePrecision]}, If[LessEqual[n, -3.7e-150], t$95$0, If[LessEqual[n, 1.5e-213], N[(0.0 / i), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(n \cdot i, \mathsf{fma}\left(0.16666666666666666, i, 0.5\right), n\right) \cdot 100\\
\mathbf{if}\;n \leq -3.7 \cdot 10^{-150}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 1.5 \cdot 10^{-213}:\\
\;\;\;\;\frac{0}{i}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -3.70000000000000001e-150 or 1.49999999999999993e-213 < n Initial program 19.5%
Taylor expanded in i around 0
+-commutativeN/A
Applied rewrites69.9%
Taylor expanded in n around inf
Applied rewrites72.0%
if -3.70000000000000001e-150 < n < 1.49999999999999993e-213Initial program 59.3%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
distribute-rgt-inN/A
lower-fma.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-*.f64N/A
Applied rewrites18.8%
Taylor expanded in i around 0
distribute-rgt-outN/A
metadata-evalN/A
mul0-rgtN/A
lower-/.f6466.7
Applied rewrites66.7%
Final simplification71.5%
(FPCore (i n) :precision binary64 (let* ((t_0 (* (fma (fma 16.666666666666668 i 50.0) i 100.0) n))) (if (<= n -3.7e-150) t_0 (if (<= n 1.5e-213) (/ 0.0 i) t_0))))
double code(double i, double n) {
double t_0 = fma(fma(16.666666666666668, i, 50.0), i, 100.0) * n;
double tmp;
if (n <= -3.7e-150) {
tmp = t_0;
} else if (n <= 1.5e-213) {
tmp = 0.0 / i;
} else {
tmp = t_0;
}
return tmp;
}
function code(i, n) t_0 = Float64(fma(fma(16.666666666666668, i, 50.0), i, 100.0) * n) tmp = 0.0 if (n <= -3.7e-150) tmp = t_0; elseif (n <= 1.5e-213) tmp = Float64(0.0 / i); else tmp = t_0; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(N[(16.666666666666668 * i + 50.0), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision]}, If[LessEqual[n, -3.7e-150], t$95$0, If[LessEqual[n, 1.5e-213], N[(0.0 / i), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(\mathsf{fma}\left(16.666666666666668, i, 50\right), i, 100\right) \cdot n\\
\mathbf{if}\;n \leq -3.7 \cdot 10^{-150}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 1.5 \cdot 10^{-213}:\\
\;\;\;\;\frac{0}{i}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -3.70000000000000001e-150 or 1.49999999999999993e-213 < n Initial program 19.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.f6485.9
Applied rewrites85.9%
Taylor expanded in i around 0
Applied rewrites72.0%
if -3.70000000000000001e-150 < n < 1.49999999999999993e-213Initial program 59.3%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
distribute-rgt-inN/A
lower-fma.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-*.f64N/A
Applied rewrites18.8%
Taylor expanded in i around 0
distribute-rgt-outN/A
metadata-evalN/A
mul0-rgtN/A
lower-/.f6466.7
Applied rewrites66.7%
(FPCore (i n) :precision binary64 (let* ((t_0 (* (fma 0.5 i 1.0) (* 100.0 n)))) (if (<= n -3.7e-150) t_0 (if (<= n 2.8e-217) (/ 0.0 i) t_0))))
double code(double i, double n) {
double t_0 = fma(0.5, i, 1.0) * (100.0 * n);
double tmp;
if (n <= -3.7e-150) {
tmp = t_0;
} else if (n <= 2.8e-217) {
tmp = 0.0 / i;
} else {
tmp = t_0;
}
return tmp;
}
function code(i, n) t_0 = Float64(fma(0.5, i, 1.0) * Float64(100.0 * n)) tmp = 0.0 if (n <= -3.7e-150) tmp = t_0; elseif (n <= 2.8e-217) tmp = Float64(0.0 / i); else tmp = t_0; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(0.5 * i + 1.0), $MachinePrecision] * N[(100.0 * n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -3.7e-150], t$95$0, If[LessEqual[n, 2.8e-217], N[(0.0 / i), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(0.5, i, 1\right) \cdot \left(100 \cdot n\right)\\
\mathbf{if}\;n \leq -3.7 \cdot 10^{-150}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 2.8 \cdot 10^{-217}:\\
\;\;\;\;\frac{0}{i}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -3.70000000000000001e-150 or 2.8e-217 < n Initial program 19.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.f6485.9
Applied rewrites85.9%
Applied rewrites85.9%
Taylor expanded in i around 0
Applied rewrites68.7%
if -3.70000000000000001e-150 < n < 2.8e-217Initial program 59.3%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
distribute-rgt-inN/A
lower-fma.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-*.f64N/A
Applied rewrites18.8%
Taylor expanded in i around 0
distribute-rgt-outN/A
metadata-evalN/A
mul0-rgtN/A
lower-/.f6466.7
Applied rewrites66.7%
Final simplification68.5%
(FPCore (i n) :precision binary64 (let* ((t_0 (* (fma 50.0 i 100.0) n))) (if (<= n -3.7e-150) t_0 (if (<= n 2.8e-217) (/ 0.0 i) t_0))))
double code(double i, double n) {
double t_0 = fma(50.0, i, 100.0) * n;
double tmp;
if (n <= -3.7e-150) {
tmp = t_0;
} else if (n <= 2.8e-217) {
tmp = 0.0 / i;
} else {
tmp = t_0;
}
return tmp;
}
function code(i, n) t_0 = Float64(fma(50.0, i, 100.0) * n) tmp = 0.0 if (n <= -3.7e-150) tmp = t_0; elseif (n <= 2.8e-217) tmp = Float64(0.0 / i); else tmp = t_0; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(50.0 * i + 100.0), $MachinePrecision] * n), $MachinePrecision]}, If[LessEqual[n, -3.7e-150], t$95$0, If[LessEqual[n, 2.8e-217], N[(0.0 / i), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(50, i, 100\right) \cdot n\\
\mathbf{if}\;n \leq -3.7 \cdot 10^{-150}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 2.8 \cdot 10^{-217}:\\
\;\;\;\;\frac{0}{i}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -3.70000000000000001e-150 or 2.8e-217 < n Initial program 19.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.f6485.9
Applied rewrites85.9%
Taylor expanded in i around 0
Applied rewrites68.7%
if -3.70000000000000001e-150 < n < 2.8e-217Initial program 59.3%
lift-*.f64N/A
lift-/.f64N/A
lift--.f64N/A
div-subN/A
lift-/.f64N/A
clear-numN/A
sub-negN/A
distribute-rgt-inN/A
lower-fma.f64N/A
lift-/.f64N/A
associate-/r/N/A
lower-*.f64N/A
lower-/.f64N/A
lift-+.f64N/A
+-commutativeN/A
lower-+.f64N/A
lower-*.f64N/A
Applied rewrites18.8%
Taylor expanded in i around 0
distribute-rgt-outN/A
metadata-evalN/A
mul0-rgtN/A
lower-/.f6466.7
Applied rewrites66.7%
(FPCore (i n) :precision binary64 (if (<= i 6.3e+49) (* 100.0 n) (* (* 50.0 i) n)))
double code(double i, double n) {
double tmp;
if (i <= 6.3e+49) {
tmp = 100.0 * n;
} else {
tmp = (50.0 * i) * n;
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if (i <= 6.3d+49) then
tmp = 100.0d0 * n
else
tmp = (50.0d0 * i) * n
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (i <= 6.3e+49) {
tmp = 100.0 * n;
} else {
tmp = (50.0 * i) * n;
}
return tmp;
}
def code(i, n): tmp = 0 if i <= 6.3e+49: tmp = 100.0 * n else: tmp = (50.0 * i) * n return tmp
function code(i, n) tmp = 0.0 if (i <= 6.3e+49) tmp = Float64(100.0 * n); else tmp = Float64(Float64(50.0 * i) * n); end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (i <= 6.3e+49) tmp = 100.0 * n; else tmp = (50.0 * i) * n; end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[i, 6.3e+49], N[(100.0 * n), $MachinePrecision], N[(N[(50.0 * i), $MachinePrecision] * n), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq 6.3 \cdot 10^{+49}:\\
\;\;\;\;100 \cdot n\\
\mathbf{else}:\\
\;\;\;\;\left(50 \cdot i\right) \cdot n\\
\end{array}
\end{array}
if i < 6.30000000000000007e49Initial program 18.8%
Taylor expanded in i around 0
lower-*.f6466.9
Applied rewrites66.9%
if 6.30000000000000007e49 < i Initial program 45.1%
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.f6463.1
Applied rewrites63.1%
Taylor expanded in i around 0
Applied rewrites42.4%
Taylor expanded in i around inf
Applied rewrites42.4%
(FPCore (i n) :precision binary64 (* (fma 50.0 i 100.0) n))
double code(double i, double n) {
return fma(50.0, i, 100.0) * n;
}
function code(i, n) return Float64(fma(50.0, i, 100.0) * n) end
code[i_, n_] := N[(N[(50.0 * i + 100.0), $MachinePrecision] * n), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(50, i, 100\right) \cdot n
\end{array}
Initial program 23.7%
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.f6480.1
Applied rewrites80.1%
Taylor expanded in i around 0
Applied rewrites63.3%
(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 23.7%
Taylor expanded in i around 0
lower-*.f6455.3
Applied rewrites55.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 2024331
(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))))