Kapitel 11 Jacobi–Verfahren

11.5. IMPLEMENTIERUNG 105
34 if (vec)
35 {
36 // update v(r,p) and v(r,q)
37 // ...
38 }
39 }
40
41 #undef a
42
43 } // end Jacobi rotation
44 } // end sweep
45 } // end iterations
46
47 sortEigenvalue(a,vec);
48 }
Aufgabe 11.1— Vervollständigen Sie die Funktion eigenvalue. Fügen Sie die
entsprechenden Befehle an den mit
1 // ...
gekennzeichneten Stellen ein. Benutzen Sie die Dateien
http://physik.uni-graz.at/~uxh/lineare-algebra/Jacobi.h
http://physik.uni-graz.at/~uxh/lineare-algebra/Jacobi.cc
Eine Testroutine ist in
http://physik.uni-graz.at/~uxh/lineare-algebra/test-Jacobi.cc
gegeben, bei der die Eigenwerte und Eigenvektoren der Matrix
bestimmt werden.
A =
⎛
⎜
⎝
17 3 2 13
3 11 11 8
2 11 8 12
13 8 12 2
⎞
⎟
⎠ , (11.16)