]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Clarify something about the matrix in step-20. 9470/head
authorWolfgang Bangerth <bangerth@colostate.edu>
Mon, 3 Feb 2020 02:37:04 +0000 (19:37 -0700)
committerWolfgang Bangerth <bangerth@colostate.edu>
Tue, 4 Feb 2020 14:59:20 +0000 (07:59 -0700)
examples/step-20/doc/intro.dox

index 22b5ea333e57b4a7f31a7868c07c7054d87dfe85..86e0b8656cfbb44def62e82f9ae6eb2b7556a3b8 100644 (file)
@@ -358,17 +358,25 @@ program. We will therefore not comment much on it below.
 <h3>Linear solvers and preconditioners</h3>
 
 After assembling the linear system we are faced with the task of solving
-it. The problem here is: the matrix has a zero block at the bottom right
-(there is no term in the bilinear form that couples the pressure $p$ with the
-pressure test function $q$), and it is indefinite. At least it is
-symmetric. In other words: the Conjugate Gradient method is not going to
-work since it is only applicable to problems in which the matrix is
-symmetric and positive definite.
+it. The problem here is that the matrix possesses two undesirable properties:
+- It is <a href="https://en.wikipedia.org/wiki/Definiteness_of_a_matrix">indefinite</a>,
+  i.e., it has both positive and negative eigenvalues.
+  We don't want to prove this property here, but note that this is true
+  for all matrices of the form
+  $\left(\begin{array}{cc} M & B \\ B^T & 0 \end{array}\right)$
+  such as the one here where $M$ is positive definite.
+- The matrix has a zero block at the bottom right (there is no term in
+  the bilinear form that couples the pressure $p$ with the
+  pressure test function $q$).
+
+At least it is symmetric, but the first issue above still means that
+the Conjugate Gradient method is not going to work since it is only
+applicable to problems in which the matrix is symmetric and positive definite.
 We would have to resort to other iterative solvers instead, such as
 MinRes, SymmLQ, or GMRES, that can deal with indefinite systems. However, then
-the next problem immediately surfaces: due to the zero block, there are zeros
-on the diagonal and none of the usual, "simple" preconditioners (Jacobi, SSOR) will work
-as they require division by diagonal elements.
+the next problem immediately surfaces: Due to the zero block, there are zeros
+on the diagonal and none of the usual, "simple" preconditioners (Jacobi, SSOR)
+will work as they require division by diagonal elements.
 
 For the matrix sizes we expect to run with this program, the by far simplest
 approach would be to just use a direct solver (in particular, the

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.