2.6M - 1. Institut für Theoretische Physik - Universität Stuttgart

4. Parallelisierung
36 t r a n s f e r (2�dmax+5)=dimag (E( j z v ( imax ) ) )
37 c a l l MPI Send ( t r a n s f e r , 2 �DMAX+5, MPI Double Precision ,
38 & imin , 4 2 , MPI Comm world , i e r r o r )
39 e l s e i f (Rank . eq . imin ) then
40 c a l l MPI Recv ( t r a n s f e r , 2 �DMAX+5, MPI Double Precision ,
41 & imax , 4 2 , MPI Comm World , status , i e r r o r )
42 do i =1,DMAX
43 r ( i , j z v ( imin)+1)= t r a n s f e r ( i )
44 F( i , j z v ( imin)+1)= t r a n s f e r ( i+dmax)
45 enddo
46 Psi ( j z v ( imin)+1)= t r a n s f e r (2�dmax+1)
47 WF( j z v ( imin)+1)=dcmplx ( t r a n s f e r (2�dmax+2) ,
48 & t r a n s f e r (2�dmax+3))
49 E( j z v ( imin )+1) =dcmplx ( t r a n s f e r (2�dmax+4) ,
50 & t r a n s f e r (2�dmax+5))
51 endif
52 j z v ( imin)= j z v ( imin)+1
53 j z v ( imax)= j z v ( imax)−1
54 goto 5000
55 5100 continue
56 j z=j z v (Rank)
57 END
Die Funktionsweise ist dabei wie folgt: die MPI-Funktion MPI AllGather vereint alle
Walker im Vektor jzv. Wenn an einem Prozessor ein relativer Überhang und an einem
der anderen Prozessoren ein relatives Defizit an Walkern besteht, so wird der überzählige
Walker verschoben. Der Prozessor mit dem Überhang schaltet dazu in den Sendemodus
(MPI Send), der mit dem Defizit in den Empfangsmodus (MPI Recv). Es müssen alle
Daten des Walkers verschoben werden. Dazu gehören alle zur aktuellen Position gespeicherten
Informationen wie die Position selbst (r), die Quantenkraft (F), der Wert der
Wellenfunktion (Psi), der Gewichtsfaktor (WF) und die lokale Energie (E). Die für diesen
Fall von MPI zur Verfügung gestellten sog. ” derived datatypes“ erwiesen sich als nicht
praktikabel, weshalb künstlich eine entsprechende Gemeinschaftsvariable (transfer) geschaffen
wurde.
4.1.3. Programmende
Die Statistikdaten werden am Ende des Blockes, also während des Programmablaufs
ständig zum Hauptprozeß (Rang 0) geschickt und separat verarbeitet, so daß am Ende
des Programms nur noch die Ausgabe der Ergebnisse steht (siehe Listing 4.3). Nach der
Angabe der Kernladungszahl, der Elektronenzahl N und der Magnetfeldstärke (β, βZ =
β/Z 2 ), folgen die Anzahl der Konfigurationen (Walker) und die Anzahl der parallel
verwendeten Rechner, sowie die gesamte Rechenzeit in Sekunden. Anschließend folgt
60