]> https://gitweb.dealii.org/ - dealii.git/commitdiff
fix some of the multigrid 64 bit indexing
authorTimo Heister <timo.heister@gmail.com>
Mon, 25 Mar 2013 14:21:11 +0000 (14:21 +0000)
committerTimo Heister <timo.heister@gmail.com>
Mon, 25 Mar 2013 14:21:11 +0000 (14:21 +0000)
git-svn-id: https://svn.dealii.org/branches/branch_bigger_global_dof_indices_4@29025 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/source/dofs/dof_handler_policy.cc

index e4ec53924140f50b606154df903bc8c9aa6d8810..7155778e9a3c947c6306b1508b4d797f7d71fc3a 100644 (file)
@@ -1316,7 +1316,7 @@ namespace internal
           const typename dealii::internal::p4est::types<dim>::quadrant &p4est_cell,
           const typename DoFHandler<dim,spacedim>::level_cell_iterator &dealii_cell,
           const typename dealii::internal::p4est::types<dim>::quadrant &quadrant,
-          unsigned int *dofs,
+          dealii::types::global_dof_index* dofs,
           unsigned int level)
         {
           if (internal::p4est::quadrant_is_equal<dim>(p4est_cell, quadrant))
@@ -1570,8 +1570,11 @@ namespace internal
               typename dealii::internal::p4est::types<dim>::quadrant *quadrant
                 =reinterpret_cast<typename dealii::internal::p4est::types<dim>::quadrant *>(ptr);
               ptr+=cells*sizeof(typename dealii::internal::p4est::types<dim>::quadrant);
-              dealii::types::global_dof_index *dofs=reinterpret_cast<dealii::types::global_dof_index *>(ptr);
+              dealii::types::global_dof_index *dofs
+               = reinterpret_cast<dealii::types::global_dof_index*>(ptr);
 
+             // the dofs pointer contains for each cell the number of dofs
+             // on that cell (dofs[0]) followed by the dof indices itself.
               for (unsigned int c=0; c<cells; ++c, dofs+=1+dofs[0])
                 {
                   typename DoFHandler<dim,spacedim>::level_cell_iterator
@@ -1809,8 +1812,11 @@ namespace internal
               typename dealii::internal::p4est::types<dim>::quadrant *quadrant
                 =reinterpret_cast<typename dealii::internal::p4est::types<dim>::quadrant *>(ptr);
               ptr+=cells*sizeof(typename dealii::internal::p4est::types<dim>::quadrant);
-              unsigned int *dofs=reinterpret_cast<unsigned int *>(ptr);
+              dealii::types::global_dof_index *dofs
+               = reinterpret_cast<dealii::types::global_dof_index*>(ptr);
 
+             // the dofs pointer contains for each cell the number of dofs
+             // on that cell (dofs[0]) followed by the dof indices itself.
               for (unsigned int c=0; c<cells; ++c, dofs+=1+dofs[0])
                 {
                   typename DoFHandler<dim,spacedim>::level_cell_iterator

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.