Nektar::LibUtilities::StdPyrData 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::StdPyrData::getNumberOfBndCoefficients ( int  Na, int  Nb, int  Nc  )

inline

Definition at line 267 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()

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

inline

Definition at line 240 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::FieldIO::CheckFieldDefinition(), Nektar::LibUtilities::GetNumberOfCoefficients(), Nektar::FieldUtils::ProcessEquiSpacedOutput::Process(), and Nektar::StdRegions::StdPyrExp::v_CalcNumberOfCoefficients().