/usr/include/deal.II/dofs/dof_levels.h is in libdeal.ii-dev 8.5.1-3.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 | // ---------------------------------------------------------------------
//
// Copyright (C) 1998 - 2015 by the deal.II authors
//
// This file is part of the deal.II library.
//
// The deal.II library is free software; you can use it, redistribute
// it, and/or modify it under the terms of the GNU Lesser General
// Public License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.
// The full text of the license can be found in the file LICENSE at
// the top level of the deal.II distribution.
//
// ---------------------------------------------------------------------
#ifndef dealii__dof_levels_h
#define dealii__dof_levels_h
#include <deal.II/base/config.h>
#include <deal.II/base/exceptions.h>
#include <deal.II/base/memory_consumption.h>
#include <deal.II/dofs/dof_objects.h>
#include <vector>
DEAL_II_NAMESPACE_OPEN
namespace internal
{
namespace DoFHandler
{
/**
* Structure for storing degree of freedom information for cells,
* organized by levels.
*
* We store are cached values for the DoF indices on each cell
* in#cell_dof_indices_cache, since this is a frequently requested
* operation. The values are set by
* DoFCellAccessor::update_cell_dof_indices_cache and are used by
* DoFCellAccessor::get_dof_indices.
*
* Note that vertices are separate from, and in fact have nothing to do
* with cells. The indices of degrees of freedom located on vertices
* therefore are not stored here, but rather in member variables of the
* dealii::DoFHandler class.
*
* The indices of degrees of freedom located on lower dimensional objects,
* i.e. on lines for 2D and on quads and lines for 3D are treated
* similarly than that on cells. However, these geometrical objects, which
* are called faces as a generalisation, are not organised in a
* hierarchical structure of levels. Therefore, the degrees of freedom
* located on these objects are stored in separate classes, namely the
* <tt>DoFFaces</tt> classes.
*
* Access to this object is usually through the
* DoFAccessor::set_dof_index() and DoFAccessor::dof_index() functions or
* similar functions of derived classes that in turn access the member
* variables using the DoFHandler::get_dof_index() and corresponding
* setter functions. Knowledge of the actual data format is therefore
* encapsulated to the present hierarchy of classes as well as the
* dealii::DoFHandler class.
*
* @author Wolfgang Bangerth, 1998, 2006, Guido Kanschat, 2012
*/
template <int dim>
class DoFLevel
{
public:
/**
* Cache for the DoF indices on cells. The size of this array equals the
* number of cells on a given level times selected_fe.dofs_per_cell.
*/
std::vector<types::global_dof_index> cell_dof_indices_cache;
/**
* The object containing dof-indices and related access-functions
*/
DoFObjects<dim> dof_object;
/**
* Return a pointer to the beginning of the DoF indices cache for a
* given cell.
*
* @param obj_index The number of the cell we are looking at.
* @param dofs_per_cell The number of DoFs per cell for this cell.
* @return A pointer to the first DoF index for the current cell. The
* next dofs_per_cell indices are for the current cell.
*/
const types::global_dof_index *
get_cell_cache_start (const unsigned int obj_index,
const unsigned int dofs_per_cell) const;
/**
* Determine an estimate for the memory consumption (in bytes) of this
* object.
*/
std::size_t memory_consumption () const;
/**
* Read or write the data of this object to or from a stream for the
* purpose of serialization
*/
template <class Archive>
void serialize(Archive &ar,
const unsigned int version);
};
template <int dim>
inline
const types::global_dof_index *
DoFLevel<dim>::get_cell_cache_start (const unsigned int obj_index,
const unsigned int dofs_per_cell) const
{
Assert (obj_index*dofs_per_cell+dofs_per_cell
<=
cell_dof_indices_cache.size(),
ExcInternalError());
return &cell_dof_indices_cache[obj_index*dofs_per_cell];
}
template <int dim>
inline
std::size_t
DoFLevel<dim>::memory_consumption () const
{
return (MemoryConsumption::memory_consumption (cell_dof_indices_cache) +
MemoryConsumption::memory_consumption (dof_object));
}
template <int dim>
template <class Archive>
inline
void
DoFLevel<dim>::serialize (Archive &ar,
const unsigned int)
{
ar &cell_dof_indices_cache;
ar &dof_object;
}
}
}
DEAL_II_NAMESPACE_CLOSE
#endif
|