DOLFIN
DOLFIN C++ interface
PointCell.h
1// Copyright (C) 2007-2007 Kristian B. Oelgaard
2//
3// This file is part of DOLFIN.
4//
5// DOLFIN is free software: you can redistribute it and/or modify
6// it under the terms of the GNU Lesser General Public License as published by
7// the Free Software Foundation, either version 3 of the License, or
8// (at your option) any later version.
9//
10// DOLFIN is distributed in the hope that it will be useful,
11// but WITHOUT ANY WARRANTY; without even the implied warranty of
12// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13// GNU Lesser General Public License for more details.
14//
15// You should have received a copy of the GNU Lesser General Public License
16// along with DOLFIN. If not, see <http://www.gnu.org/licenses/>.
17//
18// Modified by Anders Logg 2008-2013
19// Modified by Kristoffer Selim 2008
20//
21// First added: 2007-12-12
22// Last changed: 2016-05-05
23
24#ifndef __POINT_CELL_H
25#define __POINT_CELL_H
26
27#include <vector>
28#include <boost/multi_array.hpp>
29#include "CellType.h"
30
31namespace dolfin
32{
33
35
36 class PointCell : public CellType
37 {
38 public:
39
41 PointCell() : CellType(Type::point, Type::point) {}
42
44 bool is_simplex() const
45 { return true; }
46
48 std::size_t dim() const;
49
51 std::size_t num_entities(std::size_t dim) const;
52
55 std::size_t num_vertices(std::size_t dim) const;
56
58 std::size_t orientation(const Cell& cell) const;
59
61 void create_entities(boost::multi_array<unsigned int, 2>& e,
62 std::size_t dim, const unsigned int* v) const;
63
65 void order(Cell& cell, const std::vector<std::int64_t>&
66 local_to_global_vertex_indices) const;
67
69 double volume(const MeshEntity& triangle) const;
70
72 double circumradius(const MeshEntity& point) const;
73
75 double squared_distance(const Cell& cell, const Point& point) const;
76
79 double normal(const Cell& cell, std::size_t facet, std::size_t i) const;
80
82 Point normal(const Cell& cell, std::size_t facet) const;
83
85 Point cell_normal(const Cell& cell) const;
86
89 double facet_area(const Cell& cell, std::size_t facet) const;
90
92 void order(Cell& cell) const;
93
95 bool collides(const Cell& cell, const Point& point) const;
96
98 bool collides(const Cell& cell, const MeshEntity& entity) const;
99
101 std::string description(bool plural) const;
102
104 std::vector<std::int8_t> vtk_mapping() const
105 { return {0}; }
106
107 private:
108
109 // Find local index of edge i according to ordering convention
110 std::size_t find_edge(std::size_t i, const Cell& cell) const;
111
112 };
113
114}
115
116#endif
Definition: CellType.h:47
std::size_t num_vertices() const
Return number of vertices for cell.
Definition: CellType.h:92
Type
Enum for different cell types.
Definition: CellType.h:51
A Cell is a MeshEntity of topological codimension 0.
Definition: Cell.h:43
Definition: MeshEntity.h:43
This class implements functionality for point cell meshes.
Definition: PointCell.h:37
Point cell_normal(const Cell &cell) const
Compute normal to given cell (viewed as embedded in 1D)
Definition: PointCell.cpp:130
bool is_simplex() const
Check if cell is a simplex.
Definition: PointCell.h:44
std::size_t num_entities(std::size_t dim) const
Return number of entities of given topological dimension.
Definition: PointCell.cpp:42
std::size_t dim() const
Return topological dimension of cell.
Definition: PointCell.cpp:37
double volume(const MeshEntity &triangle) const
Compute (generalized) volume (area) of triangle.
Definition: PointCell.cpp:89
double circumradius(const MeshEntity &point) const
Compute circumradius of PointCell.
Definition: PointCell.cpp:98
double facet_area(const Cell &cell, std::size_t facet) const
Definition: PointCell.cpp:139
void order(Cell &cell, const std::vector< std::int64_t > &local_to_global_vertex_indices) const
Order entities locally (connectivity 1-0, 2-0, 2-1)
Definition: PointCell.cpp:147
std::string description(bool plural) const
Return description of cell type.
Definition: PointCell.cpp:166
double squared_distance(const Cell &cell, const Point &point) const
Compute squared distance to given point.
Definition: PointCell.cpp:106
void create_entities(boost::multi_array< unsigned int, 2 > &e, std::size_t dim, const unsigned int *v) const
Create entities e of given topological dimension from vertices v.
Definition: PointCell.cpp:80
void order(Cell &cell) const
Order entities locally.
double normal(const Cell &cell, std::size_t facet, std::size_t i) const
Definition: PointCell.cpp:112
std::vector< std::int8_t > vtk_mapping() const
Mapping of DOLFIN/UFC vertex ordering to VTK/XDMF ordering.
Definition: PointCell.h:104
std::size_t orientation(const Cell &cell) const
Return orientation of the cell.
Definition: PointCell.cpp:72
bool collides(const Cell &cell, const Point &point) const
Check whether given point is contained in cell.
Definition: PointCell.cpp:156
PointCell()
Specify cell type and facet type.
Definition: PointCell.h:41
Definition: Point.h:41
Definition: adapt.h:30