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

◆ getNumberOfBndCoefficients()

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

Definition at line 213 of file ShapeType.hpp.

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

References ASSERTL1, and ASSERTL2.

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

◆ getNumberOfCoefficients()

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 190 of file ShapeType.hpp.

191 {
192  ASSERTL2(Na > 1, "Order in 'a' direction must be > 1.");
193  ASSERTL2(Nb > 1, "Order in 'b' direction must be > 1.");
194  ASSERTL2(Nc > 1, "Order in 'c' direction must be > 1.");
195  ASSERTL1(Na <= Nc, "order in 'a' direction is higher "
196  "than order in 'c' direction");
197  ASSERTL1(Nb <= Nc, "order in 'b' direction is higher "
198  "than order in 'c' direction");
199  int nCoef = 0;
200  for (int a = 0; a < Na; ++a)
201  {
202  for (int b = 0; b < Nb - a; ++b)
203  {
204  for (int c = 0; c < Nc - a - b; ++c)
205  {
206  ++nCoef;
207  }
208  }
209  }
210  return nCoef;
211 }

References ASSERTL1, and ASSERTL2.

Referenced by Nektar::SpatialDomains::MeshPartition::CalculateElementWeight(), Nektar::LibUtilities::FieldIO::CheckFieldDefinition(), Nektar::LibUtilities::GetNumberOfCoefficients(), Nektar::FieldUtils::ProcessEquiSpacedOutput::Process(), Nektar::StdRegions::StdTetExp::v_CalcNumberOfCoefficients(), and Nektar::StdRegions::StdTetExp::v_GenMatrix().