Das Beispielprogramm TRACE

Problem-Stellung:

• Berechnet die Spur (Summe der Diagonalelemente) einer quadratischen Matrix

Sourcen:

• trace.c

Bemerkungen:

• Matrix blockweise auf quadratisches Prozessorgitter verteilt
• neuer Communicator mit MPI_Comm_split:

  diagonale Tasks:	color = 0
  andere:	color = MPI_UNDEFINED

• Beispiel-Matrizen:
  
  $$A = (a_{ij}), \ a_{ij} = i+j$$
  
  
  
  
  $$B = (b_{ij}), \ b_{ij} = i-j$$
  
  
  
  
  $$C = A * B \Rightarrow \mathrm{Tr}\, C = 0$$