- Compare and contrast relative and absolute error
- Categorize a cost as \(\mathcal{O}(n^p)\)
- Categorize an error \(\mathcal{O}(h^p)\)
- Identify algebraic vs exponential growth and convergence

- Numerical algorithms are distinguished by their
and*cost*, and the tradeoff between them.*error* - The algorithms or methods introduced in this course indicate their error and cost whenever possible. These might be exact expressions or asymptotic bounds like \(\mathcal{O}(h^2)\) as \(h \to 0\) or \(\mathcal{O}(n^3)\) as \(n \to \infty\). For asymptotics we always indicate the limit.

Results computed using numerical methods are inaccurate – they are approximations to the true values. We can represent an approximate result as a combination of the true value and some error:

\[\begin{eqnarray} \text{Approximate Result} = \text{True Value} + \text{Error} \end{eqnarray}\] \[\hat{x} = x_0 + \Delta x\]Given this problem setup we can define the absolute error as:

\[\begin{equation} \text{Absolute Error} = |x_0 - \hat{x}| \end{equation} .\]This tells us how close our approximate result is to the actual answer. However, absolute error can become an unsatisfactory and misleading representation of the error depending on the magnitude of \(x_0\).

Case 1 | Case 2 |
---|---|

\(x_0 = 0.1\), \(\hat{x} = 0.2\) | \(x_0 = 100.0\), \(\hat{x} = 100.1\) |

\(\mid x_0 - \hat{x} \mid = 0.1\) | \(\mid x_0 - \hat{x} \mid = 0.1\) |

In both of these cases, the absolute error is the same, 0.1. However, we would intuitively consider case 2 more accurate than case 1 since our approximation is double the true value in case 1. Because of this, we define the relative error, which will be an error estimate independent of the magnitude. To obtain this we simply divide the absolute error by the absolute value of the true value.

\[\begin{equation} \text{Relative Error} = \frac{|x_0 - \hat{x}|}{|x_0|} \end{equation}\]If we consider the two cases again, we can see that the relative error will be much lower in the second case.

Case 1 | Case 2 |
---|---|

\(x_0 = 0.1\), \(\hat{x} = 0.2\) | \(x_0 = 100.0\), \(\hat{x} = 100.1\) |

\(\frac{\mid x_0 - \hat{x} \mid}{\mid x_0 \mid} = 1\) | \(\frac{\mid x_0 - \hat{x} \mid}{\mid x_0 \mid} = 10^{-3}\) |

**Significant figures** of a number are digits that carry meaningful information. They are digits beginning from the leftmost nonzero digit and ending with the rightmost “correct” digit, including final zeros that are exact. For example:

- The number 3.14159 has six significant digits.
- The number 0.00035 has two significant digits.
- The number 0.000350 has three significant digits.

An approximate result \(\hat{x}\) has \(n\) **significant figures** of a true value \(x_0\) if the absolute error, \(\vert x_0 - \hat{x}\vert\), has zeros in the first \(n\) decimal places counting from the leftmost nonzero (leading) digit of \(x_0\), followed by a digit from 0 to 4.

**Example:** Assume \(x_0 = 3.141592653\) and suppose \(\hat{x}\) is the approximate result:

The number of accurate digits can be estimated by the relative error. If

\[\text{Relative Error} = \frac{|x_0 - \hat{x}|}{|x_0|} \leq 10^{-n + 1}\]then \(\hat{x}\) has **at most** \(n\) correct digits. Conversely, if an approximation has \(n\) correct digits, then the relative error is \(\leq 10^{-n+1}\).

If our calculated quantities are vectors then instead of using the absolute value function, we can use the norm instead. Thus, our formulas become \(\begin{equation} \text{Absolute Error} = \|\mathbf{x} - \mathbf{\hat{x}}\| \end{equation}\)

\[\begin{equation} \text{Relative Error} = \frac{\|\mathbf{x} - \mathbf{\hat{x}}\|}{\|\mathbf{x}\|} \end{equation}\]We take the norm of the difference (and not the difference of the norms), because we are interested in how far apart these two quantities are. This formula is similar to finding that difference then using the vector norm to find the length of that difference vector.

Rounding error is the error that occurs from rounding values in a computation. This occurs constantly since computers use finite precision. Approximating \(\frac{1}{3} = 0.33333\dots\) with a finite decimal expansion is an example of rounding error.

Truncation error is the error from using an approximate algorithm in place of an exact mathematical procedure or function. For example, in the case of evaluating functions, we may represent our function by a finite Taylor series up to degree \(n\). The truncation error is the error that is incurred by not using the \(n+1\) term and above.

Big-O notation is used to understand and describe asymptotic behavior. The definition in the cases of approching 0 or \(\infty\) are as follows:

Let \(f\) and \(g\) be two functions. Then \(f(x) = \mathcal{O}(g(x))\) as \(x \rightarrow \infty\) if and only if there exists a value \(M\) and some \(x_0\) such that \(|f(x)| \leq M|g(x)|\) \(\forall x\) where \(x\geq x_0\)

Let \(f\) and \(g\) be two functions. Then \(f(h) = \mathcal{O}(g(h))\) as \(h \rightarrow 0\) if and only if there exists a value \(M\) and some \(h_0\) such that \(|f(h)| \leq M|g(h)|\) \(\forall h\) where \(0 < h < h_0\)

But what if we want to consider the function approaching an arbitrary value? Then we can redefine the expression as:

Let \(f\) and \(g\) be two functions. Then \(f(x) = \mathcal{O}(g(x))\) as \(x \rightarrow a\) if and only if there exists a value \(M\) and some \(\delta\) such that \(|f(x)| \leq M|g(x)|\) \(\forall x\) where \(0 < |x − a| < \delta\)

We can use Big-O to describe the time complexity of our algorithms.

Consider the case of matrix-matrix multiplication. If the size of each of our matrices is \(n \times n\), then the time it will take to multiply the matrices is \(\mathcal{O}(n^3)\) meaning that \(\text{Run time} \approx C \cdot n^3\). Suppose we know that for \(n_1=1000\), the matrix-matrix multiplication takes 5 seconds. Estimate how much time it would take if we double the size of our matrices to \(2n \times 2n\).

We know that:

\[\begin{align*} \text{Time}(2n_1) &\approx C \cdot (2n_1)^3 \\ &= C \cdot 2^3 \cdot n_1^3\\ &= 8 \cdot (C \cdot n_1^3) \\ &= 8 \cdot \text{Time}(n_1) \\ &= 40 \text{ seconds} \end{align*}\]So, when we double the size of our our matrices to \(2n \times 2n\), the time becomes \((2n)^3 = 8n^3\). Thus, the runtime will be roughly 8 times as long.

We can also use Big-O notation to describe the truncation error. A numerical method is called \(n\)-th order accurate if its truncation error \(E(h)\) obeys \(E(h) = \mathcal{O}(h^n)\).

Consider solving an interpolation problem. We have an interval of length \(h\) where our interpolant is valid and we know that our approximation is order \(\mathcal{O}(h^2)\). What this means is that as we decrease h (the interval length), our error will decrease quadratically. Using the definition of Big-O, we know that \(\text{Error} = C \cdot h^2\) where \(C\) is some constant.

In some cases, we may not know the exponent in \(E(h) = \mathcal{O}(h^n)\). We can estimate it using by computing the error at two different values of \(h\). Suppose we have two quantities, \(h_1 = 0.5\) and \(h_2 = 0.25\). We compute the corresponding errors as \(E(h_1) = 0.125\) and \(E(h_2) = 0.015625\). Then, since \(E(h) = \mathcal{O}(h^n)\), we have:

\[\begin{eqnarray}\frac{0.125}{0.015625} &=\frac{E(h_1)}{E(h_2)} \\ &\approx\frac{Ch_1^n}{Ch_2^n}\\ &=\left(\frac{h_1}{h_2}\right)^n\\ \end{eqnarray}\] \[\begin{eqnarray} \implies\log\left(\frac{0.125}{0.015625}\right)= n\log\left(\frac{h_1}{h_2}\right) =n\log\left(\frac{0.5}{0.25}\right) \end{eqnarray}\]Solving this equation for \(n\), we obtain \(n = 3\).

It is important that one does not place too much importance on the constant \(M\) in the definition of Big-O notation; it is essentially arbitrary.

Suppose \(f_1(n) = 10^{-20}n^2\) and \(f_2(n) = 10^{20}n^2\). While \(f_2\) is much larger than \(f_1\) for all values of \(n\), **both** are \(\mathcal{O}(n^2)\); this is obvious if we choose any constants \(M_1 \geq 10^{-20}\) and \(M_2 \geq 10^{20}\).

However, it is also true that \(f_2(n) = \mathcal{O}(10^{-20}n^2)\) for any constant \(M \geq 10^{40}\)

\[\begin{eqnarray} f_2(n) = 10^{20}n^2 = 10^{40} \times 10^{-20}n^2 \leq M\times 10^{-20}n^2. \end{eqnarray}\]So including a constant inside the \(\mathcal{O}\) is basically meaningless.

**Question:** What is the function \(g(n)\) that gives the tightest bound on \(f_2(n) = \mathcal{O}(g(n))\)?

**Solution:** the answer is \(g(n) = n^2\). For any \(r < 2\), there is **no** constant \(M\) such that \(|f_2(n)| \leq Mn^r\) for all
\(n\) sufficiently large. So \(n^r\) for \(r < 2\) is not a bound on \(f_2\). For any \(q > 2\), there exist a pair of constants \(M_1\) and \(M_2\) such that for all \(n\) sufficiently large:

However, we **cannot** find a pair of constants \(M_3\) and \(M_4\) such that:

Thus, we cannot “fit” another function in between \(f_2(n)\) and \(n^2\), so \(n^2\) is the tightest bound.

One may be tempted to think the correct answer should actually be \(g(n) = 10^{20}n^2\); however, this does not actually provide any additional information about the growth of \(f_2\). Notice that we didn’t specify what \(M_1\) and \(M_2\) were in the inequality above. Big-O notation says **nothing** about the size of the constant. The statements

are all equivalent, in that they all give the same amount of information on the growth of \(f_2\), since the constants are not specified. Since \(10^{-20}\) is very small, it may be tempting to conclude that it is “tighter” than the other two, which is not true. Therefore, it is always best practice to avoid placing unnecessary constants inside the \(\mathcal{O}\), and we expect you do refrain from doing so in this course.

Algebraic growth/convergence is when the coefficients \(a_n\) in the sequence we are interested in behave like \(\mathcal{O}(n^{\alpha})\) for growth and \(\mathcal{O}(1/n^{\alpha})\) for convergence, where \(\alpha\) is called the algebraic index of convergence. A sequence that grows or converges algebraically is a straight line in a log-log plot.

Exponential growth/convergence is when the coefficients \(a_n\) of the sequence we are interested in behave like \(\mathcal{O}(e^{qn^{\beta}})\) for growth and \(\mathcal{O}(e^{-qn^{\beta}})\) for convergence, where \(q\) is a constant for some \(\beta > 0\). Exponential growth is much faster than algebraic growth. Exponential growth/convergence is also sometimes called spectral growth/convergence. A sequence that grows exponentially is a straight line in a log-linear plot. Exponential convergence is often further classified as supergeometric, geometric, or subgeometric convergence.

Figures from J. P. Boyd, *Chebyshev and Fourier Spectral Methods*, 2nd ed., Dover, New York, 2001.

- See this review link

- 2020-04-25 Mariana Silva mfsilva@illinois.edu: small text revisions
- 2020-02-19 Peter Sentz sentz2@illinois.edu: Add section on role of constants, change Big-Oh’s to “mathcal”
- 2020-01-26 Wanjun Jiang wjiang24@illinois.edu: add scientific notations, digits and figures
- 2018-01-31 Aming Ni amingni2@illinois.edu: changed three graphs
- 2018-01-16 Yu Meng yumeng5@illinois.edu: minor fixes throughout
- 2017-11-02 Erin Carrier ecarrie2@illinois.edu: adds changelog
- 2017-10-26 Erin Carrier ecarrie2@illinois.edu: adds review questions, minor changes throughout to better match termiology in class notes
- 2017-10-23 John Doherty jjdoher2@illinois.edu: first complete draft
- 2017-10-17 Luke Olson lukeo@illinois.edu: outline