Jacobi Iteration: MPI Program Design
Partitioning is simple
- Every point is a micro-task
Communication is simple
- 4 nearest neighbors in Cartesian mesh
- Reduction for convergence test
Agglomeration works along dimensions
- 1-D packing for high-latency machines (as minimizes number of messages)
- 2-D packing for others (most general as minimizes information sent)
- One process per processor practically required