Hello, I'm working with MultiGrid on a parallel::distributet::Triangulation and want to reorder the DoFs for the individual levels. After doing e.g. a Cuthill_McKee reordering I realised that the MGTools::make_sparsity_pattern() function throws an assertion in line 602:
Assert(sparsity.n_rows() == n_dofs, ExcDimensionMismatch(sparsity.n_rows(), n_dofs)); A simple test shows that indeed the number of level DoFs gets somehow changed during the reordering. I have attached a simple test, based on a deal.II test. The problem doesn't occur If one uses a normal (not parallel) triangulation. I'm still looking for the root cause of the problem, maybe someone has an idea here? Greetings Mathias -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/30146053-d7fd-4d09-8d42-df975f56fd06n%40googlegroups.com.
<<attachment: mg_reorder.zip>>