continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp File Reference

#include <MoFEM.hpp>

Go to the source code of this file.

Classes
struct	CommonData
struct	SkeletonFE
struct	SkeletonFE::OpFaceSide

Typedefs
using	FaceEleOnSide = MoFEM::FaceElementForcesAndSourcesCoreOnSide
using	EdgeEle = MoFEM::EdgeElementForcesAndSourcesCore
using	FaceEleOnSideOp = FaceEleOnSide::UserDataOperator
using	EdgeEleOp = EdgeEle::UserDataOperator

Functions
int	main (int argc, char *argv[])

Variables
static char	help [] = "...\n\n"

Typedef Documentation

EdgeEle

using EdgeEle = MoFEM::EdgeElementForcesAndSourcesCore

Definition at line 18 of file continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp.

EdgeEleOp

using EdgeEleOp = EdgeEle::UserDataOperator

Definition at line 20 of file continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp.

FaceEleOnSide

using FaceEleOnSide = MoFEM::FaceElementForcesAndSourcesCoreOnSide

Definition at line 17 of file continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp.

FaceEleOnSideOp

using FaceEleOnSideOp = FaceEleOnSide::UserDataOperator

Definition at line 19 of file continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp.

Function Documentation

main()

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 152 of file continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp.

                                 {
 
  // initialize petsc
  MoFEM::Core::Initialize(&argc, &argv, (char *)0, help);
 
  try {
 
    // Create MoAB database
    moab::Core moab_core;
    moab::Interface &moab = moab_core;
 
    // Create MoFEM database and link it to MoAB
    MoFEM::Core mofem_core(moab);
    MoFEM::Interface &m_field = mofem_core;
 
    // Register DM Manager
    DMType dm_name = "DMMOFEM";
    CHKERR DMRegister_MoFEM(dm_name);
 
    // Simple interface
    Simple *simple_interface;
    CHKERR m_field.getInterface(simple_interface);
    {
      // get options from command line
      CHKERR simple_interface->getOptions();
      // load mesh file
      CHKERR simple_interface->loadFile("");
 
      auto get_base = []() -> FieldApproximationBase {
        enum bases { AINSWORTH, DEMKOWICZ, LASTBASEOP };
        const char *list_bases[] = {"ainsworth", "demkowicz"};
        PetscBool flg;
        PetscInt choice_base_value = AINSWORTH;
        CHKERR PetscOptionsGetEList(PETSC_NULLPTR, NULL, "-base", list_bases,
                                    LASTBASEOP, &choice_base_value, &flg);
        if (flg == PETSC_TRUE) {
          FieldApproximationBase base = AINSWORTH_LEGENDRE_BASE;
          if (choice_base_value == AINSWORTH)
            base = AINSWORTH_LEGENDRE_BASE;
          else if (choice_base_value == DEMKOWICZ)
            base = DEMKOWICZ_JACOBI_BASE;
          return base;
        }
        return LASTBASE;
      };
 
      // add fields
      auto base = get_base();
      CHKERR simple_interface->addDomainField("FIELD", HCURL, base, 1);
      CHKERR simple_interface->addSkeletonField("FIELD", HCURL, base, 1);
      // set fields order
      CHKERR simple_interface->setFieldOrder("FIELD", 2);
      // setup problem
      CHKERR simple_interface->setUp();
      // get dm
      auto dm = simple_interface->getDM();
 
      // create elements
      CommonData elem_data;
      boost::shared_ptr<EdgeEle> skeleton_fe =
          boost::shared_ptr<EdgeEle>(new EdgeEle(m_field));
 
      skeleton_fe->getOpPtrVector().push_back(
          new OpSetContravariantPiolaTransformOnEdge2D());
      skeleton_fe->getOpPtrVector().push_back(
          new SkeletonFE(m_field, elem_data));
 
      // iterate skeleton finite elements
      CHKERR DMoFEMLoopFiniteElements(dm, simple_interface->getSkeletonFEName(),
                                      skeleton_fe);
    }
  }
  CATCH_ERRORS;
 
  // finish work cleaning memory, getting statistics, etc.
  MoFEM::Core::Finalize();
 
  return 0;
}

Variable Documentation

help

char help[] = "...\n\n"

static

Definition at line 15 of file continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp.