#include <NodalPrismSPI.h>

Inheritance diagram for Nektar::LibUtilities::NodalPrismSPI:

Public Member Functions
	~NodalPrismSPI () override

	NodalPrismSPI (const PointsKey &key)

Public Member Functions inherited from Nektar::LibUtilities::Points< NekDouble >
virtual	~Points ()

void	Initialize (void)

size_t	GetPointsDim () const

size_t	GetNumPoints () const

size_t	GetTotNumPoints () const

PointsType	GetPointsType () const

const Array< OneD, const DataType > &	GetZ () const

const Array< OneD, const DataType > &	GetW () const

void	GetZW (Array< OneD, const DataType > &z, Array< OneD, const DataType > &w) const

const Array< OneD, const NekDouble > &	GetBaryWeights () const

void	GetPoints (Array< OneD, const DataType > &x) const

void	GetPoints (Array< OneD, const DataType > &x, Array< OneD, const DataType > &y) const

void	GetPoints (Array< OneD, const DataType > &x, Array< OneD, const DataType > &y, Array< OneD, const DataType > &z) const

const MatrixSharedPtrType &	GetD (Direction dir=xDir) const

const MatrixSharedPtrType	GetI (const PointsKey &key)

const MatrixSharedPtrType	GetI (const Array< OneD, const DataType > &x)

const MatrixSharedPtrType	GetI (size_t uint, const Array< OneD, const DataType > &x)

const MatrixSharedPtrType	GetI (const Array< OneD, const DataType > &x, const Array< OneD, const DataType > &y)

const MatrixSharedPtrType	GetI (const Array< OneD, const DataType > &x, const Array< OneD, const DataType > &y, const Array< OneD, const DataType > &z)

const MatrixSharedPtrType	GetGalerkinProjection (const PointsKey &pkey)

Static Public Member Functions
static std::shared_ptr< PointsBaseType >	Create (const PointsKey &key)

Private Member Functions
	NodalPrismSPI ()=delete

	NodalPrismSPI (const NodalPrismSPI &points)=delete

void	v_CalculatePoints () final

void	v_CalculateWeights () final

void	v_CalculateDerivMatrix () final

Private Attributes
Array< OneD, NekDouble >	m_t0

Array< OneD, NekDouble >	m_t1

Array< OneD, NekDouble >	m_tw

Array< OneD, NekDouble >	m_e0

Array< OneD, NekDouble >	m_ew

size_t	m_numtri

Static Private Attributes
static bool	initPointsManager []

Additional Inherited Members
Public Types inherited from Nektar::LibUtilities::Points< NekDouble >
typedef NekDouble	DataType

typedef std::shared_ptr< NekMatrix< DataType > >	MatrixSharedPtrType

Protected Member Functions inherited from Nektar::LibUtilities::Points< NekDouble >
virtual void	v_Initialize (void)

virtual void	v_CalculatePoints ()

virtual void	v_CalculateWeights ()

Protected Attributes inherited from Nektar::LibUtilities::Points< NekDouble >
PointsKey	m_pointsKey
	Points type for this points distributions. More...

Array< OneD, DataType >	m_points [3]
	Storage for the point locations, allowing for up to a 3D points storage. More...

Array< OneD, DataType >	m_weights
	Quadrature weights for the weights. More...

Array< OneD, DataType >	m_bcweights
	Barycentric weights. More...

MatrixSharedPtrType	m_derivmatrix [3]
	Derivative matrices. More...

NekManager< PointsKey, NekMatrix< DataType >, PointsKey::opLess >	m_InterpManager

NekManager< PointsKey, NekMatrix< DataType >, PointsKey::opLess >	m_GalerkinProjectionManager

Detailed Description

Definition at line 45 of file NodalPrismSPI.h.

Constructor & Destructor Documentation

◆ ~NodalPrismSPI()

Nektar::LibUtilities::NodalPrismSPI::~NodalPrismSPI ( )

inlineoverride

Definition at line 48 of file NodalPrismSPI.h.

49 {

50 }

◆ NodalPrismSPI() [1/3]

Nektar::LibUtilities::NodalPrismSPI::NodalPrismSPI ( const PointsKey & key )

inline

Definition at line 55 of file NodalPrismSPI.h.

                                        : PointsBaseType(key)
    {
    }

◆ NodalPrismSPI() [2/3]

Nektar::LibUtilities::NodalPrismSPI::NodalPrismSPI ( )

privatedelete

◆ NodalPrismSPI() [3/3]

Nektar::LibUtilities::NodalPrismSPI::NodalPrismSPI ( const NodalPrismSPI & points )

privatedelete

Member Function Documentation

◆ Create()

std::shared_ptr< PointsBaseType > Nektar::LibUtilities::NodalPrismSPI::Create ( const PointsKey & key )

static

Definition at line 92 of file NodalPrismSPI.cpp.

{
    std::shared_ptr<PointsBaseType> returnval(
        MemoryManager<NodalPrismSPI>::AllocateSharedPtr(key));
    returnval->Initialize();
    return returnval;
}

◆ v_CalculateDerivMatrix()

void Nektar::LibUtilities::NodalPrismSPI::v_CalculateDerivMatrix ( )

finalprivate

Definition at line 88 of file NodalPrismSPI.cpp.

89{

90}

◆ v_CalculatePoints()

void Nektar::LibUtilities::NodalPrismSPI::v_CalculatePoints ( )

finalprivatevirtual

Reimplemented from Nektar::LibUtilities::Points< NekDouble >.

Definition at line 42 of file NodalPrismSPI.cpp.

{
    // Allocate the storage for points
    size_t numPoints = GetNumPoints();
 
    PointsKey e(numPoints, eGaussLobattoLegendre);
    PointsManager()[e]->GetPoints(m_e0);
    m_ew = PointsManager()[e]->GetW();
 
    PointsKey t(numPoints, eNodalTriSPI);
    PointsManager()[t]->GetPoints(m_t0, m_t1);
    m_tw     = PointsManager()[t]->GetW();
    m_numtri = m_tw.size();
 
    for (size_t i = 0; i < 3; i++)
    {
        m_points[i] = Array<OneD, DataType>(m_numtri * numPoints);
    }
 
    for (size_t j = 0, ct = 0; j < numPoints; j++)
    {
        for (size_t i = 0; i < m_numtri; i++, ct++)
        {
            // need to flip y and z because of quad orientation
            m_points[0][ct] = m_t0[i];
            m_points[1][ct] = m_e0[j];
            m_points[2][ct] = m_t1[i];
        }
    }
}

References Nektar::LibUtilities::eGaussLobattoLegendre, Nektar::LibUtilities::eNodalTriSPI, Nektar::LibUtilities::Points< NekDouble >::GetNumPoints(), m_e0, m_ew, m_numtri, Nektar::LibUtilities::Points< NekDouble >::m_points, m_t0, m_t1, m_tw, and Nektar::LibUtilities::PointsManager().

◆ v_CalculateWeights()

void Nektar::LibUtilities::NodalPrismSPI::v_CalculateWeights ( )

finalprivatevirtual

Reimplemented from Nektar::LibUtilities::Points< NekDouble >.

Definition at line 73 of file NodalPrismSPI.cpp.

{
    size_t numPoints = GetNumPoints();
 
    m_weights = Array<OneD, DataType>(m_numtri * numPoints);
 
    for (size_t j = 0, ct = 0; j < numPoints; j++)
    {
        for (size_t i = 0; i < m_numtri; i++, ct++)
        {
            m_weights[ct] = m_tw[i] * m_ew[j];
        }
    }
}