Nektar::Collections::PhysDeriv_IterPerExp Class Reference

Phys deriv operator using element-wise operation. More...

Inheritance diagram for Nektar::Collections::PhysDeriv_IterPerExp:

Collaboration diagram for Nektar::Collections::PhysDeriv_IterPerExp:

Public Member Functions
virtual	~PhysDeriv_IterPerExp ()
virtual void	operator() (const Array< OneD, const NekDouble > &input, Array< OneD, NekDouble > &output0, Array< OneD, NekDouble > &output1, Array< OneD, NekDouble > &output2, Array< OneD, NekDouble > &wsp)
	Perform operation. More...
virtual void	operator() (int dir, const Array< OneD, const NekDouble > &input, Array< OneD, NekDouble > &output, Array< OneD, NekDouble > &wsp)
Public Member Functions inherited from Nektar::Collections::Operator
	Operator (std::vector< StdRegions::StdExpansionSharedPtr > pCollExp, boost::shared_ptr< CoalescedGeomData > GeomData)
	Constructor. More...
virtual COLLECTIONS_EXPORT	~Operator ()
int	GetWspSize ()
	Get the size of the required workspace. More...

Protected Attributes
Array< TwoD, const NekDouble >	m_derivFac
int	m_dim
int	m_coordim
Protected Attributes inherited from Nektar::Collections::Operator
StdRegions::StdExpansionSharedPtr	m_stdExp
unsigned int	m_numElmt
unsigned int	m_wspSize

Private Member Functions
	PhysDeriv_IterPerExp (vector< StdRegions::StdExpansionSharedPtr > pCollExp, CoalescedGeomDataSharedPtr pGeomData)

Detailed Description

Phys deriv operator using element-wise operation.

Definition at line 229 of file PhysDeriv.cpp.

Constructor & Destructor Documentation

virtual Nektar::Collections::PhysDeriv_IterPerExp::~PhysDeriv_IterPerExp ( )

inlinevirtual

Definition at line 234 of file PhysDeriv.cpp.

        {
        }

Nektar::Collections::PhysDeriv_IterPerExp::PhysDeriv_IterPerExp ( vector< StdRegions::StdExpansionSharedPtr >  pCollExp, CoalescedGeomDataSharedPtr  pGeomData  )

inlineprivate

Definition at line 322 of file PhysDeriv.cpp.

            : Operator(pCollExp, pGeomData)
        {
            int nqtot = 1;
            LibUtilities::PointsKeyVector PtsKey = m_stdExp->GetPointsKeys();
            m_dim = PtsKey.size();
            m_coordim = m_stdExp->GetCoordim();

            for(int i = 0; i < m_dim; ++i)
            {
                nqtot *= PtsKey[i].GetNumPoints();
            }
            m_derivFac = pGeomData->GetDerivFactors(pCollExp);
            m_wspSize = 3*nqtot*m_numElmt;
        }

Member Function Documentation

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

inlinevirtual

Perform operation.

Implements Nektar::Collections::Operator.

Definition at line 238 of file PhysDeriv.cpp.

References Vmath::Vmul(), and Vmath::Vvtvp().

        {
            int nPhys = m_stdExp->GetTotPoints();
            int ntot = m_numElmt*nPhys;
            Array<OneD, NekDouble> tmp0,tmp1,tmp2;
            Array<OneD, Array<OneD, NekDouble> > Diff(3);
            Array<OneD, Array<OneD, NekDouble> > out(3);
            out[0] = output0;  out[1] = output1;  out[2] = output2;

            for(int i = 0; i < m_dim; ++i)
            {
                Diff[i] = wsp + i*ntot;
            }

            // calculate local derivatives
            for (int i = 0; i < m_numElmt; ++i)
            {
                m_stdExp->PhysDeriv(input + i*nPhys,
                                    tmp0 = Diff[0] + i*nPhys,
                                    tmp1 = Diff[1] + i*nPhys,
                                    tmp2 = Diff[2] + i*nPhys);
            }

            // calculate full derivative
            for(int i = 0; i < m_coordim; ++i)
            {
                Vmath::Vmul(ntot,m_derivFac[i*m_dim],1,Diff[0],1,out[i],1);
                for(int j = 1; j < m_dim; ++j)
                {
                    Vmath::Vvtvp (ntot, m_derivFac[i*m_dim+j], 1,
                                        Diff[j],               1,
                                        out[i],                1,
                                        out[i],                1);
                }
            }
        }

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

inlinevirtual

Implements Nektar::Collections::Operator.

Definition at line 280 of file PhysDeriv.cpp.

References Vmath::Vmul(), and Vmath::Vvtvp().

        {
            int nPhys = m_stdExp->GetTotPoints();
            int ntot = m_numElmt*nPhys;
            Array<OneD, NekDouble> tmp0,tmp1,tmp2;
            Array<OneD, Array<OneD, NekDouble> > Diff(3);

            for(int i = 0; i < m_dim; ++i)
            {
                Diff[i] = wsp + i*ntot;
            }

            // calculate local derivatives
            for (int i = 0; i < m_numElmt; ++i)
            {
                m_stdExp->PhysDeriv(input + i*nPhys,
                                    tmp0 = Diff[0] + i*nPhys,
                                    tmp1 = Diff[1] + i*nPhys,
                                    tmp2 = Diff[2] + i*nPhys);
            }

            // calculate full derivative
            Vmath::Vmul(ntot,m_derivFac[dir*m_dim],1,Diff[0],1,output,1);
            for(int j = 1; j < m_dim; ++j)
            {
                Vmath::Vvtvp (ntot, m_derivFac[dir*m_dim+j], 1,
                                    Diff[j],               1,
                                    output,                1,
                                    output,                1);
            }
        }

Member Data Documentation

int Nektar::Collections::PhysDeriv_IterPerExp::m_coordim

protected

Definition at line 319 of file PhysDeriv.cpp.

Array<TwoD, const NekDouble> Nektar::Collections::PhysDeriv_IterPerExp::m_derivFac

protected

Definition at line 317 of file PhysDeriv.cpp.

int Nektar::Collections::PhysDeriv_IterPerExp::m_dim

protected

Definition at line 318 of file PhysDeriv.cpp.