43 namespace LibUtilities
54 "Must set factor in BLPoints key");
56 if (fabs(r-1.0) < 1e-6)
59 for (
unsigned int i = 0; i <
npts; ++i)
69 for (
unsigned int i = 1; i <
npts; ++i)
79 vector<NekDouble> tmp(npts);
80 for (
unsigned int i = 0; i <
npts; ++i)
85 for (
unsigned int i = 0; i <
npts; ++i)
107 returnval->Initialize();
121 return GetI(numpoints, xpoints);
126 ASSERTL0(pkey.
GetPointsDim()==1,
"Fourier Points can only interp to other 1d point distributions");
136 return GetI(numpoints, x);
#define ASSERTL0(condition, msg)
virtual void CalculatePoints()
NekDouble GetFactor() const
General purpose memory allocation routines with the ability to allocate from thread specific memory p...
Array< OneD, DataType > m_points[3]
boost::shared_ptr< NekMatrix< NekDouble > > CreateMatrix(const PointsKey &pkey)
unsigned int GetNumPoints() const
unsigned int GetPointsDim() const
PointsManagerT & PointsManager(void)
Defines a specification for a set of points.
virtual void CalculateDerivMatrix()
static const NekDouble kNekUnsetDouble
void CalculateDerivMatrix()
1D power law distribution for boundary layer points
static boost::shared_ptr< PointsBaseType > Create(const PointsKey &key)
unsigned int GetNumPoints() const
void CalculateInterpMatrix(unsigned int npts, const Array< OneD, const NekDouble > &xpoints, Array< OneD, NekDouble > &interp)
const MatrixSharedPtrType GetI(const PointsKey &pkey)
NekManager< PointsKey, NekMatrix< DataType >, PointsKey::opLess > m_InterpManager
PointsType GetPointsType() const