Nektar::Collections::PhysDeriv_SumFac_Tet Class Reference

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

Inheritance diagram for Nektar::Collections::PhysDeriv_SumFac_Tet:

Public Member Functions
virtual	~PhysDeriv_SumFac_Tet ()
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, std::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_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
Array< OneD, NekDouble >	m_fac2
Array< OneD, NekDouble >	m_fac3
Protected Attributes inherited from Nektar::Collections::Operator
StdRegions::StdExpansionSharedPtr	m_stdExp
unsigned int	m_numElmt
unsigned int	m_wspSize

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

Detailed Description

Phys deriv operator using sum-factorisation (Tet)

Definition at line 1078 of file PhysDeriv.cpp.

Constructor & Destructor Documentation

~PhysDeriv_SumFac_Tet()

virtual Nektar::Collections::PhysDeriv_SumFac_Tet::~PhysDeriv_SumFac_Tet ( )

inlinevirtual

Definition at line 1083 of file PhysDeriv.cpp.

        {
        }

PhysDeriv_SumFac_Tet()

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

inlineprivate

Definition at line 1272 of file PhysDeriv.cpp.

References Nektar::Collections::ePhysDeriv, Nektar::Collections::eSumFac, Nektar::LibUtilities::eTetrahedron, Nektar::Collections::GetOperatorFactory(), and Nektar::LibUtilities::NekFactory< tKey, tBase, tParam >::RegisterCreatorFunction().

            : 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_coordim = pCollExp[0]->GetCoordim();

            m_derivFac = pGeomData->GetDerivFactors(pCollExp);

            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;

            const Array<OneD, const NekDouble>& z0
                                            = m_stdExp->GetBasis(0)->GetZ();
            const Array<OneD, const NekDouble>& z1
                                            = m_stdExp->GetBasis(1)->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);
            m_fac2 = Array<OneD, NekDouble>(m_nquad0*m_nquad1*m_nquad2);
            m_fac3 = Array<OneD, NekDouble>(m_nquad0*m_nquad1*m_nquad2);
            // calculate 2.0/((1-eta_1)(1-eta_2))
            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]
                               = 4.0/((1-z1[j])*(1-z2[k]));
                        m_fac1[i + j*m_nquad0 + k*m_nquad0*m_nquad1]
                               = 2.0*(1+z0[i])/((1-z1[j])*(1-z2[k]));
                        m_fac2[i + j*m_nquad0 + k*m_nquad0*m_nquad1]
                               = 2.0/(1-z2[k]);
                        m_fac3[i + j*m_nquad0 + k*m_nquad0*m_nquad1]
                               = (1+z1[j])/(1-z2[k]);
                    }
                }
            }

        }

Member Function Documentation

operator()() [1/2]

virtual void Nektar::Collections::PhysDeriv_SumFac_Tet::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 1087 of file PhysDeriv.cpp.

References Blas::Dgemm(), 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 < 3; ++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);

            // dEta2
            for(int  i = 0; i < m_numElmt; ++i)
            {
                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);
            }

            for(int  i = 0; i < m_numElmt; ++i)
            {

                // dEta1
                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);
                }

                // dxi2 = (1 + eta_1)/(1 -eta_2)*dEta1 + dEta2
                Vmath::Vvtvp(nPhys, m_fac3.get(),            1,
                                    Diff[1].get() + i*nPhys, 1,
                                    Diff[2].get() + i*nPhys, 1,
                                    Diff[2].get() + i*nPhys, 1);

                // dxi1 =  2/(1 - eta_2) dEta1
                Vmath::Vmul(nPhys,  m_fac2.get(),            1,
                                    Diff[1].get() + i*nPhys, 1,
                                    Diff[1].get() + i*nPhys, 1);

                // dxi1 = 2.0(1+eta_0)/((1-eta_1)(1-eta_2)) dEta0 + dxi1
                Vmath::Vvtvp(nPhys, m_fac1.get(),            1,
                                    Diff[0].get() + i*nPhys, 1,
                                    Diff[1].get() + i*nPhys, 1,
                                    Diff[1].get() + i*nPhys, 1);

                // dxi2 = 2.0(1+eta_0)/((1-eta_1)(1-eta_2)) dEta0 + dxi2
                Vmath::Vvtvp(nPhys, m_fac1.get(),            1,
                                    Diff[0].get() + i*nPhys, 1,
                                    Diff[2].get() + i*nPhys, 1,
                                    Diff[2].get() + i*nPhys, 1);

                // dxi0 = 4.0/((1-eta_1)(1-eta_2)) dEta0
                Vmath::Vmul(nPhys,  m_fac0.get(),            1,
                                    Diff[0].get() + i*nPhys, 1,
                                    Diff[0].get() + i*nPhys, 1);

            }

            // calculate full derivative
            for(int i = 0; i < m_coordim; ++i)
            {
                Vmath::Vmul(ntot,m_derivFac[i*3],1,Diff[0],1,out[i],1);
                for(int j = 1; j < 3; ++j)
                {
                    Vmath::Vvtvp (ntot, m_derivFac[i*3+j], 1,
                                        Diff[j], 1, out[i], 1, out[i], 1);
                }
            }
        }

operator()() [2/2]

virtual void Nektar::Collections::PhysDeriv_SumFac_Tet::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 1175 of file PhysDeriv.cpp.

References Blas::Dgemm(), 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 < 3; ++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);

            // dEta2
            for(int  i = 0; i < m_numElmt; ++i)
            {
                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);
            }

            for(int  i = 0; i < m_numElmt; ++i)
            {

                // dEta1
                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);
                }

                // dxi2 = (1 + eta_1)/(1 -eta_2)*dEta1 + dEta2
                Vmath::Vvtvp(nPhys, m_fac3.get(),            1,
                                    Diff[1].get() + i*nPhys, 1,
                                    Diff[2].get() + i*nPhys, 1,
                                    Diff[2].get() + i*nPhys, 1);

                // dxi1 =  2/(1 - eta_2) dEta1
                Vmath::Vmul(nPhys,  m_fac2.get(),            1,
                                    Diff[1].get() + i*nPhys, 1,
                                    Diff[1].get() + i*nPhys, 1);

                // dxi1 = 2.0(1+eta_0)/((1-eta_1)(1-eta_2)) dEta0 + dxi1
                Vmath::Vvtvp(nPhys, m_fac1.get(),            1,
                                    Diff[0].get() + i*nPhys, 1,
                                    Diff[1].get() + i*nPhys, 1,
                                    Diff[1].get() + i*nPhys, 1);

                // dxi2 = 2.0(1+eta_0)/((1-eta_1)(1-eta_2)) dEta0 + dxi2
                Vmath::Vvtvp(nPhys, m_fac1.get(),            1,
                                    Diff[0].get() + i*nPhys, 1,
                                    Diff[2].get() + i*nPhys, 1,
                                    Diff[2].get() + i*nPhys, 1);

                // dxi0 = 4.0/((1-eta_1)(1-eta_2)) dEta0
                Vmath::Vmul(nPhys,  m_fac0.get(),            1,
                                    Diff[0].get() + i*nPhys, 1,
                                    Diff[0].get() + i*nPhys, 1);

            }

            // calculate full derivative
            Vmath::Vmul(ntot,m_derivFac[dir*3],1,Diff[0],1,output,1);
            for(int j = 1; j < 3; ++j)
            {
                Vmath::Vvtvp (ntot, m_derivFac[dir*3+j], 1,
                                    Diff[j], 1, output, 1, output, 1);
            }
        }

Member Data Documentation

m_coordim

int Nektar::Collections::PhysDeriv_SumFac_Tet::m_coordim

protected

Definition at line 1259 of file PhysDeriv.cpp.

m_Deriv0

NekDouble* Nektar::Collections::PhysDeriv_SumFac_Tet::m_Deriv0

protected

Definition at line 1263 of file PhysDeriv.cpp.

m_Deriv1

NekDouble* Nektar::Collections::PhysDeriv_SumFac_Tet::m_Deriv1

protected

Definition at line 1264 of file PhysDeriv.cpp.

m_Deriv2

NekDouble* Nektar::Collections::PhysDeriv_SumFac_Tet::m_Deriv2

protected

Definition at line 1265 of file PhysDeriv.cpp.

m_derivFac

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

protected

Definition at line 1258 of file PhysDeriv.cpp.

m_fac0

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

protected

Definition at line 1266 of file PhysDeriv.cpp.

m_fac1

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

protected

Definition at line 1267 of file PhysDeriv.cpp.

m_fac2

Array<OneD, NekDouble> Nektar::Collections::PhysDeriv_SumFac_Tet::m_fac2

protected

Definition at line 1268 of file PhysDeriv.cpp.

m_fac3

Array<OneD, NekDouble> Nektar::Collections::PhysDeriv_SumFac_Tet::m_fac3

protected

Definition at line 1269 of file PhysDeriv.cpp.

m_nquad0

const int Nektar::Collections::PhysDeriv_SumFac_Tet::m_nquad0

protected

Definition at line 1260 of file PhysDeriv.cpp.

m_nquad1

const int Nektar::Collections::PhysDeriv_SumFac_Tet::m_nquad1

protected

Definition at line 1261 of file PhysDeriv.cpp.

m_nquad2

const int Nektar::Collections::PhysDeriv_SumFac_Tet::m_nquad2

protected

Definition at line 1262 of file PhysDeriv.cpp.