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

Inheritance diagram for Nektar::Collections::IProductWRTDerivBase_SumFac_Tri:

Public Member Functions
	~IProductWRTDerivBase_SumFac_Tri () final=default

void	operator() (const Array< OneD, const NekDouble > &entry0, Array< OneD, NekDouble > &entry1, Array< OneD, NekDouble > &entry2, Array< OneD, NekDouble > &entry3, Array< OneD, NekDouble > &wsp) final

void	operator() (int dir, const Array< OneD, const NekDouble > &input, Array< OneD, NekDouble > &output, Array< OneD, NekDouble > &wsp) final

Public Member Functions inherited from Nektar::Collections::Operator
	Operator (std::vector< LocalRegions::ExpansionSharedPtr > pCollExp, std::shared_ptr< CoalescedGeomData > GeomData, StdRegions::FactorMap factors)
	Constructor.

virtual	~Operator ()=default

virtual COLLECTIONS_EXPORT void	UpdateFactors (StdRegions::FactorMap factors)
	Update the supplied factor map.

virtual COLLECTIONS_EXPORT void	UpdateVarcoeffs (StdRegions::VarCoeffMap &varcoeffs)
	Update the supplied variable coefficients.

unsigned int	GetWspSize ()
	Get the size of the required workspace.

unsigned int	GetNumElmt ()
	Get number of elements.

StdRegions::StdExpansionSharedPtr	GetExpSharedPtr ()
	Get expansion pointer.

unsigned int	GetInputSize (void)

unsigned int	GetOutputSize (void)

unsigned int	GetPhysSize (void)

unsigned int	GetCoeffSize (void)

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_jacWStdW

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
bool	m_isDeformed

StdRegions::StdExpansionSharedPtr	m_stdExp

unsigned int	m_numElmt
	number of elements that the operator is applied on

unsigned int	m_nqe

unsigned int	m_wspSize

unsigned int	m_inputSize
	number of modes or quadrature points that are passed as input to an operator

unsigned int	m_outputSize
	number of modes or quadrature points that are taken as output from an operator

Private Member Functions
	IProductWRTDerivBase_SumFac_Tri (vector< LocalRegions::ExpansionSharedPtr > pCollExp, CoalescedGeomDataSharedPtr pGeomData, StdRegions::FactorMap factors)

Private Member Functions inherited from Nektar::Collections::IProductWRTDerivBase_Helper
	IProductWRTDerivBase_Helper ()

Detailed Description

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

Definition at line 948 of file IProductWRTDerivBase.cpp.

Constructor & Destructor Documentation

◆ ~IProductWRTDerivBase_SumFac_Tri()

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

finaldefault

◆ IProductWRTDerivBase_SumFac_Tri()

Nektar::Collections::IProductWRTDerivBase_SumFac_Tri::IProductWRTDerivBase_SumFac_Tri	(	vector< LocalRegions::ExpansionSharedPtr >	pCollExp,
		CoalescedGeomDataSharedPtr	pGeomData,
		StdRegions::FactorMap	factors
	)

inlineprivate

Definition at line 1096 of file IProductWRTDerivBase.cpp.

        : Operator(pCollExp, pGeomData, factors), IProductWRTDerivBase_Helper(),
          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())
    {
        m_coordim  = pCollExp[0]->GetCoordim();
        m_derivFac = pGeomData->GetDerivFactors(pCollExp);
        m_jacWStdW = 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]);
            }
        }
    }

References Nektar::LibUtilities::eModified_A, m_coordim, m_derivFac, m_fac0, m_fac1, m_jacWStdW, m_nmodes0, m_nmodes1, m_nquad0, m_nquad1, Nektar::Collections::Operator::m_numElmt, m_sortTopVertex, Nektar::Collections::Operator::m_stdExp, Nektar::Collections::Operator::m_wspSize, and tinysimd::max().

Member Function Documentation

◆ operator()() [1/2]

void Nektar::Collections::IProductWRTDerivBase_SumFac_Tri::operator()	(	const Array< OneD, const NekDouble > &	entry0,
		Array< OneD, NekDouble > &	entry1,
		Array< OneD, NekDouble > &	entry2,
		Array< OneD, NekDouble > &	entry3,
		Array< OneD, NekDouble > &	wsp
	)

inlinefinalvirtual

This method calculates:

\[ (d\phi/dx,in[0]) + (d\phi/dy,in[1]) \]

which can be represented in terms of local cartesian derivaties as:

\[ ((d\phi/d\xi_0\, d\xi_0/dx + d\phi/d\xi_1\, d\xi_1/dx),in[0]) + \]

\[ ((d\phi/d\xi_0\, d\xi_0/dy + d\phi/d\xi_1\, d\xi_1/dy),in[1]) + \]

where we note that

\[ d\phi/d\xi_0 = d\phi/d\eta_0\, d\eta_0/d\xi_0 = d\phi/d\eta_0 2/(1-\eta_1) \]

\[ d\phi/d\xi_1 = d\phi/d\eta_1\, d\eta_1/d\xi_1 + d\phi/d\eta_1\, d\eta_1/d\xi_1 = d\phi/d\eta_0 (1+\eta_0)/(1-\eta_1) + d\phi/d\eta_1 \]

and so the full inner products are

\[ (d\phi/dx,in[0]) + (dphi/dy,in[1]) = (d\phi/d\eta_0, ((2/(1-\eta_1) (d\xi_0/dx in[0] + d\xi_0/dy in[1]) + (1-\eta_0)/(1-\eta_1) (d\xi_1/dx in[0]+d\xi_1/dy in[1])) + (d\phi/d\eta_1, (d\xi_1/dx in[0] + d\xi_1/dy in[1])) \]

Implements Nektar::Collections::Operator.

Definition at line 987 of file IProductWRTDerivBase.cpp.

    {
        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;
 
        if (m_isDeformed)
        {
            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 + 2 * j], 1, in[j], 1,
                                 tmp[i], 1, tmp[i], 1);
                }
            }
        }
        else
        {
            Array<OneD, NekDouble> t;
            for (int e = 0; e < m_numElmt; ++e)
            {
                // calculate dx/dxi in[0] + dy/dxi in[1]
                for (int i = 0; i < 2; ++i)
                {
                    Vmath::Smul(m_nqe, m_derivFac[i][e], in[0] + e * m_nqe, 1,
                                t = tmp[i] + e * m_nqe, 1);
                    for (int j = 1; j < m_coordim; ++j)
                    {
                        Vmath::Svtvp(m_nqe, m_derivFac[i + 2 * j][e],
                                     in[j] + e * m_nqe, 1, tmp[i] + e * m_nqe,
                                     1, t = tmp[i] + e * m_nqe, 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].data() + i * nPhys, 1,
                        tmp[0].data() + i * nPhys, 1);
 
            // scale tmp[1] by geometric factor (1+z0)/(1-z1)
            Vmath::Vvtvp(nPhys, &m_fac1[0], 1, tmp[1].data() + i * nPhys, 1,
                         tmp[0].data() + i * nPhys, 1,
                         tmp[0].data() + 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_jacWStdW, 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_jacWStdW, tmp[1], tmp[0],
                    wsp1);
 
        Vmath::Vadd(m_numElmt * nmodes, tmp[0], 1, output, 1, output, 1);
    }

References m_base0, m_base1, m_coordim, m_derbase0, m_derbase1, m_derivFac, m_fac0, m_fac1, Nektar::Collections::Operator::m_isDeformed, m_jacWStdW, m_nmodes0, m_nmodes1, Nektar::Collections::Operator::m_nqe, m_nquad0, m_nquad1, Nektar::Collections::Operator::m_numElmt, m_sortTopVertex, Nektar::Collections::Operator::m_stdExp, tinysimd::max(), Vmath::Smul(), Vmath::Svtvp(), Nektar::Collections::TriIProduct(), Vmath::Vadd(), Vmath::Vmul(), and Vmath::Vvtvp().

◆ operator()() [2/2]

void Nektar::Collections::IProductWRTDerivBase_SumFac_Tri::operator()	(	int	dir,
		const Array< OneD, const NekDouble > &	input,
		Array< OneD, NekDouble > &	output,
		Array< OneD, NekDouble > &	wsp
	)

inlinefinalvirtual

Implements Nektar::Collections::Operator.

Definition at line 1069 of file IProductWRTDerivBase.cpp.

    {
        NEKERROR(ErrorUtil::efatal, "Not valid for this operator.");
    }

References Nektar::ErrorUtil::efatal, and NEKERROR.

Member Data Documentation

◆ m_base0

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

protected

Definition at line 1087 of file IProductWRTDerivBase.cpp.

Referenced by operator()().

◆ m_base1

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

protected

Definition at line 1088 of file IProductWRTDerivBase.cpp.

Referenced by operator()().

◆ m_colldir0

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

protected

Definition at line 1082 of file IProductWRTDerivBase.cpp.

◆ m_colldir1

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

protected

Definition at line 1083 of file IProductWRTDerivBase.cpp.

◆ m_coordim

int Nektar::Collections::IProductWRTDerivBase_SumFac_Tri::m_coordim

protected

Definition at line 1084 of file IProductWRTDerivBase.cpp.

Referenced by IProductWRTDerivBase_SumFac_Tri(), and operator()().

◆ m_derbase0

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

protected

Definition at line 1089 of file IProductWRTDerivBase.cpp.

Referenced by operator()().

◆ m_derbase1

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

protected

Definition at line 1090 of file IProductWRTDerivBase.cpp.

Referenced by operator()().

◆ m_derivFac

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

protected

Definition at line 1085 of file IProductWRTDerivBase.cpp.

Referenced by IProductWRTDerivBase_SumFac_Tri(), and operator()().

◆ m_fac0

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

protected

Definition at line 1091 of file IProductWRTDerivBase.cpp.

Referenced by IProductWRTDerivBase_SumFac_Tri(), and operator()().

◆ m_fac1

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

protected

Definition at line 1092 of file IProductWRTDerivBase.cpp.

Referenced by IProductWRTDerivBase_SumFac_Tri(), and operator()().

◆ m_jacWStdW

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

protected

Definition at line 1086 of file IProductWRTDerivBase.cpp.

Referenced by IProductWRTDerivBase_SumFac_Tri(), and operator()().

◆ m_nmodes0

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

protected

Definition at line 1080 of file IProductWRTDerivBase.cpp.

Referenced by IProductWRTDerivBase_SumFac_Tri(), and operator()().

◆ m_nmodes1

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

protected

Definition at line 1081 of file IProductWRTDerivBase.cpp.

Referenced by IProductWRTDerivBase_SumFac_Tri(), and operator()().

◆ m_nquad0

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

protected

Definition at line 1078 of file IProductWRTDerivBase.cpp.

Referenced by IProductWRTDerivBase_SumFac_Tri(), and operator()().

◆ m_nquad1

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

protected

Definition at line 1079 of file IProductWRTDerivBase.cpp.

Referenced by IProductWRTDerivBase_SumFac_Tri(), and operator()().

◆ m_sortTopVertex

bool Nektar::Collections::IProductWRTDerivBase_SumFac_Tri::m_sortTopVertex

protected

Definition at line 1093 of file IProductWRTDerivBase.cpp.

Referenced by IProductWRTDerivBase_SumFac_Tri(), and operator()().

Public Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ ~IProductWRTDerivBase_SumFac_Tri()

◆ IProductWRTDerivBase_SumFac_Tri()

Member Function Documentation

◆ operator()() [1/2]

◆ operator()() [2/2]

Member Data Documentation

◆ m_base0

◆ m_base1

◆ m_colldir0

◆ m_colldir1

◆ m_coordim

◆ m_derbase0

◆ m_derbase1

◆ m_derivFac

◆ m_fac0

◆ m_fac1

◆ m_jacWStdW

◆ m_nmodes0

◆ m_nmodes1

◆ m_nquad0

◆ m_nquad1

◆ m_sortTopVertex