Nektar::Collections::IProductWRTDerivBase_SumFac_Tet Class Reference

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

Inheritance diagram for Nektar::Collections::IProductWRTDerivBase_SumFac_Tet:

Collaboration diagram for Nektar::Collections::IProductWRTDerivBase_SumFac_Tet:

Public Member Functions
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
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
Array< OneD, NekDouble >	m_fac0
Array< OneD, NekDouble >	m_fac1
Array< OneD, NekDouble >	m_fac2
Array< OneD, NekDouble >	m_fac3
bool	m_sortTopEdge
Protected Attributes inherited from Nektar::Collections::Operator
StdRegions::StdExpansionSharedPtr	m_stdExp
unsigned int	m_numElmt
unsigned int	m_wspSize

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

Detailed Description

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

Definition at line 1015 of file IProductWRTDerivBase.cpp.

Constructor & Destructor Documentation

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

inlineprivate

Definition at line 1174 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_nquad2  (m_stdExp->GetNumPoints(2)),
              m_nmodes0 (m_stdExp->GetBasisNumModes(0)),
              m_nmodes1 (m_stdExp->GetBasisNumModes(1)),
              m_nmodes2 (m_stdExp->GetBasisNumModes(2)),
              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);


            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]
                               = (1+z0[i])*0.5;
                        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])*0.5;
                    }
                }
            }

            if(m_stdExp->GetBasis(0)->GetBasisType()
                    == LibUtilities::eModified_A)
            {
                m_sortTopEdge = true;
            }
            else
            {
                m_sortTopEdge = false;
            }

        }

Member Function Documentation

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

References Nektar::Collections::TetIProduct(), 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 (dphi/dx,in[0]) = ((dphi/dxi_0 dxi_0/dx +
            //                               dphi/dxi_1 dxi_1/dx +
            //                               dphi/dxi_2 dxi_2/dx),in[0])
            //     +     (dphi/dy,in[1]) = ((dphi/dxi_0 dxi_0/dy +
            //                               dphi/dxi_1 dxi_1/dy +
            //                               dphi/dxi_2 dxi_2/dy),in[1])
            //     +     (dphi/dz,in[2]) = ((dphi/dxi_0 dxi_0/dz +
            //                               dphi/dxi_1 dxi_1/dz +
            //                               dphi/dxi_2 dxi_2/dz),in[1])
            //
            // Note dphi/dxi_0  =
            //             dphi/deta_0 4/((1-eta_1)(1-eta2))
            //
            //      dphi/dxi_1  =
            //             dphi/deta_0 2(1+eta_0)/((1-eta_1)(1-eta_2)) +
            //             dphi/deta_1 2/(1-eta_2)
            //
            //      dphi/dxi_2  =
            //             dphi/deta_0 2(1+eta_0)/((1-eta_1)(1-eta_2)) +
            //             dphi/deta_1 (1+eta_1)/(1-eta_2)  + dphi/deta_2
            //
            // and so the full inner products are
            //
            // (dphi/dx,in[0]) + (dphi/dy,in[1]) + (dphi/dz,in[2])
            //    = (dphi/deta_0, fac0 (tmp0 + fac1(tmp1 + tmp2)))
            //    + (dphi/deta_1, fac2 (tmp1 + fac3 tmp2))
            //    + (dphi/deta_2, tmp2)
            //
            // tmp0 = (dxi_0/dx in[0] + dxi_0/dy in[1] + dxi_0/dz in[2])
            // tmp1 = (dxi_1/dx in[0] + dxi_1/dy in[1] + dxi_1/dz in[2])
            // tmp2 = (dxi_2/dx in[0] + dxi_2/dy in[1] + dxi_2/dz in[2])

            // fac0 = 4/((1-eta_1)(1-eta2))
            // fac1 = (1+eta_0)/2
            // fac2 = 2/(1-eta_2)
            // fac3 = (1+eta_1)/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;

            // Sort into eta factors
            for (int i = 0; i < m_numElmt; ++i)
            {
                // add tmp[1] + tmp[2]
                Vmath::Vadd(nPhys, tmp[1].get() + i*nPhys, 1,
                                   tmp[2].get() + i*nPhys, 1,
                                   wsp1.get(),             1);

                // scale wsp1 by fac1 and add to tmp0
                Vmath::Vvtvp(nPhys,&m_fac1[0],1,wsp1.get(),1,
                             tmp[0].get()+i*nPhys,1,tmp[0].get()+i*nPhys,1);

                // scale tmp[0] by fac0
                Vmath::Vmul(nPhys,&m_fac0[0],1,tmp[0].get()+i*nPhys,1,
                            tmp[0].get()+i*nPhys,1);

                // scale tmp[2] by fac3 and add to tmp1
                Vmath::Vvtvp(nPhys,&m_fac3[0],1,tmp[2].get()+i*nPhys,1,
                             tmp[1].get()+i*nPhys,1,tmp[1].get()+i*nPhys,1);

                // scale tmp[1] by fac2
                Vmath::Vmul(nPhys,&m_fac2[0],1,tmp[1].get()+i*nPhys,1,
                            tmp[1].get()+i*nPhys,1);
            }


            // calculate Iproduct WRT Std Deriv
            TetIProduct(m_sortTopEdge, 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);

            TetIProduct(m_sortTopEdge, 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);

            TetIProduct(m_sortTopEdge, 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_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 1143 of file IProductWRTDerivBase.cpp.

References ASSERTL0.

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

Member Data Documentation

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

protected

Definition at line 1160 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1161 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1162 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1163 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1164 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1165 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1166 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1167 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1168 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1169 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1170 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1159 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Tet::m_nmodes0

protected

Definition at line 1156 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Tet::m_nmodes1

protected

Definition at line 1157 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Tet::m_nmodes2

protected

Definition at line 1158 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Tet::m_nquad0

protected

Definition at line 1153 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Tet::m_nquad1

protected

Definition at line 1154 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Tet::m_nquad2

protected

Definition at line 1155 of file IProductWRTDerivBase.cpp.

bool Nektar::Collections::IProductWRTDerivBase_SumFac_Tet::m_sortTopEdge

protected

Definition at line 1171 of file IProductWRTDerivBase.cpp.