# designated by a single-character variable, since that can be
# reference using $D only, i.e. without the parentheses that are
# required for most other parameters, as e.g. in $(target).
-D = ../../
+# D = ../../
+D = ~/deal.II_mg/deal.II
# The last field specifies the names of data and other files that
# shall be deleted when calling `make clean'. Object and backup files,
mg_dofs_per_component.resize (nlevels);
for (unsigned int level=0; level<nlevels; ++level)
mg_dofs_per_component[level].resize (2);
-
- MGTools::count_dofs_per_component (dof_handler, mg_dofs_per_component,
+
+ MGTools::count_dofs_per_block (dof_handler, mg_dofs_per_component,
block_component);
+ for (unsigned int level=0; level<nlevels; ++level)
+ std::cout << mg_dofs_per_component[level].size () << std::endl;
+
for (unsigned int level=0; level<nlevels; ++level)
{
BlockCompressedSparsityPattern bcsp (mg_dofs_per_component[level],
mg_smoother(mg_vector_memory);
MGLevelObject<typename Smoother::AdditionalData>
- smoother_data (0, triangulation.n_levels());
+ smoother_data (0, triangulation.n_levels()-1);
for (unsigned int level=0; level<triangulation.n_levels(); ++level)
smoother_data[level].A_preconditioner = mg_A_preconditioner[level].get();
{
try
{
- deallog.depth_console (0);
+// deallog.depth_console (0);
StokesProblem<2> flow_problem(1);
flow_problem.run ();