Block orderings for tensor-product grids in two and three dimensions (Q1610719)

Two-line and two-plane orderings for a convection-diffusion model problem in two and three dimensions are considered. The strategies are aimed at introducing dense diagonal blocks, at the price of a sligh increase of the bandwidth of the matrix, compared to natural lexicographic ordering. A convergence analysis is performed for the block Jacobi scheme. Furthermore a two-step preconditioning technique is proposed. For three dimensional problems this approach yields a viable alternative to the incomplete LU approach, and may be easier to implement in parallel environments. The analysis is illustrated and validated by numerical examples.