simple_l2_only.cpp File Reference

#include <MoFEM.hpp>

Go to the source code of this file.

Classes
struct	ElementsAndOps< DIM >
struct	ElementsAndOps< 2 >

Typedefs
using	DomainEle = ElementsAndOps< SPACE_DIM >::DomainEle
using	DomainEleOp = DomainEle::UserDataOperator
using	BoundaryEle = ElementsAndOps< SPACE_DIM >::BoundaryEle
using	BoundaryEleOp = BoundaryEle::UserDataOperator
using	EntData = EntitiesFieldData::EntData

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

Variables
static char	help [] = "...\n\n"
constexpr int	SPACE_DIM = 2

Typedef Documentation

BoundaryEle

using BoundaryEle = ElementsAndOps<SPACE_DIM>::BoundaryEle

Definition at line 31 of file simple_l2_only.cpp.

BoundaryEleOp

using BoundaryEleOp = BoundaryEle::UserDataOperator

Examples

simple_l2_only.cpp.

Definition at line 32 of file simple_l2_only.cpp.

DomainEle

using DomainEle = ElementsAndOps<SPACE_DIM>::DomainEle

Definition at line 29 of file simple_l2_only.cpp.

DomainEleOp

using DomainEleOp = DomainEle::UserDataOperator

Examples

simple_l2_only.cpp.

Definition at line 30 of file simple_l2_only.cpp.

EntData

using EntData = EntitiesFieldData::EntData

Examples

simple_l2_only.cpp.

Definition at line 33 of file simple_l2_only.cpp.

Function Documentation

main()

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

Examples

simple_l2_only.cpp.

Definition at line 35 of file simple_l2_only.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();
 
      CHKERR simple_interface->addDomainField("FIELD", L2,
                                              AINSWORTH_LEGENDRE_BASE, 1);
 
      simple_interface->getAddBoundaryFE() = true;
      simple_interface->getAddSkeletonFE() = true;
 
      // set fields order
      CHKERR simple_interface->setFieldOrder("FIELD", 1);
      // setup problem
      CHKERR simple_interface->setUp();
 
      int count_fe;
      int count_side_fe;
 
      PipelineManager *pipeline_mng = m_field.getInterface<PipelineManager>();
 
      // Create OP for side FE
      auto op_side_fe = new DomainEleOp(NOSPACE, DomainEleOp::OPSPACE);
      op_side_fe->doWorkRhsHook = [&](DataOperator *op_ptr, int side,
                                      EntityType type,
                                      EntitiesFieldData::EntData &data) {
        auto domain_op = static_cast<DomainEleOp *>(op_ptr);
        MoFEMFunctionBegin;
 
        MOFEM_LOG("SELF", Sev::verbose)
            << "Side element name [ " << count_side_fe << " ] "
            << domain_op->getFEName();
 
        ++count_side_fe;
 
        MoFEMFunctionReturn(0);
      };
 
      // Create side FE
      auto side_fe = boost::make_shared<DomainEle>(m_field);
      side_fe->getOpPtrVector().push_back(op_side_fe);
 
      // Create boundary FE OP
 
      auto do_work_rhs = [&](DataOperator *op_ptr, int side, EntityType type,
                             EntitiesFieldData::EntData &data) {
        auto bdy_op = static_cast<BoundaryEleOp *>(op_ptr);
        MoFEMFunctionBegin;
 
        MOFEM_LOG("SELF", Sev::verbose) << "Element name  [ " << count_fe
                                        << " ] " << bdy_op->getFEName();
 
        CHKERR bdy_op->loopSide(simple_interface->getDomainFEName(),
                                side_fe.get(), SPACE_DIM);
 
        ++count_fe;
 
        MoFEMFunctionReturn(0);
      };
 
      auto op_bdy_fe = new BoundaryEleOp(NOSPACE, DomainEleOp::OPSPACE);
      op_bdy_fe->doWorkRhsHook = do_work_rhs;
 
      auto op_skeleton_fe = new BoundaryEleOp(NOSPACE, DomainEleOp::OPSPACE);
      op_skeleton_fe->doWorkRhsHook = do_work_rhs;
 
      // Count boundary
      count_fe = 0;
      count_side_fe = 0;
 
      pipeline_mng->getOpBoundaryRhsPipeline().push_back(op_bdy_fe);
      pipeline_mng->getOpSkeletonRhsPipeline().push_back(op_skeleton_fe);
      pipeline_mng->loopFiniteElements();
 
      MOFEM_LOG("SELF", Sev::inform) << "Number of elements " << count_fe;
      MOFEM_LOG("SELF", Sev::inform)
          << "Number of side elements " << count_side_fe;
 
      if (count_fe != 16)
        SETERRQ(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
                "Wrong numbers of FEs");
      if (count_side_fe != 24)
        SETERRQ(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
                "Wrong numbers of side FEs");
    }
  }
  CATCH_ERRORS;
 
  // finish work cleaning memory, getting statistics, etc.
  MoFEM::Core::Finalize();
 
  return 0;
}

Variable Documentation

help

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

static

Examples

simple_l2_only.cpp.

Definition at line 16 of file simple_l2_only.cpp.

SPACE_DIM

constexpr int SPACE_DIM = 2

constexpr

Examples

simple_l2_only.cpp.

Definition at line 18 of file simple_l2_only.cpp.