Hi Mathias, Awesome. This looks correct to me. Would you be interested in creating a pull request to fix this in the master version?
Thanks, Timo On Mon, Nov 2, 2020 at 1:14 PM mathias....@posteo.de <mathias.anselm...@posteo.de> wrote: > > Hello again, > > I guess I have found the problem in dof_handler_policy.cc: In line 4565 the > n_global_dofs is set to the n_dofs(), not to the n_dofs(level). > Changing this fixes the problem for me. > I have created a commit in my dealii fork, where you can see it. > > Greetings > > Mathias > mathias....@posteo.de schrieb am Montag, 2. November 2020 um 18:06:27 UTC+1: >> >> >> 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/73afa160-05bd-420b-93ac-62ffac19b8b8n%40googlegroups.com. -- Timo Heister http://www.math.clemson.edu/~heister/ -- 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/CAMRj59H%3D8u%2BcDdKtRPePThjVLV1YhobrTSQKSXTU_pb7PePZyg%40mail.gmail.com.
