Prof. Dr. Bernd Hartke, Universit¨at Kiel, hartke@pctc.uni-kiel.de
Anleitung zu den Routinen jacobi und eigsrt:
subroutine jacobi(a,n,lda,e,v,nrot):
diagonalisiert die Matrix a
a: real*8-Matrix der Gr¨oße (n,n),
Eingabe: zu diagonalisierende Matrix (ist nachher zerst¨ort)
n: integer,
Eingabe: tats¨achliche Dimension des Diagonalisierungsproblems (hier: gleich Anzahl der Basisfunktionen) (Begrenzung n < 100 kann durch offensichtliche Modifikation des Programmtextes erh¨oht werden.)
lda: integer,
Eingabe: (feste) erste Dimension, mit der die Matrix a im aufrufenden Programm tats¨achlich vereinbart wurde
e: real*8-Vektor der L¨ange (n),
Ausgabe: Eigenwerte (ungeordnet!!)
v: real*8-Matrix der Gr¨oße (n,n),
Ausgabe: enth¨alt die Eigenvektoren als Spaltenvektoren (in derselben Ordnung wie im Vektor e)
nrot: integer,
Ausgabe: Anzahl der elementaren Iterationsschritte, die die Routine bis zur L¨osung ben¨otigt hat (hier weniger interessant).
subroutine eigsrt(e,v,n,lda):
Serviceprogramm zur Sortierung der Eigenwerte und -vektoren nach der Gr¨oße der Eigenwerte. Bedeutung der Parameter dieselbe wie bei jacobi.