PhysInterp1DScaled operator using matrix free implementation. More...

Inheritance diagram for Nektar::Collections::PhysInterp1DScaled_MatrixFree:

Public Member Functions
	~PhysInterp1DScaled_MatrixFree () final=default

void	operator() (const Array< OneD, const NekDouble > &input, Array< OneD, NekDouble > &output0, Array< OneD, NekDouble > &output1, Array< OneD, NekDouble > &output2, Array< OneD, NekDouble > &wsp) final
	Perform operation. More...

void	operator() (int dir, const Array< OneD, const NekDouble > &input, Array< OneD, NekDouble > &output, Array< OneD, NekDouble > &wsp) final

void	UpdateFactors (StdRegions::FactorMap factors) override
	Update the supplied factor map. More...

Public Member Functions inherited from Nektar::Collections::Operator
	Operator (std::vector< StdRegions::StdExpansionSharedPtr > pCollExp, std::shared_ptr< CoalescedGeomData > GeomData, StdRegions::FactorMap factors)
	Constructor. More...

virtual	~Operator ()=default

virtual COLLECTIONS_EXPORT void	operator() (const Array< OneD, const NekDouble > &input, Array< OneD, NekDouble > &output0, Array< OneD, NekDouble > &output1, Array< OneD, NekDouble > &output2, Array< OneD, NekDouble > &wsp=NullNekDouble1DArray)=0
	Perform operation. More...

virtual COLLECTIONS_EXPORT void	operator() (int dir, const Array< OneD, const NekDouble > &input, Array< OneD, NekDouble > &output, Array< OneD, NekDouble > &wsp=NullNekDouble1DArray)=0

virtual COLLECTIONS_EXPORT void	UpdateFactors (StdRegions::FactorMap factors)
	Update the supplied factor map. More...

virtual COLLECTIONS_EXPORT void	UpdateVarcoeffs (StdRegions::VarCoeffMap &varcoeffs)
	Update the supplied variable coefficients. More...

unsigned int	GetWspSize ()
	Get the size of the required workspace. More...

unsigned int	GetNumElmt ()
	Get number of elements. More...

StdRegions::StdExpansionSharedPtr	GetExpSharedPtr ()
	Get expansion pointer. More...

unsigned int	GetInputSize (bool defaultIn=true)

unsigned int	GetOutputSize (bool defaultOut=true)

Private Member Functions
	PhysInterp1DScaled_MatrixFree (vector< StdRegions::StdExpansionSharedPtr > pCollExp, CoalescedGeomDataSharedPtr pGeomData, StdRegions::FactorMap factors)

Private Member Functions inherited from Nektar::Collections::MatrixFreeBase
	MatrixFreeBase ()=default

	MatrixFreeBase (const unsigned int nIn, const unsigned int nOut, const unsigned int nCollSize)

void	UpdateFactors (StdRegions::FactorMap factors) override
	Update the supplied factor map. More...

	PhysInterp1DScaled_Helper ()

Private Attributes
std::shared_ptr< MatrixFree::PhysInterp1DScaled >	m_oper

Private Attributes inherited from Nektar::Collections::MatrixFreeBase
bool	m_isPadded {false}
	flag for padding More...

unsigned int	m_nElmtPad
	size after padding More...

unsigned int	m_nIn
	actural size of input array More...

unsigned int	m_nOut
	actural size of output array More...

Private Attributes inherited from Nektar::Collections::PhysInterp1DScaled_Helper
StdRegions::FactorMap	m_factors

Additional Inherited Members
Protected Attributes inherited from Nektar::Collections::Operator
bool	m_isDeformed

StdRegions::StdExpansionSharedPtr	m_stdExp

unsigned int	m_numElmt
	number of elements that the operator is applied on More...

unsigned int	m_nqe

unsigned int	m_wspSize

unsigned int	m_inputSize
	number of modes or quadrature points that are passed as input to an operator More...

unsigned int	m_outputSize
	number of modes or quadrature points that are taken as output from an operator More...

unsigned int	m_inputSizeOther
	Number of modes or quadrature points, opposite to m_inputSize. More...

unsigned int	m_outputSizeOther
	Number of modes or quadrature points, opposite to m_outputSize. More...

Detailed Description

PhysInterp1DScaled operator using matrix free implementation.

Definition at line 121 of file PhysInterp1DScaled.cpp.

Constructor & Destructor Documentation

◆ ~PhysInterp1DScaled_MatrixFree()

Nektar::Collections::PhysInterp1DScaled_MatrixFree::~PhysInterp1DScaled_MatrixFree ( )

finaldefault

◆ PhysInterp1DScaled_MatrixFree()

Nektar::Collections::PhysInterp1DScaled_MatrixFree::PhysInterp1DScaled_MatrixFree	(	vector< StdRegions::StdExpansionSharedPtr >	pCollExp,
		CoalescedGeomDataSharedPtr	pGeomData,
		StdRegions::FactorMap	factors
	)

inlineprivate

Definition at line 181 of file PhysInterp1DScaled.cpp.

        : Operator(pCollExp, pGeomData, factors),
          MatrixFreeBase(pCollExp[0]->GetStdExp()->GetTotPoints(),
                         pCollExp[0]->GetStdExp()->GetTotPoints(),
                         pCollExp.size()),
          PhysInterp1DScaled_Helper()
    {
        const auto dim = pCollExp[0]->GetStdExp()->GetShapeDimension();
 
        // Definition of basis vector
        std::vector<LibUtilities::BasisSharedPtr> basis(dim);
        for (unsigned int i = 0; i < dim; ++i)
        {
            basis[i] = pCollExp[0]->GetBasis(i);
        }
 
        // Get shape type
        auto shapeType = pCollExp[0]->GetStdExp()->DetShapeType();
 
        // Generate operator string and create operator.
        std::string op_string = "PhysInterp1DScaled";
        op_string += MatrixFree::GetOpstring(shapeType, false);
        auto oper = MatrixFree::GetOperatorFactory().CreateInstance(
            op_string, basis, pCollExp.size());
 
        m_oper =
            std::dynamic_pointer_cast<MatrixFree::PhysInterp1DScaled>(oper);
        ASSERTL0(m_oper, "Failed to cast pointer.");
 
        NekDouble scale; // declaration of the scaling factor to be used for the
                         // output size
        if (factors[StdRegions::eFactorConst] != 0)
        {
            m_factors = factors;
            scale     = m_factors[StdRegions::eFactorConst];
        }
        else
        {
            scale                               = 1.5;
            m_factors[StdRegions::eFactorConst] = 1.5;
        }
        // Update the scaling factor inside MatrixFreeOps using the SetLamda
        // routine
        m_oper->SetScalingFactor(scale);
        m_oper->SetUpInterp1D(basis, scale);
        m_nOut = m_nElmtPad * m_outputSize / m_numElmt;
    }

References ASSERTL0, Nektar::LibUtilities::NekFactory< tKey, tBase, tParam >::CreateInstance(), Nektar::StdRegions::eFactorConst, Nektar::VarcoeffHashingTest::factors, Nektar::Collections::GetOperatorFactory(), Nektar::Collections::PhysInterp1DScaled_Helper::m_factors, Nektar::Collections::MatrixFreeBase::m_nElmtPad, Nektar::Collections::MatrixFreeBase::m_nOut, Nektar::Collections::Operator::m_numElmt, m_oper, and Nektar::Collections::Operator::m_outputSize.

Member Function Documentation

◆ operator()() [1/2]

void Nektar::Collections::PhysInterp1DScaled_MatrixFree::operator()	(	const Array< OneD, const NekDouble > &	input,
		Array< OneD, NekDouble > &	output0,
		Array< OneD, NekDouble > &	output1,
		Array< OneD, NekDouble > &	output2,
		Array< OneD, NekDouble > &	wsp
	)

inlinefinalvirtual

Perform operation.

Implements Nektar::Collections::Operator.

Definition at line 130 of file PhysInterp1DScaled.cpp.

    {
        (*m_oper)(input, output0);
    }

◆ operator()() [2/2]

void Nektar::Collections::PhysInterp1DScaled_MatrixFree::operator()	(	int	dir,
		const Array< OneD, const NekDouble > &	input,
		Array< OneD, NekDouble > &	output,
		Array< OneD, NekDouble > &	wsp
	)

inlinefinalvirtual

Implements Nektar::Collections::Operator.

Definition at line 139 of file PhysInterp1DScaled.cpp.

    {
        NEKERROR(ErrorUtil::efatal,
                 "PhysInterp1DScaled_MatrixFree: Not valid for this operator.");
    }

References Nektar::ErrorUtil::efatal, and NEKERROR.

◆ UpdateFactors()

void Nektar::Collections::PhysInterp1DScaled_MatrixFree::UpdateFactors ( StdRegions::FactorMap factors )

inlineoverridevirtual

Update the supplied factor map.

Reimplemented from Nektar::Collections::Operator.

Definition at line 148 of file PhysInterp1DScaled.cpp.

    {
 
        // Set lambda for this call
        auto x = factors.find(StdRegions::eFactorConst);
        ASSERTL1(
            x != factors.end(),
            "Constant factor not defined: " +
                std::string(
                    StdRegions::ConstFactorTypeMap[StdRegions::eFactorConst]));
        // Update the factors member inside PhysInterp1DScaled_MatrixFree
        // class
        m_factors = factors;
 
        const auto dim = m_stdExp->GetShapeDimension();
 
        // Definition of basis vector
        std::vector<LibUtilities::BasisSharedPtr> basis(dim);
        for (unsigned int i = 0; i < dim; ++i)
        {
            basis[i] = m_stdExp->GetBasis(i);
        }
 
        // Update the scaling factor inside MatrixFreeOps using the SetLamda
        // routine
        m_oper->SetScalingFactor(x->second);
        m_oper->SetUpInterp1D(basis, m_factors[StdRegions::eFactorConst]);
        m_nOut = m_nElmtPad * m_outputSize / m_numElmt;
    }

References ASSERTL1, Nektar::StdRegions::ConstFactorTypeMap, Nektar::StdRegions::eFactorConst, Nektar::VarcoeffHashingTest::factors, Nektar::Collections::PhysInterp1DScaled_Helper::m_factors, Nektar::Collections::MatrixFreeBase::m_nElmtPad, Nektar::Collections::MatrixFreeBase::m_nOut, Nektar::Collections::Operator::m_numElmt, m_oper, Nektar::Collections::Operator::m_outputSize, and Nektar::Collections::Operator::m_stdExp.

Member Data Documentation

◆ m_oper

std::shared_ptr<MatrixFree::PhysInterp1DScaled> Nektar::Collections::PhysInterp1DScaled_MatrixFree::m_oper

private

Definition at line 179 of file PhysInterp1DScaled.cpp.

Referenced by PhysInterp1DScaled_MatrixFree(), and UpdateFactors().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members