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

#include <GaussPoints.h>

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

Public Member Functions

virtual ~GaussPoints ()
 
std::shared_ptr< NekMatrix< NekDouble > > CreateMatrix (const PointsKey &pkey)
 
std::shared_ptr< NekMatrix< NekDouble > > CreateGPMatrix (const PointsKey &pkey)
 
 GaussPoints (const PointsKey &pkey)
 
- Public Member Functions inherited from Nektar::LibUtilities::Points< NekDouble >
virtual ~Points ()
 
void Initialize (void)
 
unsigned int GetPointsDim () const
 
unsigned int GetNumPoints () const
 
unsigned int 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 MatrixSharedPtrType & GetD (Direction dir=xDir) const
 
const MatrixSharedPtrType GetI (const PointsKey &key)
 
const MatrixSharedPtrType GetI (const Array< OneD, const DataType > &x)
 
const MatrixSharedPtrType GetI (unsigned int 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< Points< NekDouble > > Create (const PointsKey &pkey)
 

Protected Member Functions

virtual const std::shared_ptr< NekMatrix< NekDouble > > v_GetI (const PointsKey &pkey) override
 
virtual const std::shared_ptr< NekMatrix< NekDouble > > v_GetI (const Array< OneD, const NekDouble > &x) override
 
virtual const std::shared_ptr< NekMatrix< NekDouble > > v_GetI (unsigned int numpoints, const Array< OneD, const NekDouble > &x) override
 
virtual const std::shared_ptr< NekMatrix< NekDouble > > v_GetGalerkinProjection (const PointsKey &pkey) override
 
- Protected Member Functions inherited from Nektar::LibUtilities::Points< NekDouble >
virtual void v_Initialize (void)
 
virtual void v_CalculateBaryWeights ()
 This function calculates the barycentric weights used for enhanced interpolation speed. More...
 
 Points (const PointsKey &key)
 
virtual const MatrixSharedPtrType v_GetI (const Array< OneD, const DataType > &x, const Array< OneD, const DataType > &y)
 
virtual const MatrixSharedPtrType v_GetI (const Array< OneD, const DataType > &x, const Array< OneD, const DataType > &y, const Array< OneD, const DataType > &z)
 

Private Member Functions

 GaussPoints ()
 These should not be called. All creation is done using the constructor requiring the key, declared above. More...
 
 GaussPoints (const GaussPoints &points)
 
virtual void v_CalculatePoints () override
 
virtual void v_CalculateWeights () override
 
virtual void v_CalculateDerivMatrix () override
 
void CalculateInterpMatrix (unsigned int npts, const Array< OneD, const NekDouble > &xpoints, Array< OneD, NekDouble > &interp)
 
std::shared_ptr< NekMatrix< NekDouble > > CalculateGalerkinProjectionMatrix (const PointsKey &pkey)
 
NekDouble LagrangeInterpolant (NekDouble x, int npts, const Array< OneD, const NekDouble > &xpts, const Array< OneD, const NekDouble > &funcvals)
 functions used by the Kronrod points More...
 
NekDouble LagrangePoly (NekDouble x, int pt, int npts, const Array< OneD, const NekDouble > &xpts)
 
NekDouble LagrangePolyDeriv (NekDouble x, int pt, int npts, const Array< OneD, const NekDouble > &xpts)
 

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, DataType > m_points [3]
 Storage for the point locations, allowing for up to a 3D points storage. More...
 
Array< OneD, DataType > m_weights
 Quadrature weights for the weights. More...
 
Array< OneD, DataType > m_bcweights
 Barycentric weights. More...
 
MatrixSharedPtrType m_derivmatrix [3]
 Derivative matrices. More...
 
NekManager< PointsKey, NekMatrix< DataType >, PointsKey::opLess > m_InterpManager
 
NekManager< PointsKey, NekMatrix< DataType >, PointsKey::opLess > m_GalerkinProjectionManager
 

Detailed Description

Definition at line 47 of file GaussPoints.h.

Constructor & Destructor Documentation

◆ ~GaussPoints()

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

Definition at line 50 of file GaussPoints.h.

51  {
52  }

◆ GaussPoints() [1/3]

Nektar::LibUtilities::GaussPoints::GaussPoints ( const PointsKey &  pkey)
inline

Definition at line 63 of file GaussPoints.h.

63  : PointsBaseType(pkey)
64  {
65  namespace pl = std::placeholders;
66  m_InterpManager.RegisterCreator(
67  PointsKey(0, eGaussGaussLegendre),
68  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
69  m_InterpManager.RegisterCreator(
70  PointsKey(0, eGaussRadauMLegendre),
71  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
72  m_InterpManager.RegisterCreator(
73  PointsKey(0, eGaussRadauPLegendre),
74  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
75  m_InterpManager.RegisterCreator(
76  PointsKey(0, eGaussLobattoLegendre),
77  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
78  m_InterpManager.RegisterCreator(
79  PointsKey(0, eGaussGaussChebyshev),
80  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
81  m_InterpManager.RegisterCreator(
82  PointsKey(0, eGaussRadauMChebyshev),
83  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
84  m_InterpManager.RegisterCreator(
85  PointsKey(0, eGaussRadauPChebyshev),
86  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
87  m_InterpManager.RegisterCreator(
88  PointsKey(0, eGaussLobattoChebyshev),
89  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
90  m_InterpManager.RegisterCreator(
91  PointsKey(0, eGaussRadauMAlpha0Beta1),
92  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
93  m_InterpManager.RegisterCreator(
94  PointsKey(0, eGaussRadauMAlpha0Beta2),
95  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
96  m_InterpManager.RegisterCreator(
97  PointsKey(0, eGaussRadauMAlpha1Beta0),
98  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
99  m_InterpManager.RegisterCreator(
100  PointsKey(0, eGaussRadauMAlpha2Beta0),
101  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
102  m_InterpManager.RegisterCreator(
103  PointsKey(0, eGaussKronrodLegendre),
104  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
105  m_InterpManager.RegisterCreator(
106  PointsKey(0, eGaussRadauKronrodMLegendre),
107  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
108  m_InterpManager.RegisterCreator(
109  PointsKey(0, eGaussRadauKronrodMAlpha1Beta0),
110  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
111  m_InterpManager.RegisterCreator(
112  PointsKey(0, eGaussLobattoKronrodLegendre),
113  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
114  m_InterpManager.RegisterCreator(
115  PointsKey(0, eFourierEvenlySpaced),
116  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
117  m_InterpManager.RegisterCreator(
118  PointsKey(0, ePolyEvenlySpaced),
119  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
120  m_InterpManager.RegisterCreator(
121  PointsKey(0, eBoundaryLayerPoints),
122  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
123  m_InterpManager.RegisterCreator(
124  PointsKey(0, eBoundaryLayerPointsRev),
125  std::bind(&GaussPoints::CreateMatrix, this, pl::_1));
126  m_GalerkinProjectionManager.RegisterCreator(
127  PointsKey(0, eGaussGaussLegendre),
128  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
129  m_GalerkinProjectionManager.RegisterCreator(
130  PointsKey(0, eGaussRadauMLegendre),
131  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
132  m_GalerkinProjectionManager.RegisterCreator(
133  PointsKey(0, eGaussRadauPLegendre),
134  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
135  m_GalerkinProjectionManager.RegisterCreator(
136  PointsKey(0, eGaussLobattoLegendre),
137  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
138  m_GalerkinProjectionManager.RegisterCreator(
139  PointsKey(0, eGaussGaussChebyshev),
140  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
141  m_GalerkinProjectionManager.RegisterCreator(
142  PointsKey(0, eGaussRadauMChebyshev),
143  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
144  m_GalerkinProjectionManager.RegisterCreator(
145  PointsKey(0, eGaussRadauPChebyshev),
146  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
147  m_GalerkinProjectionManager.RegisterCreator(
148  PointsKey(0, eGaussLobattoChebyshev),
149  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
150  m_GalerkinProjectionManager.RegisterCreator(
151  PointsKey(0, eGaussRadauMAlpha0Beta1),
152  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
153  m_GalerkinProjectionManager.RegisterCreator(
154  PointsKey(0, eGaussRadauMAlpha0Beta2),
155  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
156  m_GalerkinProjectionManager.RegisterCreator(
157  PointsKey(0, eGaussRadauMAlpha1Beta0),
158  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
159  m_GalerkinProjectionManager.RegisterCreator(
160  PointsKey(0, eGaussRadauMAlpha2Beta0),
161  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
162  m_GalerkinProjectionManager.RegisterCreator(
163  PointsKey(0, eGaussKronrodLegendre),
164  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
165  m_GalerkinProjectionManager.RegisterCreator(
166  PointsKey(0, eGaussRadauKronrodMLegendre),
167  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
168  m_GalerkinProjectionManager.RegisterCreator(
169  PointsKey(0, eGaussRadauKronrodMAlpha1Beta0),
170  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
171  m_GalerkinProjectionManager.RegisterCreator(
172  PointsKey(0, eGaussLobattoKronrodLegendre),
173  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
174  m_GalerkinProjectionManager.RegisterCreator(
175  PointsKey(0, eFourierEvenlySpaced),
176  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
177  m_GalerkinProjectionManager.RegisterCreator(
178  PointsKey(0, ePolyEvenlySpaced),
179  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
180  m_GalerkinProjectionManager.RegisterCreator(
181  PointsKey(0, eBoundaryLayerPoints),
182  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
183  m_GalerkinProjectionManager.RegisterCreator(
184  PointsKey(0, eBoundaryLayerPointsRev),
185  std::bind(&GaussPoints::CreateGPMatrix, this, pl::_1));
186  }
Nektar::LibUtilities::GaussPoints::CreateGPMatrix
std::shared_ptr< NekMatrix< NekDouble > > CreateGPMatrix(const PointsKey &pkey)
Definition: GaussPoints.cpp:481
Nektar::LibUtilities::GaussPoints::CreateMatrix
std::shared_ptr< NekMatrix< NekDouble > > CreateMatrix(const PointsKey &pkey)
Definition: GaussPoints.cpp:370
Nektar::LibUtilities::Points< NekDouble >::m_InterpManager
NekManager< PointsKey, NekMatrix< DataType >, PointsKey::opLess > m_InterpManager
Definition: Points.h:383
Nektar::LibUtilities::Points< NekDouble >::m_GalerkinProjectionManager
NekManager< PointsKey, NekMatrix< DataType >, PointsKey::opLess > m_GalerkinProjectionManager
Definition: Points.h:385
Nektar::LibUtilities::PointsBaseType
Points< NekDouble > PointsBaseType
Definition: FoundationsFwd.hpp:74
Nektar::LibUtilities::eGaussRadauMLegendre
@ eGaussRadauMLegendre
1D Gauss-Radau-Legendre quadrature points, pinned at x=-1
Definition: PointsType.h:49
Nektar::LibUtilities::eBoundaryLayerPoints
@ eBoundaryLayerPoints
1D power law distribution for boundary layer points
Definition: PointsType.h:79
Nektar::LibUtilities::eGaussRadauMAlpha0Beta1
@ eGaussRadauMAlpha0Beta1
Gauss Radau pinned at x=-1,.
Definition: PointsType.h:60
Nektar::LibUtilities::eGaussRadauKronrodMLegendre
@ eGaussRadauKronrodMLegendre
1D Gauss-Radau-Kronrod-Legendre quadrature points, pinned at x=-1
Definition: PointsType.h:69
Nektar::LibUtilities::eGaussLobattoChebyshev
@ eGaussLobattoChebyshev
1D Gauss-Lobatto-Legendre quadrature points
Definition: PointsType.h:59
Nektar::LibUtilities::eFourierEvenlySpaced
@ eFourierEvenlySpaced
1D Evenly-spaced points using Fourier Fit
Definition: PointsType.h:76
Nektar::LibUtilities::eGaussRadauPChebyshev
@ eGaussRadauPChebyshev
1D Gauss-Radau-Chebyshev quadrature points, pinned at x=1
Definition: PointsType.h:57
Nektar::LibUtilities::eGaussRadauMChebyshev
@ eGaussRadauMChebyshev
1D Gauss-Radau-Chebyshev quadrature points, pinned at x=-1
Definition: PointsType.h:55
Nektar::LibUtilities::eGaussRadauKronrodMAlpha1Beta0
@ eGaussRadauKronrodMAlpha1Beta0
1D Gauss-Radau-Kronrod-Legendre pinned at x=-1,
Definition: PointsType.h:71
Nektar::LibUtilities::eGaussLobattoLegendre
@ eGaussLobattoLegendre
1D Gauss-Lobatto-Legendre quadrature points
Definition: PointsType.h:53
Nektar::LibUtilities::eGaussGaussChebyshev
@ eGaussGaussChebyshev
1D Gauss-Gauss-Chebyshev quadrature points
Definition: PointsType.h:54
Nektar::LibUtilities::eBoundaryLayerPointsRev
@ eBoundaryLayerPointsRev
1D power law distribution for boundary layer points
Definition: PointsType.h:81
Nektar::LibUtilities::ePolyEvenlySpaced
@ ePolyEvenlySpaced
1D Evenly-spaced points using Lagrange polynomial
Definition: PointsType.h:75
Nektar::LibUtilities::eGaussLobattoKronrodLegendre
@ eGaussLobattoKronrodLegendre
1D Lobatto Kronrod quadrature points
Definition: PointsType.h:74
Nektar::LibUtilities::eGaussGaussLegendre
@ eGaussGaussLegendre
1D Gauss-Gauss-Legendre quadrature points
Definition: PointsType.h:48
Nektar::LibUtilities::eGaussRadauPLegendre
@ eGaussRadauPLegendre
1D Gauss-Radau-Legendre quadrature points, pinned at x=1
Definition: PointsType.h:51

References CreateGPMatrix(), CreateMatrix(), Nektar::LibUtilities::eBoundaryLayerPoints, Nektar::LibUtilities::eBoundaryLayerPointsRev, Nektar::LibUtilities::eFourierEvenlySpaced, Nektar::LibUtilities::eGaussGaussChebyshev, Nektar::LibUtilities::eGaussGaussLegendre, Nektar::LibUtilities::eGaussLobattoChebyshev, Nektar::LibUtilities::eGaussLobattoKronrodLegendre, Nektar::LibUtilities::eGaussLobattoLegendre, Nektar::LibUtilities::eGaussRadauKronrodMAlpha1Beta0, Nektar::LibUtilities::eGaussRadauKronrodMLegendre, Nektar::LibUtilities::eGaussRadauMAlpha0Beta1, Nektar::LibUtilities::eGaussRadauMChebyshev, Nektar::LibUtilities::eGaussRadauMLegendre, Nektar::LibUtilities::eGaussRadauPChebyshev, Nektar::LibUtilities::eGaussRadauPLegendre, Nektar::LibUtilities::ePolyEvenlySpaced, Nektar::LibUtilities::Points< NekDouble >::m_GalerkinProjectionManager, and Nektar::LibUtilities::Points< NekDouble >::m_InterpManager.

◆ GaussPoints() [2/3]

Nektar::LibUtilities::GaussPoints::GaussPoints ( )
private

These should not be called. All creation is done using the constructor requiring the key, declared above.

◆ GaussPoints() [3/3]

Nektar::LibUtilities::GaussPoints::GaussPoints ( const GaussPoints &  points)
private

Member Function Documentation

◆ CalculateGalerkinProjectionMatrix()

std::shared_ptr< NekMatrix< NekDouble > > Nektar::LibUtilities::GaussPoints::CalculateGalerkinProjectionMatrix ( const PointsKey &  pkey)
private

Definition at line 491 of file GaussPoints.cpp.

493 {
494  int numpointsfrom = pkey.GetNumPoints();
495  int numpointsto = GetNumPoints();
496 
497  Array<OneD, const NekDouble> weightsfrom;
498 
499  weightsfrom = PointsManager()[pkey]->GetW();
500 
501  std::shared_ptr<NekMatrix<NekDouble>> Interp = GetI(pkey);
502 
503  Array<OneD, NekDouble> GalProj(numpointsfrom * numpointsto);
504 
505  // set up inner product matrix and multiply by inverse of
506  // diagaonal mass matrix
507  for (int i = 0; i < numpointsto; ++i)
508  {
509  Vmath::Vmul(numpointsfrom, Interp->GetPtr().get() + i * numpointsfrom,
510  1, &weightsfrom[0], 1, &GalProj[0] + i, numpointsto);
511  Vmath::Smul(numpointsfrom, 1.0 / m_weights[i], &GalProj[0] + i,
512  numpointsto, &GalProj[0] + i, numpointsto);
513  }
514 
515  NekDouble *t = GalProj.data();
516  std::shared_ptr<NekMatrix<NekDouble>> returnval(
517  MemoryManager<NekMatrix<NekDouble>>::AllocateSharedPtr(
518  numpointsto, numpointsfrom, t));
519 
520  return returnval;
521 }
Nektar::LibUtilities::Points< NekDouble >::GetNumPoints
unsigned int GetNumPoints() const
Definition: Points.h:273
Nektar::LibUtilities::Points< NekDouble >::m_weights
Array< OneD, DataType > m_weights
Quadrature weights for the weights.
Definition: Points.h:377
Nektar::LibUtilities::Points< NekDouble >::GetI
const MatrixSharedPtrType GetI(const PointsKey &key)
Definition: Points.h:336
Nektar::LibUtilities::PointsManager
PointsManagerT & PointsManager(void)
Definition: ManagerAccess.cpp:42
Nektar::NekDouble
double NekDouble
Definition: NektarUnivTypeDefs.hpp:43
Vmath::Vmul
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.cpp:209
Vmath::Smul
void Smul(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha*x.
Definition: Vmath.cpp:248

References Nektar::LibUtilities::Points< NekDouble >::GetI(), Nektar::LibUtilities::PointsKey::GetNumPoints(), Nektar::LibUtilities::Points< NekDouble >::GetNumPoints(), Nektar::LibUtilities::Points< NekDouble >::m_weights, Nektar::LibUtilities::PointsManager(), Vmath::Smul(), and Vmath::Vmul().

Referenced by CreateGPMatrix().

◆ CalculateInterpMatrix()

void Nektar::LibUtilities::GaussPoints::CalculateInterpMatrix ( unsigned int  npts,
const Array< OneD, const NekDouble > &  xpoints,
Array< OneD, NekDouble > &  interp 
)
private

Definition at line 272 of file GaussPoints.cpp.

275 {
276  switch (m_pointsKey.GetPointsType())
277  {
278  case eGaussGaussLegendre:
279  Polylib::Imgj(interp.data(), m_points[0].data(), xpoints.data(),
280  GetNumPoints(), npts, 0.0, 0.0);
281  break;
282 
283  case eGaussRadauMLegendre:
284  Polylib::Imgrjm(interp.data(), m_points[0].data(), xpoints.data(),
285  GetNumPoints(), npts, 0.0, 0.0);
286  break;
287 
288  case eGaussRadauPLegendre:
289  Polylib::Imgrjp(interp.data(), m_points[0].data(), xpoints.data(),
290  GetNumPoints(), npts, 0.0, 0.0);
291  break;
292 
293  case eGaussLobattoLegendre:
294  Polylib::Imglj(interp.data(), m_points[0].data(), xpoints.data(),
295  GetNumPoints(), npts, 0.0, 0.0);
296  break;
297 
298  case eGaussGaussChebyshev:
299  Polylib::Imgj(interp.data(), m_points[0].data(), xpoints.data(),
300  GetNumPoints(), npts, -0.5, -0.5);
301  break;
302 
303  case eGaussRadauMChebyshev:
304  Polylib::Imgrjm(interp.data(), m_points[0].data(), xpoints.data(),
305  GetNumPoints(), npts, -0.5, -0.5);
306  break;
307 
308  case eGaussRadauPChebyshev:
309  Polylib::Imgrjp(interp.data(), m_points[0].data(), xpoints.data(),
310  GetNumPoints(), npts, -0.5, -0.5);
311  break;
312 
313  case eGaussLobattoChebyshev:
314  Polylib::Imglj(interp.data(), m_points[0].data(), xpoints.data(),
315  GetNumPoints(), npts, -0.5, -0.5);
316  break;
317 
318  case eGaussRadauMAlpha0Beta1:
319  Polylib::Imgrjm(interp.data(), m_points[0].data(), xpoints.data(),
320  GetNumPoints(), npts, 0.0, 1.0);
321  break;
322 
323  case eGaussRadauMAlpha0Beta2:
324  Polylib::Imgrjm(interp.data(), m_points[0].data(), xpoints.data(),
325  GetNumPoints(), npts, 0.0, 2.0);
326  break;
327 
328  case eGaussRadauMAlpha1Beta0:
329  Polylib::Imgrjm(interp.data(), m_points[0].data(), xpoints.data(),
330  GetNumPoints(), npts, 1.0, 0.0);
331  break;
332 
333  case eGaussRadauMAlpha2Beta0:
334  Polylib::Imgrjm(interp.data(), m_points[0].data(), xpoints.data(),
335  GetNumPoints(), npts, 2.0, 0.0);
336  break;
337 
338  case eGaussKronrodLegendre:
339  case eGaussRadauKronrodMLegendre:
340  case eGaussRadauKronrodMAlpha1Beta0:
341  case eGaussLobattoKronrodLegendre:
342  {
343  for (unsigned int i = 0; i < npts; ++i)
344  {
345  for (unsigned int j = 0; j < m_pointsKey.GetNumPoints(); ++j)
346  {
347  interp[i + j * npts] = LagrangePoly(
348  xpoints[i], j, m_pointsKey.GetNumPoints(), m_points[0]);
349  }
350  }
351  }
352  break;
353 
354  default:
355  NEKERROR(ErrorUtil::efatal,
356  "Unknown Gauss quadrature point distribution requested");
357  }
358 }
NEKERROR
#define NEKERROR(type, msg)
Assert Level 0 – Fundamental assert which is used whether in FULLDEBUG, DEBUG or OPT compilation mode...
Definition: ErrorUtil.hpp:209
Nektar::LibUtilities::GaussPoints::LagrangePoly
NekDouble LagrangePoly(NekDouble x, int pt, int npts, const Array< OneD, const NekDouble > &xpts)
Definition: GaussPoints.cpp:429
Nektar::LibUtilities::Points< NekDouble >::m_points
Array< OneD, DataType > m_points[3]
Storage for the point locations, allowing for up to a 3D points storage.
Definition: Points.h:375
Nektar::LibUtilities::Points< NekDouble >::m_pointsKey
PointsKey m_pointsKey
Points type for this points distributions.
Definition: Points.h:372
Nektar::LibUtilities::PointsKey::GetPointsType
PointsType GetPointsType() const
Definition: Points.h:109
Nektar::LibUtilities::PointsKey::GetNumPoints
unsigned int GetNumPoints() const
Definition: Points.h:104
Polylib::Imgrjp
void Imgrjp(double *im, const double *zgrj, const double *zm, const int nz, const int mz, const double alpha, const double beta)
Interpolation Operator from Gauss-Radau-Jacobi points (including z=1) to an arbitrary distrubtion at ...
Definition: Polylib.cpp:1054
Polylib::Imgj
void Imgj(double *im, const double *zgj, const double *zm, const int nz, const int mz, const double alpha, const double beta)
Interpolation Operator from Gauss-Jacobi points to an arbitrary distribution at points zm.
Definition: Polylib.cpp:992
Polylib::Imgrjm
void Imgrjm(double *im, const double *zgrj, const double *zm, const int nz, const int mz, const double alpha, const double beta)
Interpolation Operator from Gauss-Radau-Jacobi points (including z=-1) to an arbitrary distrubtion at...
Definition: Polylib.cpp:1023
Polylib::Imglj
void Imglj(double *im, const double *zglj, const double *zm, const int nz, const int mz, const double alpha, const double beta)
Interpolation Operator from Gauss-Lobatto-Jacobi points to an arbitrary distrubtion at points zm.
Definition: Polylib.cpp:1085

References Nektar::Array< OneD, DataType >::data(), Nektar::ErrorUtil::efatal, Nektar::LibUtilities::eGaussGaussChebyshev, Nektar::LibUtilities::eGaussGaussLegendre, Nektar::LibUtilities::eGaussLobattoChebyshev, Nektar::LibUtilities::eGaussLobattoKronrodLegendre, Nektar::LibUtilities::eGaussLobattoLegendre, Nektar::LibUtilities::eGaussRadauKronrodMAlpha1Beta0, Nektar::LibUtilities::eGaussRadauKronrodMLegendre, Nektar::LibUtilities::eGaussRadauMAlpha0Beta1, Nektar::LibUtilities::eGaussRadauMChebyshev, Nektar::LibUtilities::eGaussRadauMLegendre, Nektar::LibUtilities::eGaussRadauPChebyshev, Nektar::LibUtilities::eGaussRadauPLegendre, Nektar::LibUtilities::PointsKey::GetNumPoints(), Nektar::LibUtilities::Points< NekDouble >::GetNumPoints(), Nektar::LibUtilities::PointsKey::GetPointsType(), Polylib::Imgj(), Polylib::Imglj(), Polylib::Imgrjm(), Polylib::Imgrjp(), LagrangePoly(), Nektar::LibUtilities::Points< NekDouble >::m_points, Nektar::LibUtilities::Points< NekDouble >::m_pointsKey, and NEKERROR.

Referenced by v_GetI().

◆ Create()

std::shared_ptr< Points< NekDouble > > Nektar::LibUtilities::GaussPoints::Create ( const PointsKey &  pkey)
static

Definition at line 360 of file GaussPoints.cpp.

361 {
362  std::shared_ptr<Points<NekDouble>> returnval(
363  MemoryManager<GaussPoints>::AllocateSharedPtr(pkey));
364 
365  returnval->Initialize();
366 
367  return returnval;
368 }
Nektar::MemoryManager::AllocateSharedPtr
static std::shared_ptr< DataType > AllocateSharedPtr(const Args &...args)
Allocate a shared pointer from the memory pool.
Definition: NekMemoryManager.hpp:168

◆ CreateGPMatrix()

std::shared_ptr< NekMatrix< NekDouble > > Nektar::LibUtilities::GaussPoints::CreateGPMatrix ( const PointsKey &  pkey)

Definition at line 481 of file GaussPoints.cpp.

483 {
484  std::shared_ptr<NekMatrix<NekDouble>> returnval =
485  CalculateGalerkinProjectionMatrix(pkey);
486 
487  // Delegate to function below
488  return returnval;
489 }
Nektar::LibUtilities::GaussPoints::CalculateGalerkinProjectionMatrix
std::shared_ptr< NekMatrix< NekDouble > > CalculateGalerkinProjectionMatrix(const PointsKey &pkey)
Definition: GaussPoints.cpp:492

References CalculateGalerkinProjectionMatrix().

Referenced by GaussPoints().

◆ CreateMatrix()

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

Definition at line 370 of file GaussPoints.cpp.

372 {
373  int numpoints = pkey.GetNumPoints();
374  Array<OneD, const NekDouble> xpoints;
375 
376  PointsManager()[pkey]->GetPoints(xpoints);
377 
378  // Delegate to function below
379  return GetI(numpoints, xpoints);
380 }

References Nektar::LibUtilities::Points< NekDouble >::GetI(), Nektar::LibUtilities::PointsKey::GetNumPoints(), and Nektar::LibUtilities::PointsManager().

Referenced by GaussPoints().

◆ LagrangeInterpolant()

NekDouble Nektar::LibUtilities::GaussPoints::LagrangeInterpolant ( NekDouble  x,
int  npts,
const Array< OneD, const NekDouble > &  xpts,
const Array< OneD, const NekDouble > &  funcvals 
)
private

functions used by the Kronrod points

Definition at line 416 of file GaussPoints.cpp.

419 {
420  NekDouble sum = 0.0;
421 
422  for (int i = 0; i < npts; ++i)
423  {
424  sum += funcvals[i] * LagrangePoly(x, i, npts, xpts);
425  }
426  return sum;
427 }

References LagrangePoly().

◆ LagrangePoly()

NekDouble Nektar::LibUtilities::GaussPoints::LagrangePoly ( NekDouble  x,
int  pt,
int  npts,
const Array< OneD, const NekDouble > &  xpts 
)
private

Definition at line 429 of file GaussPoints.cpp.

431 {
432  NekDouble h = 1.0;
433 
434  for (int i = 0; i < pt; ++i)
435  {
436  h = h * (x - xpts[i]) / (xpts[pt] - xpts[i]);
437  }
438 
439  for (int i = pt + 1; i < npts; ++i)
440  {
441  h = h * (x - xpts[i]) / (xpts[pt] - xpts[i]);
442  }
443 
444  return h;
445 }

Referenced by CalculateInterpMatrix(), and LagrangeInterpolant().

◆ LagrangePolyDeriv()

NekDouble Nektar::LibUtilities::GaussPoints::LagrangePolyDeriv ( NekDouble  x,
int  pt,
int  npts,
const Array< OneD, const NekDouble > &  xpts 
)
private

Definition at line 447 of file GaussPoints.cpp.

449 {
450  NekDouble h;
451  NekDouble y = 0.0;
452 
453  for (int j = 0; j < npts; ++j)
454  {
455  if (j != pt)
456  {
457  h = 1.0;
458  for (int i = 0; i < npts; ++i)
459  {
460  if (i != pt)
461  {
462  if (i != j)
463  {
464  h = h * (x - xpts[i]);
465  }
466  h = h / (xpts[pt] - xpts[i]);
467  }
468  }
469  y = y + h;
470  }
471  }
472  return y;
473 }

Referenced by v_CalculateDerivMatrix().

◆ v_CalculateDerivMatrix()

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

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

Definition at line 184 of file GaussPoints.cpp.

185 {
186  // Allocate the derivative matrix
187  PointsBaseType::v_CalculateDerivMatrix();
188 
189  int numpoints = m_pointsKey.GetNumPoints();
190  int totpoints = m_pointsKey.GetTotNumPoints();
191  Array<OneD, NekDouble> dmtempSharedArray =
192  Array<OneD, NekDouble>(totpoints * totpoints);
193  NekDouble *dmtemp = dmtempSharedArray.data();
194 
195  switch (m_pointsKey.GetPointsType())
196  {
197  case eGaussGaussLegendre:
198  Polylib::Dgj(dmtemp, m_points[0].data(), numpoints, 0.0, 0.0);
199  break;
200 
201  case eGaussRadauMLegendre:
202  Polylib::Dgrjm(dmtemp, m_points[0].data(), numpoints, 0.0, 0.0);
203  break;
204 
205  case eGaussRadauPLegendre:
206  Polylib::Dgrjp(dmtemp, m_points[0].data(), numpoints, 0.0, 0.0);
207  break;
208 
209  case eGaussLobattoLegendre:
210  Polylib::Dglj(dmtemp, m_points[0].data(), numpoints, 0.0, 0.0);
211  break;
212 
213  case eGaussGaussChebyshev:
214  Polylib::Dgj(dmtemp, m_points[0].data(), numpoints, -0.5, -0.5);
215  break;
216 
217  case eGaussRadauMChebyshev:
218  Polylib::Dgrjm(dmtemp, m_points[0].data(), numpoints, -0.5, -0.5);
219  break;
220 
221  case eGaussRadauPChebyshev:
222  Polylib::Dgrjp(dmtemp, m_points[0].data(), numpoints, -0.5, -0.5);
223  break;
224 
225  case eGaussLobattoChebyshev:
226  Polylib::Dglj(dmtemp, m_points[0].data(), numpoints, -0.5, -0.5);
227  break;
228 
229  case eGaussRadauMAlpha0Beta1:
230  Polylib::Dgrjm(dmtemp, m_points[0].data(), numpoints, 0.0, 1.0);
231  break;
232 
233  case eGaussRadauMAlpha0Beta2:
234  Polylib::Dgrjm(dmtemp, m_points[0].data(), numpoints, 0.0, 2.0);
235  break;
236 
237  case eGaussRadauMAlpha1Beta0:
238  Polylib::Dgrjm(dmtemp, m_points[0].data(), numpoints, 1.0, 0.0);
239  break;
240 
241  case eGaussRadauMAlpha2Beta0:
242  Polylib::Dgrjm(dmtemp, m_points[0].data(), numpoints, 2.0, 0.0);
243  break;
244 
245  case eGaussKronrodLegendre:
246  case eGaussRadauKronrodMLegendre:
247  case eGaussRadauKronrodMAlpha1Beta0:
248  case eGaussLobattoKronrodLegendre:
249  {
250  for (unsigned int i = 0; i < m_pointsKey.GetNumPoints(); ++i)
251  {
252  for (unsigned int j = 0; j < m_pointsKey.GetNumPoints(); ++j)
253  {
254  (*m_derivmatrix[0])(i, j) = LagrangePolyDeriv(
255  m_points[0][i], j, m_pointsKey.GetNumPoints(),
256  m_points[0]);
257  }
258  }
259  return;
260  }
261  break;
262 
263  default:
264  NEKERROR(ErrorUtil::efatal,
265  "Unknown Gauss quadrature point distribution requested");
266  }
267 
268  std::copy(dmtemp, dmtemp + totpoints * totpoints,
269  m_derivmatrix[0]->begin());
270 }
Nektar::LibUtilities::GaussPoints::LagrangePolyDeriv
NekDouble LagrangePolyDeriv(NekDouble x, int pt, int npts, const Array< OneD, const NekDouble > &xpts)
Definition: GaussPoints.cpp:447
Nektar::LibUtilities::Points< NekDouble >::m_derivmatrix
MatrixSharedPtrType m_derivmatrix[3]
Derivative matrices.
Definition: Points.h:381
Nektar::LibUtilities::PointsKey::GetTotNumPoints
unsigned int GetTotNumPoints() const
Definition: Points.h:177
CellMLToNektar.pycml.copy
def copy(self)
Definition: pycml.py:2663
Polylib::Dgj
void Dgj(double *D, const double *z, const int np, const double alpha, const double beta)
Compute the Derivative Matrix and its transpose associated with the Gauss-Jacobi zeros.
Definition: Polylib.cpp:627
Polylib::Dglj
void Dglj(double *D, const double *z, const int np, const double alpha, const double beta)
Compute the Derivative Matrix associated with the Gauss-Lobatto-Jacobi zeros.
Definition: Polylib.cpp:787
Polylib::Dgrjm
void Dgrjm(double *D, const double *z, const int np, const double alpha, const double beta)
Compute the Derivative Matrix and its transpose associated with the Gauss-Radau-Jacobi zeros with a z...
Definition: Polylib.cpp:674
Polylib::Dgrjp
void Dgrjp(double *D, const double *z, const int np, const double alpha, const double beta)
Compute the Derivative Matrix associated with the Gauss-Radau-Jacobi zeros with a zero at z=1.
Definition: Polylib.cpp:730

References CellMLToNektar.pycml::copy(), Polylib::Dgj(), Polylib::Dglj(), Polylib::Dgrjm(), Polylib::Dgrjp(), Nektar::ErrorUtil::efatal, Nektar::LibUtilities::eGaussGaussChebyshev, Nektar::LibUtilities::eGaussGaussLegendre, Nektar::LibUtilities::eGaussLobattoChebyshev, Nektar::LibUtilities::eGaussLobattoKronrodLegendre, Nektar::LibUtilities::eGaussLobattoLegendre, Nektar::LibUtilities::eGaussRadauKronrodMAlpha1Beta0, Nektar::LibUtilities::eGaussRadauKronrodMLegendre, Nektar::LibUtilities::eGaussRadauMAlpha0Beta1, Nektar::LibUtilities::eGaussRadauMChebyshev, Nektar::LibUtilities::eGaussRadauMLegendre, Nektar::LibUtilities::eGaussRadauPChebyshev, Nektar::LibUtilities::eGaussRadauPLegendre, Nektar::LibUtilities::PointsKey::GetNumPoints(), Nektar::LibUtilities::PointsKey::GetPointsType(), Nektar::LibUtilities::PointsKey::GetTotNumPoints(), LagrangePolyDeriv(), Nektar::LibUtilities::Points< NekDouble >::m_derivmatrix, Nektar::LibUtilities::Points< NekDouble >::m_points, Nektar::LibUtilities::Points< NekDouble >::m_pointsKey, NEKERROR, and Nektar::LibUtilities::Points< NekDouble >::v_CalculateDerivMatrix().

◆ v_CalculatePoints()

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

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

Definition at line 83 of file GaussPoints.cpp.

84 {
85  // Allocate the storage for points and weights
86  PointsBaseType::v_CalculatePoints();
87  PointsBaseType::v_CalculateWeights();
88 
89  int numpoints = m_pointsKey.GetNumPoints();
90 
91  switch (m_pointsKey.GetPointsType())
92  {
93  case eGaussGaussLegendre:
94  Polylib::zwgj(m_points[0].data(), m_weights.data(), numpoints, 0.0,
95  0.0);
96  break;
97 
98  case eGaussRadauMLegendre:
99  Polylib::zwgrjm(m_points[0].data(), m_weights.data(), numpoints,
100  0.0, 0.0);
101  break;
102 
103  case eGaussRadauPLegendre:
104  Polylib::zwgrjp(m_points[0].data(), m_weights.data(), numpoints,
105  0.0, 0.0);
106  break;
107 
108  case eGaussLobattoLegendre:
109  Polylib::zwglj(m_points[0].data(), m_weights.data(), numpoints, 0.0,
110  0.0);
111  break;
112 
113  case eGaussGaussChebyshev:
114  Polylib::zwgj(m_points[0].data(), m_weights.data(), numpoints, -0.5,
115  -0.5);
116  break;
117 
118  case eGaussRadauMChebyshev:
119  Polylib::zwgrjm(m_points[0].data(), m_weights.data(), numpoints,
120  -0.5, -0.5);
121  break;
122 
123  case eGaussRadauPChebyshev:
124  Polylib::zwgrjp(m_points[0].data(), m_weights.data(), numpoints,
125  -0.5, -0.5);
126  break;
127 
128  case eGaussLobattoChebyshev:
129  Polylib::zwglj(m_points[0].data(), m_weights.data(), numpoints,
130  -0.5, -0.5);
131  break;
132 
133  case eGaussRadauMAlpha0Beta1:
134  Polylib::zwgrjm(m_points[0].data(), m_weights.data(), numpoints,
135  0.0, 1.0);
136  break;
137 
138  case eGaussRadauMAlpha0Beta2:
139  Polylib::zwgrjm(m_points[0].data(), m_weights.data(), numpoints,
140  0.0, 2.0);
141  break;
142 
143  case eGaussRadauMAlpha1Beta0:
144  Polylib::zwgrjm(m_points[0].data(), m_weights.data(), numpoints,
145  1.0, 0.0);
146  break;
147 
148  case eGaussRadauMAlpha2Beta0:
149  Polylib::zwgrjm(m_points[0].data(), m_weights.data(), numpoints,
150  2.0, 0.0);
151  break;
152 
153  case eGaussKronrodLegendre:
154  Polylib::zwgk(m_points[0].data(), m_weights.data(), numpoints, 0.0,
155  0.0);
156  break;
157 
158  case eGaussRadauKronrodMLegendre:
159  Polylib::zwrk(m_points[0].data(), m_weights.data(), numpoints, 0.0,
160  0.0);
161  break;
162 
163  case eGaussRadauKronrodMAlpha1Beta0:
164  Polylib::zwrk(m_points[0].data(), m_weights.data(), numpoints, 1.0,
165  0.0);
166  break;
167 
168  case eGaussLobattoKronrodLegendre:
169  Polylib::zwlk(m_points[0].data(), m_weights.data(), numpoints, 0.0,
170  0.0);
171  break;
172 
173  default:
174  NEKERROR(ErrorUtil::efatal,
175  "Unknown Gauss quadrature point distribution requested");
176  }
177 }
Polylib::zwgrjm
void zwgrjm(double *z, double *w, const int np, const double alpha, const double beta)
Gauss-Radau-Jacobi zeros and weights with end point at z=-1.
Definition: Polylib.cpp:161
Polylib::zwgrjp
void zwgrjp(double *z, double *w, const int np, const double alpha, const double beta)
Gauss-Radau-Jacobi zeros and weights with end point at z=1.
Definition: Polylib.cpp:202
Polylib::zwglj
void zwglj(double *z, double *w, const int np, const double alpha, const double beta)
Gauss-Lobatto-Jacobi zeros and weights with end point at z=-1,1.
Definition: Polylib.cpp:241
Polylib::zwgj
void zwgj(double *z, double *w, const int np, const double alpha, const double beta)
Gauss-Jacobi zeros and weights.
Definition: Polylib.cpp:133
Polylib::zwgk
void zwgk(double *z, double *w, const int npt, const double alpha, const double beta)
Gauss-Kronrod-Jacobi zeros and weights.
Definition: Polylib.cpp:290
Polylib::zwrk
void zwrk(double *z, double *w, const int npt, const double alpha, const double beta)
Gauss-Radau-Kronrod-Jacobi zeros and weights.
Definition: Polylib.cpp:362
Polylib::zwlk
void zwlk(double *z, double *w, const int npt, const double alpha, const double beta)
Gauss-Lobatto-Kronrod-Jacobi zeros and weights.
Definition: Polylib.cpp:468

References Nektar::Array< OneD, DataType >::data(), Nektar::ErrorUtil::efatal, Nektar::LibUtilities::eGaussGaussChebyshev, Nektar::LibUtilities::eGaussGaussLegendre, Nektar::LibUtilities::eGaussLobattoChebyshev, Nektar::LibUtilities::eGaussLobattoKronrodLegendre, Nektar::LibUtilities::eGaussLobattoLegendre, Nektar::LibUtilities::eGaussRadauKronrodMAlpha1Beta0, Nektar::LibUtilities::eGaussRadauKronrodMLegendre, Nektar::LibUtilities::eGaussRadauMAlpha0Beta1, Nektar::LibUtilities::eGaussRadauMChebyshev, Nektar::LibUtilities::eGaussRadauMLegendre, Nektar::LibUtilities::eGaussRadauPChebyshev, Nektar::LibUtilities::eGaussRadauPLegendre, Nektar::LibUtilities::PointsKey::GetNumPoints(), Nektar::LibUtilities::PointsKey::GetPointsType(), Nektar::LibUtilities::Points< NekDouble >::m_points, Nektar::LibUtilities::Points< NekDouble >::m_pointsKey, Nektar::LibUtilities::Points< NekDouble >::m_weights, NEKERROR, Nektar::LibUtilities::Points< NekDouble >::v_CalculatePoints(), Nektar::LibUtilities::Points< NekDouble >::v_CalculateWeights(), Polylib::zwgj(), Polylib::zwgk(), Polylib::zwglj(), Polylib::zwgrjm(), Polylib::zwgrjp(), Polylib::zwlk(), and Polylib::zwrk().

◆ v_CalculateWeights()

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

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

Definition at line 179 of file GaussPoints.cpp.

180 {
181  // For Gauss Quadrature, this is done as part of the points computation
182 }

◆ v_GetGalerkinProjection()

const std::shared_ptr< NekMatrix< NekDouble > > Nektar::LibUtilities::GaussPoints::v_GetGalerkinProjection ( const PointsKey &  pkey)
overrideprotectedvirtual

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

Definition at line 475 of file GaussPoints.cpp.

477 {
478  return m_GalerkinProjectionManager[pkey];
479 }

References Nektar::LibUtilities::Points< NekDouble >::m_GalerkinProjectionManager.

◆ v_GetI() [1/3]

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

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

Definition at line 391 of file GaussPoints.cpp.

393 {
394  int numpoints = 1;
395 
396  // Delegate to function below
397  return GetI(numpoints, x);
398 }

References Nektar::LibUtilities::Points< NekDouble >::GetI().

◆ v_GetI() [2/3]

const std::shared_ptr< NekMatrix< NekDouble > > Nektar::LibUtilities::GaussPoints::v_GetI ( const PointsKey &  pkey)
overrideprotectedvirtual

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

Definition at line 382 of file GaussPoints.cpp.

384 {
385  ASSERTL0(pkey.GetPointsDim() == 1,
386  "Gauss Points can only interp to other 1d point distributions");
387 
388  return m_InterpManager[pkey];
389 }
ASSERTL0
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:215

References ASSERTL0, Nektar::LibUtilities::PointsKey::GetPointsDim(), and Nektar::LibUtilities::Points< NekDouble >::m_InterpManager.

◆ v_GetI() [3/3]

const std::shared_ptr< NekMatrix< NekDouble > > Nektar::LibUtilities::GaussPoints::v_GetI ( unsigned int  numpoints,
const Array< OneD, const NekDouble > &  x 
)
overrideprotectedvirtual

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

Definition at line 400 of file GaussPoints.cpp.

402 {
403  Array<OneD, NekDouble> interp(GetNumPoints() * numpoints);
404 
405  CalculateInterpMatrix(numpoints, x, interp);
406 
407  NekDouble *t = interp.data();
408  unsigned int np = GetNumPoints();
409  std::shared_ptr<NekMatrix<NekDouble>> returnval(
410  MemoryManager<NekMatrix<NekDouble>>::AllocateSharedPtr(numpoints, np,
411  t));
412 
413  return returnval;
414 }
Nektar::LibUtilities::GaussPoints::CalculateInterpMatrix
void CalculateInterpMatrix(unsigned int npts, const Array< OneD, const NekDouble > &xpoints, Array< OneD, NekDouble > &interp)
Definition: GaussPoints.cpp:272

References CalculateInterpMatrix(), and Nektar::LibUtilities::Points< NekDouble >::GetNumPoints().

Member Data Documentation

◆ initPointsManager

bool Nektar::LibUtilities::GaussPoints::initPointsManager
staticprivate

Definition at line 199 of file GaussPoints.h.

Generated on Sat Mar 25 2023 12:35:58 for Nektar++ by doxygen 1.9.1