From: Abner J. Salgado Date: Tue, 6 Oct 2009 20:21:07 +0000 (+0000) Subject: Changed the description of the initialization of the pressure space X-Git-Tag: v8.0.0~6970 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=75fa523ac6a9005bd0e9be82650b9acbfb589d73;p=dealii.git Changed the description of the initialization of the pressure space matrices. We do not have homogeneous Neumann bc's anymore, and so that comment was not needed git-svn-id: https://svn.dealii.org/trunk@19735 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/examples/step-35/step-35.cc b/deal.II/examples/step-35/step-35.cc index eb0e107c62..9792ee2ae5 100644 --- a/deal.II/examples/step-35/step-35.cc +++ b/deal.II/examples/step-35/step-35.cc @@ -874,23 +874,9 @@ NavierStokesProjection::initialize_velocity_matrices() vel_Laplace); } - // For the initialization of the matrices - // that act on the pressure space it is worth - // noticing one small detail. Since the - // projection step involves the solution of a - // Poisson equation with homogeneous Neumann - // boundary conditions, we need somehow to - // regularize this problem, that is to pick a - // solution. The way we do it is by setting - // the value of the pressure solution at the - // first node (wherever it is) to zero. This - // regularizes the problem and does not - // increase the number of entries in the - // sparsity pattern we use as the solution - // shown in step-11 would do. In the end it - // achieves the solution: it makes sure that - // the pressure is not left undetermined up - // to a constant. + // The initialization of the matrices + // that act on the pressure space is similar + // to the ones that act on the velocity space. template void NavierStokesProjection::initialize_pressure_matrices()