From d37a35bc29b7945e78eb75698bc738bc9cb1bf05 Mon Sep 17 00:00:00 2001 From: Denis Davydov Date: Thu, 29 Mar 2018 17:14:45 +0200 Subject: [PATCH] doc: minor fixes in step-37 --- examples/step-37/doc/results.dox | 6 ++---- examples/step-37/step-37.cc | 1 + 2 files changed, 3 insertions(+), 4 deletions(-) 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); -- 2.39.5