Inner product WRT deriv base operator using sum-factorisation (Quad) More...

Inheritance diagram for Nektar::Collections::IProductWRTDerivBase_SumFac_Quad:

Collaboration diagram for Nektar::Collections::IProductWRTDerivBase_SumFac_Quad:

Public Member Functions
virtual	~IProductWRTDerivBase_SumFac_Quad ()

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
const int	m_nquad0

const int	m_nquad1

const int	m_nmodes0

const int	m_nmodes1

const bool	m_colldir0

const bool	m_colldir1

int	m_coordim

Array< TwoD, const NekDouble >	m_derivFac

Array< OneD, const NekDouble >	m_jac

Array< OneD, const NekDouble >	m_base0

Array< OneD, const NekDouble >	m_base1

Array< OneD, const NekDouble >	m_derbase0

Array< OneD, const NekDouble >	m_derbase1

Protected Attributes inherited from Nektar::Collections::Operator
StdRegions::StdExpansionSharedPtr	m_stdExp

unsigned int	m_numElmt

unsigned int	m_wspSize

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

Detailed Description

Inner product WRT deriv base operator using sum-factorisation (Quad)

Definition at line 570 of file IProductWRTDerivBase.cpp.

Constructor & Destructor Documentation

virtual Nektar::Collections::IProductWRTDerivBase_SumFac_Quad::~IProductWRTDerivBase_SumFac_Quad ( )

inlinevirtual

Definition at line 575 of file IProductWRTDerivBase.cpp.

576 {

577 }

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

inlineprivate

Definition at line 655 of file IProductWRTDerivBase.cpp.

             : Operator(pCollExp, pGeomData),
               m_nquad0  (m_stdExp->GetNumPoints(0)),
               m_nquad1  (m_stdExp->GetNumPoints(1)),
               m_nmodes0 (m_stdExp->GetBasisNumModes(0)),
               m_nmodes1 (m_stdExp->GetBasisNumModes(1)),
               m_colldir0(m_stdExp->GetBasis(0)->Collocation()),
               m_colldir1(m_stdExp->GetBasis(1)->Collocation()),
               m_base0   (m_stdExp->GetBasis(0)->GetBdata()),
               m_base1   (m_stdExp->GetBasis(1)->GetBdata()),
               m_derbase0(m_stdExp->GetBasis(0)->GetDbdata()),
               m_derbase1(m_stdExp->GetBasis(1)->GetDbdata())
         {
             LibUtilities::PointsKeyVector PtsKey = m_stdExp->GetPointsKeys();
             m_coordim = m_stdExp->GetCoordim();
 
             m_derivFac = pGeomData->GetDerivFactors(pCollExp);
             m_jac      = pGeomData->GetJacWithStdWeights(pCollExp);
             m_wspSize  = 4 * m_numElmt * (max(m_nquad0*m_nquad1,
                                               m_nmodes0*m_nmodes1));
         }

Member Function Documentation

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

References Nektar::Collections::QuadIProduct(), 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, wsp1;
             Array<OneD, Array<OneD, NekDouble> > tmp(2);
 
             in[0] = entry0; in[1] = entry1; in[2] = entry2;
 
             output = (m_coordim == 2)? entry2: entry3;
 
             tmp[0] = wsp; tmp[1] = wsp + nmax;
             wsp1   = wsp + 2*nmax;
 
             // calculate dx/dxi in[0] + dy/dxi in[1]
             for(int i = 0; i < 2; ++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*2],1,
                                   in[j],1, tmp[i], 1, tmp[i],1);
                 }
             }
 
             // Iproduct wrt derivative of base 0
             QuadIProduct(false, m_colldir1,m_numElmt,
                          m_nquad0,   m_nquad1,
                          m_nmodes0,  m_nmodes1,
                          m_derbase0, m_base1,
                          m_jac, tmp[0], output, wsp1);
 
             // Iproduct wrt derivative of base 1
             QuadIProduct(m_colldir0, false, m_numElmt,
                          m_nquad0,   m_nquad1,
                          m_nmodes0,  m_nmodes1,
                          m_base0, m_derbase1,
                          m_jac, tmp[1],  tmp[0], wsp1);
 
             Vmath::Vadd(m_numElmt*nmodes,tmp[0],1,output,1,output,1);
         }

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

References ASSERTL0.

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

Member Data Documentation

Array<OneD, const NekDouble> Nektar::Collections::IProductWRTDerivBase_SumFac_Quad::m_base0

protected

Definition at line 649 of file IProductWRTDerivBase.cpp.

Array<OneD, const NekDouble> Nektar::Collections::IProductWRTDerivBase_SumFac_Quad::m_base1

protected

Definition at line 650 of file IProductWRTDerivBase.cpp.

const bool Nektar::Collections::IProductWRTDerivBase_SumFac_Quad::m_colldir0

protected

Definition at line 644 of file IProductWRTDerivBase.cpp.

const bool Nektar::Collections::IProductWRTDerivBase_SumFac_Quad::m_colldir1

protected

Definition at line 645 of file IProductWRTDerivBase.cpp.

int Nektar::Collections::IProductWRTDerivBase_SumFac_Quad::m_coordim

protected

Definition at line 646 of file IProductWRTDerivBase.cpp.

Array<OneD, const NekDouble> Nektar::Collections::IProductWRTDerivBase_SumFac_Quad::m_derbase0

protected

Definition at line 651 of file IProductWRTDerivBase.cpp.

Array<OneD, const NekDouble> Nektar::Collections::IProductWRTDerivBase_SumFac_Quad::m_derbase1

protected

Definition at line 652 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 647 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 648 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Quad::m_nmodes0

protected

Definition at line 642 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Quad::m_nmodes1

protected

Definition at line 643 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Quad::m_nquad0

protected

Definition at line 640 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Quad::m_nquad1

protected

Definition at line 641 of file IProductWRTDerivBase.cpp.

Public Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation