Inner product operator using standard matrix approach. More...

Inheritance diagram for Nektar::Collections::IProductWRTBase_StdMat:

Public Member Functions
	~IProductWRTBase_StdMat () final=default

void	operator() (const Array< OneD, const NekDouble > &input, Array< OneD, NekDouble > &output, 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	CheckFactors (StdRegions::FactorMap factors, int coll_phys_offset) override
	Check the validity of 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	CheckFactors (StdRegions::FactorMap factors, int coll_phys_offset)=0
	Check the validity of the supplied factor map. 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 ()

unsigned int	GetOutputSize ()

Protected Attributes
DNekMatSharedPtr	m_mat

Array< OneD, const NekDouble >	m_jac

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

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

Additional Inherited Members
Protected Member Functions inherited from Nektar::Collections::Operator
virtual int	v_GetInputSize ()
	This purely virtual function needs to be set-up for every operator inside Collections. It is responsible for returning the size of input collection, that the operator is applied on either in physical or coefficient space. More...

virtual int	v_GetOutputSize ()
	This purely virtual function needs to be set-up for every operator inside Collections. It is responsible for returning the output size either in physical or coefficient space of an operator inside Collections. More...

Protected Member Functions inherited from Nektar::Collections::IProductWRTBase_Helper
	IProductWRTBase_Helper ()

Detailed Description

Inner product operator using standard matrix approach.

Definition at line 75 of file IProductWRTBase.cpp.

Constructor & Destructor Documentation

◆ ~IProductWRTBase_StdMat()

Nektar::Collections::IProductWRTBase_StdMat::~IProductWRTBase_StdMat ( )

finaldefault

◆ IProductWRTBase_StdMat()

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

inlineprivate

Definition at line 130 of file IProductWRTBase.cpp.

        : Operator(pCollExp, pGeomData, factors), IProductWRTBase_Helper()
    {
        m_jac = pGeomData->GetJac(pCollExp);
        StdRegions::StdMatrixKey key(StdRegions::eIProductWRTBase,
                                     m_stdExp->DetShapeType(), *m_stdExp);
        m_mat     = m_stdExp->GetStdMatrix(key);
        m_nqe     = m_stdExp->GetTotPoints();
        m_wspSize = m_nqe * m_numElmt;
    }

References Nektar::StdRegions::eIProductWRTBase, m_jac, m_mat, Nektar::Collections::Operator::m_nqe, Nektar::Collections::Operator::m_numElmt, Nektar::Collections::Operator::m_stdExp, and Nektar::Collections::Operator::m_wspSize.

Member Function Documentation

◆ CheckFactors()

void Nektar::Collections::IProductWRTBase_StdMat::CheckFactors	(	StdRegions::FactorMap	factors,
		int	coll_phys_offset
	)

inlineoverridevirtual

Check the validity of the supplied factor map.

Implements Nektar::Collections::Operator.

Definition at line 119 of file IProductWRTBase.cpp.

    {
        ASSERTL0(false, "Not valid for this operator.");
    }

References ASSERTL0.

◆ operator()() [1/2]

void Nektar::Collections::IProductWRTBase_StdMat::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 83 of file IProductWRTBase.cpp.

    {
        ASSERTL1(wsp.size() == m_wspSize, "Incorrect workspace size");
 
        if (m_isDeformed)
        {
            Vmath::Vmul(m_jac.size(), m_jac, 1, input, 1, wsp, 1);
        }
        else
        {
            Array<OneD, NekDouble> tmp;
            for (int e = 0; e < m_numElmt; ++e)
            {
                Vmath::Smul(m_nqe, m_jac[e], input + e * m_nqe, 1,
                            tmp = wsp + e * m_nqe, 1);
            }
        }
 
        Blas::Dgemm('N', 'N', m_mat->GetRows(), m_numElmt, m_mat->GetColumns(),
                    1.0, m_mat->GetRawPtr(), m_mat->GetRows(), wsp.get(),
                    m_stdExp->GetTotPoints(), 0.0, output.get(),
                    m_stdExp->GetNcoeffs());
    }

References ASSERTL1, Blas::Dgemm(), Nektar::Collections::Operator::m_isDeformed, m_jac, m_mat, Nektar::Collections::Operator::m_nqe, Nektar::Collections::Operator::m_numElmt, Nektar::Collections::Operator::m_stdExp, Nektar::Collections::Operator::m_wspSize, Vmath::Smul(), and Vmath::Vmul().

◆ operator()() [2/2]

void Nektar::Collections::IProductWRTBase_StdMat::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 111 of file IProductWRTBase.cpp.

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

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

Member Data Documentation

◆ m_jac

Array<OneD, const NekDouble> Nektar::Collections::IProductWRTBase_StdMat::m_jac

protected

Definition at line 127 of file IProductWRTBase.cpp.

Referenced by IProductWRTBase_StdMat(), and operator()().

◆ m_mat

DNekMatSharedPtr Nektar::Collections::IProductWRTBase_StdMat::m_mat

protected

Definition at line 126 of file IProductWRTBase.cpp.

Referenced by IProductWRTBase_StdMat(), and operator()().

Public Member Functions

Protected Attributes

Private Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ~IProductWRTBase_StdMat()

◆ IProductWRTBase_StdMat()

Member Function Documentation

◆ CheckFactors()

◆ operator()() [1/2]

◆ operator()() [2/2]

Member Data Documentation

◆ m_jac

◆ m_mat