From 75fa523ac6a9005bd0e9be82650b9acbfb589d73 Mon Sep 17 00:00:00 2001 From: "Abner J. Salgado" Date: Tue, 6 Oct 2009 20:21:07 +0000 Subject: [PATCH] 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 --- deal.II/examples/step-35/step-35.cc | 20 +++----------------- 1 file changed, 3 insertions(+), 17 deletions(-) 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() -- 2.39.5