From ee6351c7d3972c9bc00405981b31c84baa452ce4 Mon Sep 17 00:00:00 2001 From: bangerth Date: Tue, 5 Jan 2010 15:52:30 +0000 Subject: [PATCH] Fix one problem. git-svn-id: https://svn.dealii.org/trunk@20293 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-42/step-42.cc | 31 +++++++++++++++++++---------- 1 file changed, 21 insertions(+), 10 deletions(-) diff --git a/deal.II/examples/step-42/step-42.cc b/deal.II/examples/step-42/step-42.cc index 45f674a223..a90bc97f44 100644 --- a/deal.II/examples/step-42/step-42.cc +++ b/deal.II/examples/step-42/step-42.cc @@ -394,7 +394,8 @@ void StokesProblem::setup_dofs () std::vector dofs_per_block (2); - DoFTools::count_dofs_per_block (dof_handler, dofs_per_block, block_component); + DoFTools::count_dofs_per_block (dof_handler, dofs_per_block, + block_component); const unsigned int n_u = dofs_per_block[0], n_p = dofs_per_block[1]; @@ -447,17 +448,21 @@ void StokesProblem::setup_dofs () mg_interface_matrices.clear (); mg_sparsity.resize(0, nlevels-1); - MGTools::count_dofs_per_block (dof_handler, mg_dofs_per_component); + mg_dofs_per_component.resize (nlevels); + for (unsigned int level=0; level0) - mg_interface_matrices[level].reinit (mg_sparsity[level]); + mg_interface_matrices[level].reinit (mg_sparsity[level]); } } @@ -654,7 +659,13 @@ void StokesProblem::assemble_multigrid () } cell->get_mg_dof_indices (local_dof_indices); - boundary_constraints[level] +// for (unsigned int i=0; i::run () { { std::vector subdivisions (dim, 1); - subdivisions[0] = 4; + subdivisions[0] = 1; const Point bottom_left = (dim == 2 ? Point(0,0) : -- 2.39.5