]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Add a paragraph on the suboptimality of the solver.
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Sun, 20 Jun 2010 18:47:11 +0000 (18:47 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Sun, 20 Jun 2010 18:47:11 +0000 (18:47 +0000)
git-svn-id: https://svn.dealii.org/trunk@21237 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-20/doc/intro.dox

index f96eae3f878c087bbeeb431dd4b302e6afabb680..0f992b73621218bcae5d838ae918ae0e57ec6bc0 100644 (file)
@@ -342,10 +342,19 @@ the next problem immediately surfaces: due to the zero block, there are zeros
 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*}

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.