Nektar::Collections::IProductWRTDerivBase_IterPerExp Class Reference

Inner product WRT deriv base operator using element-wise operation. More...

Inheritance diagram for Nektar::Collections::IProductWRTDerivBase_IterPerExp:

Collaboration diagram for Nektar::Collections::IProductWRTDerivBase_IterPerExp:

Public Member Functions
virtual	~IProductWRTDerivBase_IterPerExp ()
virtual void	operator() (const Array< OneD, const NekDouble > &entry0, Array< OneD, NekDouble > &entry1, Array< OneD, NekDouble > &entry2, Array< OneD, NekDouble > &entry3, 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
Array< OneD, const NekDouble >	m_jac
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
	IProductWRTDerivBase_IterPerExp (vector< StdRegions::StdExpansionSharedPtr > pCollExp, CoalescedGeomDataSharedPtr pGeomData)

Detailed Description

Inner product WRT deriv base operator using element-wise operation.

Definition at line 227 of file IProductWRTDerivBase.cpp.

Constructor & Destructor Documentation

virtual Nektar::Collections::IProductWRTDerivBase_IterPerExp::~IProductWRTDerivBase_IterPerExp ( )

inlinevirtual

Definition at line 232 of file IProductWRTDerivBase.cpp.

        {
        }

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

inlineprivate

Definition at line 312 of file IProductWRTDerivBase.cpp.

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

            int nqtot  = m_stdExp->GetTotPoints();

            m_derivFac = pGeomData->GetDerivFactors(pCollExp);
            m_jac      = pGeomData->GetJac(pCollExp);
            m_wspSize  = m_dim*nqtot*m_numElmt;
        }

Member Function Documentation

virtual void Nektar::Collections::IProductWRTDerivBase_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 236 of file IProductWRTDerivBase.cpp.

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

        {
            unsigned int nPhys  = m_stdExp->GetTotPoints();
            unsigned int ntot   = m_numElmt*nPhys;
            unsigned int nmodes = m_stdExp->GetNcoeffs();
            unsigned int nmax   = max(ntot,m_numElmt*nmodes);
            Array<OneD, Array<OneD, const NekDouble> > in(3);
            Array<OneD, NekDouble> output, tmp1;
            Array<OneD, Array<OneD, NekDouble> > tmp(3);

            in[0] = entry0; in[1] = entry1; in[2] = entry2;

            output = (m_coordim == 3)? entry3: (m_coordim == 2)?
                entry2: entry1;

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

            // calculate dx/dxi in[0] + dy/dxi in[2] + dz/dxi in[3]
            for(int i = 0; i < m_dim; ++i)
            {
                Vmath::Vmul (ntot,m_derivFac[i],1, in[0],1,
                             tmp[i],1);
                for(int j = 1; j < m_coordim; ++j)
                {
                    Vmath::Vvtvp (ntot,m_derivFac[i +j*m_dim],1,
                                  in[j],1, tmp[i], 1, tmp[i],1);
                }
            }

            // calculate Iproduct WRT Std Deriv
            // first component
            Vmath::Vmul(ntot,m_jac,1,tmp[0],1,tmp[0],1);
            for(int n = 0; n < m_numElmt; ++n)
            {
                m_stdExp->IProductWRTDerivBase(0,tmp[0]+n*nPhys,
                                               tmp1 = output + n*nmodes);
            }

            // other components
            for(int i = 1; i < m_dim; ++i)
            {
                // multiply by Jacobian
                Vmath::Vmul(ntot,m_jac,1,tmp[i],1,tmp[i],1);
                for(int n = 0; n < m_numElmt; ++n)
                {
                    m_stdExp->IProductWRTDerivBase(i,tmp[i]+n*nPhys,tmp[0]);
                    Vmath::Vadd(nmodes,tmp[0],1,output+n*nmodes,1,
                                tmp1 = output+n*nmodes,1);
                }
            }
        }

virtual void Nektar::Collections::IProductWRTDerivBase_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 296 of file IProductWRTDerivBase.cpp.

References ASSERTL0.

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

Member Data Documentation

int Nektar::Collections::IProductWRTDerivBase_IterPerExp::m_coordim

protected

Definition at line 309 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 306 of file IProductWRTDerivBase.cpp.

int Nektar::Collections::IProductWRTDerivBase_IterPerExp::m_dim

protected

Definition at line 308 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 307 of file IProductWRTDerivBase.cpp.