Nektar++
Public Member Functions | Static Public Member Functions | Protected Member Functions | Private Member Functions | Static Private Attributes | List of all members
Nektar::LibUtilities::FourierSingleModePoints Class Reference

#include <FourierSingleModePoints.h>

Inheritance diagram for Nektar::LibUtilities::FourierSingleModePoints:
[legend]

Public Member Functions

virtual ~FourierSingleModePoints ()
 
std::shared_ptr< NekMatrix< NekDouble > > CreateMatrix (const PointsKey &pkey)
 
 FourierSingleModePoints (const PointsKey &key)
 
- Public Member Functions inherited from Nektar::LibUtilities::Points< NekDouble >
virtual ~Points ()
 
void Initialize (void)
 
size_t GetPointsDim () const
 
size_t GetNumPoints () const
 
size_t GetTotNumPoints () const
 
PointsType GetPointsType () const
 
const Array< OneD, const DataType > & GetZ () const
 
const Array< OneD, const DataType > & GetW () const
 
void GetZW (Array< OneD, const DataType > &z, Array< OneD, const DataType > &w) const
 
const Array< OneD, const NekDouble > & GetBaryWeights () const
 
void GetPoints (Array< OneD, const DataType > &x) const
 
void GetPoints (Array< OneD, const DataType > &x, Array< OneD, const DataType > &y) const
 
void GetPoints (Array< OneD, const DataType > &x, Array< OneD, const DataType > &y, Array< OneD, const DataType > &z) const
 
const MatrixSharedPtrTypeGetD (Direction dir=xDir) const
 
const MatrixSharedPtrType GetI (const PointsKey &key)
 
const MatrixSharedPtrType GetI (const Array< OneD, const DataType > &x)
 
const MatrixSharedPtrType GetI (size_t uint, const Array< OneD, const DataType > &x)
 
const MatrixSharedPtrType GetI (const Array< OneD, const DataType > &x, const Array< OneD, const DataType > &y)
 
const MatrixSharedPtrType GetI (const Array< OneD, const DataType > &x, const Array< OneD, const DataType > &y, const Array< OneD, const DataType > &z)
 
const MatrixSharedPtrType GetGalerkinProjection (const PointsKey &pkey)
 

Static Public Member Functions

static std::shared_ptr< PointsBaseTypeCreate (const PointsKey &key)
 

Protected Member Functions

virtual const MatrixSharedPtrType v_GetI (const PointsKey &pkey) override final
 
virtual const MatrixSharedPtrType v_GetI (const Array< OneD, const NekDouble > &x) override final
 
virtual const MatrixSharedPtrType v_GetI (size_t numpoints, const Array< OneD, const NekDouble > &x) override final
 
- Protected Member Functions inherited from Nektar::LibUtilities::Points< NekDouble >
virtual void v_Initialize (void)
 
virtual void v_CalculatePoints ()
 
virtual void v_CalculateWeights ()
 

Private Member Functions

 FourierSingleModePoints ()=delete
 
 FourierSingleModePoints (const FourierSingleModePoints &points)=delete
 
virtual void v_CalculatePoints () override
 
virtual void v_CalculateWeights () override
 
virtual void v_CalculateDerivMatrix () override
 
void CalculateInterpMatrix (size_t npts, const Array< OneD, const NekDouble > &xpoints, Array< OneD, NekDouble > &interp)
 
NekDouble PeriodicSincFunction (const NekDouble x, const NekDouble h)
 

Static Private Attributes

static bool initPointsManager []
 

Additional Inherited Members

- Public Types inherited from Nektar::LibUtilities::Points< NekDouble >
typedef NekDouble DataType
 
typedef std::shared_ptr< NekMatrix< DataType > > MatrixSharedPtrType
 
- Protected Attributes inherited from Nektar::LibUtilities::Points< NekDouble >
PointsKey m_pointsKey
 Points type for this points distributions. More...
 
Array< OneD, DataTypem_points [3]
 Storage for the point locations, allowing for up to a 3D points storage. More...
 
Array< OneD, DataTypem_weights
 Quadrature weights for the weights. More...
 
Array< OneD, DataTypem_bcweights
 Barycentric weights. More...
 
MatrixSharedPtrType m_derivmatrix [3]
 Derivative matrices. More...
 
NekManager< PointsKey, NekMatrix< DataType >, PointsKey::opLessm_InterpManager
 
NekManager< PointsKey, NekMatrix< DataType >, PointsKey::opLessm_GalerkinProjectionManager
 

Detailed Description

Definition at line 44 of file FourierSingleModePoints.h.

Constructor & Destructor Documentation

◆ ~FourierSingleModePoints()

virtual Nektar::LibUtilities::FourierSingleModePoints::~FourierSingleModePoints ( )
inlinevirtual

Definition at line 47 of file FourierSingleModePoints.h.

48 {
49 }

◆ FourierSingleModePoints() [1/3]

Nektar::LibUtilities::FourierSingleModePoints::FourierSingleModePoints ( const PointsKey key)
inline

Definition at line 56 of file FourierSingleModePoints.h.

56 : PointsBaseType(key)
57 {
58 namespace pl = std::placeholders;
59 m_InterpManager.RegisterCreator(
60 PointsKey(0, eGaussGaussLegendre),
61 std::bind(&FourierSingleModePoints::CreateMatrix, this, pl::_1));
62 m_InterpManager.RegisterCreator(
63 PointsKey(0, eGaussRadauMLegendre),
64 std::bind(&FourierSingleModePoints::CreateMatrix, this, pl::_1));
65 m_InterpManager.RegisterCreator(
66 PointsKey(0, eGaussRadauPLegendre),
67 std::bind(&FourierSingleModePoints::CreateMatrix, this, pl::_1));
68 m_InterpManager.RegisterCreator(
69 PointsKey(0, eGaussLobattoLegendre),
70 std::bind(&FourierSingleModePoints::CreateMatrix, this, pl::_1));
71 m_InterpManager.RegisterCreator(
72 PointsKey(0, eGaussGaussChebyshev),
73 std::bind(&FourierSingleModePoints::CreateMatrix, this, pl::_1));
74 m_InterpManager.RegisterCreator(
75 PointsKey(0, eGaussRadauMChebyshev),
76 std::bind(&FourierSingleModePoints::CreateMatrix, this, pl::_1));
77 m_InterpManager.RegisterCreator(
78 PointsKey(0, eGaussRadauPChebyshev),
79 std::bind(&FourierSingleModePoints::CreateMatrix, this, pl::_1));
80 m_InterpManager.RegisterCreator(
81 PointsKey(0, eGaussLobattoChebyshev),
82 std::bind(&FourierSingleModePoints::CreateMatrix, this, pl::_1));
83 m_InterpManager.RegisterCreator(
84 PointsKey(0, eGaussRadauMAlpha0Beta1),
85 std::bind(&FourierSingleModePoints::CreateMatrix, this, pl::_1));
86 m_InterpManager.RegisterCreator(
87 PointsKey(0, eGaussRadauMAlpha0Beta2),
88 std::bind(&FourierSingleModePoints::CreateMatrix, this, pl::_1));
89 m_InterpManager.RegisterCreator(
90 PointsKey(0, eGaussRadauMAlpha1Beta0),
91 std::bind(&FourierSingleModePoints::CreateMatrix, this, pl::_1));
92 m_InterpManager.RegisterCreator(
93 PointsKey(0, eGaussRadauMAlpha2Beta0),
94 std::bind(&FourierSingleModePoints::CreateMatrix, this, pl::_1));
95 m_InterpManager.RegisterCreator(
96 PointsKey(0, ePolyEvenlySpaced),
97 std::bind(&FourierSingleModePoints::CreateMatrix, this, pl::_1));
98 m_InterpManager.RegisterCreator(
99 PointsKey(0, eFourierEvenlySpaced),
100 std::bind(&FourierSingleModePoints::CreateMatrix, this, pl::_1));
101 m_InterpManager.RegisterCreator(
102 PointsKey(0, eFourierSingleModeSpaced),
103 std::bind(&FourierSingleModePoints::CreateMatrix, this, pl::_1));
104 }
std::shared_ptr< NekMatrix< NekDouble > > CreateMatrix(const PointsKey &pkey)
NekManager< PointsKey, NekMatrix< DataType >, PointsKey::opLess > m_InterpManager
Definition: Points.h:369
Points< NekDouble > PointsBaseType
@ eGaussRadauMLegendre
1D Gauss-Radau-Legendre quadrature points, pinned at x=-1
Definition: PointsType.h:49
@ eGaussRadauMAlpha0Beta1
Gauss Radau pinned at x=-1,.
Definition: PointsType.h:60
@ eGaussLobattoChebyshev
1D Gauss-Lobatto-Legendre quadrature points
Definition: PointsType.h:59
@ eFourierEvenlySpaced
1D Evenly-spaced points using Fourier Fit
Definition: PointsType.h:76
@ eGaussRadauPChebyshev
1D Gauss-Radau-Chebyshev quadrature points, pinned at x=1
Definition: PointsType.h:57
@ eGaussRadauMChebyshev
1D Gauss-Radau-Chebyshev quadrature points, pinned at x=-1
Definition: PointsType.h:55
@ eGaussLobattoLegendre
1D Gauss-Lobatto-Legendre quadrature points
Definition: PointsType.h:53
@ eGaussGaussChebyshev
1D Gauss-Gauss-Chebyshev quadrature points
Definition: PointsType.h:54
@ ePolyEvenlySpaced
1D Evenly-spaced points using Lagrange polynomial
Definition: PointsType.h:75
@ eGaussGaussLegendre
1D Gauss-Gauss-Legendre quadrature points
Definition: PointsType.h:48
@ eFourierSingleModeSpaced
1D Non Evenly-spaced points for Single Mode analysis
Definition: PointsType.h:77
@ eGaussRadauPLegendre
1D Gauss-Radau-Legendre quadrature points, pinned at x=1
Definition: PointsType.h:51

References CreateMatrix(), Nektar::LibUtilities::eFourierEvenlySpaced, Nektar::LibUtilities::eFourierSingleModeSpaced, Nektar::LibUtilities::eGaussGaussChebyshev, Nektar::LibUtilities::eGaussGaussLegendre, Nektar::LibUtilities::eGaussLobattoChebyshev, Nektar::LibUtilities::eGaussLobattoLegendre, Nektar::LibUtilities::eGaussRadauMAlpha0Beta1, Nektar::LibUtilities::eGaussRadauMChebyshev, Nektar::LibUtilities::eGaussRadauMLegendre, Nektar::LibUtilities::eGaussRadauPChebyshev, Nektar::LibUtilities::eGaussRadauPLegendre, Nektar::LibUtilities::ePolyEvenlySpaced, and Nektar::LibUtilities::Points< NekDouble >::m_InterpManager.

◆ FourierSingleModePoints() [2/3]

Nektar::LibUtilities::FourierSingleModePoints::FourierSingleModePoints ( )
privatedelete

◆ FourierSingleModePoints() [3/3]

Nektar::LibUtilities::FourierSingleModePoints::FourierSingleModePoints ( const FourierSingleModePoints points)
privatedelete

Member Function Documentation

◆ CalculateInterpMatrix()

void Nektar::LibUtilities::FourierSingleModePoints::CalculateInterpMatrix ( size_t  npts,
const Array< OneD, const NekDouble > &  xpoints,
Array< OneD, NekDouble > &  interp 
)
private

Definition at line 138 of file FourierSingleModePoints.cpp.

141{
142 boost::ignore_unused(npts, xpoints, interp);
143
144 ASSERTL0(
145 false,
146 "CalculateInterpMatrix not available for Fourier Single Mode Points");
147}
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:215

References ASSERTL0.

◆ Create()

std::shared_ptr< Points< NekDouble > > Nektar::LibUtilities::FourierSingleModePoints::Create ( const PointsKey key)
static

Definition at line 86 of file FourierSingleModePoints.cpp.

88{
89 std::shared_ptr<Points<NekDouble>> returnval(
91
92 returnval->Initialize();
93
94 return returnval;
95}
static std::shared_ptr< DataType > AllocateSharedPtr(const Args &...args)
Allocate a shared pointer from the memory pool.

◆ CreateMatrix()

std::shared_ptr< NekMatrix< NekDouble > > Nektar::LibUtilities::FourierSingleModePoints::CreateMatrix ( const PointsKey pkey)

Definition at line 97 of file FourierSingleModePoints.cpp.

99{
100 boost::ignore_unused(pkey);
101
102 ASSERTL0(false,
103 "CreateMatrix not available for Fourier Single Mode Points");
104
105 return nullptr;
106}

References ASSERTL0.

Referenced by FourierSingleModePoints().

◆ PeriodicSincFunction()

NekDouble Nektar::LibUtilities::FourierSingleModePoints::PeriodicSincFunction ( const NekDouble  x,
const NekDouble  h 
)
private

◆ v_CalculateDerivMatrix()

void Nektar::LibUtilities::FourierSingleModePoints::v_CalculateDerivMatrix ( )
overrideprivatevirtual

Definition at line 80 of file FourierSingleModePoints.cpp.

81{
82
83 Points<NekDouble>::v_CalculateDerivMatrix();
84}

◆ v_CalculatePoints()

void Nektar::LibUtilities::FourierSingleModePoints::v_CalculatePoints ( )
overrideprivatevirtual

Reimplemented from Nektar::LibUtilities::Points< NekDouble >.

Definition at line 48 of file FourierSingleModePoints.cpp.

49{
50 // Allocate the storage for points
52 size_t npts = m_pointsKey.GetNumPoints();
53
54 if (npts == 1)
55 {
56 m_points[0][0] = 0.25;
57 }
58 else
59 {
60 ASSERTL0(npts == 2,
61 "Fourier points for single mode analysis should be 2");
62
63 m_points[0][0] = 0.0;
64 m_points[0][1] = 0.5;
65 }
66}
Array< OneD, DataType > m_points[3]
Storage for the point locations, allowing for up to a 3D points storage.
Definition: Points.h:361
PointsKey m_pointsKey
Points type for this points distributions.
Definition: Points.h:358
size_t GetNumPoints() const
Definition: Points.h:90

References ASSERTL0, Nektar::LibUtilities::PointsKey::GetNumPoints(), Nektar::LibUtilities::Points< NekDouble >::m_points, Nektar::LibUtilities::Points< NekDouble >::m_pointsKey, and Nektar::LibUtilities::Points< NekDouble >::v_CalculatePoints().

◆ v_CalculateWeights()

void Nektar::LibUtilities::FourierSingleModePoints::v_CalculateWeights ( )
overrideprivatevirtual

Reimplemented from Nektar::LibUtilities::Points< NekDouble >.

Definition at line 68 of file FourierSingleModePoints.cpp.

69{
70 // Allocate the storage for points
72
73 size_t npts = m_pointsKey.GetNumPoints();
74 for (size_t i = 0; i < npts; ++i)
75 {
76 m_weights[i] = 1.0;
77 }
78}
Array< OneD, DataType > m_weights
Quadrature weights for the weights.
Definition: Points.h:363

References Nektar::LibUtilities::PointsKey::GetNumPoints(), Nektar::LibUtilities::Points< NekDouble >::m_pointsKey, Nektar::LibUtilities::Points< NekDouble >::m_weights, and Nektar::LibUtilities::Points< NekDouble >::v_CalculateWeights().

◆ v_GetI() [1/3]

const std::shared_ptr< NekMatrix< NekDouble > > Nektar::LibUtilities::FourierSingleModePoints::v_GetI ( const Array< OneD, const NekDouble > &  x)
finaloverrideprotectedvirtual

Definition at line 118 of file FourierSingleModePoints.cpp.

120{
121 boost::ignore_unused(x);
122
123 ASSERTL0(false, "Interp not available for Fourier Single Mode Points");
124
125 return nullptr;
126}

References ASSERTL0.

◆ v_GetI() [2/3]

const std::shared_ptr< NekMatrix< NekDouble > > Nektar::LibUtilities::FourierSingleModePoints::v_GetI ( const PointsKey pkey)
finaloverrideprotectedvirtual

Definition at line 108 of file FourierSingleModePoints.cpp.

110{
111 boost::ignore_unused(pkey);
112
113 ASSERTL0(false, "Interp not available for Fourier Single Mode Points");
114
115 return nullptr;
116}

References ASSERTL0.

◆ v_GetI() [3/3]

const std::shared_ptr< NekMatrix< NekDouble > > Nektar::LibUtilities::FourierSingleModePoints::v_GetI ( size_t  numpoints,
const Array< OneD, const NekDouble > &  x 
)
finaloverrideprotectedvirtual

Definition at line 128 of file FourierSingleModePoints.cpp.

130{
131 boost::ignore_unused(numpoints, x);
132
133 ASSERTL0(false, "Interp not available for Fourier Single Mode Points");
134
135 return nullptr;
136}

References ASSERTL0.

Member Data Documentation

◆ initPointsManager

bool Nektar::LibUtilities::FourierSingleModePoints::initPointsManager
staticprivate
Initial value:
= {
static std::shared_ptr< PointsBaseType > Create(const PointsKey &key)
bool RegisterCreator(const KeyType &key, const CreateFuncType &createFunc)
Register the given function and associate it with the key. The return value is just to facilitate cal...
Definition: NekManager.hpp:169
PointsManagerT & PointsManager(void)

Definition at line 115 of file FourierSingleModePoints.h.