Nektar::LibUtilities::GaussPoints Class Reference

#include <GaussPoints.h>

Inheritance diagram for Nektar::LibUtilities::GaussPoints:

Public Member Functions
	~GaussPoints () override
std::shared_ptr< NekMatrix< NekDouble > >	CreateMatrix (const PointsKey &pkey)
std::shared_ptr< NekMatrix< NekDouble > >	CreateGPMatrix (const PointsKey &pkey)
	GaussPoints (const PointsKey &pkey)
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< Points< NekDouble > >	Create (const PointsKey &pkey)

Protected Member Functions
const std::shared_ptr< NekMatrix< NekDouble > >	v_GetI (const PointsKey &pkey) override
const std::shared_ptr< NekMatrix< NekDouble > >	v_GetI (const Array< OneD, const NekDouble > &x) override
const std::shared_ptr< NekMatrix< NekDouble > >	v_GetI (size_t numpoints, const Array< OneD, const NekDouble > &x) override
const std::shared_ptr< NekMatrix< NekDouble > >	v_GetGalerkinProjection (const PointsKey &pkey) override
Protected Member Functions inherited from Nektar::LibUtilities::Points< NekDouble >
virtual void	v_Initialize (void)

Private Member Functions
	GaussPoints ()=delete
	GaussPoints (const GaussPoints &points)=delete
void	v_CalculatePoints () final
void	v_CalculateWeights () final
void	v_CalculateDerivMatrix () final
void	CalculateInterpMatrix (size_t npts, const Array< OneD, const NekDouble > &xpoints, Array< OneD, NekDouble > &interp)
std::shared_ptr< NekMatrix< NekDouble > >	CalculateGalerkinProjectionMatrix (const PointsKey &pkey)

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.
Array< OneD, DataType >	m_points [3]
	Storage for the point locations, allowing for up to a 3D points storage.
Array< OneD, DataType >	m_weights
	Quadrature weights for the weights.
Array< OneD, DataType >	m_bcweights
	Barycentric weights.
MatrixSharedPtrType	m_derivmatrix [3]
	Derivative matrices.
NekManager< PointsKey, NekMatrix< DataType >, PointsKey::opLess >	m_InterpManager
NekManager< PointsKey, NekMatrix< DataType >, PointsKey::opLess >	m_GalerkinProjectionManager

Detailed Description

Definition at line 42 of file GaussPoints.h.

Constructor & Destructor Documentation

~GaussPoints()

Nektar::LibUtilities::GaussPoints::~GaussPoints ( )

inlineoverride

Definition at line 45 of file GaussPoints.h.

    {
    }

GaussPoints() [1/3]

Nektar::LibUtilities::GaussPoints::GaussPoints ( const PointsKey &  pkey )

inline

Definition at line 58 of file GaussPoints.h.

                                       : PointsBaseType(pkey)
    {
        namespace pl = std::placeholders;
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussGaussLegendre),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussRadauMLegendre),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussRadauPLegendre),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussLobattoLegendre),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussGaussChebyshev),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussRadauMChebyshev),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussRadauPChebyshev),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussLobattoChebyshev),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussRadauMAlpha0Beta1),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussRadauMAlpha0Beta2),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussRadauMAlpha1Beta0),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussRadauMAlpha2Beta0),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussKronrodLegendre),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussRadauKronrodMLegendre),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussRadauKronrodMAlpha1Beta0),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussLobattoKronrodLegendre),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eFourierEvenlySpaced),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, ePolyEvenlySpaced),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eBoundaryLayerPoints),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eBoundaryLayerPointsRev),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussLegendreWithMP),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_InterpManager.RegisterCreator(
            PointsKey(0, eGaussLegendreWithM),
            std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eGaussGaussLegendre),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eGaussRadauMLegendre),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eGaussRadauPLegendre),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eGaussLobattoLegendre),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eGaussGaussChebyshev),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eGaussRadauMChebyshev),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eGaussRadauPChebyshev),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eGaussLobattoChebyshev),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eGaussRadauMAlpha0Beta1),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eGaussRadauMAlpha0Beta2),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eGaussRadauMAlpha1Beta0),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eGaussRadauMAlpha2Beta0),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eGaussKronrodLegendre),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eGaussRadauKronrodMLegendre),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eGaussRadauKronrodMAlpha1Beta0),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eGaussLobattoKronrodLegendre),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eFourierEvenlySpaced),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, ePolyEvenlySpaced),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eBoundaryLayerPoints),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eBoundaryLayerPointsRev),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eGaussLegendreWithMP),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
        m_GalerkinProjectionManager.RegisterCreator(
            PointsKey(0, eGaussLegendreWithM),
            std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
    }

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

GaussPoints() [2/3]

Nektar::LibUtilities::GaussPoints::GaussPoints ( )

privatedelete

GaussPoints() [3/3]

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

privatedelete

Member Function Documentation

CalculateGalerkinProjectionMatrix()

std::shared_ptr< NekMatrix< NekDouble > > Nektar::LibUtilities::GaussPoints::CalculateGalerkinProjectionMatrix ( const PointsKey &  pkey )

private

Definition at line 466 of file GaussPoints.cpp.

{
    size_t numpointsfrom = pkey.GetNumPoints();
    size_t numpointsto   = GetNumPoints();
 
    Array<OneD, const NekDouble> weightsfrom;
 
    weightsfrom = PointsManager()[pkey]->GetW();
 
    std::shared_ptr<NekMatrix<NekDouble>> Interp = GetI(pkey);
 
    Array<OneD, NekDouble> GalProj(numpointsfrom * numpointsto);
 
    // set up inner product matrix and multiply by inverse of
    // diagaonal mass matrix
    for (size_t i = 0; i < numpointsto; ++i)
    {
        Vmath::Vmul(numpointsfrom, Interp->GetPtr().data() + i * numpointsfrom,
                    1, &weightsfrom[0], 1, &GalProj[0] + i, numpointsto);
        if (m_weights[i] == 0.0)
        {
            Vmath::Zero(numpointsto, &GalProj[0] + i, numpointsto);
        }
        else
        {
            Vmath::Smul(numpointsfrom, 1.0 / m_weights[i], &GalProj[0] + i,
                        numpointsto, &GalProj[0] + i, numpointsto);
        }
    }
 
    NekDouble *t = GalProj.data();
    std::shared_ptr<NekMatrix<NekDouble>> returnval(
        MemoryManager<NekMatrix<NekDouble>>::AllocateSharedPtr(
            numpointsto, numpointsfrom, t));
 
    return returnval;
}

References Nektar::LibUtilities::Points< NekDouble >::GetI(), Nektar::LibUtilities::PointsKey::GetNumPoints(), Nektar::LibUtilities::Points< NekDouble >::GetNumPoints(), Nektar::LibUtilities::Points< NekDouble >::m_weights, Nektar::LibUtilities::PointsManager(), Vmath::Smul(), Vmath::Vmul(), and Vmath::Zero().

Referenced by CreateGPMatrix().

CalculateInterpMatrix()

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

private

Definition at line 283 of file GaussPoints.cpp.

{
    switch (m_pointsKey.GetPointsType())
    {
        case eGaussGaussLegendre:
            Polylib::Imgj(interp.data(), m_points[0].data(), xpoints.data(),
                          GetNumPoints(), npts, 0.0, 0.0);
            break;
 
        case eGaussLegendreWithMP:
            // first / last line of interp are zeros:
            for (int i = 0; i < interp.size(); ++i)
            {
                interp[i] = 0.0;
            }
            // only interpolate from interior points
            Polylib::Imgj(interp.data() + npts, m_points[0].data() + 1,
                          xpoints.data(), GetNumPoints() - 2, npts, 0.0, 0.0);
            break;
 
        case eGaussLegendreWithM:
            // first / last line of interp are zeros:
            for (int i = 0; i < interp.size(); ++i)
            {
                interp[i] = 0.0;
            }
            // only interpolate from interior points
            Polylib::Imgj(interp.data() + npts, m_points[0].data() + 1,
                          xpoints.data(), GetNumPoints() - 1, npts, 0.0, 0.0);
            break;
 
        case eGaussRadauMLegendre:
            Polylib::Imgrjm(interp.data(), m_points[0].data(), xpoints.data(),
                            GetNumPoints(), npts, 0.0, 0.0);
            break;
 
        case eGaussRadauPLegendre:
            Polylib::Imgrjp(interp.data(), m_points[0].data(), xpoints.data(),
                            GetNumPoints(), npts, 0.0, 0.0);
            break;
 
        case eGaussLobattoLegendre:
            Polylib::Imglj(interp.data(), m_points[0].data(), xpoints.data(),
                           GetNumPoints(), npts, 0.0, 0.0);
            break;
 
        case eGaussGaussChebyshev:
            Polylib::Imgj(interp.data(), m_points[0].data(), xpoints.data(),
                          GetNumPoints(), npts, -0.5, -0.5);
            break;
 
        case eGaussRadauMChebyshev:
            Polylib::Imgrjm(interp.data(), m_points[0].data(), xpoints.data(),
                            GetNumPoints(), npts, -0.5, -0.5);
            break;
 
        case eGaussRadauPChebyshev:
            Polylib::Imgrjp(interp.data(), m_points[0].data(), xpoints.data(),
                            GetNumPoints(), npts, -0.5, -0.5);
            break;
 
        case eGaussLobattoChebyshev:
            Polylib::Imglj(interp.data(), m_points[0].data(), xpoints.data(),
                           GetNumPoints(), npts, -0.5, -0.5);
            break;
 
        case eGaussRadauMAlpha0Beta1:
            Polylib::Imgrjm(interp.data(), m_points[0].data(), xpoints.data(),
                            GetNumPoints(), npts, 0.0, 1.0);
            break;
 
        case eGaussRadauMAlpha0Beta2:
            Polylib::Imgrjm(interp.data(), m_points[0].data(), xpoints.data(),
                            GetNumPoints(), npts, 0.0, 2.0);
            break;
 
        case eGaussRadauMAlpha1Beta0:
            Polylib::Imgrjm(interp.data(), m_points[0].data(), xpoints.data(),
                            GetNumPoints(), npts, 1.0, 0.0);
            break;
 
        case eGaussRadauMAlpha2Beta0:
            Polylib::Imgrjm(interp.data(), m_points[0].data(), xpoints.data(),
                            GetNumPoints(), npts, 2.0, 0.0);
            break;
 
        case eGaussKronrodLegendre:
        case eGaussRadauKronrodMLegendre:
        case eGaussRadauKronrodMAlpha1Beta0:
        case eGaussLobattoKronrodLegendre:
        {
            for (size_t i = 0; i < npts; ++i)
            {
                for (size_t j = 0; j < m_pointsKey.GetNumPoints(); ++j)
                {
                    interp[i + j * npts] =
                        Polylib::laginterp(xpoints[i], j, &m_points[0][0],
                                           m_pointsKey.GetNumPoints());
                }
            }
        }
        break;
 
        default:
            NEKERROR(ErrorUtil::efatal,
                     "Unknown Gauss quadrature point distribution requested");
    }
}

References Nektar::Array< OneD, DataType >::data(), Nektar::ErrorUtil::efatal, Nektar::LibUtilities::eGaussGaussChebyshev, Nektar::LibUtilities::eGaussGaussLegendre, Nektar::LibUtilities::eGaussLegendreWithM, Nektar::LibUtilities::eGaussLegendreWithMP, Nektar::LibUtilities::eGaussLobattoChebyshev, Nektar::LibUtilities::eGaussLobattoKronrodLegendre, Nektar::LibUtilities::eGaussLobattoLegendre, Nektar::LibUtilities::eGaussRadauKronrodMAlpha1Beta0, Nektar::LibUtilities::eGaussRadauKronrodMLegendre, Nektar::LibUtilities::eGaussRadauMAlpha0Beta1, Nektar::LibUtilities::eGaussRadauMChebyshev, Nektar::LibUtilities::eGaussRadauMLegendre, Nektar::LibUtilities::eGaussRadauPChebyshev, Nektar::LibUtilities::eGaussRadauPLegendre, Nektar::LibUtilities::PointsKey::GetNumPoints(), Nektar::LibUtilities::Points< NekDouble >::GetNumPoints(), Nektar::LibUtilities::PointsKey::GetPointsType(), Polylib::Imgj(), Polylib::Imglj(), Polylib::Imgrjm(), Polylib::Imgrjp(), Polylib::laginterp(), Nektar::LibUtilities::Points< NekDouble >::m_points, Nektar::LibUtilities::Points< NekDouble >::m_pointsKey, and NEKERROR.

Referenced by v_GetI().

Create()

std::shared_ptr< Points< NekDouble > > Nektar::LibUtilities::GaussPoints::Create ( const PointsKey &  pkey )

static

Definition at line 394 of file GaussPoints.cpp.

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

CreateGPMatrix()

std::shared_ptr< NekMatrix< NekDouble > > Nektar::LibUtilities::GaussPoints::CreateGPMatrix

(

const PointsKey &

pkey

)

Definition at line 456 of file GaussPoints.cpp.

{
    std::shared_ptr<NekMatrix<NekDouble>> returnval =
        CalculateGalerkinProjectionMatrix(pkey);
 
    // Delegate to function below
    return returnval;
}

References CalculateGalerkinProjectionMatrix().

Referenced by GaussPoints().

CreateMatrix()

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

(

const PointsKey &

pkey

)

Definition at line 404 of file GaussPoints.cpp.

{
    size_t numpoints = pkey.GetNumPoints();
    Array<OneD, const NekDouble> xpoints;
 
    PointsManager()[pkey]->GetPoints(xpoints);
 
    // Delegate to function below
    return GetI(numpoints, xpoints);
}

References Nektar::LibUtilities::Points< NekDouble >::GetI(), Nektar::LibUtilities::PointsKey::GetNumPoints(), and Nektar::LibUtilities::PointsManager().

Referenced by GaussPoints().

v_CalculateDerivMatrix()

void Nektar::LibUtilities::GaussPoints::v_CalculateDerivMatrix ( )

finalprivate

Definition at line 193 of file GaussPoints.cpp.

{
    // Allocate the derivative matrix
    PointsBaseType::v_CalculateDerivMatrix();
 
    size_t numpoints = m_pointsKey.GetNumPoints();
    size_t totpoints = m_pointsKey.GetTotNumPoints();
    Array<OneD, NekDouble> dmtempSharedArray =
        Array<OneD, NekDouble>(totpoints * totpoints);
    NekDouble *dmtemp = dmtempSharedArray.data();
 
    switch (m_pointsKey.GetPointsType())
    {
        case eGaussGaussLegendre:
            Polylib::Dgj(dmtemp, m_points[0].data(), numpoints, 0.0, 0.0);
            break;
 
        case eGaussRadauMLegendre:
            Polylib::Dgrjm(dmtemp, m_points[0].data(), numpoints, 0.0, 0.0);
            break;
 
        case eGaussRadauPLegendre:
            Polylib::Dgrjp(dmtemp, m_points[0].data(), numpoints, 0.0, 0.0);
            break;
 
        case eGaussLobattoLegendre:
            Polylib::Dglj(dmtemp, m_points[0].data(), numpoints, 0.0, 0.0);
            break;
 
        case eGaussGaussChebyshev:
            Polylib::Dgj(dmtemp, m_points[0].data(), numpoints, -0.5, -0.5);
            break;
 
        case eGaussRadauMChebyshev:
            Polylib::Dgrjm(dmtemp, m_points[0].data(), numpoints, -0.5, -0.5);
            break;
 
        case eGaussRadauPChebyshev:
            Polylib::Dgrjp(dmtemp, m_points[0].data(), numpoints, -0.5, -0.5);
            break;
 
        case eGaussLobattoChebyshev:
            Polylib::Dglj(dmtemp, m_points[0].data(), numpoints, -0.5, -0.5);
            break;
 
        case eGaussRadauMAlpha0Beta1:
            Polylib::Dgrjm(dmtemp, m_points[0].data(), numpoints, 0.0, 1.0);
            break;
 
        case eGaussRadauMAlpha0Beta2:
            Polylib::Dgrjm(dmtemp, m_points[0].data(), numpoints, 0.0, 2.0);
            break;
 
        case eGaussRadauMAlpha1Beta0:
            Polylib::Dgrjm(dmtemp, m_points[0].data(), numpoints, 1.0, 0.0);
            break;
 
        case eGaussRadauMAlpha2Beta0:
            Polylib::Dgrjm(dmtemp, m_points[0].data(), numpoints, 2.0, 0.0);
            break;
 
        case eGaussLegendreWithMP:
        case eGaussLegendreWithM:
        case eGaussKronrodLegendre:
        case eGaussRadauKronrodMLegendre:
        case eGaussRadauKronrodMAlpha1Beta0:
        case eGaussLobattoKronrodLegendre:
        {
            for (size_t i = 0; i < m_pointsKey.GetNumPoints(); ++i)
            {
                for (size_t j = 0; j < m_pointsKey.GetNumPoints(); ++j)
                {
                    (*m_derivmatrix[0])(i, j) = Polylib::laginterpderiv(
                        m_points[0][i], j, &m_points[0][0],
                        m_pointsKey.GetNumPoints());
                }
            }
            return;
        }
        break;
 
        default:
            NEKERROR(ErrorUtil::efatal,
                     "Unknown Gauss quadrature point distribution requested");
    }
 
    std::copy(dmtemp, dmtemp + totpoints * totpoints,
              m_derivmatrix[0]->begin());
}

References Polylib::Dgj(), Polylib::Dglj(), Polylib::Dgrjm(), Polylib::Dgrjp(), Nektar::ErrorUtil::efatal, Nektar::LibUtilities::eGaussGaussChebyshev, Nektar::LibUtilities::eGaussGaussLegendre, Nektar::LibUtilities::eGaussLegendreWithM, Nektar::LibUtilities::eGaussLegendreWithMP, Nektar::LibUtilities::eGaussLobattoChebyshev, Nektar::LibUtilities::eGaussLobattoKronrodLegendre, Nektar::LibUtilities::eGaussLobattoLegendre, Nektar::LibUtilities::eGaussRadauKronrodMAlpha1Beta0, Nektar::LibUtilities::eGaussRadauKronrodMLegendre, Nektar::LibUtilities::eGaussRadauMAlpha0Beta1, Nektar::LibUtilities::eGaussRadauMChebyshev, Nektar::LibUtilities::eGaussRadauMLegendre, Nektar::LibUtilities::eGaussRadauPChebyshev, Nektar::LibUtilities::eGaussRadauPLegendre, Nektar::LibUtilities::PointsKey::GetNumPoints(), Nektar::LibUtilities::PointsKey::GetPointsType(), Nektar::LibUtilities::PointsKey::GetTotNumPoints(), Polylib::laginterpderiv(), Nektar::LibUtilities::Points< NekDouble >::m_derivmatrix, Nektar::LibUtilities::Points< NekDouble >::m_points, Nektar::LibUtilities::Points< NekDouble >::m_pointsKey, and NEKERROR.

v_CalculatePoints()

void Nektar::LibUtilities::GaussPoints::v_CalculatePoints ( )

finalprivatevirtual

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

Definition at line 79 of file GaussPoints.cpp.

{
    // Allocate the storage for points and weights
    PointsBaseType::v_CalculatePoints();
    PointsBaseType::v_CalculateWeights();
 
    size_t numpoints = m_pointsKey.GetNumPoints();
 
    switch (m_pointsKey.GetPointsType())
    {
        case eGaussGaussLegendre:
            Polylib::zwgj(m_points[0].data(), m_weights.data(), numpoints, 0.0,
                          0.0);
            break;
        case eGaussLegendreWithMP:
            Polylib::zwgj(m_points[0].data() + 1, m_weights.data() + 1,
                          numpoints - 2, 0.0, 0.0);
            m_points[0][0]             = -1.0;
            m_points[0][numpoints - 1] = 1.0;
            m_weights[0]               = 0.0;
            m_weights[numpoints - 1]   = 0.0;
            break;
        case eGaussLegendreWithM:
            Polylib::zwgj(m_points[0].data() + 1, m_weights.data() + 1,
                          numpoints - 1, 0.0, 0.0);
            m_points[0][0] = -1.0;
            m_weights[0]   = 0.0;
            break;
        case eGaussRadauMLegendre:
            Polylib::zwgrjm(m_points[0].data(), m_weights.data(), numpoints,
                            0.0, 0.0);
            break;
 
        case eGaussRadauPLegendre:
            Polylib::zwgrjp(m_points[0].data(), m_weights.data(), numpoints,
                            0.0, 0.0);
            break;
 
        case eGaussLobattoLegendre:
            Polylib::zwglj(m_points[0].data(), m_weights.data(), numpoints, 0.0,
                           0.0);
            break;
 
        case eGaussGaussChebyshev:
            Polylib::zwgj(m_points[0].data(), m_weights.data(), numpoints, -0.5,
                          -0.5);
            break;
 
        case eGaussRadauMChebyshev:
            Polylib::zwgrjm(m_points[0].data(), m_weights.data(), numpoints,
                            -0.5, -0.5);
            break;
 
        case eGaussRadauPChebyshev:
            Polylib::zwgrjp(m_points[0].data(), m_weights.data(), numpoints,
                            -0.5, -0.5);
            break;
 
        case eGaussLobattoChebyshev:
            Polylib::zwglj(m_points[0].data(), m_weights.data(), numpoints,
                           -0.5, -0.5);
            break;
 
        case eGaussRadauMAlpha0Beta1:
            Polylib::zwgrjm(m_points[0].data(), m_weights.data(), numpoints,
                            0.0, 1.0);
            break;
 
        case eGaussRadauMAlpha0Beta2:
            Polylib::zwgrjm(m_points[0].data(), m_weights.data(), numpoints,
                            0.0, 2.0);
            break;
 
        case eGaussRadauMAlpha1Beta0:
            Polylib::zwgrjm(m_points[0].data(), m_weights.data(), numpoints,
                            1.0, 0.0);
            break;
 
        case eGaussRadauMAlpha2Beta0:
            Polylib::zwgrjm(m_points[0].data(), m_weights.data(), numpoints,
                            2.0, 0.0);
            break;
 
        case eGaussKronrodLegendre:
            Polylib::zwgk(m_points[0].data(), m_weights.data(), numpoints, 0.0,
                          0.0);
            break;
 
        case eGaussRadauKronrodMLegendre:
            Polylib::zwrk(m_points[0].data(), m_weights.data(), numpoints, 0.0,
                          0.0);
            break;
 
        case eGaussRadauKronrodMAlpha1Beta0:
            Polylib::zwrk(m_points[0].data(), m_weights.data(), numpoints, 1.0,
                          0.0);
            break;
 
        case eGaussLobattoKronrodLegendre:
            Polylib::zwlk(m_points[0].data(), m_weights.data(), numpoints, 0.0,
                          0.0);
            break;
 
        default:
            NEKERROR(ErrorUtil::efatal,
                     "Unknown Gauss quadrature point distribution requested");
    }
}

References Nektar::Array< OneD, DataType >::data(), Nektar::ErrorUtil::efatal, Nektar::LibUtilities::eGaussGaussChebyshev, Nektar::LibUtilities::eGaussGaussLegendre, Nektar::LibUtilities::eGaussLegendreWithM, Nektar::LibUtilities::eGaussLegendreWithMP, Nektar::LibUtilities::eGaussLobattoChebyshev, Nektar::LibUtilities::eGaussLobattoKronrodLegendre, Nektar::LibUtilities::eGaussLobattoLegendre, Nektar::LibUtilities::eGaussRadauKronrodMAlpha1Beta0, Nektar::LibUtilities::eGaussRadauKronrodMLegendre, Nektar::LibUtilities::eGaussRadauMAlpha0Beta1, Nektar::LibUtilities::eGaussRadauMChebyshev, Nektar::LibUtilities::eGaussRadauMLegendre, Nektar::LibUtilities::eGaussRadauPChebyshev, Nektar::LibUtilities::eGaussRadauPLegendre, Nektar::LibUtilities::PointsKey::GetNumPoints(), Nektar::LibUtilities::PointsKey::GetPointsType(), Nektar::LibUtilities::Points< NekDouble >::m_points, Nektar::LibUtilities::Points< NekDouble >::m_pointsKey, Nektar::LibUtilities::Points< NekDouble >::m_weights, NEKERROR, Nektar::LibUtilities::Points< NekDouble >::v_CalculatePoints(), Nektar::LibUtilities::Points< NekDouble >::v_CalculateWeights(), Polylib::zwgj(), Polylib::zwgk(), Polylib::zwglj(), Polylib::zwgrjm(), Polylib::zwgrjp(), Polylib::zwlk(), and Polylib::zwrk().

v_CalculateWeights()

void Nektar::LibUtilities::GaussPoints::v_CalculateWeights ( )

finalprivatevirtual

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

Definition at line 188 of file GaussPoints.cpp.

{
    // For Gauss Quadrature, this is done as part of the points computation
}

v_GetGalerkinProjection()

const std::shared_ptr< NekMatrix< NekDouble > > Nektar::LibUtilities::GaussPoints::v_GetGalerkinProjection ( const PointsKey &  pkey )

overrideprotected

Definition at line 450 of file GaussPoints.cpp.

{
    return m_GalerkinProjectionManager[pkey];
}

References Nektar::LibUtilities::Points< NekDouble >::m_GalerkinProjectionManager.

v_GetI() [1/3]

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

overrideprotected

Definition at line 425 of file GaussPoints.cpp.

{
    size_t numpoints = 1;
 
    // Delegate to function below
    return GetI(numpoints, x);
}

References Nektar::LibUtilities::Points< NekDouble >::GetI().

v_GetI() [2/3]

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

overrideprotected

Definition at line 416 of file GaussPoints.cpp.

{
    ASSERTL0(pkey.GetPointsDim() == 1,
             "Gauss Points can only interp to other 1d point distributions");
 
    return m_InterpManager[pkey];
}

References ASSERTL0, Nektar::LibUtilities::PointsKey::GetPointsDim(), and Nektar::LibUtilities::Points< NekDouble >::m_InterpManager.

v_GetI() [3/3]

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

overrideprotected

Definition at line 434 of file GaussPoints.cpp.

{
    Array<OneD, NekDouble> interp(GetNumPoints() * numpoints);
 
    CalculateInterpMatrix(numpoints, x, interp);
 
    NekDouble *t = interp.data();
    size_t np    = GetNumPoints();
    std::shared_ptr<NekMatrix<NekDouble>> returnval(
        MemoryManager<NekMatrix<NekDouble>>::AllocateSharedPtr(numpoints, np,
                                                               t));
 
    return returnval;
}

References CalculateInterpMatrix(), and Nektar::LibUtilities::Points< NekDouble >::GetNumPoints().

Member Data Documentation

initPointsManager

bool Nektar::LibUtilities::GaussPoints::initPointsManager

staticprivate

Definition at line 206 of file GaussPoints.h.