49 namespace LibUtilities
66 for(
unsigned int i=0; i < numPoints-2; ++i)
87 for(
unsigned int j=0; j < 3; ++j)
100 for(
unsigned int j=0; j < 3; ++j)
114 for(
unsigned int j=0; j < 6; ++j)
157 for(
int i = 0; i < rows; ++i ) {
158 for(
int j = 0; j < cols; ++j ) {
159 interp[j + i*cols] = interMat(i,j);
184 returnval->Initialize();
194 const int nVerts = 3;
196 const int nBoundaryPoints = 3*nEdgeInteriorPoints + 3;
198 if(nEdgeInteriorPoints==0)
205 for(i = cnt = istart; i < nBoundaryPoints; i++)
216 iend = istart + nEdgeInteriorPoints;
217 for(i = istart; i < iend; i++)
219 for(j = istart+1; j < iend; j++)
231 for(i = cnt = istart; i < nBoundaryPoints; i++)
242 iend = istart + nEdgeInteriorPoints;
243 for(i = istart; i < iend; i++)
245 for(j = istart+1; j < iend; j++)
257 for(i = cnt = istart; i < nBoundaryPoints; i++)
267 iend = istart + nEdgeInteriorPoints;
268 for(i = istart; i < iend; i++)
270 for(j = istart+1; j < iend; j++)
static const unsigned int perm3B_2d[3][3]
virtual void CalculatePoints()
General purpose memory allocation routines with the ability to allocate from thread specific memory p...
MatrixSharedPtrType m_derivmatrix[3]
Array< OneD, DataType > m_points[3]
static const unsigned int NodalTriFeketeNPTS[NodalTriFeketeAvailable]
NekVector< NekDouble > MakeQuadratureWeights(const NekVector< NekDouble > &x, const NekVector< NekDouble > &y)
static const NekDouble kNekZeroTol
void CalculateInterpMatrix(const Array< OneD, const NekDouble > &xi, const Array< OneD, const NekDouble > &yi, Array< OneD, NekDouble > &interp)
unsigned int GetNumPoints() const
unsigned int GetTotNumPoints() const
Points< NekDouble >::MatrixSharedPtrType GetXDerivativeMatrix(const NekVector< NekDouble > &x, const NekVector< NekDouble > &y, const NekVector< NekDouble > &xi, const NekVector< NekDouble > &yi)
Array< OneD, DataType > m_weights
Defines a specification for a set of points.
virtual void CalculateDerivMatrix()
static const unsigned int perm3A_2d[3][3]
unsigned int GetRows() const
NekMatrix< NekDouble > GetInterpolationMatrix(const NekVector< NekDouble > &x, const NekVector< NekDouble > &y, const NekVector< NekDouble > &xi, const NekVector< NekDouble > &yi)
void CalculateDerivMatrix()
void NodalPointReorder2d()
static boost::shared_ptr< PointsBaseType > Create(const PointsKey &key)
static const unsigned int perm6_2d[6][3]
unsigned int GetTotNumPoints() const
static const NekDouble NodalTriFeketeData[][6]
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Points< NekDouble >::MatrixSharedPtrType GetYDerivativeMatrix(const NekVector< NekDouble > &x, const NekVector< NekDouble > &y, const NekVector< NekDouble > &xi, const NekVector< NekDouble > &yi)
Array< OneD, DataType > & GetPtr()
virtual void CalculateWeights()