]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Fix a few problems with the logic...
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 1 Sep 2006 21:33:27 +0000 (21:33 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 1 Sep 2006 21:33:27 +0000 (21:33 +0000)
git-svn-id: https://svn.dealii.org/trunk@13799 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/dofs/dof_tools.cc

index 0e90ade63971f5d68b480fac42c23aae2a15b633..74b9358d4c7543e352db9326054735b5a8522ebd 100644 (file)
@@ -2144,10 +2144,14 @@ namespace internal
              switch (mother_face_dominates)
                {
                  case FiniteElementDomination::this_element_dominates:
+                 case FiniteElementDomination::either_element_can_dominate:
                  {
-                                                    // Case 1 (the simple case):
-                                                    // The coarse element dominates
-                                                    // the elements on the subfaces
+                                                    // Case 1 (the simple
+                                                    // case): The coarse
+                                                    // element dominates the
+                                                    // elements on the
+                                                    // subfaces (or they are
+                                                    // all the same)
                    const unsigned int n_dofs_on_mother = cell->get_fe().dofs_per_face;
                    dofs_on_mother.resize (n_dofs_on_mother);
 
@@ -2492,19 +2496,6 @@ namespace internal
                    break;
                  }
 
-                 case FiniteElementDomination::either_element_can_dominate:
-                 {
-                                                    // hm, it isn't quite
-                                                    // clear what exactly we
-                                                    // would have to do
-                                                    // here. sit tight until
-                                                    // someone trips over the
-                                                    // following statement
-                                                    // and see what exactly
-                                                    // is going on
-                   Assert (false, ExcNotImplemented());
-                 }
-
                  default:
                                                         // we shouldn't get here
                        Assert (false, ExcInternalError());
@@ -2595,6 +2586,16 @@ namespace internal
                                                       // neither element has
                                                       // any constraints on
                                                       // its neighbor.
+                                                      //
+                                                      // the only case we know
+                                                      // for sure how to deal
+                                                      // with is if neither
+                                                      // element has any DoFs
+                                                      // on faces at all
+                     Assert (cell->get_fe().dofs_per_face == 0, ExcNotImplemented());
+                     Assert (cell->neighbor(face)->get_fe().dofs_per_face == 0,
+                             ExcNotImplemented());
+                     
                      break;
                    }
                    

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.