Conjugate Gradient: HPF Issues
Each step classically data-parallel
- Reductions for dot products
- Array assignments to p and r
- Sparse matrix-vector multiply is data-parallel
Communication is probably well-structured
- Sparse matrix-vector multiply depends on PDE
Simple mapping suffices for regular grids
- Local PDE operators suggest BLOCK for all