next up previous contents
Next: BP_EIV Up: BP_MAT Previous: I/O Files   Contents

MPI implementation

The MPI implementation computes all columns in an interleaved fashion. Each node computes only the nprocs + myid, where jb, nprocs, myid are correspondingly the current column, the number of processes (nodes), and ID number of the node.

 For column = 1,ncfg,nprocs
    For row = column to ncfg
       compute the matrix elements for all 
       three types HNR, HZETA, HSPIN

Each node reads input angular data for the current column and writes the contributions to a set of files. The filenames contain the processor ID as a four digit number, Figure  8.30.

Figure 8.30: Subroutine brevalr().
\begin{figure}\begin{center}
\centerline{\psfig{figure=tex/fig/mpi_io_bp_mat.eps}}\end{center}\end{figure}



2001-10-11