debug_crack_mesh.cpp File Reference

#include <MoFEM.hpp>
#include <cmath>
#include <cholesky.hpp>
#include <EshelbianPlasticity.hpp>

Go to the source code of this file.

Macros
#define	SINGULARITY

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

Variables
constexpr int	adolc_tag = 1
static char	help [] = "...\n\n"

Macro Definition Documentation

SINGULARITY

#define SINGULARITY

Definition at line 27 of file debug_crack_mesh.cpp.

Function Documentation

main()

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

Definition at line 50 of file debug_crack_mesh.cpp.

                                 {
 
  // initialize petsc
  const char param_file[] = "param_file.petsc";
  MoFEM::Core::Initialize(&argc, &argv, param_file, help);
 
  // Add logging channel for example
  auto core_log = logging::core::get();
  core_log->add_sink(LogManager::createSink(LogManager::getStrmWorld(), "EP"));
  LogManager::setLog("EP");
  MOFEM_LOG_TAG("EP", "ep");
  core_log->add_sink(
      LogManager::createSink(LogManager::getStrmSelf(), "EPSELF"));
  LogManager::setLog("EPSELF");
  MOFEM_LOG_TAG("EPSELF", "ep");
 
#ifdef ENABLE_PYTHON_BINDING
  Py_Initialize();
  np::initialize();
  MOFEM_LOG("EP", Sev::inform) << "Python initialised";
#else
  MOFEM_LOG("EP", Sev::inform) << "Python NOT initialised";
#endif
 
  core_log->add_sink(
      LogManager::createSink(LogManager::getStrmSync(), "EPSYNC"));
  LogManager::setLog("EPSYNC");
  MOFEM_LOG_TAG("EPSYNC", "ep");
 
  try {
 
    // Get mesh file
    PetscBool flg = PETSC_TRUE;
    char mesh_file_name[255];
    CHKERR PetscOptionsGetString(PETSC_NULLPTR, "", "-file_name", mesh_file_name,
                                 255, &flg);
    // Register DM Manager
    DMType dm_name = "DMMOFEM";
    CHKERR DMRegister_MoFEM(dm_name);
    DMType dm_name_mg = "DMMOFEM_MG";
    CHKERR DMRegister_MGViaApproxOrders(dm_name_mg);
 
    // Create MoAB database
    moab::Core moab_core;
    moab::Interface &moab = moab_core;
 
    ParallelComm *pcomm = ParallelComm::get_pcomm(&moab, MYPCOMM_INDEX);
    if (pcomm == NULL)
      pcomm = new ParallelComm(&moab, PETSC_COMM_WORLD);
    // Read mesh to MOAB
    PetscBool fully_distributed = PETSC_FALSE;
    CHKERR PetscOptionsGetBool(PETSC_NULLPTR, "", "-fully_distributed",
                               &fully_distributed, PETSC_NULLPTR);
    if (fully_distributed) {
      const char *option;
      if (pcomm->proc_config().proc_size() == 1)
        option = "";
      else
        option = "PARALLEL=READ_PART;"
                 "PARALLEL_RESOLVE_SHARED_ENTS;"
                 "PARTITION=PARALLEL_PARTITION";
      CHKERR moab.load_file(mesh_file_name, 0, option);
    } else {
      CHKERR CommInterface::loadFileRootProcAllRestDistributed(
          moab, mesh_file_name, SPACE_DIM);
    }
 
    // Create MoFEM database and link it to MoAB
    MoFEM::Core mofem_core(moab);
    MoFEM::Interface &m_field = mofem_core;
 
    BitRefLevel bit_level0 = BitRefLevel().set(0);
    CHKERR m_field.getInterface<BitRefManager>()->setBitRefLevelByDim(
        0, 3, bit_level0);
    
    // Data structures
    EshelbianCore ep(m_field);
 
    auto meshset_ptr = get_temp_meshset_ptr(moab);
    CHKERR moab.add_entities(
        *meshset_ptr, CommInterface::getPartEntities(moab, pcomm->rank()));
    auto get_adj = [&](Range ents, int dim) {
      Range adj;
      CHKERR moab.get_adjacencies(ents, dim, false, adj,
                                  moab::Interface::UNION);
      return adj;
    };
    CHKERR moab.add_entities(
        *meshset_ptr,
        get_adj(CommInterface::getPartEntities(moab, pcomm->rank())
                    .subset_by_dimension(SPACE_DIM),
                2));
    CHKERR moab.add_entities(
        *meshset_ptr,
        get_adj(CommInterface::getPartEntities(moab, pcomm->rank())
                    .subset_by_dimension(SPACE_DIM),
                1));
    CHKERR moab.add_entities(
        *meshset_ptr,
        get_adj(CommInterface::getPartEntities(moab, pcomm->rank())
                    .subset_by_dimension(SPACE_DIM),
                0));
 
    // create growing crack surface meshset
    CHKERR ep.createCrackSurfaceMeshset();
    CHKERR ep.createExchangeVectors(Sev::inform);
 
    CHKERR ep.projectGeometry();
    CHKERR ep.addCrackSurfaces(true);
    CHKERR ep.projectGeometry();
  }
  CATCH_ERRORS;
 
  // finish work cleaning memory, getting statistics, etc.
  MoFEM::Core::Finalize();
 
#ifdef ENABLE_PYTHON_BINDING
  if (Py_FinalizeEx() < 0) {
    exit(120);
  }
#endif
}

Variable Documentation

adolc_tag

int adolc_tag = 1

constexpr

Examples

ep.cpp.

Definition at line 25 of file debug_crack_mesh.cpp.

help

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

static

Definition at line 48 of file debug_crack_mesh.cpp.