Nektar::Collections::BwdTrans_SumFac_Tet Class Reference

Backward transform operator using sum-factorisation (Tet) More...

Inheritance diagram for Nektar::Collections::BwdTrans_SumFac_Tet:

Public Member Functions
	~BwdTrans_SumFac_Tet () final=default
void	operator() (const Array< OneD, const NekDouble > &input, Array< OneD, NekDouble > &output0, Array< OneD, NekDouble > &output1, Array< OneD, NekDouble > &output2, Array< OneD, NekDouble > &wsp) final
	Perform operation. More...
void	operator() (int dir, const Array< OneD, const NekDouble > &input, Array< OneD, NekDouble > &output, Array< OneD, NekDouble > &wsp) final
void	CheckFactors (StdRegions::FactorMap factors, int coll_phys_offset) override
	Check the validity of the supplied factor map. More...
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	CheckFactors (StdRegions::FactorMap factors, int coll_phys_offset)=0
	Check the validity of the supplied factor map. 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 ()
unsigned int	GetOutputSize ()

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_base0
Array< OneD, const NekDouble >	m_base1
Array< OneD, const NekDouble >	m_base2
bool	m_sortTopEdge
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...

Private Member Functions
	BwdTrans_SumFac_Tet (vector< StdRegions::StdExpansionSharedPtr > pCollExp, CoalescedGeomDataSharedPtr pGeomData, StdRegions::FactorMap factors)

Additional Inherited Members
Protected Member Functions inherited from Nektar::Collections::Operator
virtual int	v_GetInputSize ()
	This purely virtual function needs to be set-up for every operator inside Collections. It is responsible for returning the size of input collection, that the operator is applied on either in physical or coefficient space. More...
virtual int	v_GetOutputSize ()
	This purely virtual function needs to be set-up for every operator inside Collections. It is responsible for returning the output size either in physical or coefficient space of an operator inside Collections. More...
Protected Member Functions inherited from Nektar::Collections::BwdTrans_Helper
	BwdTrans_Helper ()

Detailed Description

Backward transform operator using sum-factorisation (Tet)

Definition at line 830 of file BwdTrans.cpp.

Constructor & Destructor Documentation

~BwdTrans_SumFac_Tet()

Nektar::Collections::BwdTrans_SumFac_Tet::~BwdTrans_SumFac_Tet ( )

finaldefault

BwdTrans_SumFac_Tet()

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

inlineprivate

Definition at line 965 of file BwdTrans.cpp.

        : Operator(pCollExp, pGeomData, factors), BwdTrans_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_wspSize = m_numElmt * (m_nquad2 * m_nmodes0 *
                                     (2 * m_nmodes1 - m_nmodes0 + 1) / 2 +
                                 m_nquad2 * m_nquad1 * m_nmodes0);
 
        if (m_stdExp->GetBasis(0)->GetBasisType() == LibUtilities::eModified_A)
        {
            m_sortTopEdge = true;
        }
        else
        {
            m_sortTopEdge = false;
        }
    }

References Nektar::LibUtilities::eModified_A, m_nmodes0, m_nmodes1, m_nquad1, m_nquad2, Nektar::Collections::Operator::m_numElmt, m_sortTopEdge, Nektar::Collections::Operator::m_stdExp, and Nektar::Collections::Operator::m_wspSize.

Member Function Documentation

CheckFactors()

void Nektar::Collections::BwdTrans_SumFac_Tet::CheckFactors ( StdRegions::FactorMap  factors, int  coll_phys_offset  )

inlineoverridevirtual

Check the validity of the supplied factor map.

Implements Nektar::Collections::Operator.

Definition at line 946 of file BwdTrans.cpp.

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

References ASSERTL0.

operator()() [1/2]

void Nektar::Collections::BwdTrans_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  )

inlinefinalvirtual

Perform operation.

Implements Nektar::Collections::Operator.

Definition at line 838 of file BwdTrans.cpp.

    {
        ASSERTL1(wsp.size() == m_wspSize, "Incorrect workspace size");
 
        Array<OneD, NekDouble> tmp = wsp;
        Array<OneD, NekDouble> tmp1 =
            tmp + m_numElmt * m_nquad2 * m_nmodes0 *
                      (2 * m_nmodes1 - m_nmodes0 + 1) / 2;
 
        int mode    = 0;
        int mode1   = 0;
        int cnt     = 0;
        int ncoeffs = m_stdExp->GetNcoeffs();
 
        // Perform summation over '2' direction
        for (int i = 0; i < m_nmodes0; ++i)
        {
            for (int j = 0; j < m_nmodes1 - i; ++j, ++cnt)
            {
                Blas::Dgemm('N', 'N', m_nquad2, m_numElmt, m_nmodes2 - i - j,
                            1.0, m_base2.get() + mode * m_nquad2, m_nquad2,
                            input.get() + mode1, ncoeffs, 0.0,
                            tmp.get() + cnt * m_nquad2 * m_numElmt, m_nquad2);
                mode += m_nmodes2 - i - j;
                mode1 += m_nmodes2 - i - j;
            }
 
            // increment mode in case m_nmodes1!=m_nmodes2
            mode += (m_nmodes2 - m_nmodes1) * (m_nmodes2 - m_nmodes1 + 1) / 2;
        }
 
        // vertex mode - currently (1+c)/2 x (1-b)/2 x (1-a)/2
        // component is evaluated
        if (m_sortTopEdge)
        {
            for (int i = 0; i < m_numElmt; ++i)
            {
                // top singular vertex
                // (1+c)/2 x (1+b)/2 x (1-a)/2 component
                Blas::Daxpy(m_nquad2, input[1 + i * ncoeffs],
                            m_base2.get() + m_nquad2, 1,
                            &tmp[m_nquad2 * m_numElmt] + i * m_nquad2, 1);
 
                // top singular vertex
                // (1+c)/2 x (1-b)/2 x (1+a)/2 component
                Blas::Daxpy(
                    m_nquad2, input[1 + i * ncoeffs], m_base2.get() + m_nquad2,
                    1, &tmp[m_nmodes1 * m_nquad2 * m_numElmt] + i * m_nquad2,
                    1);
            }
        }
 
        // Perform summation over '1' direction
        mode = 0;
        for (int i = 0; i < m_nmodes0; ++i)
        {
            Blas::Dgemm('N', 'T', m_nquad1, m_nquad2 * m_numElmt, m_nmodes1 - i,
                        1.0, m_base1.get() + mode * m_nquad1, m_nquad1,
                        tmp.get() + mode * m_nquad2 * m_numElmt,
                        m_nquad2 * m_numElmt, 0.0,
                        tmp1.get() + i * m_nquad1 * m_nquad2 * m_numElmt,
                        m_nquad1);
            mode += m_nmodes1 - i;
        }
 
        // fix for modified basis by adding additional split of
        // top and base singular vertex modes as well as singular
        // edge
        if (m_sortTopEdge)
        {
            // this could probably be a dgemv or higher if we
            // made a specialised m_base1[m_nuqad1] array
            // containing multiply copies
            for (int i = 0; i < m_numElmt; ++i)
            {
                // sort out singular vertices and singular
                // edge components with (1+b)/2 (1+a)/2 form
                for (int j = 0; j < m_nquad2; ++j)
                {
                    Blas::Daxpy(m_nquad1,
                                tmp[m_nquad2 * m_numElmt + i * m_nquad2 + j],
                                m_base1.get() + m_nquad1, 1,
                                &tmp1[m_nquad1 * m_nquad2 * m_numElmt] +
                                    i * m_nquad1 * m_nquad2 + j * m_nquad1,
                                1);
                }
            }
        }
 
        // Perform summation over '0' direction
        Blas::Dgemm('N', 'T', m_nquad0, m_nquad1 * m_nquad2 * m_numElmt,
                    m_nmodes0, 1.0, m_base0.get(), m_nquad0, tmp1.get(),
                    m_nquad1 * m_nquad2 * m_numElmt, 0.0, output0.get(),
                    m_nquad0);
    }

References ASSERTL1, Blas::Daxpy(), Blas::Dgemm(), m_base0, m_base1, m_base2, m_nmodes0, m_nmodes1, m_nmodes2, m_nquad0, m_nquad1, m_nquad2, Nektar::Collections::Operator::m_numElmt, m_sortTopEdge, Nektar::Collections::Operator::m_stdExp, and Nektar::Collections::Operator::m_wspSize.

operator()() [2/2]

void Nektar::Collections::BwdTrans_SumFac_Tet::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 938 of file BwdTrans.cpp.

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

References ASSERTL0.

Member Data Documentation

m_base0

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

protected

Definition at line 959 of file BwdTrans.cpp.

Referenced by operator()().

m_base1

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

protected

Definition at line 960 of file BwdTrans.cpp.

Referenced by operator()().

m_base2

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

protected

Definition at line 961 of file BwdTrans.cpp.

Referenced by operator()().

m_nmodes0

const int Nektar::Collections::BwdTrans_SumFac_Tet::m_nmodes0

protected

Definition at line 956 of file BwdTrans.cpp.

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

m_nmodes1

const int Nektar::Collections::BwdTrans_SumFac_Tet::m_nmodes1

protected

Definition at line 957 of file BwdTrans.cpp.

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

m_nmodes2

const int Nektar::Collections::BwdTrans_SumFac_Tet::m_nmodes2

protected

Definition at line 958 of file BwdTrans.cpp.

Referenced by operator()().

m_nquad0

const int Nektar::Collections::BwdTrans_SumFac_Tet::m_nquad0

protected

Definition at line 953 of file BwdTrans.cpp.

Referenced by operator()().

m_nquad1

const int Nektar::Collections::BwdTrans_SumFac_Tet::m_nquad1

protected

Definition at line 954 of file BwdTrans.cpp.

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

m_nquad2

const int Nektar::Collections::BwdTrans_SumFac_Tet::m_nquad2

protected

Definition at line 955 of file BwdTrans.cpp.

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

m_sortTopEdge

bool Nektar::Collections::BwdTrans_SumFac_Tet::m_sortTopEdge

protected

Definition at line 962 of file BwdTrans.cpp.

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