Nektar++
Functions
Nektar::LibUtilities::StdTetData Namespace Reference

Functions

int getNumberOfCoefficients (int Na, int Nb, int Nc)
 
int getNumberOfBndCoefficients (int Na, int Nb, int Nc)
 

Function Documentation

int Nektar::LibUtilities::StdTetData::getNumberOfBndCoefficients ( int  Na,
int  Nb,
int  Nc 
)
inline

Definition at line 209 of file ShapeType.hpp.

References ASSERTL1, and ASSERTL2.

Referenced by Nektar::LibUtilities::MeshPartition::CalculateElementWeight(), and Nektar::StdRegions::StdTetExp::v_NumBndryCoeffs().

210  {
211  ASSERTL2(Na > 1, "Order in 'a' direction must be > 1.");
212  ASSERTL2(Nb > 1, "Order in 'b' direction must be > 1.");
213  ASSERTL2(Nc > 1, "Order in 'c' direction must be > 1.");
214  ASSERTL1(Na <= Nc, "order in 'a' direction is higher "
215  "than order in 'c' direction");
216  ASSERTL1(Nb <= Nc, "order in 'b' direction is higher "
217  "than order in 'c' direction");
218 
219  int nCoef = Na*(Na+1)/2 + (Nb-Na)*Na // base
220  + Na*(Na+1)/2 + (Nc-Na)*Na // front
221  + 2*(Nb*(Nb+1)/2 + (Nc-Nb)*Nb)// 2 other sides
222  - Na - 2*Nb - 3*Nc // less edges
223  + 4; // plus vertices
224 
225  return nCoef;
226  }
#define ASSERTL2(condition, msg)
Assert Level 2 – Debugging which is used FULLDEBUG compilation mode. This level assert is designed t...
Definition: ErrorUtil.hpp:213
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:191
int Nektar::LibUtilities::StdTetData::getNumberOfCoefficients ( int  Na,
int  Nb,
int  Nc 
)
inline

Adds up the number of cells in a truncated Nc by Nc by Nc pyramid, where the longest Na rows and longest Nb columns are kept. Example: (Na, Nb, Nc) = (3, 4, 5); The number of coefficients is the sum of the elements of the following matrix:

|5 4 3 2 0| |4 3 2 0 | |3 2 0 | |0 0 | |0 |

Sum = 28 = number of tet coefficients.

Definition at line 186 of file ShapeType.hpp.

References ASSERTL1, and ASSERTL2.

Referenced by Nektar::LibUtilities::MeshPartition::CalculateElementWeight(), Nektar::LibUtilities::FieldIO::CheckFieldDefinition(), Nektar::MultiRegions::ExpList3D::ExpList3D(), Nektar::LibUtilities::GetNumberOfCoefficients(), Nektar::Utilities::ProcessEquiSpacedOutput::SetupEquiSpacedField(), Nektar::StdRegions::StdTetExp::v_CalcNumberOfCoefficients(), and Nektar::StdRegions::StdTetExp::v_GenMatrix().

187  {
188  ASSERTL2(Na > 1, "Order in 'a' direction must be > 1.");
189  ASSERTL2(Nb > 1, "Order in 'b' direction must be > 1.");
190  ASSERTL2(Nc > 1, "Order in 'c' direction must be > 1.");
191  ASSERTL1(Na <= Nc, "order in 'a' direction is higher "
192  "than order in 'c' direction");
193  ASSERTL1(Nb <= Nc, "order in 'b' direction is higher "
194  "than order in 'c' direction");
195  int nCoef = 0;
196  for (int a = 0; a < Na; ++a)
197  {
198  for (int b = 0; b < Nb - a; ++b)
199  {
200  for (int c = 0; c < Nc - a - b; ++c)
201  {
202  ++nCoef;
203  }
204  }
205  }
206  return nCoef;
207  }
#define ASSERTL2(condition, msg)
Assert Level 2 – Debugging which is used FULLDEBUG compilation mode. This level assert is designed t...
Definition: ErrorUtil.hpp:213
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:191