Reference documentation for deal.II version Git 58efed1 2018-02-20 08:43:45 -0500
dof_handler.h
1 // ---------------------------------------------------------------------
2 //
3 // Copyright (C) 1998 - 2018 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 
410  virtual void distribute_dofs (const FiniteElement<dim,spacedim> &fe);
411 
420  DEAL_II_DEPRECATED
421  virtual void distribute_mg_dofs (const FiniteElement<dim, spacedim> &fe);
422 
428  virtual void distribute_mg_dofs ();
429 
435  bool has_level_dofs() const;
436 
451  bool has_active_dofs() const;
452 
460 
464  virtual void clear ();
465 
518  void renumber_dofs (const std::vector<types::global_dof_index> &new_numbers);
519 
524  void renumber_dofs (const unsigned int level,
525  const std::vector<types::global_dof_index> &new_numbers);
526 
555  unsigned int max_couplings_between_dofs () const;
556 
569  unsigned int max_couplings_between_boundary_dofs () const;
570 
571  /*--------------------------------------*/
572 
577  /*
578  * @{
579  */
580 
584  cell_iterator begin (const unsigned int level = 0) const;
585 
598  active_cell_iterator begin_active(const unsigned int level = 0) const;
599 
604  cell_iterator end () const;
605 
610  cell_iterator end (const unsigned int level) const;
611 
617  active_cell_iterator end_active (const unsigned int level) const;
618 
619 
624  level_cell_iterator begin_mg (const unsigned int level = 0) const;
625 
630  level_cell_iterator end_mg (const unsigned int level) const;
631 
636  level_cell_iterator end_mg () const;
637 
653 
696 
709 
725  IteratorRange<cell_iterator> cell_iterators_on_level (const unsigned int level) const;
726 
743 
760  IteratorRange<level_cell_iterator> mg_cell_iterators_on_level (const unsigned int level) const;
761 
762  /*
763  * @}
764  */
765 
766 
767  /*---------------------------------------*/
768 
769 
794 
802  types::global_dof_index n_dofs (const unsigned int level) const;
803 
808 
821  template <typename number>
823  n_boundary_dofs (const std::map<types::boundary_id, const Function<spacedim,number>*> &boundary_ids) const;
824 
830  n_boundary_dofs (const std::set<types::boundary_id> &boundary_ids) const;
831 
847  const BlockInfo &block_info() const;
848 
849 
869  unsigned int n_locally_owned_dofs() const;
870 
876  const IndexSet &locally_owned_dofs() const;
877 
882  const IndexSet &locally_owned_mg_dofs(const unsigned int level) const;
883 
897  const std::vector<IndexSet> &
899 
916  const std::vector<types::global_dof_index> &
918 
931  const std::vector<IndexSet> &
932  locally_owned_mg_dofs_per_processor (const unsigned int level) const;
933 
940  get_fe (const unsigned int index=0) const;
941 
949 
954 
963  virtual std::size_t memory_consumption () const;
964 
969  template <class Archive>
970  void save (Archive &ar, const unsigned int version) const;
971 
976  template <class Archive>
977  void load (Archive &ar, const unsigned int version);
978 
979  BOOST_SERIALIZATION_SPLIT_MEMBER()
980 
981 
996  types::global_dof_index,
997  << "The given list of new dof indices is not consecutive: "
998  << "the index " << arg1 << " does not exist.");
1004  int,
1005  << "The given level " << arg1
1006  << " is not in the valid range!");
1017  int,
1018  << "You tried to do something on level " << arg1
1019  << ", but this level is empty.");
1020 
1021 
1022 private:
1023 
1028 
1032  SmartPointer<const Triangulation<dim,spacedim>,DoFHandler<dim,spacedim> >
1034 
1035 
1040  std::unique_ptr<const hp::FECollection<dim,spacedim> > fe_collection;
1041 
1046  std::unique_ptr<::internal::DoFHandler::Policy::PolicyBase<dim,spacedim> > policy;
1047 
1056 
1060  std::vector<::internal::DoFHandler::NumberCache> mg_number_cache;
1061 
1069  {
1070  public:
1074  MGVertexDoFs ();
1075 
1081  void init (const unsigned int coarsest_level,
1082  const unsigned int finest_level,
1083  const unsigned int dofs_per_vertex);
1084 
1088  unsigned int get_coarsest_level () const;
1089 
1093  unsigned int get_finest_level () const;
1094 
1100  get_index (const unsigned int level,
1101  const unsigned int dof_number,
1102  const unsigned int dofs_per_vertex) const;
1103 
1108  void set_index (const unsigned int level,
1109  const unsigned int dof_number,
1110  const unsigned int dofs_per_vertex,
1111  const types::global_dof_index index);
1112 
1113  private:
1117  unsigned int coarsest_level;
1118 
1122  unsigned int finest_level;
1123 
1133  std::unique_ptr<types::global_dof_index[]> indices;
1134  };
1135 
1139  void clear_space ();
1140 
1144  void clear_mg_space ();
1145 
1146  template <int structdim>
1147  types::global_dof_index get_dof_index (const unsigned int obj_level,
1148  const unsigned int obj_index,
1149  const unsigned int fe_index,
1150  const unsigned int local_index) const;
1151 
1152  template <int structdim>
1153  void set_dof_index (const unsigned int obj_level,
1154  const unsigned int obj_index,
1155  const unsigned int fe_index,
1156  const unsigned int local_index,
1157  const types::global_dof_index global_index) const;
1158 
1162  std::vector<types::global_dof_index> vertex_dofs;
1163 
1168  std::vector<MGVertexDoFs> mg_vertex_dofs;
1169 
1174  std::vector<std::unique_ptr<::internal::DoFHandler::DoFLevel<dim> > > levels;
1175 
1176  std::vector<std::unique_ptr<::internal::DoFHandler::DoFLevel<dim> > > mg_levels;
1177 
1183  std::unique_ptr<::internal::DoFHandler::DoFFaces<dim> > faces;
1184 
1185  std::unique_ptr<::internal::DoFHandler::DoFFaces<dim> > mg_faces;
1186 
1190  template <int, class, bool> friend class DoFAccessor;
1191  template <class, bool> friend class DoFCellAccessor;
1192  friend struct ::internal::DoFAccessor::Implementation;
1193  friend struct ::internal::DoFCellAccessor::Implementation;
1194 
1195  friend struct ::internal::DoFHandler::Implementation;
1196  friend struct ::internal::DoFHandler::Policy::Implementation;
1197 
1198  // explicitly check for sensible template arguments, but not on windows
1199  // because MSVC creates bogus warnings during normal compilation
1200 #ifndef DEAL_II_MSVC
1201  static_assert (dim<=spacedim,
1202  "The dimension <dim> of a DoFHandler must be less than or "
1203  "equal to the space dimension <spacedim> in which it lives.");
1204 #endif
1205 };
1206 
1207 
1208 
1209 
1210 #ifndef DOXYGEN
1211 
1212 /* ----------------------- Inline functions ---------------------------------- */
1213 
1214 
1215 template <int dim, int spacedim>
1216 inline
1217 bool
1219 {
1220  return mg_number_cache.size()>0;
1221 }
1222 
1223 
1224 
1225 template <int dim, int spacedim>
1226 inline
1227 bool
1229 {
1230  return number_cache.n_global_dofs>0;
1231 }
1232 
1233 
1234 
1235 template <int dim, int spacedim>
1236 inline
1239 {
1240  return number_cache.n_global_dofs;
1241 }
1242 
1243 
1244 
1245 template <int dim, int spacedim>
1246 inline
1248 {
1249  Assert(has_level_dofs(), ExcMessage("n_dofs(level) can only be called after distribute_mg_dofs()"));
1250  Assert (level < mg_number_cache.size (), ExcInvalidLevel (level));
1251  return mg_number_cache[level].n_global_dofs;
1252 }
1253 
1254 
1255 
1256 template <int dim, int spacedim>
1257 unsigned int
1259 {
1260  return number_cache.n_locally_owned_dofs;
1261 }
1262 
1263 
1264 
1265 template <int dim, int spacedim>
1266 const IndexSet &
1268 {
1269  return number_cache.locally_owned_dofs;
1270 }
1271 
1272 
1273 
1274 template <int dim, int spacedim>
1275 const IndexSet &
1276 DoFHandler<dim, spacedim>::locally_owned_mg_dofs(const unsigned int level) const
1277 {
1278  Assert(level < this->get_triangulation().n_global_levels(),
1279  ExcMessage("invalid level in locally_owned_mg_dofs"));
1280  return mg_number_cache[level].locally_owned_dofs;
1281 }
1282 
1283 
1284 
1285 template <int dim, int spacedim>
1286 const std::vector<types::global_dof_index> &
1288 {
1289  return number_cache.n_locally_owned_dofs_per_processor;
1290 }
1291 
1292 
1293 
1294 template <int dim, int spacedim>
1295 const std::vector<IndexSet> &
1297 {
1298  return number_cache.locally_owned_dofs_per_processor;
1299 }
1300 
1301 
1302 
1303 template <int dim, int spacedim>
1304 const std::vector<IndexSet> &
1306 {
1307  Assert(level < this->get_triangulation().n_global_levels(), ExcMessage("invalid level in locally_owned_mg_dofs_per_processor"));
1308  return mg_number_cache[level].locally_owned_dofs_per_processor;
1309 }
1310 
1311 
1312 
1313 template <int dim, int spacedim>
1314 inline
1317 (const unsigned int index) const
1318 {
1319  (void) index;
1320  Assert(index == 0, ExcMessage("There is only one FiniteElement stored. The index must be zero!"));
1321  return get_fe_collection()[0];
1322 }
1323 
1324 
1325 
1326 template <int dim, int spacedim>
1327 inline
1330 {
1331  Assert(fe_collection != nullptr,
1332  ExcMessage("You are trying to access the DoFHandler's FECollection object before it has been initialized."));
1333  return *fe_collection;
1334 }
1335 
1336 
1337 
1338 template <int dim, int spacedim>
1339 inline
1342 {
1343  Assert(tria != nullptr,
1344  ExcMessage("This DoFHandler object has not been associated "
1345  "with a triangulation."));
1346  return *tria;
1347 }
1348 
1349 
1350 
1351 template <int dim, int spacedim>
1352 inline
1353 const BlockInfo &
1355 {
1356  return block_info_object;
1357 }
1358 
1359 
1360 
1361 template <int dim, int spacedim>
1362 template <typename number>
1365 {
1366  // extract the set of boundary ids and forget about the function object pointers
1367  std::set<types::boundary_id> boundary_ids_only;
1368  for (typename std::map<types::boundary_id, const Function<spacedim,number>*>::const_iterator
1369  p = boundary_ids.begin();
1370  p != boundary_ids.end(); ++p)
1371  boundary_ids_only.insert (p->first);
1372 
1373  // then just hand everything over to the other function that does the work
1374  return n_boundary_dofs(boundary_ids_only);
1375 }
1376 
1377 
1378 
1379 namespace internal
1380 {
1388  template <int dim, int spacedim>
1389  std::string policy_to_string(const ::internal::DoFHandler::Policy::PolicyBase<dim,spacedim> &policy);
1390 }
1391 
1392 
1393 
1394 template <int dim, int spacedim>
1395 template <class Archive>
1396 void DoFHandler<dim,spacedim>::save (Archive &ar,
1397  const unsigned int) const
1398 {
1399  ar &block_info_object;
1400  ar &vertex_dofs;
1401  ar &number_cache;
1402 
1403  // some versions of gcc have trouble with loading vectors of
1404  // std::unique_ptr objects because std::unique_ptr does not
1405  // have a copy constructor. do it one level at a time
1406  unsigned int n_levels = levels.size();
1407  ar &n_levels;
1408  for (unsigned int i = 0; i < levels.size(); ++i)
1409  ar &levels[i];
1410 
1411  // boost dereferences a nullptr when serializing a nullptr
1412  // at least up to 1.65.1. This causes problems with clang-5.
1413  // Therefore, work around it.
1414  bool faces_is_nullptr = (faces.get()==nullptr);
1415  ar &faces_is_nullptr;
1416  if (!faces_is_nullptr)
1417  ar &faces;
1418 
1419  // write out the number of triangulation cells and later check during
1420  // loading that this number is indeed correct; same with something that
1421  // identifies the FE and the policy
1422  unsigned int n_cells = tria->n_cells();
1423  std::string fe_name = this->get_fe(0).get_name();
1424  std::string policy_name = internal::policy_to_string(*policy);
1425 
1426  ar &n_cells &fe_name &policy_name;
1427 }
1428 
1429 
1430 
1431 template <int dim, int spacedim>
1432 template <class Archive>
1433 void DoFHandler<dim,spacedim>::load (Archive &ar,
1434  const unsigned int)
1435 {
1436  ar &block_info_object;
1437  ar &vertex_dofs;
1438  ar &number_cache;
1439 
1440  // boost::serialization can restore pointers just fine, but if the
1441  // pointer object still points to something useful, that object is not
1442  // destroyed and we end up with a memory leak. consequently, first delete
1443  // previous content before re-loading stuff
1444  levels.clear();
1445  faces.reset();
1446 
1447  // some versions of gcc have trouble with loading vectors of
1448  // std::unique_ptr objects because std::unique_ptr does not
1449  // have a copy constructor. do it one level at a time
1450  unsigned int size;
1451  ar &size;
1452  levels.resize(size);
1453  for (unsigned int i = 0; i < levels.size(); ++i)
1454  {
1455  std::unique_ptr<internal::DoFHandler::DoFLevel<dim>> level;
1456  ar &level;
1457  levels[i] = std::move(level);
1458  }
1459 
1460  //Workaround for nullptr, see in save().
1461  bool faces_is_nullptr = true;
1462  ar &faces_is_nullptr;
1463  if (!faces_is_nullptr)
1464  ar &faces;
1465 
1466  // these are the checks that correspond to the last block in the save()
1467  // function
1468  unsigned int n_cells;
1469  std::string fe_name;
1470  std::string policy_name;
1471 
1472  ar &n_cells &fe_name &policy_name;
1473 
1474  AssertThrow (n_cells == tria->n_cells(),
1475  ExcMessage ("The object being loaded into does not match the triangulation "
1476  "that has been stored previously."));
1477  AssertThrow (fe_name == this->get_fe(0).get_name(),
1478  ExcMessage ("The finite element associated with this DoFHandler does not match "
1479  "the one that was associated with the DoFHandler previously stored."));
1480  AssertThrow (policy_name == internal::policy_to_string(*policy),
1481  ExcMessage ("The policy currently associated with this DoFHandler ("
1482  + internal::policy_to_string(*policy)
1483  + ") does not match the one that was associated with the "
1484  "DoFHandler previously stored ("
1485  + policy_name
1486  + ")."));
1487 }
1488 
1489 
1490 
1491 
1492 template <int dim, int spacedim>
1493 inline
1495 DoFHandler<dim, spacedim>::MGVertexDoFs::get_index (const unsigned int level,
1496  const unsigned int dof_number,
1497  const unsigned int dofs_per_vertex) const
1498 {
1499  Assert ((level >= coarsest_level) && (level <= finest_level), ExcInvalidLevel (level));
1500  return indices[dofs_per_vertex * (level - coarsest_level) + dof_number];
1501 }
1502 
1503 
1504 
1505 template <int dim, int spacedim>
1506 inline
1507 void
1508 DoFHandler<dim, spacedim>::MGVertexDoFs::set_index (const unsigned int level,
1509  const unsigned int dof_number,
1510  const unsigned int dofs_per_vertex,
1511  const types::global_dof_index index)
1512 {
1513  Assert ((level >= coarsest_level) && (level <= finest_level), ExcInvalidLevel (level));
1514  indices[dofs_per_vertex * (level - coarsest_level) + dof_number] = index;
1515 }
1516 
1517 
1518 #endif // DOXYGEN
1519 
1520 DEAL_II_NAMESPACE_CLOSE
1521 
1522 #endif
DoFHandler & operator=(const DoFHandler &)=delete
unsigned int max_couplings_between_boundary_dofs() const
std::vector< MGVertexDoFs > mg_vertex_dofs
Definition: dof_handler.h:1168
active_cell_iterator end_active(const unsigned int level) const
Definition: dof_handler.cc:775
active_cell_iterator begin_active(const unsigned int level=0) const
Definition: dof_handler.cc:737
Iterator points to a valid object.
std::vector<::internal::DoFHandler::NumberCache > mg_number_cache
Definition: dof_handler.h:1060
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:787
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:1033
STL namespace.
bool has_level_dofs() const
#define AssertThrow(cond, exc)
Definition: exceptions.h:410
const std::vector< types::global_dof_index > & n_locally_owned_dofs_per_processor() const
::internal::DoFHandler::NumberCache number_cache
Definition: dof_handler.h:1055
unsigned int max_couplings_between_dofs() const
void initialize(const Triangulation< dim, spacedim > &tria, const FiniteElement< dim, spacedim > &fe)
Definition: dof_handler.cc:701
std::unique_ptr<::internal::DoFHandler::Policy::PolicyBase< dim, spacedim > > policy
Definition: dof_handler.h:1046
static::ExceptionBase & ExcInvalidBoundaryIndicator()
ActiveSelector::face_iterator face_iterator
Definition: dof_handler.h:275
BlockInfo block_info_object
Definition: dof_handler.h:1027
IteratorRange< active_cell_iterator > active_cell_iterators_on_level(const unsigned int level) const
Definition: dof_handler.cc:866
static::ExceptionBase & ExcFacesHaveNoLevel()
static::ExceptionBase & ExcMessage(std::string arg1)
level_cell_iterator end_mg() const
Definition: dof_handler.cc:813
ActiveSelector::cell_iterator cell_iterator
Definition: dof_handler.h:257
#define DeclException1(Exception1, type1, outsequence)
Definition: exceptions.h:605
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:1174
virtual void distribute_mg_dofs()
#define Assert(cond, exc)
Definition: exceptions.h:349
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:1040
const Triangulation< dim, spacedim > * tria
unsigned int n_locally_owned_dofs() const
static::ExceptionBase & ExcGridsDoNotMatch()
#define DeclException0(Exception0)
Definition: exceptions.h:582
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:753
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:855
const IndexSet & locally_owned_mg_dofs(const unsigned int level) const
IteratorRange< cell_iterator > cell_iterators() const
Definition: dof_handler.cc:822
void renumber_dofs(const std::vector< types::global_dof_index > &new_numbers)
std::unique_ptr< types::global_dof_index[]> indices
Definition: dof_handler.h:1133
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:725
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:843
types::global_dof_index n_boundary_dofs() const
Definition: dof_handler.cc:891
unsigned int coarsest_level
Definition: dof_handler.h:1117
void clear_space()
IteratorRange< level_cell_iterator > mg_cell_iterators_on_level(const unsigned int level) const
Definition: dof_handler.cc:877
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:995
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:966
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:832
std::unique_ptr<::internal::DoFHandler::DoFFaces< dim > > faces
Definition: dof_handler.h:1183
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:1162
int level() const