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

Inheritance diagram for Nektar::Collections::IProductWRTDerivBase_SumFac_Hex:

Collaboration diagram for Nektar::Collections::IProductWRTDerivBase_SumFac_Hex:

Public Member Functions
virtual	~IProductWRTDerivBase_SumFac_Hex ()

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_nquad2

const int	m_nmodes0

const int	m_nmodes1

const int	m_nmodes2

const bool	m_colldir0

const bool	m_colldir1

const bool	m_colldir2

Array< OneD, const NekDouble >	m_jac

Array< OneD, const NekDouble >	m_base0

Array< OneD, const NekDouble >	m_base1

Array< OneD, const NekDouble >	m_base2

Array< OneD, const NekDouble >	m_derbase0

Array< OneD, const NekDouble >	m_derbase1

Array< OneD, const NekDouble >	m_derbase2

Array< TwoD, const NekDouble >	m_derivFac

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

unsigned int	m_numElmt

unsigned int	m_wspSize

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

Detailed Description

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

Definition at line 877 of file IProductWRTDerivBase.cpp.

Constructor & Destructor Documentation

virtual Nektar::Collections::IProductWRTDerivBase_SumFac_Hex::~IProductWRTDerivBase_SumFac_Hex ( )

inlinevirtual

Definition at line 882 of file IProductWRTDerivBase.cpp.

883 {

884 }

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

inlineprivate

Definition at line 976 of file IProductWRTDerivBase.cpp.

             : Operator  (pCollExp, pGeomData),
               m_nquad0  (m_stdExp->GetNumPoints(0)),
               m_nquad1  (m_stdExp->GetNumPoints(1)),
               m_nquad2  (m_stdExp->GetNumPoints(2)),
               m_nmodes0 (m_stdExp->GetBasisNumModes(0)),
               m_nmodes1 (m_stdExp->GetBasisNumModes(1)),
               m_nmodes2 (m_stdExp->GetBasisNumModes(2)),
               m_colldir0(m_stdExp->GetBasis(0)->Collocation()),
               m_colldir1(m_stdExp->GetBasis(1)->Collocation()),
               m_colldir2(m_stdExp->GetBasis(2)->Collocation()),
               m_base0   (m_stdExp->GetBasis(0)->GetBdata()),
               m_base1   (m_stdExp->GetBasis(1)->GetBdata()),
               m_base2   (m_stdExp->GetBasis(2)->GetBdata()),
               m_derbase0(m_stdExp->GetBasis(0)->GetDbdata()),
               m_derbase1(m_stdExp->GetBasis(1)->GetDbdata()),
               m_derbase2(m_stdExp->GetBasis(2)->GetDbdata())
 
         {
             m_jac      = pGeomData->GetJacWithStdWeights(pCollExp);
             m_wspSize  = 6 * m_numElmt * (max(m_nquad0*m_nquad1*m_nquad2,
                                               m_nmodes0*m_nmodes1*m_nmodes2));
             m_derivFac = pGeomData->GetDerivFactors(pCollExp);
         }

Member Function Documentation

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

References Nektar::Collections::HexIProduct(), 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(3);
 
             in[0] = entry0; in[1] = entry1;
             in[2] = entry2;
 
             output =  entry3;
 
             for(int i = 0; i < 3; ++i)
             {
                 tmp[i] = wsp + i*nmax;
             }
 
             // calculate dx/dxi in[0] + dy/dxi in[1] + dz/dxi in[2]
             for(int i = 0; i < 3; ++i)
             {
                 Vmath::Vmul (ntot,m_derivFac[i],1, in[0],1,
                              tmp[i],1);
                 for(int j = 1; j < 3; ++j)
                 {
                     Vmath::Vvtvp (ntot,m_derivFac[i+3*j],1,
                                   in[j],1, tmp[i], 1, tmp[i],1);
                 }
             }
 
             wsp1   = wsp + 3*nmax;
 
             // calculate Iproduct WRT Std Deriv
             HexIProduct(false,m_colldir1,m_colldir2, m_numElmt,
                         m_nquad0,   m_nquad1,  m_nquad2,
                         m_nmodes0,  m_nmodes1, m_nmodes2,
                         m_derbase0, m_base1,   m_base2,
                         m_jac,tmp[0],output,wsp1);
 
             HexIProduct(m_colldir0,false,m_colldir2, m_numElmt,
                         m_nquad0,  m_nquad1,   m_nquad2,
                         m_nmodes0, m_nmodes1,  m_nmodes2,
                         m_base0,   m_derbase1, m_base2,
                         m_jac,tmp[1],tmp[0],wsp1);
             Vmath::Vadd(m_numElmt*nmodes,tmp[0],1,output,1,output,1);
 
             HexIProduct(m_colldir0,m_colldir1,false, m_numElmt,
                         m_nquad0,  m_nquad1,  m_nquad2,
                         m_nmodes0, m_nmodes1, m_nmodes2,
                         m_base0,   m_base1,   m_derbase2,
                         m_jac,tmp[2],tmp[0],wsp1);
             Vmath::Vadd(m_numElmt*nmodes,tmp[0],1,output,1,output,1);
         }

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

References ASSERTL0.

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

Member Data Documentation

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

protected

Definition at line 967 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 968 of file IProductWRTDerivBase.cpp.

Array<OneD, const NekDouble> Nektar::Collections::IProductWRTDerivBase_SumFac_Hex::m_base2

protected

Definition at line 969 of file IProductWRTDerivBase.cpp.

const bool Nektar::Collections::IProductWRTDerivBase_SumFac_Hex::m_colldir0

protected

Definition at line 963 of file IProductWRTDerivBase.cpp.

const bool Nektar::Collections::IProductWRTDerivBase_SumFac_Hex::m_colldir1

protected

Definition at line 964 of file IProductWRTDerivBase.cpp.

const bool Nektar::Collections::IProductWRTDerivBase_SumFac_Hex::m_colldir2

protected

Definition at line 965 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 970 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 971 of file IProductWRTDerivBase.cpp.

Array<OneD, const NekDouble> Nektar::Collections::IProductWRTDerivBase_SumFac_Hex::m_derbase2

protected

Definition at line 972 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 973 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 966 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Hex::m_nmodes0

protected

Definition at line 960 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Hex::m_nmodes1

protected

Definition at line 961 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Hex::m_nmodes2

protected

Definition at line 962 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Hex::m_nquad0

protected

Definition at line 957 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Hex::m_nquad1

protected

Definition at line 958 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Hex::m_nquad2

protected

Definition at line 959 of file IProductWRTDerivBase.cpp.

Public Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation