45 namespace LibUtilities
62 for(
unsigned int i=0; i < numPoints-2; ++i)
87 for(
unsigned int j=0; j < 4; ++j)
105 for(
unsigned int j=0; j < 6; ++j)
123 for(
unsigned int j=0; j < 12; ++j)
141 for(
unsigned int j=0; j < 12; ++j)
159 for(
unsigned int j=0; j < 12; ++j)
176 for(
unsigned int j=0; j < 24; ++j)
224 boost::shared_ptr<NekMatrix<NekDouble> > mat =
225 m_util->GetInterpolationMatrix(xi);
226 Vmath::Vcopy(mat->GetRows() * mat->GetColumns(), mat->GetRawPtr(),
243 returnval->Initialize();
252 const int nVerts = 4;
257 if(nEdgeInteriorPoints==0)
264 for(
int i = cnt = istart; i < nAllPoints; i++)
277 iend = istart + nEdgeInteriorPoints;
278 for(
int i = istart; i < iend; i++)
280 for(
int j = istart+1; j < iend; j++)
293 for(
int i = cnt = istart; i < nAllPoints; i++)
306 iend = istart + nEdgeInteriorPoints;
307 for(
int i = istart; i < iend; i++)
309 for(
int j = istart+1; j < iend; j++)
322 for(
int i = cnt = istart; i < nAllPoints; i++)
335 iend = istart + nEdgeInteriorPoints;
336 for(
int i = istart; i < iend; i++)
338 for(
int j = istart+1; j < iend; j++)
351 for(
int i = cnt = istart; i < nAllPoints; i++)
364 iend = istart + nEdgeInteriorPoints;
365 for(
int i = istart; i < iend; i++)
367 for(
int j = istart+1; j < iend; j++)
380 for(
int i = cnt = istart; i < nAllPoints; i++)
393 iend = istart + nEdgeInteriorPoints;
394 for(
int i = istart; i < iend; i++)
396 for(
int j = istart+1; j < iend; j++)
409 for(
int i = cnt = istart; i < nAllPoints; i++)
422 iend = istart + nEdgeInteriorPoints;
423 for(
int i = istart; i < iend; i++)
425 for(
int j = istart+1; j < iend; j++)
444 for(
int i = cnt = istart; i < nAllPoints; i++)
456 iend = istart + nFaceInteriorPoints;
461 for(
int i = istart; i < iend - 1; i++)
480 for(
int i = offset+istart; i < offset+istart + npl - 1; i++)
497 for(
int i = cnt = istart; i < nAllPoints; i++)
509 iend = istart + nFaceInteriorPoints;
514 for(
int i = istart; i < iend - 1; i++)
533 for(
int i = offset+istart; i < offset+istart + npl - 1; i++)
550 for(
int i = cnt = istart; i < nAllPoints; i++)
562 iend = istart + nFaceInteriorPoints;
567 for(
int i = istart; i < iend - 1; i++)
586 for(
int i = offset+istart; i < offset+istart + npl - 1; i++)
603 for(
int i = cnt = istart; i < nAllPoints; i++)
615 iend = istart + nFaceInteriorPoints;
620 for(
int i = istart; i < iend - 1; i++)
639 for(
int i = offset+istart; i < offset+istart + npl - 1; i++)
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
virtual void CalculatePoints()
General purpose memory allocation routines with the ability to allocate from thread specific memory p...
void CalculateInterpMatrix(const Array< OneD, const NekDouble > &xia, const Array< OneD, const NekDouble > &yia, const Array< OneD, const NekDouble > &zia, Array< OneD, NekDouble > &interp)
MatrixSharedPtrType m_derivmatrix[3]
Array< OneD, DataType > m_points[3]
void CalculateDerivMatrix()
static const unsigned int perm12C_3d[12][4]
static const NekDouble kNekZeroTol
static const unsigned int perm6_3d[6][4]
unsigned int GetNumPoints() const
unsigned int GetTotNumPoints() const
Array< OneD, DataType > m_weights
Defines a specification for a set of points.
virtual void CalculateDerivMatrix()
unsigned int GetRows() const
static const unsigned int perm4_3d[4][4]
static const unsigned int NodalTetElecNPTS[NodalTetElecAvailable]
static const unsigned int perm24_3d[24][4]
static const NekDouble NodalTetElecData[][9]
#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)
static boost::shared_ptr< PointsBaseType > Create(const PointsKey &key)
static const unsigned int perm12A_3d[12][4]
void NodalPointReorder3d()
Array< OneD, DataType > & GetPtr()
boost::shared_ptr< NodalUtilTetrahedron > m_util
static const unsigned int perm12B_3d[12][4]
virtual void CalculateWeights()