From: Denis Davydov <davydden@gmail.com>
Date: Thu, 29 Mar 2018 15:14:45 +0000 (+0200)
Subject: doc: minor fixes in step-37
X-Git-Tag: v9.0.0-rc1~238^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=d37a35bc29b7945e78eb75698bc738bc9cb1bf05;p=dealii.git

doc: minor fixes in step-37
---

diff --git a/examples/step-37/doc/results.dox b/examples/step-37/doc/results.dox
index 789b843e44..5a7789d33c 100644
--- a/examples/step-37/doc/results.dox
+++ b/examples/step-37/doc/results.dox
@@ -432,9 +432,7 @@ the set of all nodes. The set $\mathcal N_D\subset \mathcal N$ is the subset
 of the nodes that are subject to Dirichlet boundary conditions where the
 solution is forced to equal $u_i = g_i = g(\mathbf{x}_i)$ as the interpolation
 of boundary values on the Dirichlet-constrained node points $i\in \mathcal
-N_D$. In this formula, where we have assumed that the degrees of freedom are
-numbered such that the ones subject to a Dirichlet condition are all appended
-to the end of the list of degrees of freedom. We then insert this solution
+N_D$. We then insert this solution
 representation into the weak form, e.g. the Laplacian shown above, and move
 the known quantities to the right hand side:
 @f{eqnarray*}
@@ -447,7 +445,7 @@ In this formula, the equations are tested for all basis functions $\varphi_i$
 with $i\in N \setminus \mathcal N_D$ that are not related to the nodes
 constrained by Dirichlet conditions.
 
-In the implementation in deal.II, the integrals $(\varphi_i,\varphi_j)_\Omega$
+In the implementation in deal.II, the integrals $(\nabla \varphi_i,\nabla \varphi_j)_\Omega$
 on the right hand side are already contained in the local matrix contributions
 we assemble on each cell. When using
 ConstraintMatrix::distributed_local_to_global() as first described in the
diff --git a/examples/step-37/step-37.cc b/examples/step-37/step-37.cc
index f19600b3c8..23d44253d1 100644
--- a/examples/step-37/step-37.cc
+++ b/examples/step-37/step-37.cc
@@ -1094,6 +1094,7 @@ namespace Step37
 
     time.reset();
     time.start();
+    constraints.set_zero(solution);
     cg.solve (system_matrix, solution, system_rhs,
               preconditioner);