Test of the power_method C-function

JakobCHR.com

Quick Navigation:

Personal:

Go to Home
MS Research
PhD Research
Curriculum Vitae

General:

Linux

Soon to come:

Matlab
On-line Stores
Cycling
Medicine & Health
LaTeX
OOP & C++
Sony PCM-R500 DAT

Next: Test of the fractional_method Up: C-functions Previous: C-functions Contents Index

Test of the power_method C-function

To verify the correctness of the C-implementation of the power method we performed a test on the super-critical reactor and the sequence of iterates of the exact eigenvalue error, $E_{\lambda,k}$ , defined by (4.10) is seen in Figure 4.18. When compared to the sequence obtained by the Matlab-function (see Figure 4.9) it is concluded that the C-implementation gives the same results (except that in the C-function case the first iteration is not plotted). The power method was called with an accuracy criterion of $\epsilon = 10^{-8}$ and the resulting eigenvector, $\hspace{0.2ex}\underline{\phi}{}\hspace{0.15ex}_{C}$ , was compared to the eigenvector, $\hspace{0.2ex}\underline{\phi}{}\hspace{0.15ex}_{\mbox{\protect\scriptsize Matlab}}$ , obtained in section 4.3.1 by the Matlab function eig with the following result

$\begin{displaymath} \Vert \hspace{0.2ex}\underline{\phi}{}\hspace{0.15ex}_C - \... ...scriptsize Matlab}} \Vert _\infty \approx 6.0552 \cdot 10^{-8} \end{displaymath}$

(4.15)

Therefore, we conclude that the C-function works correctly.

$\begin{figure} % latex2html id marker 17323\rule{\textwidth}{0.2mm} \rule{0cm}... ...{0.2ex}\underline{\phi}{}\hspace{0.15ex}^{(0)} = [1,1,\ldots,1]^T$.}\end{figure}$

Top Quality
Developed with
Danish
Brain Power