Compute coarse grid matrices in one loop only. This makes
the code a tiny bit cleaner and implies that a
coarse matrix is completely setup before being
used to setup the next coarser matrix.

This allows to construct the coarse matrix occupation
pattern depending on the fine matrix entries.

[[Imported from SVN: r6119]]