Nektar::Collections::PhysDeriv_MatrixFree Class Reference

Phys deriv operator using matrix free operators. More...

Inheritance diagram for Nektar::Collections::PhysDeriv_MatrixFree:

Public Member Functions
	~PhysDeriv_MatrixFree () final
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
virtual void	CheckFactors (StdRegions::FactorMap factors, int coll_phys_offset)
	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 COLLECTIONS_EXPORT	~Operator ()
unsigned int	GetWspSize ()
	Get the size of the required workspace. More...
unsigned int	GetNumElmt ()
	Get expansion pointer. More...
StdRegions::StdExpansionSharedPtr	GetExpSharedPtr ()
	Get expansion pointer. More...

Private Member Functions
	PhysDeriv_MatrixFree (vector< StdRegions::StdExpansionSharedPtr > pCollExp, CoalescedGeomDataSharedPtr pGeomData, StdRegions::FactorMap factors)
Private Member Functions inherited from Nektar::Collections::MatrixFreeOneInMultiOut
	MatrixFreeOneInMultiOut (const unsigned int coordim, const unsigned int nIn, const unsigned int nOut, const unsigned int nCollSize)
	Constructor. More...

Private Attributes
std::shared_ptr< MatrixFree::PhysDeriv >	m_oper
Private Attributes inherited from Nektar::Collections::MatrixFreeOneInMultiOut
unsigned short	m_coordim
	coordinates dimension More...
Array< OneD, NekDouble >	m_input
	padded input/output vectors More...
Array< OneD, Array< OneD, NekDouble > >	m_output

Additional Inherited Members
Protected Attributes inherited from Nektar::Collections::Operator
bool	m_isDeformed
StdRegions::StdExpansionSharedPtr	m_stdExp
unsigned int	m_numElmt
unsigned int	m_nqe
unsigned int	m_wspSize

Detailed Description

Phys deriv operator using matrix free operators.

Definition at line 277 of file PhysDeriv.cpp.

Constructor & Destructor Documentation

~PhysDeriv_MatrixFree()

Nektar::Collections::PhysDeriv_MatrixFree::~PhysDeriv_MatrixFree ( )

inlinefinal

Definition at line 282 of file PhysDeriv.cpp.

    {
    }

PhysDeriv_MatrixFree()

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

inlineprivate

Definition at line 359 of file PhysDeriv.cpp.

        : Operator(pCollExp, pGeomData, factors),
          MatrixFreeOneInMultiOut(pCollExp[0]->GetCoordim(),
                                  pCollExp[0]->GetStdExp()->GetTotPoints(),
                                  pCollExp[0]->GetStdExp()->GetTotPoints(),
                                  pCollExp.size())
    {
        // Check if deformed
        bool deformed{pGeomData->IsDeformed(pCollExp)};
        const auto dim = pCollExp[0]->GetStdExp()->GetShapeDimension();
        
        if(m_isPadded == false) // declare local space non-padded case
        {
            int nOut = pCollExp[0]->GetStdExp()->GetTotPoints();
            m_output = Array<OneD, Array<OneD, NekDouble>> (m_coordim);
            m_output[0] = Array<OneD, NekDouble>{nOut * m_nElmtPad, 0.0};
            if(m_coordim == 2)
            {
                m_output[1] = Array<OneD, NekDouble>{nOut * m_nElmtPad, 0.0};
            }
            else if (m_coordim == 3)
            {
                m_output[1] = Array<OneD, NekDouble>{nOut * m_nElmtPad, 0.0};
                m_output[2] = Array<OneD, NekDouble>{nOut * m_nElmtPad, 0.0};
            }
        }
 
        // 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 = "PhysDeriv";
        op_string += MatrixFree::GetOpstring(shapeType, deformed);
        auto oper = MatrixFree::GetOperatorFactory().
            CreateInstance(op_string, basis, m_nElmtPad);
        
        // Set derivative factors
        oper->SetDF(pGeomData->GetDerivFactorsInterLeave
                    (pCollExp,m_nElmtPad));
        
        m_oper = std::dynamic_pointer_cast<MatrixFree::PhysDeriv>(oper);
        ASSERTL0(m_oper, "Failed to cast pointer.");
        
    }

References ASSERTL0, and Nektar::Collections::GetOperatorFactory().

Member Function Documentation

CheckFactors()

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

inlinevirtual

Check the validity of the supplied factor map.

Implements Nektar::Collections::Operator.

Definition at line 349 of file PhysDeriv.cpp.

    {
        boost::ignore_unused(factors, coll_phys_offset);
        ASSERTL0(false, "Not valid for this operator.");
    }

References ASSERTL0.

operator()() [1/2]

void Nektar::Collections::PhysDeriv_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 286 of file PhysDeriv.cpp.

    {
        boost::ignore_unused(wsp);
 
        if (m_isPadded)
        {
            // copy into padded vector
            Vmath::Vcopy(m_nIn, input, 1, m_input, 1);
            (*m_oper)(m_input, m_output);
        }
        else
        {
            (*m_oper)(input, m_output);
        }
 
        // currently using temporary local temporary space for output
        // to allow for other operator call below which is
        // directionally dependent
        switch(m_coordim)
        {
        case 1:
            Vmath::Vcopy(m_nOut, m_output[0], 1, output0, 1);
            break;
        case 2:
            Vmath::Vcopy(m_nOut, m_output[0], 1, output0, 1);
            Vmath::Vcopy(m_nOut, m_output[1], 1, output1, 1);
            break;
        case 3:
            Vmath::Vcopy(m_nOut, m_output[0], 1, output0, 1);
            Vmath::Vcopy(m_nOut, m_output[1], 1, output1, 1);
            Vmath::Vcopy(m_nOut, m_output[2], 1, output2, 1);
            break;
        default:
            NEKERROR(ErrorUtil::efatal,
                     "Unknown coordinate dimension");
            break;
        }
    }

References NEKERROR, and Vmath::Vcopy().

operator()() [2/2]

void Nektar::Collections::PhysDeriv_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 330 of file PhysDeriv.cpp.

    {
        boost::ignore_unused(wsp);
        if (m_isPadded)
        {
            // copy into padded vector
            Vmath::Vcopy(m_nIn, input, 1, m_input, 1);
            (*m_oper)(m_input, m_output);
        }
        else
        {
            (*m_oper)(input, m_output);
        }
        Vmath::Vcopy(m_nOut, m_output[dir], 1, output, 1);
    }

References Vmath::Vcopy().

Member Data Documentation

m_oper

std::shared_ptr<MatrixFree::PhysDeriv> Nektar::Collections::PhysDeriv_MatrixFree::m_oper

private

Definition at line 357 of file PhysDeriv.cpp.