Nektar::Collections::IProductWRTDerivBase_SumFac_Prism Class Reference

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

Inheritance diagram for Nektar::Collections::IProductWRTDerivBase_SumFac_Prism:

Collaboration diagram for Nektar::Collections::IProductWRTDerivBase_SumFac_Prism:

Public Member Functions
virtual	~IProductWRTDerivBase_SumFac_Prism ()
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
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_Prism (vector< StdRegions::StdExpansionSharedPtr > pCollExp, CoalescedGeomDataSharedPtr pGeomData)

Detailed Description

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

Definition at line 1253 of file IProductWRTDerivBase.cpp.

Constructor & Destructor Documentation

virtual Nektar::Collections::IProductWRTDerivBase_SumFac_Prism::~IProductWRTDerivBase_SumFac_Prism ( )

inlinevirtual

Definition at line 1258 of file IProductWRTDerivBase.cpp.

        {
        }

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

inlineprivate

Definition at line 1389 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);

            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>& 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)
                    {
                        // set up geometric factor: 2/(1-z1)
                        m_fac0[i + j*m_nquad0 + k*m_nquad0*m_nquad1]
                               = 2.0/(1-z2[k]);
                        // set up geometric factor: (1+z0)/(1-z1)
                        m_fac1[i + j*m_nquad0 + k*m_nquad0*m_nquad1]
                               = (1+z0[i])/(1-z2[k]);

                    }
                }
            }
        }

Member Function Documentation

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

References Nektar::Collections::PrismIProduct(), 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),in[0])
            //     +     (dphi/dy,in[1]) = ((dphi/dxi_0 dxi_0/dy +
            //                               dphi/dxi_1 dxi_1/dy),in[1])
            //     +     (dphi/dz,in[2]) = ((dphi/dxi_0 dxi_0/dz +
            //                               dphi/dxi_1 dxi_1/dz),in[2])
            //
            // Note dphi/dxi_0  =
            //             dphi/deta_0 deta_0/dxi_0 = dphi/deta_0 2/(1-eta_2)
            //
            //      dphi/dxi_2  =
            //             dphi/deta_0 deta_0/dxi_2 + dphi/deta_2 deta_2/dxi_2 =
            //             dphi/deta_0 (1+eta_0)/(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, ((2/(1-eta_2) (dxi_0/dx in[0] + dxi_0/dy in[1]
            //                                               + dxi_0/dz in[2])
            //            + (1_eta_0)/(1-eta_2) (dxi_2/dx in[0] + dxi_2/dy in[1]
            //                                               + dxi_2/dz in[2] ))
            //    + (dphi/deta_1, (dxi_1/dx in[0] + dxi_1/dy in[1]
            //                                    + dxi_1/dz in[2]))
            //    + (dphi/deta_2, (dxi_2/dx in[0] + dxi_2/dy in[1]
            //                                    + dxi_2/dz 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;

            // Sort into eta factors
            for (int i = 0; i < m_numElmt; ++i)
            {
                // 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 fac1 and add to tmp0
                Vmath::Vvtvp(nPhys,&m_fac1[0],1,tmp[2].get()+i*nPhys,1,
                             tmp[0].get()+i*nPhys,1,tmp[0].get()+i*nPhys,1);
            }

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

            PrismIProduct(m_sortTopVertex, 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);

            PrismIProduct(m_sortTopVertex, 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_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 1359 of file IProductWRTDerivBase.cpp.

References ASSERTL0.

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

Member Data Documentation

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

protected

Definition at line 1377 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1378 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1379 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1380 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1381 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1382 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1383 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1384 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1385 of file IProductWRTDerivBase.cpp.

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

protected

Definition at line 1376 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Prism::m_nmodes0

protected

Definition at line 1373 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Prism::m_nmodes1

protected

Definition at line 1374 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Prism::m_nmodes2

protected

Definition at line 1375 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Prism::m_nquad0

protected

Definition at line 1370 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Prism::m_nquad1

protected

Definition at line 1371 of file IProductWRTDerivBase.cpp.

const int Nektar::Collections::IProductWRTDerivBase_SumFac_Prism::m_nquad2

protected

Definition at line 1372 of file IProductWRTDerivBase.cpp.

bool Nektar::Collections::IProductWRTDerivBase_SumFac_Prism::m_sortTopVertex

protected

Definition at line 1386 of file IProductWRTDerivBase.cpp.