Hello,
I am wondering if the DoFRenumbering algorithms operate on
the dof_handler associated with distributed grids (triangulation).
in dofs/dof_renumbering.cc there is little mention of distributed dofs.
In particular I add to setup_system in step-40 of the examples:
dof_handler.distribute_dofs (fe);
dealii::DoFRenumbering::Cuthill_McKee (dof_handler); // added this
and get this error message for mpirun -np x for x > 1.
Cycle 0:
--------------------------------------------------------
An error occurred in line <1718> of file
</home/tyoung/research.code/deal.II/source/dofs/dof_handler.cc> in
function void dealii::DoFHandler<dim, spacedim>::renumber_dofs(const
std::vector<unsigned int, std::allocator<unsigned int> >&) [with int
dim = 2, int spacedim = 2] The violated condition was:
new_numbers.size() == n_locally_owned_dofs() The name and call sequence
of the exception was: ExcRenumberingIncomplete() Additional
Information: (none)
Stacktrace:
-----------
#0 /home/tyoung/research.code/deal.II/lib/libdeal_II.g.so.7.1.pre:
dealii::DoFHandler<2, 2>::renumber_dofs(std::vector<unsigned int,
std::allocator<unsigned int> > const&)
#1 /home/tyoung/research.code/deal.II/lib/libdeal_II.g.so.7.1.pre:
void dealii::DoFRenumbering::Cuthill_McKee<dealii::DoFHandler<2, 2>
>(dealii::DoFHandler<2, 2>&, bool, bool, std::vector<unsigned int,
>std::allocator<unsigned int> > const&) #2 ./step-40:
>LaplaceProblem<2>::setup_system() #3 ./step-40:
>LaplaceProblem<2>::run() #4 ./step-40: main
>--------------------------------------------------------
Which makes some sense.
I am interested to know how to go about this: As I am writing this I
wonder if a loop over all subdomains is required?
Any ideas? Thanks.
Best,
Toby
--
Toby D. Young
Assistant Professor Philosophy & Physics
Polish Academy of Sciences
www: http://www.ippt.gov.pl/~tyoung
skype: stenografia
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
