20.05.2014 Views

link to my thesis

link to my thesis

link to my thesis

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

92 CHAPTER 6. ITERATIVE EIGENVALUE SOLVERS<br />

this matrix W as follows:<br />

P = W (W T W ) −1 W T =<br />

⎛<br />

⎞<br />

0 0 0 0 0 0 0 0 0<br />

5 5<br />

0<br />

− 1 1 7<br />

− 1 − 1 1<br />

48 24 16 24 48 8 48 12<br />

5 5<br />

0<br />

− 1 1 7<br />

− 1 − 1 1<br />

24 12 8 12 24 4 24 6<br />

0 − 1 − 1 5 1<br />

− 1 5 7 1<br />

16 8 48 24 48 24 48 12<br />

1 1 1 1 1 1 1 1<br />

0 24 12 24 12 8 12 8 6<br />

.<br />

7 7<br />

0<br />

− 1 1 13<br />

− 1 5 1<br />

48 24 48 8 48 24 48 4<br />

0 − 1 − 1 5 1<br />

− 1 5 7 1<br />

8 4 24 12 24 12 24 6<br />

⎜<br />

⎝ 0 − 1 − 1 7 1 5 7 13 1 ⎟<br />

48 24 48 8 48 24 48 4 ⎠<br />

0<br />

1<br />

12<br />

1<br />

6<br />

1<br />

12<br />

1<br />

6<br />

1<br />

4<br />

1<br />

6<br />

1<br />

4<br />

1<br />

3<br />

(6.32)<br />

Note that the explicit matrix P is constructed here for analytical purposes, but is<br />

avoided in numerical computations.<br />

The vec<strong>to</strong>rs log(ˆυ) and ˆυ are computed as:<br />

⎛<br />

log(ˆυ) =P · log(υ) =<br />

⎜<br />

⎝<br />

0.000<br />

1.675<br />

3.349<br />

2.532<br />

4.207<br />

5.881<br />

5.064<br />

6.738<br />

8.413<br />

⎞<br />

⎟<br />

⎠<br />

⎛<br />

and ˆυ = e log(ˆυ) =<br />

⎜<br />

⎝<br />

1.000<br />

5.337<br />

28.499<br />

12.576<br />

67.120<br />

358.234<br />

158.149<br />

844.083<br />

4505.076<br />

⎞<br />

. (6.33)<br />

⎟<br />

⎠<br />

For a particular choice of k =(k 1 ,...,k N )=(0, 0, 1, 0, 0, 1, −1, 0, 0) the projected<br />

vec<strong>to</strong>r ˆψ is computed as:<br />

⎛<br />

ˆψ = Pψ+2π (P − I) k =<br />

⎜<br />

⎝<br />

0.000<br />

3.903<br />

1.523<br />

−2.424<br />

1.479<br />

−0.901<br />

1.435<br />

−0.945<br />

2.958<br />

⎞<br />

. (6.34)<br />

⎟<br />

⎠<br />

The entries in the projected vec<strong>to</strong>r log(ˆυ) satisfy the following relations in order<br />

<strong>to</strong> exhibit Stetter structure with respect <strong>to</strong> addition (because it is in the logarithmic

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!