Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #ifdef MOOSE_MFEM_ENABLED
      11             : 
      12             : #include "MFEMScalarFESpace.h"
      13             : 
      14             : registerMooseObject("MooseApp", MFEMScalarFESpace);
      15             : 
      16             : InputParameters
      17        9872 : MFEMScalarFESpace::validParams()
      18             : {
      19        9872 :   InputParameters params = MFEMSimplifiedFESpace::validParams();
      20       19744 :   params.addClassDescription("Convenience class to construct scalar finite element spaces.");
      21       39488 :   MooseEnum fec_types("H1 L2", "H1", true);
      22       39488 :   params.addParam<MooseEnum>("fec_type", fec_types, "Specifies the family of FE shape functions.");
      23       29616 :   params.addParam<std::string>(
      24             :       "basis", "GaussLobatto", "Specifies the quadrature basis used for scalar elements.");
      25             : 
      26       19744 :   return params;
      27        9872 : }
      28             : 
      29         381 : MFEMScalarFESpace::MFEMScalarFESpace(const InputParameters & parameters)
      30         381 :   : MFEMSimplifiedFESpace(parameters), _fec_type(parameters.get<MooseEnum>("fec_type"))
      31             : {
      32         381 : }
      33             : 
      34             : std::string
      35         381 : MFEMScalarFESpace::getFECName() const
      36             : {
      37        1143 :   const char b = mfem::BasisType::GetChar(getBasis(getParam<std::string>("basis")));
      38         381 :   std::string basis(1, b);
      39             : 
      40             :   // This is to get around an MFEM bug where if you pass the full name of the default element type,
      41             :   // it crashes
      42         381 :   if (_fec_type == "H1")
      43         984 :     basis = (basis == "G" ? "" : "@" + basis);
      44          53 :   else if (_fec_type == "L2")
      45             :   {
      46          53 :     if (basis != "g")
      47         212 :       mooseInfo("L2 finite element space only supports GaussLegendre basis. Ignoring " +
      48         159 :                 getParam<std::string>("basis") +
      49             :                 " basis choice and using GaussLegendre instead.\n");
      50             : 
      51          53 :     basis = "";
      52             :   }
      53             : 
      54         762 :   return _fec_type + basis + "_" + std::to_string(getProblemDim()) + "D_P" +
      55        1524 :          std::to_string(_fec_order);
      56         381 : }
      57             : 
      58             : int
      59         381 : MFEMScalarFESpace::getVDim() const
      60             : {
      61         381 :   return 1;
      62             : }
      63             : 
      64             : #endif