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)
156 boost::shared_ptr<NekMatrix<NekDouble> > mat =
157 m_util->GetInterpolationMatrix(xi);
158 Vmath::Vcopy(mat->GetRows() * mat->GetColumns(), mat->GetRawPtr(),
176 returnval->Initialize();
186 const int nVerts = 3;
188 const int nBoundaryPoints = 3*nEdgeInteriorPoints + 3;
190 if(nEdgeInteriorPoints==0)
197 for(i = cnt = istart; i < nBoundaryPoints; i++)
208 iend = istart + nEdgeInteriorPoints;
209 for(i = istart; i < iend; i++)
211 for(j = istart+1; j < iend; j++)
223 for(i = cnt = istart; i < nBoundaryPoints; i++)
234 iend = istart + nEdgeInteriorPoints;
235 for(i = istart; i < iend; i++)
237 for(j = istart+1; j < iend; j++)
249 for(i = cnt = istart; i < nBoundaryPoints; i++)
259 iend = istart + nEdgeInteriorPoints;
260 for(i = istart; i < iend; i++)
262 for(j = istart+1; j < iend; j++)
301 for(i = 0; i < xc.num_elements() - 1; i++)
305 std::swap(xc[i],xc[i+1]);
306 std::swap(yc[i],yc[i+1]);
320 for(i = offset; i < offset + npl - 1; i++)
324 std::swap(xc[i],xc[i+1]);
325 std::swap(yc[i],yc[i+1]);
boost::shared_ptr< NodalUtilTriangle > m_util
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
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]
size_type num_elements() const
Returns the array's size.
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
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
void CalculateDerivMatrix()
void NodalPointReorder2d()
static boost::shared_ptr< PointsBaseType > Create(const PointsKey &key)
static const unsigned int perm6_2d[6][3]
static const NekDouble NodalTriFeketeData[][6]
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Array< OneD, DataType > & GetPtr()
virtual void CalculateWeights()