Nektar::Collections::BwdTrans_SumFac_Quad Class Reference

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

Inheritance diagram for Nektar::Collections::BwdTrans_SumFac_Quad:

Public Member Functions
virtual	~BwdTrans_SumFac_Quad ()
virtual void	operator() (const Array< OneD, const NekDouble > &input, Array< OneD, NekDouble > &output, Array< OneD, NekDouble > &output1, Array< OneD, NekDouble > &output2, 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, std::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_nmodes0
const int	m_nmodes1
const bool	m_colldir0
const bool	m_colldir1
Array< OneD, const NekDouble >	m_base0
Array< OneD, const NekDouble >	m_base1
Protected Attributes inherited from Nektar::Collections::Operator
StdRegions::StdExpansionSharedPtr	m_stdExp
unsigned int	m_numElmt
unsigned int	m_wspSize

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

Detailed Description

Backward transform operator using sum-factorisation (Quad)

Definition at line 394 of file BwdTrans.cpp.

Constructor & Destructor Documentation

~BwdTrans_SumFac_Quad()

virtual Nektar::Collections::BwdTrans_SumFac_Quad::~BwdTrans_SumFac_Quad ( )

inlinevirtual

Definition at line 399 of file BwdTrans.cpp.

        {
        }

BwdTrans_SumFac_Quad()

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

inlineprivate

Definition at line 481 of file BwdTrans.cpp.

References Nektar::StdRegions::eBwdTrans, Nektar::LibUtilities::eQuadrilateral, Nektar::Collections::eSumFac, and Nektar::Collections::GetOperatorFactory().

            : Operator  (pCollExp, pGeomData),
              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_wspSize = m_nquad0*m_nmodes1*m_numElmt;
        }

Member Function Documentation

operator()() [1/2]

virtual void Nektar::Collections::BwdTrans_SumFac_Quad::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 403 of file BwdTrans.cpp.

References ASSERTL1, Blas::Dgemm(), and Vmath::Vcopy().

        {
            boost::ignore_unused(output1, output2);

            int i = 0;
            if(m_colldir0 && m_colldir1)
            {
                Vmath::Vcopy(m_numElmt*m_nmodes0*m_nmodes1,
                             input.get(), 1, output.get(), 1);
            }
            else if(m_colldir0)
            {
                Array<OneD, const NekDouble> base1
                        = m_stdExp->GetBasis(1)->GetBdata();

                for(i = 0; i < m_numElmt; ++i)
                {
                    Blas::Dgemm('N', 'T', m_nquad0, m_nquad1, m_nmodes1,
                                1.0, &input[i*m_nquad0*m_nmodes1], m_nquad0,
                                     base1.get(), m_nquad1, 0.0,
                                     &output[i*m_nquad0*m_nquad1], m_nquad0);
                }
            }
            else if(m_colldir1)
            {
                Blas::Dgemm('N', 'N', m_nquad0, m_nmodes1*m_numElmt, m_nmodes0,
                            1.0, m_base0.get(), m_nquad0,
                            &input[0],  m_nmodes0, 0.0,
                            &output[0], m_nquad0);
            }
            else
            {
                ASSERTL1(wsp.num_elements() == m_wspSize,
                         "Incorrect workspace size");

                // Those two calls correpsond to the operation
                // out = B0*in*Transpose(B1);
                Blas::Dgemm('N', 'N', m_nquad0, m_nmodes1*m_numElmt, m_nmodes0,
                            1.0, m_base0.get(), m_nquad0,
                            &input[0], m_nmodes0, 0.0,
                            &wsp[0],   m_nquad0);

                for(i = 0; i < m_numElmt; ++i)
                {
                    Blas::Dgemm('N','T', m_nquad0, m_nquad1, m_nmodes1,
                                1.0, &wsp[i*m_nquad0*m_nmodes1], m_nquad0,
                                m_base1.get(), m_nquad1, 0.0,
                                &output[i*m_nquad0*m_nquad1],    m_nquad0);
                }
            }
        }

operator()() [2/2]

virtual void Nektar::Collections::BwdTrans_SumFac_Quad::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 460 of file BwdTrans.cpp.

References ASSERTL0.

        {
            boost::ignore_unused(dir, input, output, wsp);
            ASSERTL0(false, "Not valid for this operator.");
        }

Member Data Documentation

m_base0

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

protected

Definition at line 477 of file BwdTrans.cpp.

m_base1

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

protected

Definition at line 478 of file BwdTrans.cpp.

m_colldir0

const bool Nektar::Collections::BwdTrans_SumFac_Quad::m_colldir0

protected

Definition at line 475 of file BwdTrans.cpp.

m_colldir1

const bool Nektar::Collections::BwdTrans_SumFac_Quad::m_colldir1

protected

Definition at line 476 of file BwdTrans.cpp.

m_nmodes0

const int Nektar::Collections::BwdTrans_SumFac_Quad::m_nmodes0

protected

Definition at line 473 of file BwdTrans.cpp.

m_nmodes1

const int Nektar::Collections::BwdTrans_SumFac_Quad::m_nmodes1

protected

Definition at line 474 of file BwdTrans.cpp.

m_nquad0

const int Nektar::Collections::BwdTrans_SumFac_Quad::m_nquad0

protected

Definition at line 471 of file BwdTrans.cpp.

m_nquad1

const int Nektar::Collections::BwdTrans_SumFac_Quad::m_nquad1

protected

Definition at line 472 of file BwdTrans.cpp.