continuity_check_on_skeleton_with_simple_2d_for_h1.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

Examples

continuity_check_on_skeleton_with_simple_2d_for_h1.cpp.

Definition at line 18 of file continuity_check_on_skeleton_with_simple_2d_for_h1.cpp.

EdgeEleOp

using EdgeEleOp = EdgeEle::UserDataOperator

Definition at line 20 of file continuity_check_on_skeleton_with_simple_2d_for_h1.cpp.

FaceEleOnSide

using FaceEleOnSide = MoFEM::FaceElementForcesAndSourcesCoreOnSide

Definition at line 17 of file continuity_check_on_skeleton_with_simple_2d_for_h1.cpp.

FaceEleOnSideOp

using FaceEleOnSideOp = FaceEleOnSide::UserDataOperator

Definition at line 19 of file continuity_check_on_skeleton_with_simple_2d_for_h1.cpp.

Function Documentation

main()

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

Examples

continuity_check_on_skeleton_with_simple_2d_for_h1.cpp.

Definition at line 151 of file continuity_check_on_skeleton_with_simple_2d_for_h1.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;
 
    auto core_log = logging::core::get();
    core_log->add_sink(
        LogManager::createSink(LogManager::getStrmSelf(), "ATOM"));
    LogManager::setLog("ATOM");
 
    // 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();
      // Add ghost cells needed for evaluation on skeleton
      simple_interface->getAddSkeletonFE() = true;
 
      // 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_NULL, 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", H1, base, 1);
      // CHKERR simple_interface->addDomainField("TEST_FIELD", L2,
      //                                         AINSWORTH_LEGENDRE_BASE, 1);
      CHKERR simple_interface->addSkeletonField("FIELD", H1, base, 1);
      // set fields order
      CHKERR simple_interface->setFieldOrder("FIELD", 2);
      // CHKERR simple_interface->setFieldOrder("TEST_FIELD", 1);
      // 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 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

Examples

continuity_check_on_skeleton_with_simple_2d_for_h1.cpp.

Definition at line 15 of file continuity_check_on_skeleton_with_simple_2d_for_h1.cpp.