VEM_MCC_3D_ReferenceElement.hpp

Go to the documentation of this file.

// _LICENSE_HEADER_
//
// Copyright (C) 2019 - 2025.
// Terms register on the GPL-3.0 license.
//
// This file can be redistributed and/or modified under the license terms.
//
// See top level LICENSE file for more details.
//
// This file can be used citing references in CITATION.cff file.
 
#ifndef __VEM_MCC_3D_ReferenceElement_HPP
#define __VEM_MCC_3D_ReferenceElement_HPP
 
#include "I_VEM_MCC_3D_ReferenceElement.hpp"
#include "Monomials_2D.hpp"
#include "Monomials_3D.hpp"
#include "VEM_Quadrature_3D.hpp"
 
namespace Polydim
{
namespace VEM
{
namespace MCC
{
 
class VEM_MCC_3D_Pressure_ReferenceElement final : public I_VEM_MCC_3D_Pressure_ReferenceElement
{
  public:
    VEM_MCC_3D_Pressure_ReferenceElement_Data Create(const unsigned int order) const
    {
        Polydim::Utilities::Monomials_3D monomials;
        Polydim::VEM::Quadrature::VEM_Quadrature_3D quadrature;
 
        Polydim::VEM::MCC::VEM_MCC_3D_Pressure_ReferenceElement_Data result;
 
        result.Monomials = monomials.Compute(order);
        result.Quadrature = quadrature.Compute_MCC_3D(order);
 
        result.Dimension = 3;
        result.Order = order;
        result.NumDofs0D = 0;
        result.NumDofs1D = 0;
        result.NumDofs2D = 0;
        result.NumDofs3D = (order + 1) * (order + 2) * (order + 3) / 6;
 
        return result;
    }
};
 
class VEM_MCC_3D_Velocity_ReferenceElement final : public I_VEM_MCC_3D_Velocity_ReferenceElement
{
  public:
    VEM_MCC_3D_Velocity_ReferenceElement_Data Create(const unsigned int order) const
    {
        Polydim::Utilities::Monomials_2D monomials2D;
        Polydim::Utilities::Monomials_3D monomials3D;
        Polydim::VEM::Quadrature::VEM_Quadrature_3D quadrature;
 
        Polydim::VEM::MCC::VEM_MCC_3D_Velocity_ReferenceElement_Data result;
 
        result.Monomials2D = monomials2D.Compute(order);
        result.MonomialsKp1 = monomials3D.Compute(order + 1);
        result.Quadrature = quadrature.Compute_MCC_3D(order);
 
        result.Dimension = 3;
        result.Order = order;
        result.NumDofs0D = 0;
        result.NumDofs1D = 0;
        result.NumDofs2D = (order + 1) * (order + 2) / 2;
        result.NumDofs3D = order * (order + 2) * (order + 3) / 2;
 
        return result;
    }
};
} // namespace MCC
} // namespace VEM
} // namespace Polydim
 
#endif