bc_displacement_blocks.cpp File Reference

Check displacement types in JSON names containing other BC tokens. More...

#include <MoFEM.hpp>

Go to the source code of this file.

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

Detailed Description

Check displacement types in JSON names containing other BC tokens.

Definition in file bc_displacement_blocks.cpp.

Function Documentation

main()

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

Definition at line 7 of file bc_displacement_blocks.cpp.

                                 {
  MoFEM::Core::Initialize(&argc, &argv, nullptr,
                          "Displacement block name parsing\n");
  try {
    moab::Core moab_core;
    MoFEM::Core core(moab_core);
    MoFEM::Interface &m_field = core;
    auto *meshsets = m_field.getInterface<MeshsetsManager>();
    CHKERR meshsets->addMeshset(BLOCKSET, 5, "FIX_Z");
    CHKERR meshsets->addMeshset(BLOCKSET, 6, "FIX_ALL");
    auto *config = m_field.getInterface<JsonConfigManager>();
    CHKERR config->getSubInterfaceOptions();
    CHKERR meshsets->setMeshsetFromFile();
    CHKERR config->applyMeshsets();
 
    auto *bc_manager = m_field.getInterface<BcManager>();
    // An empty problem name exercises the BC data parser without a FE problem.
    CHKERR bc_manager
        ->pushMarkDOFsOnEntities<BcDisplacementMeshsetType<BLOCKSET>>("", "P",
                                                                      false);
    const auto &blocks = bc_manager->getBcMapByBlockName();
    const auto all_it = blocks.find("_P_FIX_ALL_FIX_Z_5");
    const auto z_it = blocks.find("_P_FIX_Z_FIX_ALL_6");
    if (all_it == blocks.end() || z_it == blocks.end() ||
        !all_it->second->dispBcPtr || !z_it->second->dispBcPtr)
      SETERRQ(PETSC_COMM_WORLD, MOFEM_ATOM_TEST_INVALID,
              "Missing displacement BC data for generated JSON names");
    const auto &all = all_it->second->dispBcPtr->data;
    const auto &z = z_it->second->dispBcPtr->data;
    if (all.flag1 != 1 || all.flag2 != 1 || all.flag3 != 1 ||
        all.value1 != 1.0 || all.value2 != 2.0 || all.value3 != 3.0)
      SETERRQ(
          PETSC_COMM_WORLD, MOFEM_ATOM_TEST_INVALID,
          "FIX_ALL must constrain x, y and z despite the FIX_Z source name");
    if (z.flag1 != 0 || z.flag2 != 0 || z.flag3 != 1 || z.value3 != 6.0)
      SETERRQ(PETSC_COMM_WORLD, MOFEM_ATOM_TEST_INVALID,
              "FIX_Z must constrain only z despite the FIX_ALL source name");
  }
  CATCH_ERRORS;
  CHKERR MoFEM::Core::Finalize();
  return 0;
}