Nektar++
Public Member Functions | Protected Attributes | Private Member Functions | List of all members
Nektar::Collections::IProductWRTBase_SumFac_Seg Class Referencefinal

Inner product operator using sum-factorisation (Segment) More...

Inheritance diagram for Nektar::Collections::IProductWRTBase_SumFac_Seg:
[legend]

Public Member Functions

 ~IProductWRTBase_SumFac_Seg () final=default
 
void operator() (const Array< OneD, const NekDouble > &input, Array< OneD, NekDouble > &output, 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
 
- 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_nmodes0
 
const bool m_colldir0
 
Array< OneD, const NekDoublem_jacWStdW
 
Array< OneD, const NekDoublem_base0
 
- 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

 IProductWRTBase_SumFac_Seg (vector< StdRegions::StdExpansionSharedPtr > pCollExp, CoalescedGeomDataSharedPtr pGeomData, StdRegions::FactorMap factors)
 

Additional Inherited Members

- Protected Member Functions inherited from Nektar::Collections::IProductWRTBase_Helper
 IProductWRTBase_Helper ()
 

Detailed Description

Inner product operator using sum-factorisation (Segment)

Definition at line 456 of file IProductWRTBase.cpp.

Constructor & Destructor Documentation

◆ ~IProductWRTBase_SumFac_Seg()

Nektar::Collections::IProductWRTBase_SumFac_Seg::~IProductWRTBase_SumFac_Seg ( )
finaldefault

◆ IProductWRTBase_SumFac_Seg()

Nektar::Collections::IProductWRTBase_SumFac_Seg::IProductWRTBase_SumFac_Seg ( vector< StdRegions::StdExpansionSharedPtr pCollExp,
CoalescedGeomDataSharedPtr  pGeomData,
StdRegions::FactorMap  factors 
)
inlineprivate

Definition at line 502 of file IProductWRTBase.cpp.

505 : Operator(pCollExp, pGeomData, factors), IProductWRTBase_Helper(),
506 m_nquad0(m_stdExp->GetNumPoints(0)),
507 m_nmodes0(m_stdExp->GetBasisNumModes(0)),
508 m_colldir0(m_stdExp->GetBasis(0)->Collocation()),
509 m_base0(m_stdExp->GetBasis(0)->GetBdata())
510 {
512 m_jacWStdW = pGeomData->GetJacWithStdWeights(pCollExp);
513 }
StdRegions::StdExpansionSharedPtr m_stdExp
Definition: Operator.h:217
unsigned int m_numElmt
number of elements that the operator is applied on
Definition: Operator.h:219
Operator(std::vector< StdRegions::StdExpansionSharedPtr > pCollExp, std::shared_ptr< CoalescedGeomData > GeomData, StdRegions::FactorMap factors)
Constructor.
Definition: Operator.cpp:66
StdRegions::ConstFactorMap factors

References m_jacWStdW, m_nquad0, Nektar::Collections::Operator::m_numElmt, and Nektar::Collections::Operator::m_wspSize.

Member Function Documentation

◆ operator()() [1/2]

void Nektar::Collections::IProductWRTBase_SumFac_Seg::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 464 of file IProductWRTBase.cpp.

469 {
470 if (m_colldir0)
471 {
472 Vmath::Vmul(m_numElmt * m_nquad0, m_jacWStdW, 1, input, 1, output,
473 1);
474 }
475 else
476 {
477 Vmath::Vmul(m_numElmt * m_nquad0, m_jacWStdW, 1, input, 1, wsp, 1);
478
479 // out = B0*in;
480 Blas::Dgemm('T', 'N', m_nmodes0, m_numElmt, m_nquad0, 1.0,
481 m_base0.data(), m_nquad0, &wsp[0], m_nquad0, 0.0,
482 &output[0], m_nmodes0);
483 }
484 }
static void Dgemm(const char &transa, const char &transb, const int &m, const int &n, const int &k, const double &alpha, const double *a, const int &lda, const double *b, const int &ldb, const double &beta, double *c, const int &ldc)
BLAS level 3: Matrix-matrix multiply C = A x B where op(A)[m x k], op(B)[k x n], C[m x n] DGEMM perfo...
Definition: Blas.hpp:383
void Vmul(int n, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
Multiply vector z = x*y.
Definition: Vmath.hpp:72

References Blas::Dgemm(), m_base0, m_colldir0, m_jacWStdW, m_nmodes0, m_nquad0, Nektar::Collections::Operator::m_numElmt, and Vmath::Vmul().

◆ operator()() [2/2]

void Nektar::Collections::IProductWRTBase_SumFac_Seg::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 486 of file IProductWRTBase.cpp.

490 {
491 NEKERROR(ErrorUtil::efatal, "Not valid for this operator.");
492 }
#define NEKERROR(type, msg)
Assert Level 0 – Fundamental assert which is used whether in FULLDEBUG, DEBUG or OPT compilation mode...
Definition: ErrorUtil.hpp:202

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

Member Data Documentation

◆ m_base0

Array<OneD, const NekDouble> Nektar::Collections::IProductWRTBase_SumFac_Seg::m_base0
protected

Definition at line 499 of file IProductWRTBase.cpp.

Referenced by operator()().

◆ m_colldir0

const bool Nektar::Collections::IProductWRTBase_SumFac_Seg::m_colldir0
protected

Definition at line 497 of file IProductWRTBase.cpp.

Referenced by operator()().

◆ m_jacWStdW

Array<OneD, const NekDouble> Nektar::Collections::IProductWRTBase_SumFac_Seg::m_jacWStdW
protected

Definition at line 498 of file IProductWRTBase.cpp.

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

◆ m_nmodes0

const int Nektar::Collections::IProductWRTBase_SumFac_Seg::m_nmodes0
protected

Definition at line 496 of file IProductWRTBase.cpp.

Referenced by operator()().

◆ m_nquad0

const int Nektar::Collections::IProductWRTBase_SumFac_Seg::m_nquad0
protected

Definition at line 495 of file IProductWRTBase.cpp.

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