]> https://gitweb.dealii.org/ - dealii.git/commitdiff
number of blocks in FESystem should be just the number of base elements times their...
authorGuido Kanschat <dr.guido.kanschat@gmail.com>
Thu, 4 Oct 2007 00:19:32 +0000 (00:19 +0000)
committerGuido Kanschat <dr.guido.kanschat@gmail.com>
Thu, 4 Oct 2007 00:19:32 +0000 (00:19 +0000)
git-svn-id: https://svn.dealii.org/trunk@15263 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/fe/fe_system.cc
deal.II/deal.II/source/fe/fe_tools.cc

index ac3339edc13f7e027d791c2093d5bf1876b964e7..4c01c2037ddb7bfb87f36f1ae59f7c96bba6d05e 100644 (file)
@@ -1750,7 +1750,7 @@ FESystem<dim>::multiply_dof_numbers (const FiniteElementData<dim> &fe_data,
                                 fe_data.n_components() * N,
                                 fe_data.tensor_degree(),
                                 fe_data.conforming_space,
-                                fe_data.n_blocks() * N);
+                                N);
 }
 
 
@@ -2414,7 +2414,7 @@ FESystem<dim>::multiply_dof_numbers (const FiniteElementData<dim> &fe1,
     degree,
     typename FiniteElementData<dim>::Conformity(fe1.conforming_space
                                                & fe2.conforming_space),
-    fe1.n_blocks() * N1 + fe2.n_blocks() * N2);
+    N1 + N2);
 }
 
 
@@ -2453,7 +2453,7 @@ FESystem<dim>::multiply_dof_numbers (const FiniteElementData<dim> &fe1,
     typename FiniteElementData<dim>::Conformity(fe1.conforming_space
                                                & fe2.conforming_space
                                                & fe3.conforming_space),
-    fe1.n_blocks() * N1 + fe2.n_blocks() * N2 + fe3.n_blocks() * N3);
+    N1 + N2 + N3);
 }
 
 
index 7408578ba6beaa9e95867f81a32740180e069ab2..668b0b1ed473000b879aaf1accef7761d45d5503 100644 (file)
@@ -328,7 +328,7 @@ void FETools::compute_block_renumbering (
        k += element.base_element(b).n_dofs_per_cell();
       }
   Assert (i == element.n_blocks(), ExcInternalError());
-  
+//TODO:[GK] This does not work for a single RT  
   for (unsigned int i=0;i<element.dofs_per_cell;++i)
     {
       std::pair<unsigned int, unsigned int>

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.