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 "MFEMFESpace.h"
      13             : #include "MFEMProblem.h"
      14             : 
      15             : InputParameters
      16       11462 : MFEMFESpace::validParams()
      17             : {
      18       11462 :   InputParameters params = MFEMObject::validParams();
      19       22924 :   params.registerBase("MFEMFESpace");
      20       22924 :   params.registerSystemAttributeName("MFEMFESpace");
      21       45848 :   MooseEnum ordering("NODES VDIM", "VDIM", false);
      22       45848 :   params.addParam<MooseEnum>("ordering", ordering, "Ordering style to use for vector DoFs.");
      23       45848 :   params.addParam<int>("vdim", 1, "The number of degrees of freedom per basis function.");
      24       34386 :   params.addParam<std::string>("submesh",
      25             :                                "Submesh to define the FESpace on. Leave blank to use base mesh.");
      26       22924 :   return params;
      27       11462 : }
      28             : 
      29        2572 : MFEMFESpace::MFEMFESpace(const InputParameters & parameters)
      30             :   : MFEMObject(parameters),
      31        2572 :     _ordering(parameters.get<MooseEnum>("ordering")),
      32        2572 :     _pmesh(
      33        2572 :         parameters.isParamValid("submesh")
      34        5411 :             ? getMFEMProblem().getProblemData().submeshes.GetRef(getParam<std::string>("submesh"))
      35        7449 :             : const_cast<mfem::ParMesh &>(getMFEMProblem().mesh().getMFEMParMesh()))
      36             : {
      37        2572 : }
      38             : 
      39             : void
      40        2572 : MFEMFESpace::buildFEC() const
      41             : {
      42        5112 :   _fec = std::shared_ptr<mfem::FiniteElementCollection>(
      43        7684 :       mfem::FiniteElementCollection::New(getFECName().c_str()));
      44        2556 : }
      45             : 
      46             : void
      47        2572 : MFEMFESpace::buildFESpace() const
      48             : {
      49             :   _fespace =
      50        2572 :       std::make_shared<mfem::ParFiniteElementSpace>(&_pmesh, getFEC().get(), getVDim(), _ordering);
      51        2556 : }
      52             : 
      53             : #endif