In many applications, such as equilibrium considerations in mechanical or electrical networks or in the discretisation of boundary value problems in ordinary and partial differential equations, one obtains very large equation systems, sometimes with many millions of lines. The coefficient matrices of these equation systems are typically sparse, i.e., most matrix entries are zero.