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

Inheritance diagram for Nektar::Collections::IProductWRTDerivBase_SumFac_Tri:

Collaboration diagram for Nektar::Collections::IProductWRTDerivBase_SumFac_Tri:

Public Member Functions
virtual	~IProductWRTDerivBase_SumFac_Tri ()

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

Array< OneD, NekDouble >	m_fac0

Array< OneD, NekDouble >	m_fac1

bool	m_sortTopVertex

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

unsigned int	m_numElmt

unsigned int	m_wspSize

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

Detailed Description

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

Definition at line 691 of file IProductWRTDerivBase.cpp.

Constructor & Destructor Documentation

virtual Nektar::Collections::IProductWRTDerivBase_SumFac_Tri::~IProductWRTDerivBase_SumFac_Tri ( )

inlinevirtual

Definition at line 696 of file IProductWRTDerivBase.cpp.

697 {

698 }

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

inlineprivate

Definition at line 806 of file IProductWRTDerivBase.cpp.

References Nektar::LibUtilities::eModified_A.

             : 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));
 
             if(m_stdExp->GetBasis(0)->GetBasisType()
                     == LibUtilities::eModified_A)
             {
                 m_sortTopVertex = true;
             }
             else
             {
                 m_sortTopVertex = false;
             }
 
             const Array<OneD, const NekDouble>& z0
                                             = m_stdExp->GetBasis(0)->GetZ();
             const Array<OneD, const NekDouble>& z1
                                             = m_stdExp->GetBasis(1)->GetZ();
 
             m_fac0 = Array<OneD, NekDouble>(m_nquad0*m_nquad1);
             // set up geometric factor: 2/(1-z1)
             for (int i = 0; i < m_nquad0; ++i)
             {
                 for(int j = 0; j < m_nquad1; ++j)
                 {
                     m_fac0[i+j*m_nquad0] = 2.0/(1-z1[j]);
                 }
             }
 
             m_fac1 = Array<OneD, NekDouble>(m_nquad0*m_nquad1);
             // set up geometric factor: (1+z0)/(1-z1)
             for (int i = 0; i < m_nquad0; ++i)
             {
                 for(int j = 0; j < m_nquad1; ++j)
                 {
                     m_fac1[i+j*m_nquad0] = (1+z0[i])/(1-z1[j]);
                 }
             }
         }

Member Function Documentation

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

References Nektar::Collections::TriIProduct(), 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 (dphi/dx,in[0]) = ((dphi/dxi_0 dxi_0/dx +
             //                               dphi/dxi_1 dxi_1/dx),in[0])
             //     +     (dphi/dy,in[1]) = ((dphi/dxi_0 dxi_0/dy +
             //                               dphi/dxi_1 dxi_1/dy),in[1])
             //
             // Note dphi/dxi_0  =
             //             dphi/deta_0 deta_0/dxi_0 = dphi/deta_0 2/(1-eta_1)
             //
             //      dphi/dxi_1  =
             //             dphi/deta_1 deta_1/dxi_1 + dphi/deta_1 deta_1/dxi_1 =
             //             dphi/deta_0 (1+eta_0)/(1-eta_1) + dphi/deta_1
             //
             // and so the full inner products are
             //
             // (dphi/dx,in[0]) + (dphi/dy,in[1])
             //    = (dphi/deta_0, ((2/(1-eta_1) (dxi_0/dx in[0]+dxi_0/dy in[1])
             //            + (1_eta_0)/(1-eta_1) (dxi_1/dx in[0]+dxi_1/dy in[1]))
             //    + (dphi/deta_1, (dxi_1/dx in[0] + dxi_1/dy 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);
                 }
             }
 
             // Multiply by factor: 2/(1-z1)
             for (int i = 0; i < m_numElmt; ++i)
             {
                 // scale tmp[0] by geometric factor: 2/(1-z1)
                 Vmath::Vmul(nPhys,&m_fac0[0],1,tmp[0].get()+i*nPhys,1,
                             tmp[0].get()+i*nPhys,1);
 
                 // scale tmp[1] by geometric factor (1+z0)/(1-z1)
                 Vmath::Vvtvp(nPhys,&m_fac1[0],1,tmp[1].get()+i*nPhys,1,
                              tmp[0].get()+i*nPhys,1,tmp[0].get()+i*nPhys,1);
             }
 
             // Iproduct wrt derivative of base 0
             TriIProduct(m_sortTopVertex, 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
             TriIProduct(m_sortTopVertex, 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_Tri::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 778 of file IProductWRTDerivBase.cpp.

References ASSERTL0.

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

Member Data Documentation

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

protected

Definition at line 797 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 798 of file IProductWRTDerivBase.cpp.

const bool Nektar::Collections::IProductWRTDerivBase_SumFac_Tri::m_colldir0

protected

Definition at line 792 of file IProductWRTDerivBase.cpp.

const bool Nektar::Collections::IProductWRTDerivBase_SumFac_Tri::m_colldir1

protected

Definition at line 793 of file IProductWRTDerivBase.cpp.

int Nektar::Collections::IProductWRTDerivBase_SumFac_Tri::m_coordim

protected

Definition at line 794 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 799 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 800 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 795 of file IProductWRTDerivBase.cpp.

Array<OneD, NekDouble> Nektar::Collections::IProductWRTDerivBase_SumFac_Tri::m_fac0

protected

Definition at line 801 of file IProductWRTDerivBase.cpp.

Array<OneD, NekDouble> Nektar::Collections::IProductWRTDerivBase_SumFac_Tri::m_fac1

protected

Definition at line 802 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 796 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Tri::m_nmodes0

protected

Definition at line 790 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Tri::m_nmodes1

protected

Definition at line 791 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Tri::m_nquad0

protected

Definition at line 788 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Tri::m_nquad1

protected

Definition at line 789 of file IProductWRTDerivBase.cpp.

bool Nektar::Collections::IProductWRTDerivBase_SumFac_Tri::m_sortTopVertex

protected

Definition at line 803 of file IProductWRTDerivBase.cpp.

Public Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation