hcurl_divergence_operator_2d.cpp File Reference

#include <MoFEM.hpp>

Go to the source code of this file.

Classes
struct	OpDivergence
struct	OpFlux

Typedefs
using	FaceEle = MoFEM::FaceElementForcesAndSourcesCore
using	EdgeEle = MoFEM::EdgeElementForcesAndSourcesCore
using	FaceEleOp = FaceEle::UserDataOperator
using	EdgeEleOp = EdgeEle::UserDataOperator

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

Variables
static char	help [] = "...\n\n"
constexpr int	SPACE_DIM = 2
FTensor::Index< 'i', SPACE_DIM >	i

Typedef Documentation

EdgeEle

using EdgeEle = MoFEM::EdgeElementForcesAndSourcesCore

Definition at line 21 of file hcurl_divergence_operator_2d.cpp.

EdgeEleOp

using EdgeEleOp = EdgeEle::UserDataOperator

Definition at line 24 of file hcurl_divergence_operator_2d.cpp.

FaceEle

using FaceEle = MoFEM::FaceElementForcesAndSourcesCore

Definition at line 20 of file hcurl_divergence_operator_2d.cpp.

FaceEleOp

using FaceEleOp = FaceEle::UserDataOperator

Definition at line 23 of file hcurl_divergence_operator_2d.cpp.

Function Documentation

main()

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

Definition at line 83 of file hcurl_divergence_operator_2d.cpp.

                                 {
 
  MoFEM::Core::Initialize(&argc, &argv, (char *)0, help);
 
  try {
 
    moab::Core mb_instance;
    moab::Interface &moab = mb_instance;
 
    PetscBool flg_file = PETSC_TRUE;
    char mesh_file_name[255];
    CHKERR PetscOptionsGetString(PETSC_NULL, "", "-my_file", mesh_file_name,
                                 255, &flg_file);
    if (flg_file != PETSC_TRUE)
      SETERRQ(PETSC_COMM_SELF, MOFEM_INVALID_DATA,
              "*** ERROR -my_file (MESH FILE NEEDED)");
 
    // Read mesh to MOAB
    CHKERR moab.load_file(mesh_file_name, 0, "");
 
    // Create MoFEM instance
    MoFEM::Core core(moab);
    MoFEM::Interface &m_field = core;
 
    // set entities bit level
    BitRefLevel bit_level0 = BitRefLevel().set(0);
    CHKERR m_field.getInterface<BitRefManager>()->setBitRefLevelByDim(
        0, 2, bit_level0);
 
    // Declare elements
    enum bases { AINSWORTH, DEMKOWICZ, LASBASETOP };
    const char *list_bases[] = {"ainsworth", "demkowicz"};
    PetscBool flg;
    PetscInt choice_base_value = AINSWORTH;
    CHKERR PetscOptionsGetEList(PETSC_NULL, NULL, "-base", list_bases,
                                LASBASETOP, &choice_base_value, &flg);
    if (flg != PETSC_TRUE)
      SETERRQ(PETSC_COMM_SELF, MOFEM_IMPOSSIBLE_CASE, "base not set");
    FieldApproximationBase base = AINSWORTH_LEGENDRE_BASE;
    if (choice_base_value == AINSWORTH)
      base = AINSWORTH_LEGENDRE_BASE;
    else if (choice_base_value == DEMKOWICZ)
      base = DEMKOWICZ_JACOBI_BASE;
    int order = 5;
    CHKERR PetscOptionsGetInt(PETSC_NULL, "", "-order", &order, PETSC_NULL);
 
    CHKERR m_field.add_field("FIELD1", HCURL, base, 1);
    CHKERR m_field.add_finite_element("FACE_FE");
    // Define rows/cols and element data
    CHKERR m_field.modify_finite_element_add_field_row("FACE_FE", "FIELD1");
    CHKERR m_field.modify_finite_element_add_field_col("FACE_FE", "FIELD1");
    CHKERR m_field.modify_finite_element_add_field_data("FACE_FE", "FIELD1");
 
    CHKERR m_field.add_finite_element("EDGE_FE");
    // Define rows/cols and element data
    CHKERR m_field.modify_finite_element_add_field_row("EDGE_FE", "FIELD1");
    CHKERR m_field.modify_finite_element_add_field_col("EDGE_FE", "FIELD1");
    CHKERR m_field.modify_finite_element_add_field_data("EDGE_FE", "FIELD1");
 
    // Problem
    CHKERR m_field.add_problem("TEST_PROBLEM");
    // set finite elements for problem
    CHKERR m_field.modify_problem_add_finite_element("TEST_PROBLEM", "FACE_FE");
    CHKERR m_field.modify_problem_add_finite_element("TEST_PROBLEM", "EDGE_FE");
    // set refinement level for problem
    CHKERR m_field.modify_problem_ref_level_add_bit("TEST_PROBLEM", bit_level0);
 
    // Add entities
 
    CHKERR m_field.add_ents_to_field_by_dim(0, SPACE_DIM, "FIELD1");
    // Set order
    for (auto t : {MBEDGE, MBTRI, MBQUAD})
      CHKERR m_field.set_field_order(0, t, "FIELD1", order);
 
    // Add entities to elements
    CHKERR m_field.add_ents_to_finite_element_by_dim(0, SPACE_DIM, "FACE_FE");
    auto set_edge_elements_entities_on_mesh_skin = [&]() {
      MoFEMFunctionBegin;
      Range faces;
      CHKERR moab.get_entities_by_dimension(0, 2, faces, false);
      Skinner skin(&m_field.get_moab());
      Range faces_skin;
      CHKERR skin.find_skin(0, faces, false, faces_skin);
      CHKERR m_field.add_ents_to_finite_element_by_dim(
          faces_skin, SPACE_DIM - 1, "EDGE_FE");
      MoFEMFunctionReturn(0);
    };
    CHKERR set_edge_elements_entities_on_mesh_skin();
 
    // Build database
    CHKERR m_field.build_fields();
    // build finite elemnts
    CHKERR m_field.build_finite_elements();
    // build adjacencies
    CHKERR m_field.build_adjacencies(bit_level0);
 
    // build problem
    ProblemsManager *prb_mng_ptr;
    CHKERR m_field.getInterface(prb_mng_ptr);
    CHKERR prb_mng_ptr->buildProblem("TEST_PROBLEM", true);
    // Partition
    CHKERR prb_mng_ptr->partitionSimpleProblem("TEST_PROBLEM");
    CHKERR prb_mng_ptr->partitionFiniteElements("TEST_PROBLEM");
    CHKERR prb_mng_ptr->partitionGhostDofs("TEST_PROBLEM");
 
    // integration rule
    auto rule = [&](int, int, int p) { return 2 * p; };
 
    auto calculate_divergence = [&]() {
      double div = 0;
      FaceEle fe_face(m_field);
      fe_face.getRuleHook = rule;
      CHKERR AddHOOps<2, 2, 2>::add(fe_face.getOpPtrVector(), {HDIV});
      fe_face.getOpPtrVector().push_back(new OpDivergence(div));
      CHKERR m_field.loop_finite_elements("TEST_PROBLEM", "FACE_FE", fe_face);
      return div;
    };
 
    auto calculate_flux = [&]() {
      double flux = 0;
      EdgeEle fe_edge(m_field);
      fe_edge.getRuleHook = rule;
      CHKERR AddHOOps<1, 2, 2>::add(fe_edge.getOpPtrVector(), {HDIV});
      fe_edge.getOpPtrVector().push_back(new OpFlux(flux));
      CHKERR m_field.loop_finite_elements("TEST_PROBLEM", "EDGE_FE", fe_edge);
      return flux;
    };
 
    const double div = calculate_divergence();
    const double flux = calculate_flux();
 
    MOFEM_LOG_CHANNEL("WOLD"); // reset channel
    MOFEM_LOG_C("WORLD", Sev::inform,
                "Div = %4.3e Flux = %3.4e Error = %4.3e\n", div, flux,
                div - flux);
 
    constexpr double tol = 1e-8;
    if (std::abs(div - flux) > tol)
      SETERRQ2(PETSC_COMM_WORLD, MOFEM_ATOM_TEST_INVALID,
               "Test failed (div != flux) %3.4e != %3.4e", div, flux);
  }
  CATCH_ERRORS;
 
  CHKERR MoFEM::Core::Finalize();
}

Variable Documentation

help

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

static

Definition at line 18 of file hcurl_divergence_operator_2d.cpp.

i

FTensor::Index<'i', SPACE_DIM> i

Examples

ContactOps.hpp, ElasticityMixedFormulation.hpp, EshelbianOperators.cpp, HookeElement.cpp, HookeElement.hpp, HookeInternalStressElement.hpp, MagneticElement.hpp, NavierStokesElement.cpp, PlasticOps.hpp, PlasticOpsGeneric.hpp, PlasticOpsLargeStrains.hpp, PlasticOpsSmallStrains.hpp, Remodeling.cpp, ThermoElasticOps.hpp, UnsaturatedFlow.hpp, analytical_poisson_field_split.cpp, approx_sphere.cpp, bernstein_bezier_generate_base.cpp, cell_forces.cpp, child_and_parent.cpp, continuity_check_on_skeleton_3d.cpp, continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, dynamic_first_order_con_law.cpp, edge_and_bubble_shape_functions_on_quad.cpp, elasticity.cpp, field_blas.cpp, free_surface.cpp, gauss_points_on_outer_product.cpp, hanging_node_approx.cpp, hcurl_curl_operator.cpp, hcurl_divergence_operator_2d.cpp, hdiv_divergence_operator.cpp, heat_equation.cpp, higher_derivatives.cpp, level_set.cpp, lorentz_force.cpp, matrix_function.cpp, mixed_poisson.cpp, mortar_contact_thermal.cpp, navier_stokes.cpp, phase.cpp, plastic.cpp, plate.cpp, plot_base.cpp, poisson_2d_dis_galerkin.cpp, prism_elements_from_surface.cpp, prism_polynomial_approximation.cpp, quad_polynomial_approximation.cpp, reaction_diffusion.cpp, scalar_check_approximation.cpp, shallow_wave.cpp, simple_contact.cpp, simple_contact_thermal.cpp, simple_elasticity.cpp, simple_interface.cpp, tensor_divergence_operator.cpp, testing_jacobian_of_hook_scaled_with_density_element.cpp, thermo_elastic.cpp, and wave_equation.cpp.

Definition at line 27 of file hcurl_divergence_operator_2d.cpp.

SPACE_DIM

constexpr int SPACE_DIM = 2

constexpr

Definition at line 26 of file hcurl_divergence_operator_2d.cpp.