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

Inheritance diagram for Nektar::Collections::IProductWRTDerivBase_SumFac_Prism:

Public Member Functions
	~IProductWRTDerivBase_SumFac_Prism () 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< StdRegions::StdExpansionSharedPtr > pCollExp, std::shared_ptr< CoalescedGeomData > GeomData, StdRegions::FactorMap factors)
	Constructor. More...

virtual	~Operator ()=default

virtual COLLECTIONS_EXPORT void	operator() (const Array< OneD, const NekDouble > &input, Array< OneD, NekDouble > &output0, Array< OneD, NekDouble > &output1, Array< OneD, NekDouble > &output2, Array< OneD, NekDouble > &wsp=NullNekDouble1DArray)=0
	Perform operation. More...

virtual COLLECTIONS_EXPORT void	operator() (int dir, const Array< OneD, const NekDouble > &input, Array< OneD, NekDouble > &output, Array< OneD, NekDouble > &wsp=NullNekDouble1DArray)=0

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

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

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

unsigned int	GetNumElmt ()
	Get number of elements. More...

StdRegions::StdExpansionSharedPtr	GetExpSharedPtr ()
	Get expansion pointer. More...

unsigned int	GetInputSize (bool defaultIn=true)

unsigned int	GetOutputSize (bool defaultOut=true)

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_jacWStdW

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

StdRegions::StdExpansionSharedPtr	m_stdExp

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

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 More...

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

unsigned int	m_inputSizeOther
	Number of modes or quadrature points, opposite to m_inputSize. More...

unsigned int	m_outputSizeOther
	Number of modes or quadrature points, opposite to m_outputSize. More...

Private Member Functions
	IProductWRTDerivBase_SumFac_Prism (vector< StdRegions::StdExpansionSharedPtr > 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 (Prism)

Definition at line 1576 of file IProductWRTDerivBase.cpp.

Constructor & Destructor Documentation

◆ ~IProductWRTDerivBase_SumFac_Prism()

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

finaldefault

◆ IProductWRTDerivBase_SumFac_Prism()

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

inlineprivate

Definition at line 1747 of file IProductWRTDerivBase.cpp.

        : Operator(pCollExp, pGeomData, factors), IProductWRTDerivBase_Helper(),
          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_jacWStdW = 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]);
                }
            }
        }
    }

References Nektar::LibUtilities::eModified_A, m_derivFac, m_fac0, m_fac1, m_jacWStdW, m_nmodes0, m_nmodes1, m_nmodes2, m_nquad0, m_nquad1, m_nquad2, Nektar::Collections::Operator::m_numElmt, m_sortTopVertex, Nektar::Collections::Operator::m_stdExp, and Nektar::Collections::Operator::m_wspSize.

Member Function Documentation

◆ operator()() [1/2]

void Nektar::Collections::IProductWRTDerivBase_SumFac_Prism::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]) + (d\phi/dz,in[2]) \]

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 + d\phi/d\xi_2\, d\xi_2/dx),in[0]) + \]

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

\[ ((d\phi/d\xi_0\, d\xi_0/dz + d\phi/d\xi_1\, d\xi_1/dz + d\phi/d\xi_2\, d\xi_2/dz),in[2]) \, \]

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_2) \]

\[ d\phi/d\xi_2 = d\phi/d\eta_0 d\eta_0/d\xi_2 + d\phi/d\eta_2 d\eta_2/d\xi_2 = d\phi/d\eta_0 (1+\eta_0)/(1-\eta_2) + d\phi/d\eta_2 \]

and so the full inner products are

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

\[ (d\phi/d\eta_0, ((2/(1-\eta_2) (d\xi_0/dx in[0] + d\xi_0/dy in[1] + d\xi_0/dz in[2]) + (1-\eta_0)/(1-\eta_2) (d\xi_2/dx in[0] + d\xi_2/dy in[1] + d\xi_2/dz in[2] )) + \]

\[ (d\phi/d\eta_1, (d\xi_1/dx in[0] + d\xi_1/dy in[1] + d\xi_1/dz in[2])) + \]

\[ (d\phi/d\eta_2, (d\xi_2/dx in[0] + d\xi_2/dy in[1] + d\xi_2/dz in[2])) \]

Implements Nektar::Collections::Operator.

Definition at line 1630 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(3);
 
        in[0] = entry0;
        in[1] = entry1;
        in[2] = entry2;
 
        output = entry3;
 
        for (int i = 0; i < 3; ++i)
        {
            tmp[i] = wsp + i * nmax;
        }
 
        if (m_isDeformed)
        {
            // calculate dx/dxi in[0] + dy/dxi in[1] + dz/dxi 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);
                }
            }
        }
        else
        {
            Array<OneD, NekDouble> t;
            for (int e = 0; e < m_numElmt; ++e)
            {
                // calculate dx/dxi in[0] + dy/dxi in[1] + dz/dxi in[2]
                for (int i = 0; i < 3; ++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 < 3; ++j)
                    {
                        Vmath::Svtvp(m_nqe, m_derivFac[i + 3 * j][e],
                                     in[j] + e * m_nqe, 1, tmp[i] + e * m_nqe,
                                     1, t = tmp[i] + e * m_nqe, 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].data() + i * nPhys, 1,
                        tmp[0].data() + i * nPhys, 1);
 
            // scale tmp[2] by fac1 and add to tmp0
            Vmath::Vvtvp(nPhys, &m_fac1[0], 1, tmp[2].data() + i * nPhys, 1,
                         tmp[0].data() + i * nPhys, 1,
                         tmp[0].data() + 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_jacWStdW, 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_jacWStdW, 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_jacWStdW, tmp[2], tmp[0], wsp1);
        Vmath::Vadd(m_numElmt * nmodes, tmp[0], 1, output, 1, output, 1);
    }

References m_base0, m_base1, m_base2, m_derbase0, m_derbase1, m_derbase2, m_derivFac, m_fac0, m_fac1, Nektar::Collections::Operator::m_isDeformed, m_jacWStdW, m_nmodes0, m_nmodes1, m_nmodes2, Nektar::Collections::Operator::m_nqe, m_nquad0, m_nquad1, m_nquad2, Nektar::Collections::Operator::m_numElmt, m_sortTopVertex, Nektar::Collections::Operator::m_stdExp, Nektar::Collections::PrismIProduct(), Vmath::Smul(), Vmath::Svtvp(), Vmath::Vadd(), Vmath::Vmul(), and Vmath::Vvtvp().

◆ operator()() [2/2]

void Nektar::Collections::IProductWRTDerivBase_SumFac_Prism::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 1719 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_Prism::m_base0

protected

Definition at line 1735 of file IProductWRTDerivBase.cpp.

Referenced by operator()().

◆ m_base1

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

protected

Definition at line 1736 of file IProductWRTDerivBase.cpp.

Referenced by operator()().

◆ m_base2

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

protected

Definition at line 1737 of file IProductWRTDerivBase.cpp.

Referenced by operator()().

◆ m_derbase0

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

protected

Definition at line 1738 of file IProductWRTDerivBase.cpp.

Referenced by operator()().

◆ m_derbase1

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

protected

Definition at line 1739 of file IProductWRTDerivBase.cpp.

Referenced by operator()().

◆ m_derbase2

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

protected

Definition at line 1740 of file IProductWRTDerivBase.cpp.

Referenced by operator()().

◆ m_derivFac

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

protected

Definition at line 1741 of file IProductWRTDerivBase.cpp.

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

◆ m_fac0

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

protected

Definition at line 1742 of file IProductWRTDerivBase.cpp.

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

◆ m_fac1

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

protected

Definition at line 1743 of file IProductWRTDerivBase.cpp.

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

◆ m_jacWStdW

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

protected

Definition at line 1734 of file IProductWRTDerivBase.cpp.

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

◆ m_nmodes0

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

protected

Definition at line 1731 of file IProductWRTDerivBase.cpp.

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

◆ m_nmodes1

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

protected

Definition at line 1732 of file IProductWRTDerivBase.cpp.

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

◆ m_nmodes2

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

protected

Definition at line 1733 of file IProductWRTDerivBase.cpp.

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

◆ m_nquad0

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

protected

Definition at line 1728 of file IProductWRTDerivBase.cpp.

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

◆ m_nquad1

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

protected

Definition at line 1729 of file IProductWRTDerivBase.cpp.

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

◆ m_nquad2

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

protected

Definition at line 1730 of file IProductWRTDerivBase.cpp.

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

◆ m_sortTopVertex

bool Nektar::Collections::IProductWRTDerivBase_SumFac_Prism::m_sortTopVertex

protected

Definition at line 1744 of file IProductWRTDerivBase.cpp.

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

Public Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ ~IProductWRTDerivBase_SumFac_Prism()

◆ IProductWRTDerivBase_SumFac_Prism()

Member Function Documentation

◆ operator()() [1/2]

◆ operator()() [2/2]

Member Data Documentation

◆ m_base0

◆ m_base1

◆ m_base2

◆ m_derbase0

◆ m_derbase1

◆ m_derbase2

◆ m_derivFac

◆ m_fac0

◆ m_fac1

◆ m_jacWStdW

◆ m_nmodes0

◆ m_nmodes1

◆ m_nmodes2

◆ m_nquad0

◆ m_nquad1

◆ m_nquad2

◆ m_sortTopVertex