Dear All, What exactly do I need to do to properly partition and renumber dofs from a parallel::distributed::Triangulation<dim>? Before starting to use parallel::distributed::Triangulation<dim>, I partitioned it with GridTools::partition_triangulation (which seems to make the program crash now) and renumbered using DoFRenumbering::subdomain_wise() (which no longer works since it wouldn't make sense).
Step 40 seems to simply use dof_handler.distribute_dofs (fe); to partition, with no renumbering. Step 32 states in the comments that the mesh is partitioned with GridTools::partition_triangulation, but I find this nowhere in the actual code. It seems to only do dof_handler.distribute_dofs (fe); It then renumbers with DoFRenumbering::component_wise (stokes_dof_handler, stokes_sub_blocks); but since I'm working with a scalar, I don't think this is relevant. So, does dof_handler.distribute_dofs (fe); now take care of mesh partitioning? Many thanks, Pearl Flath
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
