/usr/include/dolfin/mesh/Facet.h is in libdolfin-dev 1.4.0+dfsg-4.
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 | // Copyright (C) 2006-2010 Anders Logg
//
// This file is part of DOLFIN.
//
// DOLFIN is free software: you can 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 3 of the License, or
// (at your option) any later version.
//
// DOLFIN is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with DOLFIN. If not, see <http://www.gnu.org/licenses/>.
//
// Modified by Garth N. Wells, 2009-2011.
//
// First added: 2006-06-02
// Last changed: 2011-04-07
#ifndef __FACET_H
#define __FACET_H
#include <memory>
#include <utility>
#include <vector>
#include "Cell.h"
#include "Mesh.h"
#include "MeshEntity.h"
#include "MeshEntityIteratorBase.h"
#include "MeshFunction.h"
namespace dolfin
{
/// A Facet is a MeshEntity of topological codimension 1.
class Facet : public MeshEntity
{
public:
/// Constructor
Facet(const Mesh& mesh, std::size_t index)
: MeshEntity(mesh, mesh.topology().dim() - 1, index) {}
/// Destructor
~Facet() {}
/// Compute component i of the normal to the facet
double normal(std::size_t i) const;
/// Compute normal to the facet
Point normal() const;
/// Return true if facet is an exterior facet (relative to global mesh,
/// so this function will return false for facets on partition
/// boundaries). Facet connectivity must be initialized before
/// calling this function.
bool exterior() const;
// FIXME: This function should take care of the case where adjacent cells
// FIXME: live on different processes
/// Return adjacent cells. An optional argument that lists for
/// each facet the index of the first cell may be given to specify
/// the ordering of the two cells. If not specified, the ordering
/// will depend on the (arbitrary) ordering of the mesh
/// connectivity.
std::pair<const Cell, const Cell>
adjacent_cells(const std::vector<std::size_t>* facet_orientation) const;
};
/// A FacetIterator is a MeshEntityIterator of topological
/// codimension 1.
typedef MeshEntityIteratorBase<Facet> FacetIterator;
/// A FacetFunction is a MeshFunction of topological codimension 1.
template <typename T> class FacetFunction : public MeshFunction<T>
{
public:
FacetFunction(const Mesh& mesh)
: MeshFunction<T>(mesh, mesh.topology().dim() - 1) {}
FacetFunction(std::shared_ptr<const Mesh> mesh)
: MeshFunction<T>(mesh, mesh->topology().dim() - 1) {}
FacetFunction(const Mesh& mesh, const T& value)
: MeshFunction<T>(mesh, mesh.topology().dim() - 1, value) {}
//FacetFunction(boost:shared_ptr<const Mesh> mesh, const T& value)
// : MeshFunction<T>(mesh, mesh->topology().dim() - 1, value) {}
};
}
#endif
|