]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Enable hp::DoFHandler on parallel meshes.
authorWolfgang Bangerth <bangerth@colostate.edu>
Sat, 22 Jul 2017 00:56:07 +0000 (18:56 -0600)
committerWolfgang Bangerth <bangerth@colostate.edu>
Mon, 21 Aug 2017 23:14:09 +0000 (17:14 -0600)
source/hp/dof_handler.cc

index d0e8ab3eab18a7171da46f4250d5cfcc407516dc..92d44ebbb0c6f4071cc2a360cfa01b14cf9b5ef7 100644 (file)
@@ -969,20 +969,11 @@ namespace hp
     tria(&tria, typeid(*this).name()),
     faces (nullptr)
   {
-    Assert ((dynamic_cast<const parallel::distributed::Triangulation< dim, spacedim >*>
-             (&tria)
-             == nullptr),
-            ExcMessage ("The given triangulation is parallel distributed but "
-                        "this class does not currently support this."));
-
-    // decide whether we need a sequential or a parallel distributed policy
+    // decide whether we need a sequential or a parallel shared/distributed policy
     if (dynamic_cast<const parallel::shared::Triangulation< dim, spacedim>*> (&*this->tria) != nullptr)
       policy.reset (new internal::DoFHandler::Policy::ParallelShared<DoFHandler<dim,spacedim> >(*this));
     else if (dynamic_cast<const parallel::distributed::Triangulation< dim, spacedim >*> (&*this->tria) != nullptr)
-      //policy.reset (new internal::DoFHandler::Policy::ParallelDistributed<DoFHandler<dim,spacedim> >(*this));
-      {
-        Assert (false, ExcNotImplemented());
-      }
+      policy.reset (new internal::DoFHandler::Policy::ParallelDistributed<DoFHandler<dim,spacedim> >(*this));
     else
       policy.reset (new internal::DoFHandler::Policy::Sequential<DoFHandler<dim,spacedim> >(*this));
 

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.