Nektar++
RefRegionParallelogram.cpp
Go to the documentation of this file.
1////////////////////////////////////////////////////////////////////////////////
2//
3// File: RefRegionParallelogram.cpp
4//
5// For more information, please see: http://www.nektar.info/
6//
7// The MIT License
8//
9// Copyright (c) 2006 Division of Applied Mathematics, Brown University (USA),
10// Department of Aeronautics, Imperial College London (UK), and Scientific
11// Computing and Imaging Institute, University of Utah (USA).
12//
13// Permission is hereby granted, free of charge, to any person obtaining a
14// copy of this software and associated documentation files (the "Software"),
15// to deal in the Software without restriction, including without limitation
16// the rights to use, copy, modify, merge, publish, distribute, sublicense,
17// and/or sell copies of the Software, and to permit persons to whom the
18// Software is furnished to do so, subject to the following conditions:
19//
20// The above copyright notice and this permission notice shall be included
21// in all copies or substantial portions of the Software.
22//
23// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
24// OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
25// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
26// THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
27// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
28// FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
29// DEALINGS IN THE SOFTWARE.
30//
31// Description: Parallelogram surface for the refinement region.
32//
33////////////////////////////////////////////////////////////////////////////////
34
36
38{
39
41 const unsigned int coordim, NekDouble radius, std::vector<NekDouble> coord1,
42 std::vector<NekDouble> coord2, std::vector<unsigned int> numModes,
43 std::vector<unsigned int> numPoints)
44 : RefRegion(coordim, radius, coord1, coord2, numModes, numPoints)
45{
46}
47
48/**
49 * @brief Check if vertex is inside the Parallelogram.
50 *
51 * @param coords coordinates of the vertex
52 * @return true or false depending on if the vertex is inside
53 * or not of the surface defined by the user.
54 */
56{
57 // This is simplification for a two-dimenion domain of the algorithm in the
58 // RefRegionCylinder::v_Contains method.
59
60 const size_t dim = coords.size();
61 Array<OneD, NekDouble> e(dim, 0.0); // direction: rb - ra
62 Array<OneD, NekDouble> m(dim - 1, 0.0); // momentum: ra x rb
63 NekDouble d = 0.0; // distance
64
65 // Direction
66 e[0] = m_coord2[0] - m_coord1[0];
67 e[1] = m_coord2[1] - m_coord1[1];
68
69 // Cross product of vectors 'coord1'(ra) and 'coord2' (rb)
70 m[0] = m_coord1[0] * m_coord2[1] - m_coord1[1] * m_coord2[0];
71
72 // Distance of P (coords) to line AB (coord1coord2) is equal or less than R
73 // d = || e x (rp - ra) || / || e ||
74
75 // rA - rP
76 Array<OneD, NekDouble> rpa(dim, 0.0);
77 rpa[0] = coords[0] - m_coord1[0];
78 rpa[1] = coords[1] - m_coord1[1];
79
80 // || e ||
81 NekDouble e_mod = sqrt(e[0] * e[0] + e[1] * e[1]);
82
83 // || e x (rp - ra) ||
84 Array<OneD, NekDouble> exrpa(dim - 1, 0.0);
85 exrpa[0] = e[0] * rpa[1] - e[1] * rpa[0];
86
87 NekDouble exrpa_mod = sqrt(exrpa[0] * exrpa[0]);
88
89 d = exrpa_mod / e_mod;
90 if (d >= m_radius)
91 {
92 return false;
93 }
94
95 // Is P between the regions below?
96 // xa, xb plus the radius
97 // ya, yb plus the radius
98 Array<OneD, bool> insideFlag(dim, false);
99 for (int i = 0; i < dim; ++i)
100 {
101 if (m_coord1[i] < m_coord2[i])
102 {
103 if ((m_coord1[i] - m_radius < coords[i]) &&
104 (m_coord2[i] + m_radius > coords[i]))
105 {
106 insideFlag[i] = true;
107 }
108 }
109 else
110 {
111 if ((m_coord1[i] + m_radius > coords[i]) &&
112 (m_coord2[i] - m_radius < coords[i]))
113 {
114 insideFlag[i] = true;
115 }
116 }
117 }
118
119 if ((insideFlag[0] == true) && (insideFlag[1] == true))
120 {
121 return true;
122 }
123 else
124 {
125 return false;
126 }
127}
128
129} // namespace Nektar::SpatialDomains
Abstract base class for the refinement surface region.
Definition: RefRegion.h:51
std::vector< NekDouble > m_coord2
Coordinate 2.
Definition: RefRegion.h:83
std::vector< NekDouble > m_coord1
Coordinate 1.
Definition: RefRegion.h:81
NekDouble m_radius
Radius of the surface region.
Definition: RefRegion.h:89
bool v_Contains(const Array< OneD, NekDouble > &coords) override
Check if vertex is inside the surface region.
RefRegionParallelogram(const unsigned int coordim, NekDouble radius, std::vector< NekDouble > coord1, std::vector< NekDouble > coord2, std::vector< unsigned int > numModes, std::vector< unsigned int > numPoints)
Constructor.
std::vector< double > d(NPUPPER *NPUPPER)
double NekDouble
scalarT< T > sqrt(scalarT< T > in)
Definition: scalar.hpp:285