DOLFIN
DOLFIN C++ interface
Loading...
Searching...
No Matches
TriangleCell.h
1// Copyright (C) 2006-2017 Anders Logg
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 Kristoffer Selim, 2008.
19// Modified by Jan Blechta 2013
20//
21// First added: 2006-06-05
22// Last changed: 2017-09-26
23
24#ifndef __TRIANGLE_CELL_H
25#define __TRIANGLE_CELL_H
26
27#include <vector>
28#include <boost/multi_array.hpp>
29#include "CellType.h"
30
31namespace dolfin
32{
33
35
36 class TriangleCell : public CellType
37 {
38 public:
39
41 TriangleCell() : CellType(Type::triangle, Type::interval) {}
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
54 std::size_t num_vertices(std::size_t dim) const;
55
57 std::size_t orientation(const Cell& cell) const;
58
60 void create_entities(boost::multi_array<unsigned int, 2>& e,
61 std::size_t dim,
62 const unsigned int* v) const;
63
65 double volume(const MeshEntity& triangle) const;
66
68 double circumradius(const MeshEntity& triangle) const;
69
71 double squared_distance(const Cell& cell, const Point& point) const;
72
77 static double squared_distance(const Point& point,
78 const Point& a,
79 const Point& b,
80 const Point& c);
81
83 double normal(const Cell& cell, std::size_t facet, std::size_t i) const;
84
86 Point normal(const Cell& cell, std::size_t facet) const;
87
89 Point cell_normal(const Cell& cell) const;
90
92 double facet_area(const Cell& cell, std::size_t facet) const;
93
95 void order(Cell& cell,
96 const std::vector<std::int64_t>& local_to_global_vertex_indices) const;
97
99 bool collides(const Cell& cell, const Point& point) const;
100
102 bool collides(const Cell& cell, const MeshEntity& entity) const;
103
105 std::string description(bool plural) const;
106
108 std::vector<std::int8_t> vtk_mapping() const
109 { return {0, 1, 2}; }
110
111 private:
112
113 // Find local index of edge i according to ordering convention
114 std::size_t find_edge(std::size_t i, const Cell& cell) const;
115 };
116
117}
118
119#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
Definition Point.h:41
This class implements functionality for triangular meshes.
Definition TriangleCell.h:37
Point cell_normal(const Cell &cell) const
Compute normal to given cell (viewed as embedded in 3D)
Definition TriangleCell.cpp:338
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 TriangleCell.cpp:91
double squared_distance(const Cell &cell, const Point &point) const
Compute squared distance to given point (3D enabled)
Definition TriangleCell.cpp:200
double circumradius(const MeshEntity &triangle) const
Compute diameter of triangle.
Definition TriangleCell.cpp:162
void order(Cell &cell, const std::vector< std::int64_t > &local_to_global_vertex_indices) const
Order entities locally.
Definition TriangleCell.cpp:388
std::vector< std::int8_t > vtk_mapping() const
Mapping of DOLFIN/UFC vertex ordering to VTK/XDMF ordering.
Definition TriangleCell.h:108
std::string description(bool plural) const
Return description of cell type.
Definition TriangleCell.cpp:461
bool is_simplex() const
Check if cell is a simplex.
Definition TriangleCell.h:44
double facet_area(const Cell &cell, std::size_t facet) const
Compute the area/length of given facet with respect to the cell.
Definition TriangleCell.cpp:369
std::size_t dim() const
Return topological dimension of cell.
Definition TriangleCell.cpp:42
TriangleCell()
Specify cell type and facet type.
Definition TriangleCell.h:41
std::size_t orientation(const Cell &cell) const
Return orientation of the cell.
Definition TriangleCell.cpp:85
bool collides(const Cell &cell, const Point &point) const
Check whether given point collides with cell.
Definition TriangleCell.cpp:451
double volume(const MeshEntity &triangle) const
Compute (generalized) volume (area) of triangle.
Definition TriangleCell.cpp:111
std::size_t num_entities(std::size_t dim) const
Return number of entities of given topological dimension.
Definition TriangleCell.cpp:47
double normal(const Cell &cell, std::size_t facet, std::size_t i) const
Compute component i of normal of given facet with respect to the cell.
Definition TriangleCell.cpp:288
Definition adapt.h:30