Phys deriv operator using sum-factorisation (Prism) More...

Inheritance diagram for Nektar::Collections::PhysDeriv_SumFac_Prism:

Collaboration diagram for Nektar::Collections::PhysDeriv_SumFac_Prism:

Public Member Functions
virtual	~PhysDeriv_SumFac_Prism ()

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

const int	m_nquad0

const int	m_nquad1

const int	m_nquad2

NekDouble *	m_Deriv0

NekDouble *	m_Deriv1

NekDouble *	m_Deriv2

Array< OneD, NekDouble >	m_fac0

Array< OneD, NekDouble >	m_fac1

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

unsigned int	m_numElmt

unsigned int	m_wspSize

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

Detailed Description

Phys deriv operator using sum-factorisation (Prism)

Definition at line 1338 of file PhysDeriv.cpp.

Constructor & Destructor Documentation

virtual Nektar::Collections::PhysDeriv_SumFac_Prism::~PhysDeriv_SumFac_Prism ( )

inlinevirtual

Definition at line 1343 of file PhysDeriv.cpp.

1344 {

1345 }

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

inlineprivate

Definition at line 1487 of file PhysDeriv.cpp.

             : Operator(pCollExp, pGeomData),
               m_nquad0  (m_stdExp->GetNumPoints(0)),
               m_nquad1  (m_stdExp->GetNumPoints(1)),
               m_nquad2  (m_stdExp->GetNumPoints(2))
         {
             LibUtilities::PointsKeyVector PtsKey = m_stdExp->GetPointsKeys();
 
             m_dim = PtsKey.size();
             m_coordim = m_stdExp->GetCoordim();
 
             m_derivFac = pGeomData->GetDerivFactors(pCollExp);
 
             const Array<OneD, const NekDouble>& z0
                                             = m_stdExp->GetBasis(0)->GetZ();
             const Array<OneD, const NekDouble>& z2
                                             = m_stdExp->GetBasis(2)->GetZ();
             m_fac0 = Array<OneD, NekDouble>(m_nquad0*m_nquad1*m_nquad2);
             m_fac1 = Array<OneD, NekDouble>(m_nquad0*m_nquad1*m_nquad2);
             for (int i = 0; i < m_nquad0; ++i)
             {
                 for(int j = 0; j < m_nquad1; ++j)
                 {
                     for(int k = 0; k < m_nquad2; ++k)
                     {
                         m_fac0[i+j*m_nquad0 + k*m_nquad0*m_nquad1] =
                             2.0/(1-z2[k]);
                         m_fac1[i+j*m_nquad0 + k*m_nquad0*m_nquad1] =
                             0.5*(1+z0[i]);
                     }
                 }
             }
 
 
 
             m_Deriv0 = &((m_stdExp->GetBasis(0)->GetD())->GetPtr())[0];
             m_Deriv1 = &((m_stdExp->GetBasis(1)->GetD())->GetPtr())[0];
             m_Deriv2 = &((m_stdExp->GetBasis(2)->GetD())->GetPtr())[0];
 
             m_wspSize = 3*m_nquad0*m_nquad1*m_nquad2*m_numElmt;
         }

Member Function Documentation

virtual void Nektar::Collections::PhysDeriv_SumFac_Prism::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 1347 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;
             }
 
             // dEta0
             Blas::Dgemm('N','N', m_nquad0,m_nquad1*m_nquad2*m_numElmt,
                         m_nquad0,1.0, m_Deriv0,m_nquad0,&input[0],
                         m_nquad0,0.0,&Diff[0][0],m_nquad0);
 
             int cnt = 0;
             for(int  i = 0; i < m_numElmt; ++i)
             {
 
                 // dEta 1
                 for (int j = 0; j < m_nquad2; ++j)
                 {
                     Blas::Dgemm('N', 'T', m_nquad0, m_nquad1, m_nquad1,
                                 1.0, &input[i*nPhys+j*m_nquad0*m_nquad1],
                                 m_nquad0, m_Deriv1, m_nquad1, 0.0,
                                 &Diff[1][i*nPhys+j*m_nquad0*m_nquad1],
                                 m_nquad0);
                 }
 
                 // dEta 2
                 Blas::Dgemm('N','T',m_nquad0*m_nquad1,m_nquad2,m_nquad2,
                             1.0, &input[i*nPhys],m_nquad0*m_nquad1,
                             m_Deriv2,m_nquad2, 0.0,&Diff[2][i*nPhys],
                             m_nquad0*m_nquad1);
 
                 // dxi0 = 2/(1-eta_2) d Eta_0
                 Vmath::Vmul(nPhys,&m_fac0[0],1,Diff[0].get()+cnt,1,
                             Diff[0].get()+cnt,1);
 
                 // dxi2 = (1+eta0)/(1-eta_2) d Eta_0 + d/dEta2;
                 Vmath::Vvtvp(nPhys,&m_fac1[0],1,Diff[0].get()+cnt,1,
                              Diff[2].get()+cnt,1,Diff[2].get()+cnt,1);
                 cnt += 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_SumFac_Prism::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 1413 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;
             }
 
             // dEta0
             Blas::Dgemm('N','N', m_nquad0,m_nquad1*m_nquad2*m_numElmt,
                         m_nquad0,1.0, m_Deriv0,m_nquad0,&input[0],
                         m_nquad0,0.0,&Diff[0][0],m_nquad0);
 
             int cnt = 0;
             for(int  i = 0; i < m_numElmt; ++i)
             {
 
                 // dEta 1
                 for (int j = 0; j < m_nquad2; ++j)
                 {
                     Blas::Dgemm('N', 'T', m_nquad0, m_nquad1, m_nquad1,
                                 1.0, &input[i*nPhys+j*m_nquad0*m_nquad1],
                                 m_nquad0, m_Deriv1, m_nquad1, 0.0,
                                 &Diff[1][i*nPhys+j*m_nquad0*m_nquad1],
                                 m_nquad0);
                 }
 
                 // dEta 2
                 Blas::Dgemm('N','T',m_nquad0*m_nquad1,m_nquad2,m_nquad2,
                             1.0, &input[i*nPhys],m_nquad0*m_nquad1,
                             m_Deriv2,m_nquad2, 0.0,&Diff[2][i*nPhys],
                             m_nquad0*m_nquad1);
 
                 // dxi0 = 2/(1-eta_2) d Eta_0
                 Vmath::Vmul(nPhys,&m_fac0[0],1,Diff[0].get()+cnt,1,
                             Diff[0].get()+cnt,1);
 
                 // dxi2 = (1+eta0)/(1-eta_2) d Eta_0 + d/dEta2;
                 Vmath::Vvtvp(nPhys,&m_fac1[0],1,Diff[0].get()+cnt,1,
                              Diff[2].get()+cnt,1,Diff[2].get()+cnt,1);
                 cnt += 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_SumFac_Prism::m_coordim

protected

Definition at line 1476 of file PhysDeriv.cpp.

NekDouble* Nektar::Collections::PhysDeriv_SumFac_Prism::m_Deriv0

protected

Definition at line 1480 of file PhysDeriv.cpp.

NekDouble* Nektar::Collections::PhysDeriv_SumFac_Prism::m_Deriv1

protected

Definition at line 1481 of file PhysDeriv.cpp.

NekDouble* Nektar::Collections::PhysDeriv_SumFac_Prism::m_Deriv2

protected

Definition at line 1482 of file PhysDeriv.cpp.

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

protected

Definition at line 1474 of file PhysDeriv.cpp.

int Nektar::Collections::PhysDeriv_SumFac_Prism::m_dim

protected

Definition at line 1475 of file PhysDeriv.cpp.

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

protected

Definition at line 1483 of file PhysDeriv.cpp.

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

protected

Definition at line 1484 of file PhysDeriv.cpp.

const int Nektar::Collections::PhysDeriv_SumFac_Prism::m_nquad0

protected

Definition at line 1477 of file PhysDeriv.cpp.

const int Nektar::Collections::PhysDeriv_SumFac_Prism::m_nquad1

protected

Definition at line 1478 of file PhysDeriv.cpp.

const int Nektar::Collections::PhysDeriv_SumFac_Prism::m_nquad2

protected

Definition at line 1479 of file PhysDeriv.cpp.

Public Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation