/usr/include/deal.II/multigrid/mg_constraints.h is in libdeal.ii-dev 6.3.1-1.1.
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 | //---------------------------------------------------------------------------
// $Id: mg_constraints.h 20853 2010-03-19 15:39:54Z bangerth $
// Version: $Name$
//
// Copyright (C) 2010 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
// to the file deal.II/doc/license.html for the text and
// further information on this license.
//
//---------------------------------------------------------------------------
#ifndef __deal2__mg_constraints_h
#define __deal2__mg_constraints_h
#include <base/config.h>
#include <base/subscriptor.h>
#include <vector>
#include <set>
DEAL_II_NAMESPACE_OPEN
template <int dim, int spacedim> class MGDoFHandler;
template <int dim> class FunctionMap;
/**
* Collection of boundary constraints and refinement edge constraints
* for level vectors.
*/
class MGConstraints : public Subscriptor
{
public:
/**
* Fill the internal data
* structures with values
* extracted from the dof
* handler.
*
* This function leaves
* #boundary_indices empty, since
* no boundary values are
* provided.
*/
template <int dim, int spacedim>
void initialize(const MGDoFHandler<dim,spacedim>& dof);
/**
* Fill the internal data
* structures with values
* extracted from the dof
* handler, applying the boundary
* values provided.
*/
template <int dim, int spacedim>
void initialize(const MGDoFHandler<dim,spacedim>& dof,
const typename FunctionMap<dim>::type& function_map,
const std::vector<bool>& component_mask = std::vector<bool>());
/**
* Reset the data structures.
*/
void clear();
/**
* Determine whether a dof index
* is subject to a boundary
* constraint.
*/
bool at_boundary(unsigned int level, unsigned int index) const;
/**
* Determine whether a dof index
* is at the refinement edge.
*/
bool at_refinement_edge(unsigned int level, unsigned int index) const;
private:
/**
* The indices of boundary dofs
* for each level.
*/
std::vector<std::set<unsigned int> > boundary_indices;
/**
* The degrees of freedom on the
* refinement edge between a
* level and coarser cells.
*/
std::vector<std::vector<bool> > refinement_edge_indices;
/**
* The degrees of freedom on the
* refinement edge between a
* level and coarser cells, which
* are also on the boundary.
*
* This is a subset of
* #refinement_edge_indices.
*/
std::vector<std::vector<bool> > refinement_edge_boundary_indices;
};
bool
MGConstraints::at_boundary(unsigned int level, unsigned int index) const
{
AssertIndexRange(level, boundary_indices.size());
AssertIndexRange(level, boundary_indices.size());
Assert (false, ExcNotImplemented());
return false;
}
DEAL_II_NAMESPACE_CLOSE
#endif
|