Nektar::SpatialDomains::RefRegionParallelogram Class Reference

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

#include <RefRegionParallelogram.h>

Inheritance diagram for Nektar::SpatialDomains::RefRegionParallelogram:

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< NekDouble >	m_coord1
	Coordinate 1. More...
std::vector< NekDouble >	m_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.

    : RefRegion(coordim, radius, coord1, coord2, numModes, numPoints)
{
}

~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

coords

coordinates 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.

{
    // This is simplification for a two-dimenion domain of the algorithm in the
    // RefRegionCylinder::v_Contains method.
 
    const size_t dim = coords.size();
    Array<OneD, NekDouble> e(dim, 0.0);     // direction: rb - ra
    Array<OneD, NekDouble> m(dim - 1, 0.0); // momentum: ra x rb
    NekDouble d = 0.0;                      // distance
 
    // Direction
    e[0] = m_coord2[0] - m_coord1[0];
    e[1] = m_coord2[1] - m_coord1[1];
 
    // Cross product of vectors 'coord1'(ra) and 'coord2' (rb)
    m[0] = m_coord1[0] * m_coord2[1] - m_coord1[1] * m_coord2[0];
 
    // Distance of P (coords) to line AB (coord1coord2) is equal or less than R
    // d = || e x (rp - ra) || / || e ||
 
    // rA - rP
    Array<OneD, NekDouble> rpa(dim, 0.0);
    rpa[0] = coords[0] - m_coord1[0];
    rpa[1] = coords[1] - m_coord1[1];
 
    // || e ||
    NekDouble e_mod = sqrt(e[0] * e[0] + e[1] * e[1]);
 
    // || e x (rp - ra) ||
    Array<OneD, NekDouble> exrpa(dim - 1, 0.0);
    exrpa[0] = e[0] * rpa[1] - e[1] * rpa[0];
 
    NekDouble exrpa_mod = sqrt(exrpa[0] * exrpa[0]);
 
    d = exrpa_mod / e_mod;
    if (d >= m_radius)
    {
        return false;
    }
 
    // Is P between the regions below?
    // xa, xb plus the radius
    // ya, yb plus the radius
    Array<OneD, bool> insideFlag(dim, false);
    for (int i = 0; i < dim; ++i)
    {
        if (m_coord1[i] < m_coord2[i])
        {
            if ((m_coord1[i] - m_radius < coords[i]) &&
                (m_coord2[i] + m_radius > coords[i]))
            {
                insideFlag[i] = true;
            }
        }
        else
        {
            if ((m_coord1[i] + m_radius > coords[i]) &&
                (m_coord2[i] - m_radius < coords[i]))
            {
                insideFlag[i] = true;
            }
        }
    }
 
    if ((insideFlag[0] == true) && (insideFlag[1] == true))
    {
        return true;
    }
    else
    {
        return false;
    }
}

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