MFEMCutTransitionSubMesh.h

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
 
#pragma once
 
#include "MFEMSubMesh.h"
#include "MFEMBlockRestrictable.h"
 
class MFEMCutTransitionSubMesh : public MFEMSubMesh, public MFEMBlockRestrictable
{
public:
  static InputParameters validParams();
  MFEMCutTransitionSubMesh(const InputParameters & parameters);
 
protected:
  virtual void buildSubMesh() override;
 
  void labelMesh(mfem::ParMesh & parent_mesh);
 
  void setAttributes(mfem::ParMesh & parent_mesh, mfem::Array<int> & transition_els);
 
  bool isInDomain(const int & el, const mfem::Array<int> & subdomains, const mfem::ParMesh & mesh);
 
  mfem::Vector findFaceNormal(const mfem::ParMesh & mesh, const int & face);
 
  bool isPositiveSideOfCut(const int & el, const int & el_vertex_on_cut, mfem::ParMesh & mesh);
 
  const BoundaryName & _cut_boundary;
  std::shared_ptr<mfem::ParSubMesh> _cut_submesh{nullptr};
  const BoundaryName & _transition_subdomain_boundary;
  const SubdomainName & _transition_subdomain;
  const SubdomainName & _closed_subdomain;
  mfem::Vector _cut_normal;
};
 
#endif