Nektar::LibUtilities::StdPyrData Namespace Reference

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

Function Documentation

getNumberOfBndCoefficients()

constexpr int Nektar::LibUtilities::StdPyrData::getNumberOfBndCoefficients ( int  Na, int  Nb, int  Nc  )

inlineconstexpr

Definition at line 318 of file ShapeType.hpp.

{
    ASSERTL1(Na > 1, "Order in 'a' direction must be > 1.");
    ASSERTL1(Nb > 1, "Order in 'b' direction must be > 1.");
    ASSERTL1(Nc > 1, "Order in 'c' direction must be > 1.");
    ASSERTL1(Na <= Nc, "Order in 'a' direction is higher "
                       "than order in 'c' direction.");
    ASSERTL1(Nb <= Nc, "Order in 'b' direction is higher "
                       "than order in 'c' direction.");
 
    return Na * Nb                                    // base
           + 2 * (Na * (Na + 1) / 2 + (Nc - Na) * Na) // front and back
           + 2 * (Nb * (Nb + 1) / 2 + (Nc - Nb) * Nb) // sides
           - 2 * Na - 2 * Nb - 4 * Nc                 // less edges
           + 5;                                       // plus vertices
}

References ASSERTL1.

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

getNumberOfCoefficients()

constexpr int Nektar::LibUtilities::StdPyrData::getNumberOfCoefficients ( int  Na, int  Nb, int  Nc  )

inlineconstexpr

Definition at line 291 of file ShapeType.hpp.

{
    ASSERTL1(Na > 1, "Order in 'a' direction must be > 1.");
    ASSERTL1(Nb > 1, "Order in 'b' direction must be > 1.");
    ASSERTL1(Nc > 1, "Order in 'c' direction must be > 1.");
    ASSERTL1(Na <= Nc, "Order in 'a' direction is higher "
                       "than order in 'c' direction.");
    ASSERTL1(Nb <= Nc, "Order in 'b' direction is higher "
                       "than order in 'c' direction.");
 
    // Count number of coefficients explicitly.
    int nCoeff = 0;
 
    // Count number of interior tet modes
    for (int a = 0; a < Na; ++a)
    {
        for (int b = 0; b < Nb; ++b)
        {
            for (int c = 0; c < Nc - std::max(a, b); ++c)
            {
                ++nCoeff;
            }
        }
    }
    return nCoeff;
}

References ASSERTL1.

Referenced by Nektar::SpatialDomains::MeshPartition::CalculateElementWeight(), Nektar::LibUtilities::GetNumberOfCoefficients(), Nektar::LibUtilities::GetNumberOfCoefficients(), Nektar::LibUtilities::GetNumberOfDataPoints(), Nektar::Collections::PyrIProduct(), Nektar::StdRegions::StdPyrExp::v_CalcNumberOfCoefficients(), and Nektar::FieldUtils::ProcessEquiSpacedOutput::v_Process().