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 
31 namespace 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
PointCell()
Specify cell type and facet type.
Definition: PointCell.h:41
std::size_t orientation(const Cell &cell) const
Return orientation of the cell.
Definition: PointCell.cpp:72
double squared_distance(const Cell &cell, const Point &point) const
Compute squared distance to given point.
Definition: PointCell.cpp:106
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::size_t num_vertices() const
Return number of vertices for cell.
Definition: CellType.h:92
double circumradius(const MeshEntity &point) const
Compute circumradius of PointCell.
Definition: PointCell.cpp:98
Definition: adapt.h:29
Definition: Point.h:40
double volume(const MeshEntity &triangle) const
Compute (generalized) volume (area) of triangle.
Definition: PointCell.cpp:89
std::size_t dim() const
Return topological dimension of cell.
Definition: PointCell.cpp:37
bool collides(const Cell &cell, const Point &point) const
Check whether given point is contained in cell.
Definition: PointCell.cpp:156
std::size_t num_entities(std::size_t dim) const
Return number of entities of given topological dimension.
Definition: PointCell.cpp:42
Type
Enum for different cell types.
Definition: CellType.h:51
A Cell is a MeshEntity of topological codimension 0.
Definition: Cell.h:42
double normal(const Cell &cell, std::size_t facet, std::size_t i) const
Definition: PointCell.cpp:112
Definition: CellType.h:46
Point cell_normal(const Cell &cell) const
Compute normal to given cell (viewed as embedded in 1D)
Definition: PointCell.cpp:130
double facet_area(const Cell &cell, std::size_t facet) const
Definition: PointCell.cpp:139
bool is_simplex() const
Check if cell is a simplex.
Definition: PointCell.h:44
Definition: MeshEntity.h:42
This class implements functionality for point cell meshes.
Definition: PointCell.h:36
std::vector< std::int8_t > vtk_mapping() const
Mapping of DOLFIN/UFC vertex ordering to VTK/XDMF ordering.
Definition: PointCell.h:104
std::string description(bool plural) const
Return description of cell type.
Definition: PointCell.cpp:166
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