
(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 19 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) 1.0) (/ i n))) (t_1 (* 100.0 t_0)))
(if (<= t_0 -40000000000.0)
t_1
(if (<= t_0 5e-188)
(/ (* (/ (expm1 (* (log1p (/ i n)) n)) i) -100.0) (/ -1.0 n))
(if (<= t_0 INFINITY)
t_1
(* (/ 1.0 (fma (fma 0.0008333333333333334 i -0.005) i 0.01)) n))))))
double code(double i, double n) {
double t_0 = (pow(((i / n) + 1.0), n) - 1.0) / (i / n);
double t_1 = 100.0 * t_0;
double tmp;
if (t_0 <= -40000000000.0) {
tmp = t_1;
} else if (t_0 <= 5e-188) {
tmp = ((expm1((log1p((i / n)) * n)) / i) * -100.0) / (-1.0 / n);
} else if (t_0 <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = (1.0 / fma(fma(0.0008333333333333334, i, -0.005), i, 0.01)) * n;
}
return tmp;
}
function code(i, n) t_0 = Float64(Float64((Float64(Float64(i / n) + 1.0) ^ n) - 1.0) / Float64(i / n)) t_1 = Float64(100.0 * t_0) tmp = 0.0 if (t_0 <= -40000000000.0) tmp = t_1; elseif (t_0 <= 5e-188) tmp = Float64(Float64(Float64(expm1(Float64(log1p(Float64(i / n)) * n)) / i) * -100.0) / Float64(-1.0 / n)); elseif (t_0 <= Inf) tmp = t_1; else tmp = Float64(Float64(1.0 / fma(fma(0.0008333333333333334, i, -0.005), i, 0.01)) * n); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(N[Power[N[(N[(i / n), $MachinePrecision] + 1.0), $MachinePrecision], n], $MachinePrecision] - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(100.0 * t$95$0), $MachinePrecision]}, If[LessEqual[t$95$0, -40000000000.0], t$95$1, If[LessEqual[t$95$0, 5e-188], N[(N[(N[(N[(Exp[N[(N[Log[1 + N[(i / n), $MachinePrecision]], $MachinePrecision] * n), $MachinePrecision]] - 1), $MachinePrecision] / i), $MachinePrecision] * -100.0), $MachinePrecision] / N[(-1.0 / n), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, Infinity], t$95$1, N[(N[(1.0 / N[(N[(0.0008333333333333334 * i + -0.005), $MachinePrecision] * i + 0.01), $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{{\left(\frac{i}{n} + 1\right)}^{n} - 1}{\frac{i}{n}}\\
t_1 := 100 \cdot t\_0\\
\mathbf{if}\;t\_0 \leq -40000000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 5 \cdot 10^{-188}:\\
\;\;\;\;\frac{\frac{\mathsf{expm1}\left(\mathsf{log1p}\left(\frac{i}{n}\right) \cdot n\right)}{i} \cdot -100}{\frac{-1}{n}}\\
\mathbf{elif}\;t\_0 \leq \infty:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(\mathsf{fma}\left(0.0008333333333333334, i, -0.005\right), i, 0.01\right)} \cdot n\\
\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)) < -4e10 or 5.0000000000000001e-188 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < +inf.0Initial program 99.9%
if -4e10 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < 5.0000000000000001e-188Initial program 24.6%
lift-*.f64N/A
lift-/.f64N/A
associate-*r/N/A
lift-/.f64N/A
frac-2negN/A
div-invN/A
associate-/r*N/A
lower-/.f64N/A
Applied rewrites98.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 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.9
Applied rewrites83.9%
Applied rewrites83.7%
Taylor expanded in i around 0
Applied rewrites100.0%
Final simplification99.0%
(FPCore (i n)
:precision binary64
(let* ((t_0 (/ (- (pow (+ (/ i n) 1.0) n) 1.0) (/ i n))) (t_1 (* 100.0 t_0)))
(if (<= t_0 -40000000000.0)
t_1
(if (<= t_0 5e-188)
(* (* (/ 100.0 i) (expm1 (* (log1p (/ i n)) n))) n)
(if (<= t_0 INFINITY)
t_1
(* (/ 1.0 (fma (fma 0.0008333333333333334 i -0.005) i 0.01)) n))))))
double code(double i, double n) {
double t_0 = (pow(((i / n) + 1.0), n) - 1.0) / (i / n);
double t_1 = 100.0 * t_0;
double tmp;
if (t_0 <= -40000000000.0) {
tmp = t_1;
} else if (t_0 <= 5e-188) {
tmp = ((100.0 / i) * expm1((log1p((i / n)) * n))) * n;
} else if (t_0 <= ((double) INFINITY)) {
tmp = t_1;
} else {
tmp = (1.0 / fma(fma(0.0008333333333333334, i, -0.005), i, 0.01)) * n;
}
return tmp;
}
function code(i, n) t_0 = Float64(Float64((Float64(Float64(i / n) + 1.0) ^ n) - 1.0) / Float64(i / n)) t_1 = Float64(100.0 * t_0) tmp = 0.0 if (t_0 <= -40000000000.0) tmp = t_1; elseif (t_0 <= 5e-188) tmp = Float64(Float64(Float64(100.0 / i) * expm1(Float64(log1p(Float64(i / n)) * n))) * n); elseif (t_0 <= Inf) tmp = t_1; else tmp = Float64(Float64(1.0 / fma(fma(0.0008333333333333334, i, -0.005), i, 0.01)) * n); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(N[Power[N[(N[(i / n), $MachinePrecision] + 1.0), $MachinePrecision], n], $MachinePrecision] - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(100.0 * t$95$0), $MachinePrecision]}, If[LessEqual[t$95$0, -40000000000.0], t$95$1, If[LessEqual[t$95$0, 5e-188], 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$0, Infinity], t$95$1, N[(N[(1.0 / N[(N[(0.0008333333333333334 * i + -0.005), $MachinePrecision] * i + 0.01), $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{{\left(\frac{i}{n} + 1\right)}^{n} - 1}{\frac{i}{n}}\\
t_1 := 100 \cdot t\_0\\
\mathbf{if}\;t\_0 \leq -40000000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 5 \cdot 10^{-188}:\\
\;\;\;\;\left(\frac{100}{i} \cdot \mathsf{expm1}\left(\mathsf{log1p}\left(\frac{i}{n}\right) \cdot n\right)\right) \cdot n\\
\mathbf{elif}\;t\_0 \leq \infty:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(\mathsf{fma}\left(0.0008333333333333334, i, -0.005\right), i, 0.01\right)} \cdot n\\
\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)) < -4e10 or 5.0000000000000001e-188 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < +inf.0Initial program 99.9%
if -4e10 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < 5.0000000000000001e-188Initial program 24.6%
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.2
Applied rewrites98.2%
if +inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) Initial program 0.0%
Taylor expanded in 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.9
Applied rewrites83.9%
Applied rewrites83.7%
Taylor expanded in i around 0
Applied rewrites100.0%
Final simplification98.8%
(FPCore (i n)
:precision binary64
(let* ((t_0 (/ (- (pow (+ (/ i n) 1.0) n) 1.0) (/ i n))))
(if (<= t_0 (- INFINITY))
(*
(fma (fma (fma 4.166666666666667 i 16.666666666666668) i 50.0) i 100.0)
n)
(if (<= t_0 INFINITY)
(* (* (/ n i) 100.0) (expm1 i))
(* (/ 1.0 (fma (fma 0.0008333333333333334 i -0.005) i 0.01)) n)))))
double code(double i, double n) {
double t_0 = (pow(((i / n) + 1.0), n) - 1.0) / (i / n);
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = fma(fma(fma(4.166666666666667, i, 16.666666666666668), i, 50.0), i, 100.0) * n;
} else if (t_0 <= ((double) INFINITY)) {
tmp = ((n / i) * 100.0) * expm1(i);
} else {
tmp = (1.0 / fma(fma(0.0008333333333333334, i, -0.005), i, 0.01)) * n;
}
return tmp;
}
function code(i, n) t_0 = Float64(Float64((Float64(Float64(i / n) + 1.0) ^ n) - 1.0) / Float64(i / n)) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = Float64(fma(fma(fma(4.166666666666667, i, 16.666666666666668), i, 50.0), i, 100.0) * n); elseif (t_0 <= Inf) tmp = Float64(Float64(Float64(n / i) * 100.0) * expm1(i)); else tmp = Float64(Float64(1.0 / fma(fma(0.0008333333333333334, i, -0.005), i, 0.01)) * n); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(N[Power[N[(N[(i / n), $MachinePrecision] + 1.0), $MachinePrecision], n], $MachinePrecision] - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[(N[(N[(4.166666666666667 * i + 16.666666666666668), $MachinePrecision] * i + 50.0), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[t$95$0, Infinity], N[(N[(N[(n / i), $MachinePrecision] * 100.0), $MachinePrecision] * N[(Exp[i] - 1), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 / N[(N[(0.0008333333333333334 * i + -0.005), $MachinePrecision] * i + 0.01), $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{{\left(\frac{i}{n} + 1\right)}^{n} - 1}{\frac{i}{n}}\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4.166666666666667, i, 16.666666666666668\right), i, 50\right), i, 100\right) \cdot n\\
\mathbf{elif}\;t\_0 \leq \infty:\\
\;\;\;\;\left(\frac{n}{i} \cdot 100\right) \cdot \mathsf{expm1}\left(i\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(\mathsf{fma}\left(0.0008333333333333334, i, -0.005\right), i, 0.01\right)} \cdot n\\
\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)) < -inf.0Initial program 100.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.f6417.5
Applied rewrites17.5%
Taylor expanded in i around 0
Applied rewrites84.3%
if -inf.0 < (/.f64 (-.f64 (pow.f64 (+.f64 #s(literal 1 binary64) (/.f64 i n)) n) #s(literal 1 binary64)) (/.f64 i n)) < +inf.0Initial program 31.2%
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.f6476.6
Applied rewrites76.6%
Applied rewrites76.5%
Applied rewrites76.5%
Applied rewrites74.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%
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.9
Applied rewrites83.9%
Applied rewrites83.7%
Taylor expanded in i around 0
Applied rewrites100.0%
Final simplification80.2%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* (* (/ (expm1 i) i) 100.0) n)))
(if (<= n -1.35e-192)
t_0
(if (<= n -2e-310)
(* 100.0 (/ (- (pow (+ (/ i n) 1.0) n) 1.0) (/ i n)))
(if (<= n 2.6e-98)
(* (/ (* (- (log i) (log n)) n) (/ i n)) 100.0)
t_0)))))
double code(double i, double n) {
double t_0 = ((expm1(i) / i) * 100.0) * n;
double tmp;
if (n <= -1.35e-192) {
tmp = t_0;
} else if (n <= -2e-310) {
tmp = 100.0 * ((pow(((i / n) + 1.0), n) - 1.0) / (i / n));
} else if (n <= 2.6e-98) {
tmp = (((log(i) - log(n)) * n) / (i / n)) * 100.0;
} else {
tmp = t_0;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = ((Math.expm1(i) / i) * 100.0) * n;
double tmp;
if (n <= -1.35e-192) {
tmp = t_0;
} else if (n <= -2e-310) {
tmp = 100.0 * ((Math.pow(((i / n) + 1.0), n) - 1.0) / (i / n));
} else if (n <= 2.6e-98) {
tmp = (((Math.log(i) - Math.log(n)) * n) / (i / n)) * 100.0;
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = ((math.expm1(i) / i) * 100.0) * n tmp = 0 if n <= -1.35e-192: tmp = t_0 elif n <= -2e-310: tmp = 100.0 * ((math.pow(((i / n) + 1.0), n) - 1.0) / (i / n)) elif n <= 2.6e-98: tmp = (((math.log(i) - math.log(n)) * n) / (i / n)) * 100.0 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 <= -1.35e-192) tmp = t_0; elseif (n <= -2e-310) tmp = Float64(100.0 * Float64(Float64((Float64(Float64(i / n) + 1.0) ^ n) - 1.0) / Float64(i / n))); elseif (n <= 2.6e-98) tmp = Float64(Float64(Float64(Float64(log(i) - log(n)) * n) / Float64(i / n)) * 100.0); 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, -1.35e-192], t$95$0, If[LessEqual[n, -2e-310], N[(100.0 * N[(N[(N[Power[N[(N[(i / n), $MachinePrecision] + 1.0), $MachinePrecision], n], $MachinePrecision] - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 2.6e-98], N[(N[(N[(N[(N[Log[i], $MachinePrecision] - N[Log[n], $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision] * 100.0), $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 -1.35 \cdot 10^{-192}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq -2 \cdot 10^{-310}:\\
\;\;\;\;100 \cdot \frac{{\left(\frac{i}{n} + 1\right)}^{n} - 1}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 2.6 \cdot 10^{-98}:\\
\;\;\;\;\frac{\left(\log i - \log n\right) \cdot n}{\frac{i}{n}} \cdot 100\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -1.34999999999999996e-192 or 2.60000000000000013e-98 < n Initial program 21.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.f6486.7
Applied rewrites86.7%
if -1.34999999999999996e-192 < n < -1.999999999999994e-310Initial program 83.5%
if -1.999999999999994e-310 < n < 2.60000000000000013e-98Initial program 25.6%
Taylor expanded in n around 0
*-commutativeN/A
lower-*.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-log.f64N/A
lower-log.f6477.2
Applied rewrites77.2%
Final simplification85.4%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* (* (/ (expm1 i) i) 100.0) n)))
(if (<= n -1.35e-192)
t_0
(if (<= n -2e-310)
(* 100.0 (/ (- (pow (+ (/ i n) 1.0) n) 1.0) (/ i n)))
(if (<= n 2.6e-98)
(* (* (* (/ n i) n) (- (log i) (log n))) 100.0)
t_0)))))
double code(double i, double n) {
double t_0 = ((expm1(i) / i) * 100.0) * n;
double tmp;
if (n <= -1.35e-192) {
tmp = t_0;
} else if (n <= -2e-310) {
tmp = 100.0 * ((pow(((i / n) + 1.0), n) - 1.0) / (i / n));
} else if (n <= 2.6e-98) {
tmp = (((n / i) * n) * (log(i) - log(n))) * 100.0;
} else {
tmp = t_0;
}
return tmp;
}
public static double code(double i, double n) {
double t_0 = ((Math.expm1(i) / i) * 100.0) * n;
double tmp;
if (n <= -1.35e-192) {
tmp = t_0;
} else if (n <= -2e-310) {
tmp = 100.0 * ((Math.pow(((i / n) + 1.0), n) - 1.0) / (i / n));
} else if (n <= 2.6e-98) {
tmp = (((n / i) * n) * (Math.log(i) - Math.log(n))) * 100.0;
} else {
tmp = t_0;
}
return tmp;
}
def code(i, n): t_0 = ((math.expm1(i) / i) * 100.0) * n tmp = 0 if n <= -1.35e-192: tmp = t_0 elif n <= -2e-310: tmp = 100.0 * ((math.pow(((i / n) + 1.0), n) - 1.0) / (i / n)) elif n <= 2.6e-98: tmp = (((n / i) * n) * (math.log(i) - math.log(n))) * 100.0 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 <= -1.35e-192) tmp = t_0; elseif (n <= -2e-310) tmp = Float64(100.0 * Float64(Float64((Float64(Float64(i / n) + 1.0) ^ n) - 1.0) / Float64(i / n))); elseif (n <= 2.6e-98) tmp = Float64(Float64(Float64(Float64(n / i) * n) * Float64(log(i) - log(n))) * 100.0); 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, -1.35e-192], t$95$0, If[LessEqual[n, -2e-310], N[(100.0 * N[(N[(N[Power[N[(N[(i / n), $MachinePrecision] + 1.0), $MachinePrecision], n], $MachinePrecision] - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[n, 2.6e-98], N[(N[(N[(N[(n / i), $MachinePrecision] * n), $MachinePrecision] * N[(N[Log[i], $MachinePrecision] - N[Log[n], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 100.0), $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 -1.35 \cdot 10^{-192}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq -2 \cdot 10^{-310}:\\
\;\;\;\;100 \cdot \frac{{\left(\frac{i}{n} + 1\right)}^{n} - 1}{\frac{i}{n}}\\
\mathbf{elif}\;n \leq 2.6 \cdot 10^{-98}:\\
\;\;\;\;\left(\left(\frac{n}{i} \cdot n\right) \cdot \left(\log i - \log n\right)\right) \cdot 100\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -1.34999999999999996e-192 or 2.60000000000000013e-98 < n Initial program 21.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.f6486.7
Applied rewrites86.7%
if -1.34999999999999996e-192 < n < -1.999999999999994e-310Initial program 83.5%
if -1.999999999999994e-310 < n < 2.60000000000000013e-98Initial program 25.6%
Taylor expanded in n around 0
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
mul-1-negN/A
unsub-negN/A
lower--.f64N/A
lower-log.f64N/A
lower-log.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6463.2
Applied rewrites63.2%
Applied rewrites77.2%
Final simplification85.4%
(FPCore (i n)
:precision binary64
(if (<= n -1.35e-192)
(* (* (/ (expm1 i) i) 100.0) n)
(if (<= n 1.44e-251)
0.0
(if (<= n 5e-48)
(*
(/
1.0
(fma
(fma
(fma (* i i) -1.388888888888889e-5 0.0008333333333333334)
i
-0.005)
i
0.01))
n)
(* (/ (* (expm1 i) n) i) 100.0)))))
double code(double i, double n) {
double tmp;
if (n <= -1.35e-192) {
tmp = ((expm1(i) / i) * 100.0) * n;
} else if (n <= 1.44e-251) {
tmp = 0.0;
} else if (n <= 5e-48) {
tmp = (1.0 / fma(fma(fma((i * i), -1.388888888888889e-5, 0.0008333333333333334), i, -0.005), i, 0.01)) * n;
} else {
tmp = ((expm1(i) * n) / i) * 100.0;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (n <= -1.35e-192) tmp = Float64(Float64(Float64(expm1(i) / i) * 100.0) * n); elseif (n <= 1.44e-251) tmp = 0.0; elseif (n <= 5e-48) tmp = Float64(Float64(1.0 / fma(fma(fma(Float64(i * i), -1.388888888888889e-5, 0.0008333333333333334), i, -0.005), i, 0.01)) * n); else tmp = Float64(Float64(Float64(expm1(i) * n) / i) * 100.0); end return tmp end
code[i_, n_] := If[LessEqual[n, -1.35e-192], N[(N[(N[(N[(Exp[i] - 1), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[n, 1.44e-251], 0.0, If[LessEqual[n, 5e-48], N[(N[(1.0 / N[(N[(N[(N[(i * i), $MachinePrecision] * -1.388888888888889e-5 + 0.0008333333333333334), $MachinePrecision] * i + -0.005), $MachinePrecision] * i + 0.01), $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision], N[(N[(N[(N[(Exp[i] - 1), $MachinePrecision] * n), $MachinePrecision] / i), $MachinePrecision] * 100.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;n \leq -1.35 \cdot 10^{-192}:\\
\;\;\;\;\left(\frac{\mathsf{expm1}\left(i\right)}{i} \cdot 100\right) \cdot n\\
\mathbf{elif}\;n \leq 1.44 \cdot 10^{-251}:\\
\;\;\;\;0\\
\mathbf{elif}\;n \leq 5 \cdot 10^{-48}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(i \cdot i, -1.388888888888889 \cdot 10^{-5}, 0.0008333333333333334\right), i, -0.005\right), i, 0.01\right)} \cdot n\\
\mathbf{else}:\\
\;\;\;\;\frac{\mathsf{expm1}\left(i\right) \cdot n}{i} \cdot 100\\
\end{array}
\end{array}
if n < -1.34999999999999996e-192Initial program 25.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.f6484.3
Applied rewrites84.3%
if -1.34999999999999996e-192 < n < 1.44000000000000009e-251Initial program 78.1%
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 rewrites30.8%
Taylor expanded in i around 0
associate-*r/N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
lower-/.f6482.6
Applied rewrites82.6%
Taylor expanded in i around 0
Applied rewrites82.6%
if 1.44000000000000009e-251 < n < 4.9999999999999999e-48Initial program 15.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.f6435.7
Applied rewrites35.7%
Applied rewrites35.7%
Taylor expanded in i around 0
Applied rewrites68.9%
if 4.9999999999999999e-48 < n Initial program 17.6%
Taylor expanded in n around inf
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-expm1.f6492.1
Applied rewrites92.1%
Final simplification84.9%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* (* (/ (expm1 i) i) 100.0) n)))
(if (<= n -1.35e-192)
t_0
(if (<= n 1.44e-251)
0.0
(if (<= n 5e-48)
(*
(/
1.0
(fma
(fma
(fma (* i i) -1.388888888888889e-5 0.0008333333333333334)
i
-0.005)
i
0.01))
n)
t_0)))))
double code(double i, double n) {
double t_0 = ((expm1(i) / i) * 100.0) * n;
double tmp;
if (n <= -1.35e-192) {
tmp = t_0;
} else if (n <= 1.44e-251) {
tmp = 0.0;
} else if (n <= 5e-48) {
tmp = (1.0 / fma(fma(fma((i * i), -1.388888888888889e-5, 0.0008333333333333334), i, -0.005), i, 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 <= -1.35e-192) tmp = t_0; elseif (n <= 1.44e-251) tmp = 0.0; elseif (n <= 5e-48) tmp = Float64(Float64(1.0 / fma(fma(fma(Float64(i * i), -1.388888888888889e-5, 0.0008333333333333334), i, -0.005), i, 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, -1.35e-192], t$95$0, If[LessEqual[n, 1.44e-251], 0.0, If[LessEqual[n, 5e-48], N[(N[(1.0 / N[(N[(N[(N[(i * i), $MachinePrecision] * -1.388888888888889e-5 + 0.0008333333333333334), $MachinePrecision] * i + -0.005), $MachinePrecision] * i + 0.01), $MachinePrecision]), $MachinePrecision] * n), $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 -1.35 \cdot 10^{-192}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 1.44 \cdot 10^{-251}:\\
\;\;\;\;0\\
\mathbf{elif}\;n \leq 5 \cdot 10^{-48}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(i \cdot i, -1.388888888888889 \cdot 10^{-5}, 0.0008333333333333334\right), i, -0.005\right), i, 0.01\right)} \cdot n\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -1.34999999999999996e-192 or 4.9999999999999999e-48 < n Initial program 22.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.f6487.7
Applied rewrites87.7%
if -1.34999999999999996e-192 < n < 1.44000000000000009e-251Initial program 78.1%
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 rewrites30.8%
Taylor expanded in i around 0
associate-*r/N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
lower-/.f6482.6
Applied rewrites82.6%
Taylor expanded in i around 0
Applied rewrites82.6%
if 1.44000000000000009e-251 < n < 4.9999999999999999e-48Initial program 15.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.f6435.7
Applied rewrites35.7%
Applied rewrites35.7%
Taylor expanded in i around 0
Applied rewrites68.9%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* (/ 1.0 (fma (fma 0.0008333333333333334 i -0.005) i 0.01)) n)))
(if (<= n -2.7e+124)
(*
(*
(fma (fma (fma 0.041666666666666664 i 0.16666666666666666) i 0.5) i 1.0)
100.0)
n)
(if (<= n -1.26e-194)
t_0
(if (<= n 1.44e-251)
0.0
(if (<= n 5.1e-14)
t_0
(*
(fma
(fma (fma 4.166666666666667 i 16.666666666666668) i 50.0)
i
100.0)
n)))))))
double code(double i, double n) {
double t_0 = (1.0 / fma(fma(0.0008333333333333334, i, -0.005), i, 0.01)) * n;
double tmp;
if (n <= -2.7e+124) {
tmp = (fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * 100.0) * n;
} else if (n <= -1.26e-194) {
tmp = t_0;
} else if (n <= 1.44e-251) {
tmp = 0.0;
} else if (n <= 5.1e-14) {
tmp = t_0;
} else {
tmp = fma(fma(fma(4.166666666666667, i, 16.666666666666668), i, 50.0), i, 100.0) * n;
}
return tmp;
}
function code(i, n) t_0 = Float64(Float64(1.0 / fma(fma(0.0008333333333333334, i, -0.005), i, 0.01)) * n) tmp = 0.0 if (n <= -2.7e+124) tmp = Float64(Float64(fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * 100.0) * n); elseif (n <= -1.26e-194) tmp = t_0; elseif (n <= 1.44e-251) tmp = 0.0; elseif (n <= 5.1e-14) tmp = t_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_] := Block[{t$95$0 = N[(N[(1.0 / N[(N[(0.0008333333333333334 * i + -0.005), $MachinePrecision] * i + 0.01), $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision]}, If[LessEqual[n, -2.7e+124], N[(N[(N[(N[(N[(0.041666666666666664 * i + 0.16666666666666666), $MachinePrecision] * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[n, -1.26e-194], t$95$0, If[LessEqual[n, 1.44e-251], 0.0, If[LessEqual[n, 5.1e-14], t$95$0, 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}
t_0 := \frac{1}{\mathsf{fma}\left(\mathsf{fma}\left(0.0008333333333333334, i, -0.005\right), i, 0.01\right)} \cdot n\\
\mathbf{if}\;n \leq -2.7 \cdot 10^{+124}:\\
\;\;\;\;\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, i, 0.16666666666666666\right), i, 0.5\right), i, 1\right) \cdot 100\right) \cdot n\\
\mathbf{elif}\;n \leq -1.26 \cdot 10^{-194}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 1.44 \cdot 10^{-251}:\\
\;\;\;\;0\\
\mathbf{elif}\;n \leq 5.1 \cdot 10^{-14}:\\
\;\;\;\;t\_0\\
\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 < -2.69999999999999978e124Initial program 16.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.f6495.5
Applied rewrites95.5%
Taylor expanded in i around 0
Applied rewrites76.4%
if -2.69999999999999978e124 < n < -1.26e-194 or 1.44000000000000009e-251 < n < 5.0999999999999997e-14Initial program 25.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.f6463.7
Applied rewrites63.7%
Applied rewrites63.7%
Taylor expanded in i around 0
Applied rewrites66.2%
if -1.26e-194 < n < 1.44000000000000009e-251Initial program 78.1%
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 rewrites30.8%
Taylor expanded in i around 0
associate-*r/N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
lower-/.f6482.6
Applied rewrites82.6%
Taylor expanded in i around 0
Applied rewrites82.6%
if 5.0999999999999997e-14 < n Initial program 18.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.9
Applied rewrites93.9%
Taylor expanded in i around 0
Applied rewrites81.4%
(FPCore (i n)
:precision binary64
(if (<= n -2.7e+124)
(*
(*
(fma (fma (fma 0.041666666666666664 i 0.16666666666666666) i 0.5) i 1.0)
100.0)
n)
(if (<= n 5e-48)
(*
(/
1.0
(fma
(fma
(fma (* i i) -1.388888888888889e-5 0.0008333333333333334)
i
-0.005)
i
0.01))
n)
(*
(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 <= -2.7e+124) {
tmp = (fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * 100.0) * n;
} else if (n <= 5e-48) {
tmp = (1.0 / fma(fma(fma((i * i), -1.388888888888889e-5, 0.0008333333333333334), i, -0.005), i, 0.01)) * n;
} 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 <= -2.7e+124) tmp = Float64(Float64(fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * 100.0) * n); elseif (n <= 5e-48) tmp = Float64(Float64(1.0 / fma(fma(fma(Float64(i * i), -1.388888888888889e-5, 0.0008333333333333334), i, -0.005), i, 0.01)) * n); 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, -2.7e+124], N[(N[(N[(N[(N[(0.041666666666666664 * i + 0.16666666666666666), $MachinePrecision] * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[n, 5e-48], N[(N[(1.0 / N[(N[(N[(N[(i * i), $MachinePrecision] * -1.388888888888889e-5 + 0.0008333333333333334), $MachinePrecision] * i + -0.005), $MachinePrecision] * i + 0.01), $MachinePrecision]), $MachinePrecision] * n), $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 -2.7 \cdot 10^{+124}:\\
\;\;\;\;\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, i, 0.16666666666666666\right), i, 0.5\right), i, 1\right) \cdot 100\right) \cdot n\\
\mathbf{elif}\;n \leq 5 \cdot 10^{-48}:\\
\;\;\;\;\frac{1}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(i \cdot i, -1.388888888888889 \cdot 10^{-5}, 0.0008333333333333334\right), i, -0.005\right), i, 0.01\right)} \cdot n\\
\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 < -2.69999999999999978e124Initial program 16.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.f6495.5
Applied rewrites95.5%
Taylor expanded in i around 0
Applied rewrites76.4%
if -2.69999999999999978e124 < n < 4.9999999999999999e-48Initial program 35.4%
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.f6459.3
Applied rewrites59.3%
Applied rewrites59.2%
Taylor expanded in i around 0
Applied rewrites64.1%
if 4.9999999999999999e-48 < n Initial program 17.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.f6492.1
Applied rewrites92.1%
Taylor expanded in i around 0
Applied rewrites80.3%
(FPCore (i n)
:precision binary64
(let* ((t_0 (/ n (fma (fma 0.0008333333333333334 i -0.005) i 0.01))))
(if (<= n -2.1e+124)
(*
(*
(fma (fma (fma 0.041666666666666664 i 0.16666666666666666) i 0.5) i 1.0)
100.0)
n)
(if (<= n -1.26e-194)
t_0
(if (<= n 1.44e-251)
0.0
(if (<= n 1.55e-15)
t_0
(*
(fma
(fma (fma 4.166666666666667 i 16.666666666666668) i 50.0)
i
100.0)
n)))))))
double code(double i, double n) {
double t_0 = n / fma(fma(0.0008333333333333334, i, -0.005), i, 0.01);
double tmp;
if (n <= -2.1e+124) {
tmp = (fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * 100.0) * n;
} else if (n <= -1.26e-194) {
tmp = t_0;
} else if (n <= 1.44e-251) {
tmp = 0.0;
} else if (n <= 1.55e-15) {
tmp = t_0;
} else {
tmp = fma(fma(fma(4.166666666666667, i, 16.666666666666668), i, 50.0), i, 100.0) * n;
}
return tmp;
}
function code(i, n) t_0 = Float64(n / fma(fma(0.0008333333333333334, i, -0.005), i, 0.01)) tmp = 0.0 if (n <= -2.1e+124) tmp = Float64(Float64(fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * 100.0) * n); elseif (n <= -1.26e-194) tmp = t_0; elseif (n <= 1.44e-251) tmp = 0.0; elseif (n <= 1.55e-15) tmp = t_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_] := Block[{t$95$0 = N[(n / N[(N[(0.0008333333333333334 * i + -0.005), $MachinePrecision] * i + 0.01), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -2.1e+124], N[(N[(N[(N[(N[(0.041666666666666664 * i + 0.16666666666666666), $MachinePrecision] * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[n, -1.26e-194], t$95$0, If[LessEqual[n, 1.44e-251], 0.0, If[LessEqual[n, 1.55e-15], t$95$0, 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}
t_0 := \frac{n}{\mathsf{fma}\left(\mathsf{fma}\left(0.0008333333333333334, i, -0.005\right), i, 0.01\right)}\\
\mathbf{if}\;n \leq -2.1 \cdot 10^{+124}:\\
\;\;\;\;\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, i, 0.16666666666666666\right), i, 0.5\right), i, 1\right) \cdot 100\right) \cdot n\\
\mathbf{elif}\;n \leq -1.26 \cdot 10^{-194}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 1.44 \cdot 10^{-251}:\\
\;\;\;\;0\\
\mathbf{elif}\;n \leq 1.55 \cdot 10^{-15}:\\
\;\;\;\;t\_0\\
\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 < -2.10000000000000011e124Initial program 16.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.f6495.5
Applied rewrites95.5%
Taylor expanded in i around 0
Applied rewrites76.4%
if -2.10000000000000011e124 < n < -1.26e-194 or 1.44000000000000009e-251 < n < 1.5499999999999999e-15Initial program 25.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.f6463.7
Applied rewrites63.7%
Applied rewrites63.7%
Applied rewrites63.7%
Taylor expanded in i around 0
Applied rewrites66.2%
if -1.26e-194 < n < 1.44000000000000009e-251Initial program 78.1%
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 rewrites30.8%
Taylor expanded in i around 0
associate-*r/N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
lower-/.f6482.6
Applied rewrites82.6%
Taylor expanded in i around 0
Applied rewrites82.6%
if 1.5499999999999999e-15 < n Initial program 18.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.9
Applied rewrites93.9%
Taylor expanded in i around 0
Applied rewrites81.4%
(FPCore (i n)
:precision binary64
(let* ((t_0 (/ n (fma (fma 0.0008333333333333334 i -0.005) i 0.01)))
(t_1
(*
(fma
(fma (fma 4.166666666666667 i 16.666666666666668) i 50.0)
i
100.0)
n)))
(if (<= n -2.1e+124)
t_1
(if (<= n -1.26e-194)
t_0
(if (<= n 1.44e-251) 0.0 (if (<= n 1.55e-15) t_0 t_1))))))
double code(double i, double n) {
double t_0 = n / fma(fma(0.0008333333333333334, i, -0.005), i, 0.01);
double t_1 = fma(fma(fma(4.166666666666667, i, 16.666666666666668), i, 50.0), i, 100.0) * n;
double tmp;
if (n <= -2.1e+124) {
tmp = t_1;
} else if (n <= -1.26e-194) {
tmp = t_0;
} else if (n <= 1.44e-251) {
tmp = 0.0;
} else if (n <= 1.55e-15) {
tmp = t_0;
} else {
tmp = t_1;
}
return tmp;
}
function code(i, n) t_0 = Float64(n / fma(fma(0.0008333333333333334, i, -0.005), i, 0.01)) t_1 = Float64(fma(fma(fma(4.166666666666667, i, 16.666666666666668), i, 50.0), i, 100.0) * n) tmp = 0.0 if (n <= -2.1e+124) tmp = t_1; elseif (n <= -1.26e-194) tmp = t_0; elseif (n <= 1.44e-251) tmp = 0.0; elseif (n <= 1.55e-15) tmp = t_0; else tmp = t_1; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(n / N[(N[(0.0008333333333333334 * i + -0.005), $MachinePrecision] * i + 0.01), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(N[(4.166666666666667 * i + 16.666666666666668), $MachinePrecision] * i + 50.0), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision]}, If[LessEqual[n, -2.1e+124], t$95$1, If[LessEqual[n, -1.26e-194], t$95$0, If[LessEqual[n, 1.44e-251], 0.0, If[LessEqual[n, 1.55e-15], t$95$0, t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{n}{\mathsf{fma}\left(\mathsf{fma}\left(0.0008333333333333334, i, -0.005\right), i, 0.01\right)}\\
t_1 := \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4.166666666666667, i, 16.666666666666668\right), i, 50\right), i, 100\right) \cdot n\\
\mathbf{if}\;n \leq -2.1 \cdot 10^{+124}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;n \leq -1.26 \cdot 10^{-194}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 1.44 \cdot 10^{-251}:\\
\;\;\;\;0\\
\mathbf{elif}\;n \leq 1.55 \cdot 10^{-15}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if n < -2.10000000000000011e124 or 1.5499999999999999e-15 < n Initial program 17.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.f6494.5
Applied rewrites94.5%
Taylor expanded in i around 0
Applied rewrites79.7%
if -2.10000000000000011e124 < n < -1.26e-194 or 1.44000000000000009e-251 < n < 1.5499999999999999e-15Initial program 25.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.f6463.7
Applied rewrites63.7%
Applied rewrites63.7%
Applied rewrites63.7%
Taylor expanded in i around 0
Applied rewrites66.2%
if -1.26e-194 < n < 1.44000000000000009e-251Initial program 78.1%
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 rewrites30.8%
Taylor expanded in i around 0
associate-*r/N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
lower-/.f6482.6
Applied rewrites82.6%
Taylor expanded in i around 0
Applied rewrites82.6%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* (/ 1.0 (fma -0.005 i 0.01)) n))
(t_1
(*
(fma
(fma (fma 4.166666666666667 i 16.666666666666668) i 50.0)
i
100.0)
n)))
(if (<= n -2.9e+124)
t_1
(if (<= n -1.35e-192)
t_0
(if (<= n 1.44e-251) 0.0 (if (<= n 5e-48) t_0 t_1))))))
double code(double i, double n) {
double t_0 = (1.0 / fma(-0.005, i, 0.01)) * n;
double t_1 = fma(fma(fma(4.166666666666667, i, 16.666666666666668), i, 50.0), i, 100.0) * n;
double tmp;
if (n <= -2.9e+124) {
tmp = t_1;
} else if (n <= -1.35e-192) {
tmp = t_0;
} else if (n <= 1.44e-251) {
tmp = 0.0;
} else if (n <= 5e-48) {
tmp = t_0;
} else {
tmp = t_1;
}
return tmp;
}
function code(i, n) t_0 = Float64(Float64(1.0 / fma(-0.005, i, 0.01)) * n) t_1 = Float64(fma(fma(fma(4.166666666666667, i, 16.666666666666668), i, 50.0), i, 100.0) * n) tmp = 0.0 if (n <= -2.9e+124) tmp = t_1; elseif (n <= -1.35e-192) tmp = t_0; elseif (n <= 1.44e-251) tmp = 0.0; elseif (n <= 5e-48) tmp = t_0; else tmp = t_1; end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(1.0 / N[(-0.005 * i + 0.01), $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(N[(4.166666666666667 * i + 16.666666666666668), $MachinePrecision] * i + 50.0), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision]}, If[LessEqual[n, -2.9e+124], t$95$1, If[LessEqual[n, -1.35e-192], t$95$0, If[LessEqual[n, 1.44e-251], 0.0, If[LessEqual[n, 5e-48], t$95$0, t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{1}{\mathsf{fma}\left(-0.005, i, 0.01\right)} \cdot n\\
t_1 := \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4.166666666666667, i, 16.666666666666668\right), i, 50\right), i, 100\right) \cdot n\\
\mathbf{if}\;n \leq -2.9 \cdot 10^{+124}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;n \leq -1.35 \cdot 10^{-192}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 1.44 \cdot 10^{-251}:\\
\;\;\;\;0\\
\mathbf{elif}\;n \leq 5 \cdot 10^{-48}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if n < -2.90000000000000021e124 or 4.9999999999999999e-48 < n Initial program 17.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.f6493.2
Applied rewrites93.2%
Taylor expanded in i around 0
Applied rewrites79.0%
if -2.90000000000000021e124 < n < -1.34999999999999996e-192 or 1.44000000000000009e-251 < n < 4.9999999999999999e-48Initial program 26.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.f6463.8
Applied rewrites63.8%
Applied rewrites63.7%
Taylor expanded in i around 0
Applied rewrites63.0%
if -1.34999999999999996e-192 < n < 1.44000000000000009e-251Initial program 78.1%
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 rewrites30.8%
Taylor expanded in i around 0
associate-*r/N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
lower-/.f6482.6
Applied rewrites82.6%
Taylor expanded in i around 0
Applied rewrites82.6%
(FPCore (i n)
:precision binary64
(if (<= n -2.1e+124)
(*
(*
(fma (fma (fma 0.041666666666666664 i 0.16666666666666666) i 0.5) i 1.0)
100.0)
n)
(if (<= n 5e-48)
(/
n
(fma
(fma (fma (* i i) -1.388888888888889e-5 0.0008333333333333334) i -0.005)
i
0.01))
(*
(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 <= -2.1e+124) {
tmp = (fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * 100.0) * n;
} else if (n <= 5e-48) {
tmp = n / fma(fma(fma((i * i), -1.388888888888889e-5, 0.0008333333333333334), i, -0.005), i, 0.01);
} 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 <= -2.1e+124) tmp = Float64(Float64(fma(fma(fma(0.041666666666666664, i, 0.16666666666666666), i, 0.5), i, 1.0) * 100.0) * n); elseif (n <= 5e-48) tmp = Float64(n / fma(fma(fma(Float64(i * i), -1.388888888888889e-5, 0.0008333333333333334), i, -0.005), i, 0.01)); 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, -2.1e+124], N[(N[(N[(N[(N[(0.041666666666666664 * i + 0.16666666666666666), $MachinePrecision] * i + 0.5), $MachinePrecision] * i + 1.0), $MachinePrecision] * 100.0), $MachinePrecision] * n), $MachinePrecision], If[LessEqual[n, 5e-48], N[(n / N[(N[(N[(N[(i * i), $MachinePrecision] * -1.388888888888889e-5 + 0.0008333333333333334), $MachinePrecision] * i + -0.005), $MachinePrecision] * i + 0.01), $MachinePrecision]), $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 -2.1 \cdot 10^{+124}:\\
\;\;\;\;\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.041666666666666664, i, 0.16666666666666666\right), i, 0.5\right), i, 1\right) \cdot 100\right) \cdot n\\
\mathbf{elif}\;n \leq 5 \cdot 10^{-48}:\\
\;\;\;\;\frac{n}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(i \cdot i, -1.388888888888889 \cdot 10^{-5}, 0.0008333333333333334\right), i, -0.005\right), i, 0.01\right)}\\
\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 < -2.10000000000000011e124Initial program 16.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.f6495.5
Applied rewrites95.5%
Taylor expanded in i around 0
Applied rewrites76.4%
if -2.10000000000000011e124 < n < 4.9999999999999999e-48Initial program 35.4%
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.f6459.3
Applied rewrites59.3%
Applied rewrites59.2%
Applied rewrites59.2%
Taylor expanded in i around 0
Applied rewrites64.0%
if 4.9999999999999999e-48 < n Initial program 17.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.f6492.1
Applied rewrites92.1%
Taylor expanded in i around 0
Applied rewrites80.3%
(FPCore (i n)
:precision binary64
(let* ((t_0 (* (/ 1.0 (fma -0.005 i 0.01)) n)))
(if (<= n -1.35e-192)
t_0
(if (<= n 1.44e-251)
0.0
(if (<= n 3.4e-48)
t_0
(* (fma (fma 16.666666666666668 i 50.0) i 100.0) n))))))
double code(double i, double n) {
double t_0 = (1.0 / fma(-0.005, i, 0.01)) * n;
double tmp;
if (n <= -1.35e-192) {
tmp = t_0;
} else if (n <= 1.44e-251) {
tmp = 0.0;
} else if (n <= 3.4e-48) {
tmp = t_0;
} else {
tmp = fma(fma(16.666666666666668, i, 50.0), i, 100.0) * n;
}
return tmp;
}
function code(i, n) t_0 = Float64(Float64(1.0 / fma(-0.005, i, 0.01)) * n) tmp = 0.0 if (n <= -1.35e-192) tmp = t_0; elseif (n <= 1.44e-251) tmp = 0.0; elseif (n <= 3.4e-48) tmp = t_0; else tmp = Float64(fma(fma(16.666666666666668, i, 50.0), i, 100.0) * n); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(N[(1.0 / N[(-0.005 * i + 0.01), $MachinePrecision]), $MachinePrecision] * n), $MachinePrecision]}, If[LessEqual[n, -1.35e-192], t$95$0, If[LessEqual[n, 1.44e-251], 0.0, If[LessEqual[n, 3.4e-48], t$95$0, N[(N[(N[(16.666666666666668 * i + 50.0), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{1}{\mathsf{fma}\left(-0.005, i, 0.01\right)} \cdot n\\
\mathbf{if}\;n \leq -1.35 \cdot 10^{-192}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 1.44 \cdot 10^{-251}:\\
\;\;\;\;0\\
\mathbf{elif}\;n \leq 3.4 \cdot 10^{-48}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(16.666666666666668, i, 50\right), i, 100\right) \cdot n\\
\end{array}
\end{array}
if n < -1.34999999999999996e-192 or 1.44000000000000009e-251 < n < 3.40000000000000028e-48Initial program 23.4%
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.f6473.3
Applied rewrites73.3%
Applied rewrites73.2%
Taylor expanded in i around 0
Applied rewrites63.9%
if -1.34999999999999996e-192 < n < 1.44000000000000009e-251Initial program 78.1%
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 rewrites30.8%
Taylor expanded in i around 0
associate-*r/N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
lower-/.f6482.6
Applied rewrites82.6%
Taylor expanded in i around 0
Applied rewrites82.6%
if 3.40000000000000028e-48 < n Initial program 17.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.f6492.1
Applied rewrites92.1%
Taylor expanded in i around 0
Applied rewrites76.3%
(FPCore (i n)
:precision binary64
(let* ((t_0 (/ n (fma -0.005 i 0.01))))
(if (<= n -1.35e-192)
t_0
(if (<= n 1.44e-251)
0.0
(if (<= n 3.4e-48)
t_0
(* (fma (fma 16.666666666666668 i 50.0) i 100.0) n))))))
double code(double i, double n) {
double t_0 = n / fma(-0.005, i, 0.01);
double tmp;
if (n <= -1.35e-192) {
tmp = t_0;
} else if (n <= 1.44e-251) {
tmp = 0.0;
} else if (n <= 3.4e-48) {
tmp = t_0;
} else {
tmp = fma(fma(16.666666666666668, i, 50.0), i, 100.0) * n;
}
return tmp;
}
function code(i, n) t_0 = Float64(n / fma(-0.005, i, 0.01)) tmp = 0.0 if (n <= -1.35e-192) tmp = t_0; elseif (n <= 1.44e-251) tmp = 0.0; elseif (n <= 3.4e-48) tmp = t_0; else tmp = Float64(fma(fma(16.666666666666668, i, 50.0), i, 100.0) * n); end return tmp end
code[i_, n_] := Block[{t$95$0 = N[(n / N[(-0.005 * i + 0.01), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[n, -1.35e-192], t$95$0, If[LessEqual[n, 1.44e-251], 0.0, If[LessEqual[n, 3.4e-48], t$95$0, N[(N[(N[(16.666666666666668 * i + 50.0), $MachinePrecision] * i + 100.0), $MachinePrecision] * n), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{n}{\mathsf{fma}\left(-0.005, i, 0.01\right)}\\
\mathbf{if}\;n \leq -1.35 \cdot 10^{-192}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 1.44 \cdot 10^{-251}:\\
\;\;\;\;0\\
\mathbf{elif}\;n \leq 3.4 \cdot 10^{-48}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(16.666666666666668, i, 50\right), i, 100\right) \cdot n\\
\end{array}
\end{array}
if n < -1.34999999999999996e-192 or 1.44000000000000009e-251 < n < 3.40000000000000028e-48Initial program 23.4%
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.f6473.3
Applied rewrites73.3%
Applied rewrites73.2%
Applied rewrites73.2%
Taylor expanded in i around 0
Applied rewrites63.8%
if -1.34999999999999996e-192 < n < 1.44000000000000009e-251Initial program 78.1%
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 rewrites30.8%
Taylor expanded in i around 0
associate-*r/N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
lower-/.f6482.6
Applied rewrites82.6%
Taylor expanded in i around 0
Applied rewrites82.6%
if 3.40000000000000028e-48 < n Initial program 17.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.f6492.1
Applied rewrites92.1%
Taylor expanded in i around 0
Applied rewrites76.3%
(FPCore (i n) :precision binary64 (let* ((t_0 (* (fma (fma 16.666666666666668 i 50.0) i 100.0) n))) (if (<= n -2.9e-188) t_0 (if (<= n 3.2e-101) 0.0 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 <= -2.9e-188) {
tmp = t_0;
} else if (n <= 3.2e-101) {
tmp = 0.0;
} 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 <= -2.9e-188) tmp = t_0; elseif (n <= 3.2e-101) tmp = 0.0; 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, -2.9e-188], t$95$0, If[LessEqual[n, 3.2e-101], 0.0, 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 -2.9 \cdot 10^{-188}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;n \leq 3.2 \cdot 10^{-101}:\\
\;\;\;\;0\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if n < -2.9000000000000001e-188 or 3.19999999999999978e-101 < n Initial program 21.4%
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.f6486.3
Applied rewrites86.3%
Taylor expanded in i around 0
Applied rewrites68.0%
if -2.9000000000000001e-188 < n < 3.19999999999999978e-101Initial program 47.0%
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.8%
Taylor expanded in i around 0
associate-*r/N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
lower-/.f6461.8
Applied rewrites61.8%
Taylor expanded in i around 0
Applied rewrites61.8%
(FPCore (i n) :precision binary64 (if (<= i -1.95) 0.0 (if (<= i 4.6e+169) (* (fma 50.0 i 100.0) n) 0.0)))
double code(double i, double n) {
double tmp;
if (i <= -1.95) {
tmp = 0.0;
} else if (i <= 4.6e+169) {
tmp = fma(50.0, i, 100.0) * n;
} else {
tmp = 0.0;
}
return tmp;
}
function code(i, n) tmp = 0.0 if (i <= -1.95) tmp = 0.0; elseif (i <= 4.6e+169) tmp = Float64(fma(50.0, i, 100.0) * n); else tmp = 0.0; end return tmp end
code[i_, n_] := If[LessEqual[i, -1.95], 0.0, If[LessEqual[i, 4.6e+169], N[(N[(50.0 * i + 100.0), $MachinePrecision] * n), $MachinePrecision], 0.0]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -1.95:\\
\;\;\;\;0\\
\mathbf{elif}\;i \leq 4.6 \cdot 10^{+169}:\\
\;\;\;\;\mathsf{fma}\left(50, i, 100\right) \cdot n\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}
\end{array}
if i < -1.94999999999999996 or 4.5999999999999999e169 < i Initial program 65.0%
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 rewrites57.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-/.f6436.3
Applied rewrites36.3%
Taylor expanded in i around 0
Applied rewrites36.3%
if -1.94999999999999996 < i < 4.5999999999999999e169Initial program 11.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.f6481.3
Applied rewrites81.3%
Taylor expanded in i around 0
Applied rewrites76.1%
(FPCore (i n) :precision binary64 (if (<= i -1.38e+14) 0.0 (if (<= i 5.8e-7) (* 100.0 n) 0.0)))
double code(double i, double n) {
double tmp;
if (i <= -1.38e+14) {
tmp = 0.0;
} else if (i <= 5.8e-7) {
tmp = 100.0 * n;
} else {
tmp = 0.0;
}
return tmp;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: tmp
if (i <= (-1.38d+14)) then
tmp = 0.0d0
else if (i <= 5.8d-7) then
tmp = 100.0d0 * n
else
tmp = 0.0d0
end if
code = tmp
end function
public static double code(double i, double n) {
double tmp;
if (i <= -1.38e+14) {
tmp = 0.0;
} else if (i <= 5.8e-7) {
tmp = 100.0 * n;
} else {
tmp = 0.0;
}
return tmp;
}
def code(i, n): tmp = 0 if i <= -1.38e+14: tmp = 0.0 elif i <= 5.8e-7: tmp = 100.0 * n else: tmp = 0.0 return tmp
function code(i, n) tmp = 0.0 if (i <= -1.38e+14) tmp = 0.0; elseif (i <= 5.8e-7) tmp = Float64(100.0 * n); else tmp = 0.0; end return tmp end
function tmp_2 = code(i, n) tmp = 0.0; if (i <= -1.38e+14) tmp = 0.0; elseif (i <= 5.8e-7) tmp = 100.0 * n; else tmp = 0.0; end tmp_2 = tmp; end
code[i_, n_] := If[LessEqual[i, -1.38e+14], 0.0, If[LessEqual[i, 5.8e-7], N[(100.0 * n), $MachinePrecision], 0.0]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;i \leq -1.38 \cdot 10^{+14}:\\
\;\;\;\;0\\
\mathbf{elif}\;i \leq 5.8 \cdot 10^{-7}:\\
\;\;\;\;100 \cdot n\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}
\end{array}
if i < -1.38e14 or 5.7999999999999995e-7 < i Initial program 54.4%
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 rewrites46.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-/.f6430.7
Applied rewrites30.7%
Taylor expanded in i around 0
Applied rewrites30.7%
if -1.38e14 < i < 5.7999999999999995e-7Initial program 8.5%
Taylor expanded in i around 0
lower-*.f6483.8
Applied rewrites83.8%
(FPCore (i n) :precision binary64 0.0)
double code(double i, double n) {
return 0.0;
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
code = 0.0d0
end function
public static double code(double i, double n) {
return 0.0;
}
def code(i, n): return 0.0
function code(i, n) return 0.0 end
function tmp = code(i, n) tmp = 0.0; end
code[i_, n_] := 0.0
\begin{array}{l}
\\
0
\end{array}
Initial program 26.1%
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 rewrites20.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-/.f6416.8
Applied rewrites16.8%
Taylor expanded in i around 0
Applied rewrites16.8%
(FPCore (i n)
:precision binary64
(let* ((t_0 (+ 1.0 (/ i n))))
(*
100.0
(/
(-
(exp
(*
n
(if (== t_0 1.0)
(/ i n)
(/ (* (/ i n) (log t_0)) (- (+ (/ i n) 1.0) 1.0)))))
1.0)
(/ i n)))))
double code(double i, double n) {
double t_0 = 1.0 + (i / n);
double tmp;
if (t_0 == 1.0) {
tmp = i / n;
} else {
tmp = ((i / n) * log(t_0)) / (((i / n) + 1.0) - 1.0);
}
return 100.0 * ((exp((n * tmp)) - 1.0) / (i / n));
}
real(8) function code(i, n)
real(8), intent (in) :: i
real(8), intent (in) :: n
real(8) :: t_0
real(8) :: tmp
t_0 = 1.0d0 + (i / n)
if (t_0 == 1.0d0) then
tmp = i / n
else
tmp = ((i / n) * log(t_0)) / (((i / n) + 1.0d0) - 1.0d0)
end if
code = 100.0d0 * ((exp((n * tmp)) - 1.0d0) / (i / n))
end function
public static double code(double i, double n) {
double t_0 = 1.0 + (i / n);
double tmp;
if (t_0 == 1.0) {
tmp = i / n;
} else {
tmp = ((i / n) * Math.log(t_0)) / (((i / n) + 1.0) - 1.0);
}
return 100.0 * ((Math.exp((n * tmp)) - 1.0) / (i / n));
}
def code(i, n): t_0 = 1.0 + (i / n) tmp = 0 if t_0 == 1.0: tmp = i / n else: tmp = ((i / n) * math.log(t_0)) / (((i / n) + 1.0) - 1.0) return 100.0 * ((math.exp((n * tmp)) - 1.0) / (i / n))
function code(i, n) t_0 = Float64(1.0 + Float64(i / n)) tmp = 0.0 if (t_0 == 1.0) tmp = Float64(i / n); else tmp = Float64(Float64(Float64(i / n) * log(t_0)) / Float64(Float64(Float64(i / n) + 1.0) - 1.0)); end return Float64(100.0 * Float64(Float64(exp(Float64(n * tmp)) - 1.0) / Float64(i / n))) end
function tmp_2 = code(i, n) t_0 = 1.0 + (i / n); tmp = 0.0; if (t_0 == 1.0) tmp = i / n; else tmp = ((i / n) * log(t_0)) / (((i / n) + 1.0) - 1.0); end tmp_2 = 100.0 * ((exp((n * tmp)) - 1.0) / (i / n)); end
code[i_, n_] := Block[{t$95$0 = N[(1.0 + N[(i / n), $MachinePrecision]), $MachinePrecision]}, N[(100.0 * N[(N[(N[Exp[N[(n * If[Equal[t$95$0, 1.0], N[(i / n), $MachinePrecision], N[(N[(N[(i / n), $MachinePrecision] * N[Log[t$95$0], $MachinePrecision]), $MachinePrecision] / N[(N[(N[(i / n), $MachinePrecision] + 1.0), $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]], $MachinePrecision] - 1.0), $MachinePrecision] / N[(i / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 1 + \frac{i}{n}\\
100 \cdot \frac{e^{n \cdot \begin{array}{l}
\mathbf{if}\;t\_0 = 1:\\
\;\;\;\;\frac{i}{n}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{i}{n} \cdot \log t\_0}{\left(\frac{i}{n} + 1\right) - 1}\\
\end{array}} - 1}{\frac{i}{n}}
\end{array}
\end{array}
herbie shell --seed 2024295
(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))))