tetCellI.H
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2011-2018 OpenFOAM Foundation
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8 License
9  This file is part of OpenFOAM.
10 
11  OpenFOAM is free software: you can redistribute it and/or modify it
12  under the terms of the GNU General Public License as published by
13  the Free Software Foundation, either version 3 of the License, or
14  (at your option) any later version.
15 
16  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
19  for more details.
20 
21  You should have received a copy of the GNU General Public License
22  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
23 
24 Description
25 
26 \*---------------------------------------------------------------------------*/
27 
28 #include "IOstreams.H"
29 
30 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
31 
33 {}
34 
35 
37 (
38  const label a,
39  const label b,
40  const label c,
41  const label d
42 )
43 {
44  operator[](0) = a;
45  operator[](1) = b;
46  operator[](2) = c;
47  operator[](3) = d;
48 }
49 
50 
52 :
53  FixedList<label, 4>(lst)
54 {}
55 
56 
58 :
59  FixedList<label, 4>(is)
60 {}
61 
62 
63 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
64 
65 inline Foam::triFace Foam::tetCell::face(const label facei) const
66 {
67  // Warning. Ordering of faces needs to be the same for a tetrahedron
68  // class, a tetrahedron cell shape model and a tetCell
69  static const label a[] = {1, 0, 0, 0};
70  static const label b[] = {2, 3, 1, 2};
71  static const label c[] = {3, 2, 3, 1};
72 
73  #ifdef FULLDEBUG
74  if (facei >= 4)
75  {
77  << "index out of range 0 -> 3. facei = " << facei
78  << abort(FatalError);
79  }
80  #endif
81 
82  return triFace
83  (
84  operator[](a[facei]),
85  operator[](b[facei]),
86  operator[](c[facei])
87  );
88 }
89 
90 
91 inline Foam::label Foam::tetCell::edgeFace(const label edgeI) const
92 {
93  // Warning. Ordering of faces needs to be the same for a tetrahedron
94  // class, a tetrahedron cell shape model and a tetCell
95  // static const label edgeFaces[6] = {2, 1, 1, 0, 0, 0};
96  static const label edgeFaces[6] = {2, 3, 1, 0, 0, 1};
97 
98  #ifdef FULLDEBUG
99  if (edgeI >= 6)
100  {
102  << "edge index out of range 0 -> 5. edgeI = " << edgeI
103  << abort(FatalError);
104  }
105  #endif
106 
107  return edgeFaces[edgeI];
108 }
109 
110 
112 (
113  const label edgeI,
114  const label facei
115 ) const
116 {
117  // Warning. Ordering of faces needs to be the same for a tetrahedron
118  // class, a tetrahedron cell shape model and a tetCell
119  static const label adjacentFace[6][4] =
120  {
121  {-1, -1, 3, 2},
122  {-1, 3, -1, 1},
123  {-1, 2, 1, -1},
124  {2, -1, 0, -1},
125  {3, -1, -1, 0},
126  {1, 0, -1, -1}
127  };
128 
129  #ifdef FULLDEBUG
130  if (facei >= 4)
131  {
133  << "face index out of range 0 -> 3. facei = " << facei
134  << abort(FatalError);
135  }
136 
137  if (edgeI >= 6)
138  {
140  << "edge index out of range 0 -> 5. edgeI = " << edgeI
141  << abort(FatalError);
142  }
143  #endif
144 
145  return adjacentFace[edgeI][facei];
146 }
147 
148 
149 inline Foam::edge Foam::tetCell::tetEdge(const label edgeI) const
150 {
151  // Warning. Ordering of edges needs to be the same for a tetrahedron
152  // class, a tetrahedron cell shape model and a tetCell
153  //
154  static const label start[] = {0, 0, 0, 3, 1, 3};
155  static const label end[] = {1, 2, 3, 1, 2, 2};
156 
157  #ifdef FULLDEBUG
158  if (edgeI >= 6)
159  {
161  << "index out of range 0 -> 5. edgeI = " << edgeI
162  << abort(FatalError);
163  }
164  #endif
165 
166  return edge(operator[](start[edgeI]), operator[](end[edgeI]));
167 }
168 
169 
171 {
172  return tetPointRef
173  (
174  points[operator[](0)],
175  points[operator[](1)],
176  points[operator[](2)],
177  points[operator[](3)]
178  );
179 }
180 
181 
182 // ************************************************************************* //
A tetrahedron primitive.
Definition: tetrahedron.H:61
tetCell()
Construct null.
Definition: tetCellI.H:32
tetrahedron< point, const point & > tetPointRef
Definition: tetPointRef.H:44
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Definition: label.H:59
label edgeAdjacentFace(const label edgeI, const label facei) const
Return face adjacent to the given face sharing the same edge.
Definition: tetCellI.H:112
tetPointRef tet(const pointField &) const
Return the tetrahedron.
Definition: tetCellI.H:170
error FatalError
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
Definition: error.H:323
An Istream is an abstract base class for all input systems (streams, files, token lists etc)...
Definition: Istream.H:57
const dimensionedScalar b
Wien displacement law constant: default SI units: [m K].
Definition: createFields.H:27
const dimensionedScalar c
Speed of light in a vacuum.
Useful combination of include files which define Sin, Sout and Serr and the use of IO streams general...
const pointField & points
An edge is a list of two point labels. The functionality it provides supports the discretisation on a...
Definition: edge.H:58
label & operator[](const label)
Return element of FixedList.
Definition: FixedListI.H:247
A triangular face using a FixedList of labels corresponding to mesh vertices.
Definition: triFace.H:68
face triFace(3)
edge tetEdge(const label edgeI) const
Return i-th edge.
Definition: tetCellI.H:149
errorManip< error > abort(error &err)
Definition: errorManip.H:131
iterator end()
Return an iterator to end traversing the FixedList.
Definition: FixedListI.H:331
label edgeFace(const label edgeI) const
Return first face adjacent to the given edge.
Definition: tetCellI.H:91
triFace face(const label facei) const
Return i-th face.
Definition: tetCellI.H:65