  Iterative methods on irregular grids have been used widely in all areas of
comptational science and engineering for solving partial differential equations
with complex geometry. They provide the flexibility to express complex shapes
with relatively low computational cost. However, the direction of the evolution
of high-performance processors in the last two decades have caused serious
degradation of the computational efficiency of iterative methods on irregular
grids, because of relatively low memory bandwidth. Data compression can in
principle reduce the necessary memory memory bandwidth of iterative methods and
thus improve the efficiency. We have implemented several data compression
algorithms on the PEZY-SC processor, using the matrix generated for the HPCG
benchmark as an example. For the SpMV (Sparse Matrix-Vector multiplication)
part of the HPCG benchmark, the best implementation without data compression
achieved 11.6Gflops/chip, close to the theoretical limit due to the memory
bandwidth. Our implementation with data compression has achieved 32.4Gflops.
This is of course rather extreme case, since the grid used in HPCG is
geometrically regular and thus its compression efficiency is very high.
However, in real applications, it is in many cases possible to make a large
part of the grid to have regular geometry, in particular when the resolution is
high. Note that we do not need to change the structure of the program, except
for the addition of the data compression/decompression subroutines. Thus, we
believe the data compression will be very useful way to improve the performance
of many applications which rely on the use of irregular grids.
