AuxiliaryLogarithmicStressTestSetup.hpp

Go to the documentation of this file.

/**
 * @file AuxiliaryLogarithmicStressTestSetup.hpp
 * @brief Shared mesh and problem setup for auxiliary formulation atoms.
 */
 
#ifndef AUXILIARY_LOGARITHMIC_STRESS_TEST_SETUP_HPP
#define AUXILIARY_LOGARITHMIC_STRESS_TEST_SETUP_HPP
 
#include <EshelbianPlasticity.hpp>
 
namespace EshelbianPlasticity {
 
inline MoFEMErrorCode runAuxiliaryLogarithmicStressTest(
    const std::function<MoFEMErrorCode(EshelbianCore &)> &test) {
  MoFEMFunctionBegin;
  CHKERR JsonConfigManager::setMeshFileFromJson();
  char mesh_file[PETSC_MAX_PATH_LEN] = "";
  PetscBool found;
  CHKERR PetscOptionsGetString(nullptr, nullptr, "-file_name", mesh_file,
                               sizeof(mesh_file), &found);
  if (!found)
    SETERRQ(PETSC_COMM_WORLD, MOFEM_ATOM_TEST_INVALID,
            "Missing JSON mesh fixture");
  CHKERR DMRegister_MoFEM("DMMOFEM");
  moab::Core moab;
  auto *parallel = new ParallelComm(&moab, PETSC_COMM_WORLD);
  CHKERR CommInterface::loadFileRootProcAllRestDistributed(moab, mesh_file,
                                                           SPACE_DIM);
  MoFEM::Core core(moab);
  MoFEM::Interface &m_field = core;
  CHKERR m_field.getInterface<MeshsetsManager>()->setMeshsetFromFile();
  CHKERR m_field.getInterface<BitRefManager>()->setBitRefLevelByDim(
      0, SPACE_DIM, BitRefLevel().set(0));
  EshelbianCore ep(m_field);
  auto meshset = get_temp_meshset_ptr(moab);
  const auto part = CommInterface::getPartEntities(moab, parallel->rank());
  CHKERR moab.add_entities(*meshset, part);
  for (const int dim : {0, 1, 2}) {
    Range adjacent;
    CHKERR moab.get_adjacencies(part.subset_by_dimension(SPACE_DIM), dim, false,
                                adjacent, moab::Interface::UNION);
    CHKERR moab.add_entities(*meshset, adjacent);
  }
  CHKERR ep.createCrackSurfaceMeshset();
  CHKERR ep.getSpatialDispBc();
  CHKERR ep.getSpatialRotationBc();
  CHKERR ep.getSpatialTractionBc();
  CHKERR ep.getSpatialTractionFreeBc();
  CHKERR ep.getExternalStrain();
  CHKERR ep.setBlockTagsOnSkin();
  CHKERR ep.createExchangeVectors(Sev::inform);
  CHKERR PhysicalEquations::create(ep);
  CHKERR ep.resolveDissipationEntities(*meshset);
  CHKERR ep.addFields(*meshset);
  CHKERR ep.projectGeometry(*meshset, 0.);
  CHKERR ep.addVolumeFiniteElement(*meshset);
  CHKERR ep.addBoundaryFiniteElement(*meshset);
  CHKERR ep.addDMs();
  CHKERR test(ep);
  MoFEMFunctionReturn(0);
}
 
} // namespace EshelbianPlasticity
 
#endif