46 namespace LibUtilities
59 for(
unsigned int i=0; i < numPoints-2; ++i)
80 for(
unsigned int j=0; j < 3; ++j)
93 for(
unsigned int j=0; j < 3; ++j)
107 for(
unsigned int j=0; j < 6; ++j)
173 for(
int i = 0; i < rows; ++i ) {
174 for(
int j = 0; j < cols; ++j ) {
175 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++)
void CalculateDerivMatrix()
static boost::shared_ptr< PointsBaseType > Create(const PointsKey &key)
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]
NekVector< NekDouble > MakeQuadratureWeights(const NekVector< NekDouble > &x, const NekVector< NekDouble > &y)
static const NekDouble kNekZeroTol
static const unsigned int NodalTriElecNPTS[NodalTriElecAvailable]
unsigned int GetNumPoints() const
void CalculateInterpMatrix(const Array< OneD, const NekDouble > &xia, const Array< OneD, const NekDouble > &yia, Array< OneD, NekDouble > &interp)
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()
void NodalPointReorder2d()
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)
static const unsigned int perm6_2d[6][3]
unsigned int GetTotNumPoints() const
#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)
static const NekDouble NodalTriElecData[][6]
Array< OneD, DataType > & GetPtr()
virtual void CalculateWeights()