MFEMDomainSubMesh.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#ifdef MOOSE_MFEM_ENABLED

#include "MFEMDomainSubMesh.h"
#include "MFEMProblem.h"

registerMooseObject("MooseApp", MFEMDomainSubMesh);

InputParameters
MFEMDomainSubMesh::validParams()
{
  InputParameters params = MFEMSubMesh::validParams();
  params += MFEMBlockRestrictable::validParams();
  params.addClassDescription("Class to construct an MFEMSubMesh formed from the subspace of the "
                             "parent mesh restricted to the set of user-specified subdomains.");
  params.addParam<BoundaryName>("submesh_boundary", "Name to assign submesh boundary.");
  return params;
}

MFEMDomainSubMesh::MFEMDomainSubMesh(const InputParameters & parameters)
  : MFEMSubMesh(parameters),
    MFEMBlockRestrictable(parameters, getMFEMProblem().mesh().getMFEMParMesh())
{
}

void
MFEMDomainSubMesh::buildSubMesh()
{
  _submesh = std::make_shared<mfem::ParSubMesh>(
      mfem::ParSubMesh::CreateFromDomain(getMesh(), getSubdomainAttributes()));
  _submesh->attribute_sets.attr_sets = getMesh().attribute_sets.attr_sets;
  _submesh->bdr_attribute_sets.attr_sets = getMesh().bdr_attribute_sets.attr_sets;

  if (isParamSetByUser("submesh_boundary"))
  {
    const BoundaryName & submesh_boundary = getParam<BoundaryName>("submesh_boundary");
    _submesh->bdr_attribute_sets.CreateAttributeSet(submesh_boundary);
    _submesh->bdr_attribute_sets.AddToAttributeSet(submesh_boundary,
                                                   getMesh().bdr_attributes.Max() + 1);
  }
}

#endif