Nektar++
Public Member Functions | Protected Member Functions | List of all members
Nektar::SpatialDomains::RefRegionParallelogram Class Reference

Derived class for the refinement surface region. More...

#include <RefRegionParallelogram.h>

Inheritance diagram for Nektar::SpatialDomains::RefRegionParallelogram:
[legend]

Public Member Functions

 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. More...
 
 ~RefRegionParallelogram () override=default
 Destructor. More...
 
- Public Member Functions inherited from Nektar::SpatialDomains::RefRegion
 RefRegion (const unsigned int coordim, NekDouble m_radius, std::vector< NekDouble > coord1, std::vector< NekDouble > coord2, std::vector< unsigned int > numModes, std::vector< unsigned int > numPoints)
 Constructor. More...
 
virtual ~RefRegion ()=default
 Destructor. More...
 
virtual bool v_Contains (const Array< OneD, NekDouble > &coords)=0
 Pure virtual fuction. More...
 
std::vector< unsigned int > GetNumModes ()
 Get the number of modes to update expansion. More...
 
std::vector< unsigned int > GetNumPoints ()
 Get the number of quadrature points to update expansion. More...
 

Protected Member Functions

bool v_Contains (const Array< OneD, NekDouble > &coords) override
 Check if vertex is inside the surface region. More...
 

Additional Inherited Members

- Protected Attributes inherited from Nektar::SpatialDomains::RefRegion
std::vector< NekDoublem_coord1
 Coordinate 1. More...
 
std::vector< NekDoublem_coord2
 Coordinate 2. More...
 
std::vector< unsigned int > m_numModes
 Number of modes. More...
 
std::vector< unsigned int > m_numPoints
 Number of quadrature points. More...
 
NekDouble m_radius
 Radius of the surface region. More...
 
unsigned int m_coordim
 Dimension of the coordinate (space dimension) More...
 

Detailed Description

Derived class for the refinement surface region.

Definition at line 48 of file RefRegionParallelogram.h.

Constructor & Destructor Documentation

◆ RefRegionParallelogram()

Nektar::SpatialDomains::RefRegionParallelogram::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.

Definition at line 40 of file RefRegionParallelogram.cpp.

44 : RefRegion(coordim, radius, coord1, coord2, numModes, numPoints)
45{
46}
RefRegion(const unsigned int coordim, NekDouble m_radius, std::vector< NekDouble > coord1, std::vector< NekDouble > coord2, std::vector< unsigned int > numModes, std::vector< unsigned int > numPoints)
Constructor.
Definition: RefRegion.cpp:45

◆ ~RefRegionParallelogram()

Nektar::SpatialDomains::RefRegionParallelogram::~RefRegionParallelogram ( )
overridedefault

Destructor.

Member Function Documentation

◆ v_Contains()

bool Nektar::SpatialDomains::RefRegionParallelogram::v_Contains ( const Array< OneD, NekDouble > &  coords)
overrideprotectedvirtual

Check if vertex is inside the surface region.

Check if vertex is inside the Parallelogram.

Parameters
coordscoordinates of the vertex
Returns
true or false depending on if the vertex is inside or not of the surface defined by the user.

Implements Nektar::SpatialDomains::RefRegion.

Definition at line 55 of file RefRegionParallelogram.cpp.

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}
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
std::vector< double > d(NPUPPER *NPUPPER)
double NekDouble
scalarT< T > sqrt(scalarT< T > in)
Definition: scalar.hpp:285

References Nektar::UnitTests::d(), Nektar::SpatialDomains::RefRegion::m_coord1, Nektar::SpatialDomains::RefRegion::m_coord2, Nektar::SpatialDomains::RefRegion::m_radius, and tinysimd::sqrt().