Reference documentation for deal.II version Git 2618e0f 2017-11-23 17:25:26 +0100
dof_handler.h
1 // ---------------------------------------------------------------------
2 //
3 // Copyright (C) 1998 - 2017 by the deal.II authors
4 //
5 // This file is part of the deal.II library.
6 //
7 // The deal.II library is free software; you can use it, redistribute
8 // it, and/or modify it under the terms of the GNU Lesser General
9 // Public License as published by the Free Software Foundation; either
10 // version 2.1 of the License, or (at your option) any later version.
11 // The full text of the license can be found in the file LICENSE at
12 // the top level of the deal.II distribution.
13 //
14 // ---------------------------------------------------------------------
15 
16 #ifndef dealii_dof_handler_h
17 #define dealii_dof_handler_h
18 
19 
20 
21 #include <deal.II/base/config.h>
22 #include <deal.II/base/exceptions.h>
23 #include <deal.II/base/smartpointer.h>
24 #include <deal.II/base/index_set.h>
25 #include <deal.II/base/iterator_range.h>
26 #include <deal.II/dofs/block_info.h>
27 #include <deal.II/dofs/dof_iterator_selector.h>
28 #include <deal.II/dofs/number_cache.h>
29 #include <deal.II/dofs/dof_faces.h>
30 #include <deal.II/dofs/dof_levels.h>
31 #include <deal.II/dofs/function_map.h>
32 #include <deal.II/hp/fe_collection.h>
33 
34 #include <boost/serialization/split_member.hpp>
35 
36 #include <vector>
37 #include <map>
38 #include <set>
39 #include <memory>
40 
41 DEAL_II_NAMESPACE_OPEN
42 
43 template <int dim, int spacedim> class FiniteElement;
44 template <int dim, int spacedim> class Triangulation;
45 
46 namespace internal
47 {
48  namespace DoFHandler
49  {
50  struct Implementation;
51 
52  namespace Policy
53  {
54  template <int dim, int spacedim> class PolicyBase;
55  struct Implementation;
56  }
57  }
58 
59  namespace DoFAccessor
60  {
61  struct Implementation;
62  }
63 
64  namespace DoFCellAccessor
65  {
66  struct Implementation;
67  }
68 }
69 
70 
191 template <int dim, int spacedim=dim>
192 class DoFHandler : public Subscriptor
193 {
194  typedef ::internal::DoFHandler::Iterators<DoFHandler<dim,spacedim>, false> ActiveSelector;
195  typedef ::internal::DoFHandler::Iterators<DoFHandler<dim,spacedim>, true> LevelSelector;
196 public:
197  typedef typename ActiveSelector::CellAccessor cell_accessor;
198  typedef typename ActiveSelector::FaceAccessor face_accessor;
199 
200  typedef typename ActiveSelector::line_iterator line_iterator;
201  typedef typename ActiveSelector::active_line_iterator active_line_iterator;
202 
203  typedef typename ActiveSelector::quad_iterator quad_iterator;
204  typedef typename ActiveSelector::active_quad_iterator active_quad_iterator;
205 
206  typedef typename ActiveSelector::hex_iterator hex_iterator;
207  typedef typename ActiveSelector::active_hex_iterator active_hex_iterator;
208 
229  typedef typename ActiveSelector::active_cell_iterator active_cell_iterator;
230 
257  typedef typename ActiveSelector::cell_iterator cell_iterator;
258 
275  typedef typename ActiveSelector::face_iterator face_iterator;
276 
288  typedef typename ActiveSelector::active_face_iterator active_face_iterator;
289 
290  typedef typename LevelSelector::CellAccessor level_cell_accessor;
291  typedef typename LevelSelector::FaceAccessor level_face_accessor;
292 
293  typedef typename LevelSelector::cell_iterator level_cell_iterator;
294  typedef typename LevelSelector::face_iterator level_face_iterator;
295 
296 
300  static const unsigned int dimension = dim;
301 
305  static const unsigned int space_dimension = spacedim;
306 
318  DEAL_II_DEPRECATED
320 
330  static const unsigned int default_fe_index = 0;
331 
337  DoFHandler ();
338 
343 
350  DoFHandler (const DoFHandler &) = delete;
351 
355  virtual ~DoFHandler ();
356 
363  DoFHandler &operator = (const DoFHandler &) = delete;
364 
370  const FiniteElement<dim,spacedim> &fe);
371 
388  virtual void distribute_dofs (const FiniteElement<dim,spacedim> &fe);
389 
398  DEAL_II_DEPRECATED
399  virtual void distribute_mg_dofs (const FiniteElement<dim, spacedim> &fe);
400 
406  virtual void distribute_mg_dofs ();
407 
413  bool has_level_dofs() const;
414 
429  bool has_active_dofs() const;
430 
438 
442  virtual void clear ();
443 
473  void renumber_dofs (const std::vector<types::global_dof_index> &new_numbers);
474 
479  void renumber_dofs (const unsigned int level,
480  const std::vector<types::global_dof_index> &new_numbers);
481 
510  unsigned int max_couplings_between_dofs () const;
511 
524  unsigned int max_couplings_between_boundary_dofs () const;
525 
526  /*--------------------------------------*/
527 
532  /*
533  * @{
534  */
535 
539  cell_iterator begin (const unsigned int level = 0) const;
540 
553  active_cell_iterator begin_active(const unsigned int level = 0) const;
554 
559  cell_iterator end () const;
560 
565  cell_iterator end (const unsigned int level) const;
566 
572  active_cell_iterator end_active (const unsigned int level) const;
573 
574 
579  level_cell_iterator begin_mg (const unsigned int level = 0) const;
580 
585  level_cell_iterator end_mg (const unsigned int level) const;
586 
591  level_cell_iterator end_mg () const;
592 
608 
651 
664 
680  IteratorRange<cell_iterator> cell_iterators_on_level (const unsigned int level) const;
681 
698 
715  IteratorRange<level_cell_iterator> mg_cell_iterators_on_level (const unsigned int level) const;
716 
717  /*
718  * @}
719  */
720 
721 
722  /*---------------------------------------*/
723 
724 
749 
757  types::global_dof_index n_dofs (const unsigned int level) const;
758 
763 
776  template <typename number>
778  n_boundary_dofs (const std::map<types::boundary_id, const Function<spacedim,number>*> &boundary_ids) const;
779 
785  n_boundary_dofs (const std::set<types::boundary_id> &boundary_ids) const;
786 
802  const BlockInfo &block_info() const;
803 
804 
818  unsigned int n_locally_owned_dofs() const;
819 
825  const IndexSet &locally_owned_dofs() const;
826 
831  const IndexSet &locally_owned_mg_dofs(const unsigned int level) const;
832 
841  const std::vector<IndexSet> &
843 
856  const std::vector<types::global_dof_index> &
858 
866  const std::vector<IndexSet> &
867  locally_owned_mg_dofs_per_processor (const unsigned int level) const;
868 
875  get_fe (const unsigned int index=0) const;
876 
884 
889 
898  virtual std::size_t memory_consumption () const;
899 
904  template <class Archive>
905  void save (Archive &ar, const unsigned int version) const;
906 
911  template <class Archive>
912  void load (Archive &ar, const unsigned int version);
913 
914  BOOST_SERIALIZATION_SPLIT_MEMBER()
915 
916 
931  types::global_dof_index,
932  << "The given list of new dof indices is not consecutive: "
933  << "the index " << arg1 << " does not exist.");
939  int,
940  << "The given level " << arg1
941  << " is not in the valid range!");
952  int,
953  << "You tried to do something on level " << arg1
954  << ", but this level is empty.");
955 
956 
957 private:
958 
963 
967  SmartPointer<const Triangulation<dim,spacedim>,DoFHandler<dim,spacedim> >
969 
970 
975  std::unique_ptr<const hp::FECollection<dim,spacedim> > fe_collection;
976 
981  std::unique_ptr<::internal::DoFHandler::Policy::PolicyBase<dim,spacedim> > policy;
982 
991 
995  std::vector<::internal::DoFHandler::NumberCache> mg_number_cache;
996 
1004  {
1005  public:
1009  MGVertexDoFs ();
1010 
1016  void init (const unsigned int coarsest_level,
1017  const unsigned int finest_level,
1018  const unsigned int dofs_per_vertex);
1019 
1023  unsigned int get_coarsest_level () const;
1024 
1028  unsigned int get_finest_level () const;
1029 
1035  get_index (const unsigned int level,
1036  const unsigned int dof_number,
1037  const unsigned int dofs_per_vertex) const;
1038 
1043  void set_index (const unsigned int level,
1044  const unsigned int dof_number,
1045  const unsigned int dofs_per_vertex,
1046  const types::global_dof_index index);
1047 
1048  private:
1052  unsigned int coarsest_level;
1053 
1057  unsigned int finest_level;
1058 
1068  std::unique_ptr<types::global_dof_index[]> indices;
1069  };
1070 
1074  void clear_space ();
1075 
1079  void clear_mg_space ();
1080 
1081  template <int structdim>
1082  types::global_dof_index get_dof_index (const unsigned int obj_level,
1083  const unsigned int obj_index,
1084  const unsigned int fe_index,
1085  const unsigned int local_index) const;
1086 
1087  template <int structdim>
1088  void set_dof_index (const unsigned int obj_level,
1089  const unsigned int obj_index,
1090  const unsigned int fe_index,
1091  const unsigned int local_index,
1092  const types::global_dof_index global_index) const;
1093 
1097  std::vector<types::global_dof_index> vertex_dofs;
1098 
1103  std::vector<MGVertexDoFs> mg_vertex_dofs;
1104 
1109  std::vector<std::unique_ptr<::internal::DoFHandler::DoFLevel<dim> > > levels;
1110 
1111  std::vector<std::unique_ptr<::internal::DoFHandler::DoFLevel<dim> > > mg_levels;
1112 
1118  std::unique_ptr<::internal::DoFHandler::DoFFaces<dim> > faces;
1119 
1120  std::unique_ptr<::internal::DoFHandler::DoFFaces<dim> > mg_faces;
1121 
1125  template <int, class, bool> friend class DoFAccessor;
1126  template <class, bool> friend class DoFCellAccessor;
1127  friend struct ::internal::DoFAccessor::Implementation;
1128  friend struct ::internal::DoFCellAccessor::Implementation;
1129 
1130  friend struct ::internal::DoFHandler::Implementation;
1131  friend struct ::internal::DoFHandler::Policy::Implementation;
1132 
1133  // explicitly check for sensible template arguments, but not on windows
1134  // because MSVC creates bogus warnings during normal compilation
1135 #ifndef DEAL_II_MSVC
1136  static_assert (dim<=spacedim,
1137  "The dimension <dim> of a DoFHandler must be less than or "
1138  "equal to the space dimension <spacedim> in which it lives.");
1139 #endif
1140 };
1141 
1142 
1143 
1144 
1145 #ifndef DOXYGEN
1146 
1147 /* ----------------------- Inline functions ---------------------------------- */
1148 
1149 
1150 template <int dim, int spacedim>
1151 inline
1152 bool
1154 {
1155  return mg_number_cache.size()>0;
1156 }
1157 
1158 
1159 
1160 template <int dim, int spacedim>
1161 inline
1162 bool
1164 {
1165  return number_cache.n_global_dofs>0;
1166 }
1167 
1168 
1169 
1170 template <int dim, int spacedim>
1171 inline
1174 {
1175  return number_cache.n_global_dofs;
1176 }
1177 
1178 
1179 
1180 template <int dim, int spacedim>
1181 inline
1183 {
1184  Assert(has_level_dofs(), ExcMessage("n_dofs(level) can only be called after distribute_mg_dofs()"));
1185  Assert (level < mg_number_cache.size (), ExcInvalidLevel (level));
1186  return mg_number_cache[level].n_global_dofs;
1187 }
1188 
1189 
1190 
1191 template <int dim, int spacedim>
1192 unsigned int
1194 {
1195  return number_cache.n_locally_owned_dofs;
1196 }
1197 
1198 
1199 
1200 template <int dim, int spacedim>
1201 const IndexSet &
1203 {
1204  return number_cache.locally_owned_dofs;
1205 }
1206 
1207 
1208 
1209 template <int dim, int spacedim>
1210 const IndexSet &
1211 DoFHandler<dim, spacedim>::locally_owned_mg_dofs(const unsigned int level) const
1212 {
1213  Assert(level < this->get_triangulation().n_global_levels(),
1214  ExcMessage("invalid level in locally_owned_mg_dofs"));
1215  return mg_number_cache[level].locally_owned_dofs;
1216 }
1217 
1218 
1219 
1220 template <int dim, int spacedim>
1221 const std::vector<types::global_dof_index> &
1223 {
1224  return number_cache.n_locally_owned_dofs_per_processor;
1225 }
1226 
1227 
1228 
1229 template <int dim, int spacedim>
1230 const std::vector<IndexSet> &
1232 {
1233  return number_cache.locally_owned_dofs_per_processor;
1234 }
1235 
1236 
1237 
1238 template <int dim, int spacedim>
1239 const std::vector<IndexSet> &
1241 {
1242  Assert(level < this->get_triangulation().n_global_levels(), ExcMessage("invalid level in locally_owned_mg_dofs_per_processor"));
1243  return mg_number_cache[level].locally_owned_dofs_per_processor;
1244 }
1245 
1246 
1247 
1248 template <int dim, int spacedim>
1249 inline
1252 (const unsigned int index) const
1253 {
1254  (void) index;
1255  Assert(index == 0, ExcMessage("There is only one FiniteElement stored. The index must be zero!"));
1256  return get_fe_collection()[0];
1257 }
1258 
1259 
1260 
1261 template <int dim, int spacedim>
1262 inline
1265 {
1266  Assert(fe_collection != nullptr,
1267  ExcMessage("You are trying to access the DoFHandler's FECollection object before it has been initialized."));
1268  return *fe_collection;
1269 }
1270 
1271 
1272 
1273 template <int dim, int spacedim>
1274 inline
1277 {
1278  Assert(tria != nullptr,
1279  ExcMessage("This DoFHandler object has not been associated "
1280  "with a triangulation."));
1281  return *tria;
1282 }
1283 
1284 
1285 
1286 template <int dim, int spacedim>
1287 inline
1288 const BlockInfo &
1290 {
1291  return block_info_object;
1292 }
1293 
1294 
1295 
1296 template <int dim, int spacedim>
1297 template <typename number>
1300 {
1301  // extract the set of boundary ids and forget about the function object pointers
1302  std::set<types::boundary_id> boundary_ids_only;
1303  for (typename std::map<types::boundary_id, const Function<spacedim,number>*>::const_iterator
1304  p = boundary_ids.begin();
1305  p != boundary_ids.end(); ++p)
1306  boundary_ids_only.insert (p->first);
1307 
1308  // then just hand everything over to the other function that does the work
1309  return n_boundary_dofs(boundary_ids_only);
1310 }
1311 
1312 
1313 
1314 namespace internal
1315 {
1323  template <int dim, int spacedim>
1324  std::string policy_to_string(const ::internal::DoFHandler::Policy::PolicyBase<dim,spacedim> &policy);
1325 }
1326 
1327 
1328 
1329 template <int dim, int spacedim>
1330 template <class Archive>
1331 void DoFHandler<dim,spacedim>::save (Archive &ar,
1332  const unsigned int) const
1333 {
1334  ar &block_info_object;
1335  ar &vertex_dofs;
1336  ar &number_cache;
1337 
1338  // some versions of gcc have trouble with loading vectors of
1339  // std::unique_ptr objects because std::unique_ptr does not
1340  // have a copy constructor. do it one level at a time
1341  unsigned int n_levels = levels.size();
1342  ar &n_levels;
1343  for (unsigned int i = 0; i < levels.size(); ++i)
1344  ar &levels[i];
1345 
1346  // boost dereferences a nullptr when serializing a nullptr
1347  // at least up to 1.65.1. This causes problems with clang-5.
1348  // Therefore, work around it.
1349  bool faces_is_nullptr = (faces.get()==nullptr);
1350  ar &faces_is_nullptr;
1351  if (!faces_is_nullptr)
1352  ar &faces;
1353 
1354  // write out the number of triangulation cells and later check during
1355  // loading that this number is indeed correct; same with something that
1356  // identifies the FE and the policy
1357  unsigned int n_cells = tria->n_cells();
1358  std::string fe_name = this->get_fe(0).get_name();
1359  std::string policy_name = internal::policy_to_string(*policy);
1360 
1361  ar &n_cells &fe_name &policy_name;
1362 }
1363 
1364 
1365 
1366 template <int dim, int spacedim>
1367 template <class Archive>
1368 void DoFHandler<dim,spacedim>::load (Archive &ar,
1369  const unsigned int)
1370 {
1371  ar &block_info_object;
1372  ar &vertex_dofs;
1373  ar &number_cache;
1374 
1375  // boost::serialization can restore pointers just fine, but if the
1376  // pointer object still points to something useful, that object is not
1377  // destroyed and we end up with a memory leak. consequently, first delete
1378  // previous content before re-loading stuff
1379  levels.clear();
1380  faces.reset();
1381 
1382  // some versions of gcc have trouble with loading vectors of
1383  // std::unique_ptr objects because std::unique_ptr does not
1384  // have a copy constructor. do it one level at a time
1385  unsigned int size;
1386  ar &size;
1387  levels.resize(size);
1388  for (unsigned int i = 0; i < levels.size(); ++i)
1389  {
1390  std::unique_ptr<internal::DoFHandler::DoFLevel<dim>> level;
1391  ar &level;
1392  levels[i] = std::move(level);
1393  }
1394 
1395  //Workaround for nullptr, see in save().
1396  bool faces_is_nullptr = true;
1397  ar &faces_is_nullptr;
1398  if (!faces_is_nullptr)
1399  ar &faces;
1400 
1401  // these are the checks that correspond to the last block in the save()
1402  // function
1403  unsigned int n_cells;
1404  std::string fe_name;
1405  std::string policy_name;
1406 
1407  ar &n_cells &fe_name &policy_name;
1408 
1409  AssertThrow (n_cells == tria->n_cells(),
1410  ExcMessage ("The object being loaded into does not match the triangulation "
1411  "that has been stored previously."));
1412  AssertThrow (fe_name == this->get_fe(0).get_name(),
1413  ExcMessage ("The finite element associated with this DoFHandler does not match "
1414  "the one that was associated with the DoFHandler previously stored."));
1415  AssertThrow (policy_name == internal::policy_to_string(*policy),
1416  ExcMessage ("The policy currently associated with this DoFHandler ("
1417  + internal::policy_to_string(*policy)
1418  + ") does not match the one that was associated with the "
1419  "DoFHandler previously stored ("
1420  + policy_name
1421  + ")."));
1422 }
1423 
1424 
1425 
1426 
1427 template <int dim, int spacedim>
1428 inline
1430 DoFHandler<dim, spacedim>::MGVertexDoFs::get_index (const unsigned int level,
1431  const unsigned int dof_number,
1432  const unsigned int dofs_per_vertex) const
1433 {
1434  Assert ((level >= coarsest_level) && (level <= finest_level), ExcInvalidLevel (level));
1435  return indices[dofs_per_vertex * (level - coarsest_level) + dof_number];
1436 }
1437 
1438 
1439 
1440 template <int dim, int spacedim>
1441 inline
1442 void
1443 DoFHandler<dim, spacedim>::MGVertexDoFs::set_index (const unsigned int level,
1444  const unsigned int dof_number,
1445  const unsigned int dofs_per_vertex,
1446  const types::global_dof_index index)
1447 {
1448  Assert ((level >= coarsest_level) && (level <= finest_level), ExcInvalidLevel (level));
1449  indices[dofs_per_vertex * (level - coarsest_level) + dof_number] = index;
1450 }
1451 
1452 
1453 #endif // DOXYGEN
1454 
1455 DEAL_II_NAMESPACE_CLOSE
1456 
1457 #endif
DoFHandler & operator=(const DoFHandler &)=delete
unsigned int max_couplings_between_boundary_dofs() const
std::vector< MGVertexDoFs > mg_vertex_dofs
Definition: dof_handler.h:1103
active_cell_iterator end_active(const unsigned int level) const
Definition: dof_handler.cc:772
active_cell_iterator begin_active(const unsigned int level=0) const
Definition: dof_handler.cc:734
Iterator points to a valid object.
std::vector<::internal::DoFHandler::NumberCache > mg_number_cache
Definition: dof_handler.h:995
static const unsigned int dimension
Definition: dof_handler.h:300
const BlockInfo & block_info() const
static const unsigned int spacedim
static::ExceptionBase & ExcInvalidLevel(int arg1)
level_cell_iterator begin_mg(const unsigned int level=0) const
Definition: dof_handler.cc:784
unsigned int boundary_id
Definition: types.h:110
virtual void clear()
const Triangulation< dim, spacedim > & get_triangulation() const
const FiniteElement< dim, spacedim > & get_fe(const unsigned int index=0) const
void clear_mg_space()
ActiveSelector::active_cell_iterator active_cell_iterator
Definition: dof_handler.h:229
SmartPointer< const Triangulation< dim, spacedim >, DoFHandler< dim, spacedim > > tria
Definition: dof_handler.h:968
STL namespace.
bool has_level_dofs() const
#define AssertThrow(cond, exc)
Definition: exceptions.h:398
const std::vector< types::global_dof_index > & n_locally_owned_dofs_per_processor() const
::internal::DoFHandler::NumberCache number_cache
Definition: dof_handler.h:990
unsigned int max_couplings_between_dofs() const
void initialize(const Triangulation< dim, spacedim > &tria, const FiniteElement< dim, spacedim > &fe)
Definition: dof_handler.cc:698
std::unique_ptr<::internal::DoFHandler::Policy::PolicyBase< dim, spacedim > > policy
Definition: dof_handler.h:981
static::ExceptionBase & ExcInvalidBoundaryIndicator()
ActiveSelector::face_iterator face_iterator
Definition: dof_handler.h:275
BlockInfo block_info_object
Definition: dof_handler.h:962
IteratorRange< active_cell_iterator > active_cell_iterators_on_level(const unsigned int level) const
Definition: dof_handler.cc:863
static::ExceptionBase & ExcFacesHaveNoLevel()
static::ExceptionBase & ExcMessage(std::string arg1)
level_cell_iterator end_mg() const
Definition: dof_handler.cc:810
ActiveSelector::cell_iterator cell_iterator
Definition: dof_handler.h:257
#define DeclException1(Exception1, type1, outsequence)
Definition: exceptions.h:593
Definition: types.h:30
unsigned int global_dof_index
Definition: types.h:88
std::vector< std::unique_ptr<::internal::DoFHandler::DoFLevel< dim > > > levels
Definition: dof_handler.h:1109
virtual void distribute_mg_dofs()
#define Assert(cond, exc)
Definition: exceptions.h:337
types::global_dof_index n_dofs() const
static const unsigned int default_fe_index
Definition: dof_handler.h:330
std::unique_ptr< const hp::FECollection< dim, spacedim > > fe_collection
Definition: dof_handler.h:975
const Triangulation< dim, spacedim > * tria
unsigned int n_locally_owned_dofs() const
static::ExceptionBase & ExcGridsDoNotMatch()
#define DeclException0(Exception0)
Definition: exceptions.h:570
void set_index(const unsigned int level, const unsigned int dof_number, const unsigned int dofs_per_vertex, const types::global_dof_index index)
cell_iterator end() const
Definition: dof_handler.cc:750
Definition: hp.h:102
virtual ~DoFHandler()
Definition: dof_handler.cc:683
IteratorRange< cell_iterator > cell_iterators_on_level(const unsigned int level) const
Definition: dof_handler.cc:852
const IndexSet & locally_owned_mg_dofs(const unsigned int level) const
IteratorRange< cell_iterator > cell_iterators() const
Definition: dof_handler.cc:819
void renumber_dofs(const std::vector< types::global_dof_index > &new_numbers)
std::unique_ptr< types::global_dof_index[]> indices
Definition: dof_handler.h:1068
void initialize_local_block_info()
bool has_active_dofs() const
const std::vector< IndexSet > & locally_owned_mg_dofs_per_processor(const unsigned int level) const
static::ExceptionBase & ExcEmptyLevel(int arg1)
static const unsigned int dim
cell_iterator begin(const unsigned int level=0) const
Definition: dof_handler.cc:722
void save(Archive &ar, const unsigned int version) const
types::global_dof_index get_index(const unsigned int level, const unsigned int dof_number, const unsigned int dofs_per_vertex) const
IteratorRange< level_cell_iterator > mg_cell_iterators() const
Definition: dof_handler.cc:840
types::global_dof_index n_boundary_dofs() const
Definition: dof_handler.cc:888
unsigned int coarsest_level
Definition: dof_handler.h:1052
void clear_space()
IteratorRange< level_cell_iterator > mg_cell_iterators_on_level(const unsigned int level) const
Definition: dof_handler.cc:874
int index() const
const hp::FECollection< dim, spacedim > & get_fe_collection() const
virtual void distribute_dofs(const FiniteElement< dim, spacedim > &fe)
Definition: dof_handler.cc:992
void load(Archive &ar, const unsigned int version)
A small class collecting the different BlockIndices involved in global, multilevel and local computat...
Definition: block_info.h:93
ActiveSelector::active_face_iterator active_face_iterator
Definition: dof_handler.h:288
static::ExceptionBase & ExcNewNumbersNotConsecutive(types::global_dof_index arg1)
static const unsigned int space_dimension
Definition: dof_handler.h:305
virtual std::size_t memory_consumption() const
Definition: dof_handler.cc:963
const types::global_dof_index invalid_dof_index
Definition: types.h:187
const IndexSet & locally_owned_dofs() const
const std::vector< IndexSet > & locally_owned_dofs_per_processor() const
const Triangulation< dim, spacedim > & get_triangulation() const
IteratorRange< active_cell_iterator > active_cell_iterators() const
Definition: dof_handler.cc:829
std::unique_ptr<::internal::DoFHandler::DoFFaces< dim > > faces
Definition: dof_handler.h:1118
static const types::global_dof_index invalid_dof_index
Definition: dof_handler.h:319
std::vector< types::global_dof_index > vertex_dofs
Definition: dof_handler.h:1097
int level() const