Cache optimization for structured and unstructured grid multigrid

Abstract. Many current computer designs employ caches and a hierarchical memory architec-ture. The speed of a code depends on how well the cache structure is exploited. The number of cache misses provides a better measure for comparing algorithms than the number of multiplies. In this paper, suitable blocking strategies for both structured and unstructured grids will be introduced. They improve the cache usage without changing the underlying algorithm. In particular, bitwise compatibility is guaranteed between the standard and the high performance implementations of the algorithms. This is illustrated by comparisons for various multigrid algorithms on a selection of dierent computers for problems in two and three dimensions. The code restructuring can yield performance improvements of factors of 2-5. This allows the modi ed codes to achieve a much higher percentage of the peak performance of the CPU than is usually observed with standard implementations. Key words. Computer architectures, iterative algorithms, multigrid, high performance com-puting, cache. AMS subject classications. 65M55, 65N55, 65F10, 68-04, 65Y99 1. Introduction. Years ag