Nektar::LibUtilities::BLPoints Class Reference

#include <BLPoints.h>

Inheritance diagram for Nektar::LibUtilities::BLPoints:

Public Member Functions
virtual	~BLPoints ()
std::shared_ptr< NekMatrix< NekDouble > >	CreateMatrix (const PointsKey &pkey)
	BLPoints (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)

Protected Member Functions
virtual const MatrixSharedPtrType	v_GetI (const PointsKey &pkey) override final
virtual const MatrixSharedPtrType	v_GetI (const Array< OneD, const NekDouble > &x) override final
virtual const MatrixSharedPtrType	v_GetI (size_t numpoints, const Array< OneD, const NekDouble > &x) override final
Protected Member Functions inherited from Nektar::LibUtilities::Points< NekDouble >
virtual void	v_Initialize (void)
virtual void	v_CalculatePoints ()
virtual void	v_CalculateWeights ()

Private Member Functions
	BLPoints ()=delete
	BLPoints (const BLPoints &points)=delete
virtual void	v_CalculatePoints () override
virtual void	v_CalculateWeights () override
virtual void	v_CalculateDerivMatrix () override
void	CalculateInterpMatrix (size_t npts, const Array< OneD, const NekDouble > &xpoints, Array< OneD, NekDouble > &interp)

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 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 44 of file BLPoints.h.

Constructor & Destructor Documentation

~BLPoints()

virtual Nektar::LibUtilities::BLPoints::~BLPoints ( )

inlinevirtual

Definition at line 47 of file BLPoints.h.

    {
    }

BLPoints() [1/3]

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

inline

Definition at line 56 of file BLPoints.h.

                                   : PointsBaseType(key)
    {
        namespace pl = std::placeholders;
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussGaussLegendre),
            std::bind(&BLPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussRadauMLegendre),
            std::bind(&BLPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussRadauPLegendre),
            std::bind(&BLPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussLobattoLegendre),
            std::bind(&BLPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussGaussChebyshev),
            std::bind(&BLPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussRadauMChebyshev),
            std::bind(&BLPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussRadauPChebyshev),
            std::bind(&BLPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussLobattoChebyshev),
            std::bind(&BLPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussRadauMAlpha0Beta1),
            std::bind(&BLPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussRadauMAlpha0Beta2),
            std::bind(&BLPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussRadauMAlpha1Beta0),
            std::bind(&BLPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussRadauMAlpha2Beta0),
            std::bind(&BLPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, ePolyEvenlySpaced),
            std::bind(&BLPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eFourierEvenlySpaced),
            std::bind(&BLPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eBoundaryLayerPoints),
            std::bind(&BLPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eBoundaryLayerPointsRev),
            std::bind(&BLPoints::CreateMatrix, this, pl::_1));
    }

References CreateMatrix(), Nektar::LibUtilities::eBoundaryLayerPoints, Nektar::LibUtilities::eBoundaryLayerPointsRev, Nektar::LibUtilities::eFourierEvenlySpaced, Nektar::LibUtilities::eGaussGaussChebyshev, Nektar::LibUtilities::eGaussGaussLegendre, Nektar::LibUtilities::eGaussLobattoChebyshev, Nektar::LibUtilities::eGaussLobattoLegendre, Nektar::LibUtilities::eGaussRadauMAlpha0Beta1, Nektar::LibUtilities::eGaussRadauMChebyshev, Nektar::LibUtilities::eGaussRadauMLegendre, Nektar::LibUtilities::eGaussRadauPChebyshev, Nektar::LibUtilities::eGaussRadauPLegendre, Nektar::LibUtilities::ePolyEvenlySpaced, and Nektar::LibUtilities::Points< NekDouble >::m_InterpManager.

BLPoints() [2/3]

Nektar::LibUtilities::BLPoints::BLPoints ( )

privatedelete

BLPoints() [3/3]

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

privatedelete

Member Function Documentation

CalculateInterpMatrix()

void Nektar::LibUtilities::BLPoints::CalculateInterpMatrix ( size_t  npts, const Array< OneD, const NekDouble > &  xpoints, Array< OneD, NekDouble > &  interp  )

private

Definition at line 157 of file BLPoints.cpp.

{
    boost::ignore_unused(npts, xpoints, interp);
 
    ASSERTL0(false,
             "CalculateInterpMatrix not available for Boundary Layer Points");
}

References ASSERTL0.

Create()

std::shared_ptr< Points< NekDouble > > Nektar::LibUtilities::BLPoints::Create ( const PointsKey &  key )

static

Definition at line 107 of file BLPoints.cpp.

{
    std::shared_ptr<Points<NekDouble>> returnval(
        MemoryManager<BLPoints>::AllocateSharedPtr(key));
 
    returnval->Initialize();
 
    return returnval;
}

CreateMatrix()

std::shared_ptr< NekMatrix< NekDouble > > Nektar::LibUtilities::BLPoints::CreateMatrix

(

const PointsKey &

pkey

)

Definition at line 117 of file BLPoints.cpp.

{
    boost::ignore_unused(pkey);
 
    ASSERTL0(false, "CreateMatrix not available for Boundary Layer Points");
 
    return nullptr;
}

References ASSERTL0.

Referenced by BLPoints().

v_CalculateDerivMatrix()

void Nektar::LibUtilities::BLPoints::v_CalculateDerivMatrix ( )

overrideprivatevirtual

Definition at line 101 of file BLPoints.cpp.

{
    //// Allocate the derivative matrix.
    Points<NekDouble>::v_CalculateDerivMatrix();
}

v_CalculatePoints()

void Nektar::LibUtilities::BLPoints::v_CalculatePoints ( )

overrideprivatevirtual

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

Definition at line 50 of file BLPoints.cpp.

{
    // Allocate the storage for points.
    PointsBaseType::v_CalculatePoints();
    size_t npts = m_pointsKey.GetNumPoints();
 
    // Derived power coefficient.
    NekDouble r = m_pointsKey.GetFactor();
    ASSERTL0(r != NekConstants::kNekUnsetDouble,
             "Must set factor in BLPoints key");
 
    if (fabs(r - 1.0) < 1e-6)
    {
        NekDouble tmp = 2.0 / (npts - 1.0);
        for (size_t i = 0; i < npts; ++i)
        {
            m_points[0][i] = -1.0 + i * tmp;
        }
    }
    else
    {
        NekDouble a    = 2.0 * (1.0 - r) / (1.0 - pow(r, (double)(npts - 1)));
        m_points[0][0] = -1.0;
 
        for (size_t i = 1; i < npts; ++i)
        {
            m_points[0][i] = m_points[0][i - 1] + a * pow(r, (double)(i - 1));
        }
 
        m_points[0][npts - 1] = 1.0;
    }
 
    if (m_pointsKey.GetPointsType() == eBoundaryLayerPointsRev)
    {
        std::vector<NekDouble> tmp(npts);
        for (size_t i = 0; i < npts; ++i)
        {
            tmp[i] = -m_points[0][npts - 1 - i];
        }
 
        for (size_t i = 0; i < npts; ++i)
        {
            m_points[0][i] = tmp[i];
        }
    }
}

References ASSERTL0, Nektar::LibUtilities::eBoundaryLayerPointsRev, Nektar::LibUtilities::PointsKey::GetFactor(), Nektar::LibUtilities::PointsKey::GetNumPoints(), Nektar::LibUtilities::PointsKey::GetPointsType(), Nektar::NekConstants::kNekUnsetDouble, Nektar::LibUtilities::Points< NekDouble >::m_points, Nektar::LibUtilities::Points< NekDouble >::m_pointsKey, and Nektar::LibUtilities::Points< NekDouble >::v_CalculatePoints().

v_CalculateWeights()

void Nektar::LibUtilities::BLPoints::v_CalculateWeights ( )

overrideprivatevirtual

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

Definition at line 97 of file BLPoints.cpp.

{
}

v_GetI() [1/3]

const std::shared_ptr< NekMatrix< NekDouble > > Nektar::LibUtilities::BLPoints::v_GetI ( const Array< OneD, const NekDouble > &  x )

finaloverrideprotectedvirtual

Definition at line 137 of file BLPoints.cpp.

{
    boost::ignore_unused(x);
 
    ASSERTL0(false, "Interp not available for Boundary Layer Points");
 
    return nullptr;
}

References ASSERTL0.

v_GetI() [2/3]

const std::shared_ptr< NekMatrix< NekDouble > > Nektar::LibUtilities::BLPoints::v_GetI ( const PointsKey &  pkey )

finaloverrideprotectedvirtual

Definition at line 127 of file BLPoints.cpp.

{
    boost::ignore_unused(pkey);
 
    ASSERTL0(false, "Interp not available for Boundary Layer Points");
 
    return nullptr;
}

References ASSERTL0.

v_GetI() [3/3]

const std::shared_ptr< NekMatrix< NekDouble > > Nektar::LibUtilities::BLPoints::v_GetI ( size_t  numpoints, const Array< OneD, const NekDouble > &  x  )

finaloverrideprotectedvirtual

Definition at line 147 of file BLPoints.cpp.

{
    boost::ignore_unused(numpoints, x);
 
    ASSERTL0(false, "Interp not available for Boundary Layer Points");
 
    return nullptr;
}

References ASSERTL0.

Member Data Documentation

initPointsManager

bool Nektar::LibUtilities::BLPoints::initPointsManager

staticprivate

Initial value:

= {
    PointsManager().RegisterCreator(PointsKey(0, eBoundaryLayerPoints),
                                    BLPoints::Create),
    PointsManager().RegisterCreator(PointsKey(0, eBoundaryLayerPointsRev),
                                    BLPoints::Create)}

Definition at line 118 of file BLPoints.h.