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 up previous contents index
Next: Verification of the discretization Up: C-functions Previous: Test of the power_method   Contents   Index

Test of the fractional_method C-function

The C-implementation of the fractional iteration method, the function fractional_method, was run with the super-critical test reactor. As previously, the eigenvalue estimate was obtained by running the power method with a convergence criterion of $\epsilon = 10^{-2}$ and with this estimate at hand the fractional iteration method was brought to convergence with an accuracy criterion of $\epsilon = 10^{-14}$. The calculated eigenvalue, $\lambda_{\mbox{\protect\scriptsize FI}}$, had the following value

\begin{displaymath}
\lambda_{0,{\mbox{\protect\scriptsize FI}}} = 1.0738684260100841
\end{displaymath} (4.16)

which is almost identical to the exact (Matlab calculated) value

\begin{displaymath}
\lambda_{0,{\mbox{\protect\scriptsize Matlab}}} = 1.073868426010085
\end{displaymath} (4.17)

The calculated eigenvector, $\hspace{0.2ex}\underline{\phi}{}\hspace{0.15ex}_{\mbox{\protect\scriptsize FI}}$, was compared to the one calculated by Matlab with the result

\begin{displaymath}
\Vert \hspace{0.2ex}\underline{\phi}{}\hspace{0.15ex}_{\mbo...
...otect\scriptsize Matlab}} \Vert _\infty = 1.021 \cdot 10^{-14}
\end{displaymath} (4.18)

Finally, the exact eigenvalue error, $E_{\lambda,k}$, defined by (4.10) is depicted in Figure 4.19. When this figure is compared to its Matlab counterpart Figure 4.15 the correctness of the C-function is confirmed.

\begin{figure}
% latex2html id marker 17370\rule{\textwidth}{0.2mm}
\rule{0cm}...
...nsolution
obtained from the power method with $\epsilon = 10^{-2}$.}\end{figure}


next up previous contents index
Next: Verification of the discretization Up: C-functions Previous: Test of the power_method   Contents   Index  
 
 
 
Revision 2.0, Copyright © 1999-2004 Jakob Christensen
http://www.JakobCHR.com
E-Mail: webmaster@JakobCHR.com
Top Quality
Developed with

Danish
Brain Power
Linux Powered!