on the diagonal and none of the usual preconditioners (Jacobi, SSOR) will work
as they require division by diagonal elements.
+In the following, we will introduce some techniques that can be used in cases
+like these. Namely, we will consider the linear system as not consisting of one
+large matrix and vectors, but we will want to decompose matrices into vectors
+into <i>blocks</i> that correspond to the individual operators that appear in
+the system. We note that the resulting solver is not optimal -- there are much
+better ways, for example those explained in the results section of step-22 --
+but that the goal is to introduce techniques rather than optimal solvers.
+
<h4>Solving using the Schur complement</h4>
-In view of this, let us take another look at the matrix. If we sort our
+In view of the difficulties using standard solvers and preconditioners
+mentioned above, let us take another look at the matrix. If we sort our
degrees of freedom so that all velocity come before all pressure variables,
then we can subdivide the linear system $AX=B$ into the following blocks:
@f{eqnarray*}