EshelbianPlasticity::solve_elastic_setup Struct Reference

Static Public Member Functions
static auto	setup (EshelbianCore *ep_ptr, TS ts, Vec x, bool set_ts_monitor)

Detailed Description

Definition at line 2808 of file EshelbianPlasticity.cpp.

Member Function Documentation

setup()

static auto EshelbianPlasticity::solve_elastic_setup::setup ( EshelbianCore *  ep_ptr, TS  ts, Vec  x, bool  set_ts_monitor  )

inlinestatic

Definition at line 2810 of file EshelbianPlasticity.cpp.

                                                {
 
#ifdef PYTHON_SDF
    auto setup_sdf = [&]() {
      boost::shared_ptr<ContactOps::SDFPython> sdf_python_ptr;
 
      auto file_exists = [](std::string myfile) {
        std::ifstream file(myfile.c_str());
        if (file) {
          return true;
        }
        return false;
      };
 
      char sdf_file_name[255] = "sdf.py";
      CHKERR PetscOptionsGetString(PETSC_NULL, PETSC_NULL, "-sdf_file",
                                   sdf_file_name, 255, PETSC_NULL);
 
      if (file_exists(sdf_file_name)) {
        MOFEM_LOG("EP", Sev::inform) << sdf_file_name << " file found";
        sdf_python_ptr = boost::make_shared<ContactOps::SDFPython>();
        CHKERR sdf_python_ptr->sdfInit(sdf_file_name);
        ContactOps::sdfPythonWeakPtr = sdf_python_ptr;
      } else {
        MOFEM_LOG("EP", Sev::warning) << sdf_file_name << " file NOT found";
      }
      return sdf_python_ptr;
    };
    auto sdf_python_ptr = setup_sdf();
#endif
 
    auto setup_ts_monitor = [&]() {
      boost::shared_ptr<TsCtx> ts_ctx;
      DMMoFEMGetTsCtx(ep_ptr->dmElastic, ts_ctx);
      if (set_ts_monitor) {
        CHKERR TSMonitorSet(ts, TsMonitorSet, ts_ctx.get(), PETSC_NULL);
        auto monitor_ptr = boost::make_shared<EshelbianMonitor>(*ep_ptr);
        ts_ctx->getLoopsMonitor().push_back(
            TsCtx::PairNameFEMethodPtr(ep_ptr->elementVolumeName, monitor_ptr));
      }
      return std::make_tuple(ts_ctx);
    };
 
    auto setup_snes_monitor = [&]() {
      MoFEMFunctionBeginHot;
      SNES snes;
      CHKERR TSGetSNES(ts, &snes);
      PetscViewerAndFormat *vf;
      CHKERR PetscViewerAndFormatCreate(PETSC_VIEWER_STDOUT_WORLD,
                                        PETSC_VIEWER_DEFAULT, &vf);
      CHKERR SNESMonitorSet(
          snes,
          (MoFEMErrorCode(*)(SNES, PetscInt, PetscReal,
                             void *))SNESMonitorFields,
          vf, (MoFEMErrorCode(*)(void **))PetscViewerAndFormatDestroy);
      MoFEMFunctionReturnHot(0);
    };
 
    auto setup_section = [&]() {
      PetscSection section_raw;
      CHKERR DMGetSection(ep_ptr->dmElastic, &section_raw);
      int num_fields;
      CHKERR PetscSectionGetNumFields(section_raw, &num_fields);
      for (int ff = 0; ff != num_fields; ff++) {
        const char *field_name;
        CHKERR PetscSectionGetFieldName(section_raw, ff, &field_name);
        MOFEM_LOG_C("EP", Sev::inform, "Field %d name %s", ff, field_name);
      }
      return section_raw;
    };
 
    auto set_vector_on_mesh = [&]() {
      MoFEMFunctionBeginHot;
      CHKERR DMoFEMMeshToLocalVector(ep_ptr->dmElastic, x, INSERT_VALUES,
                                     SCATTER_FORWARD);
      CHKERR VecGhostUpdateBegin(x, INSERT_VALUES, SCATTER_FORWARD);
      CHKERR VecGhostUpdateEnd(x, INSERT_VALUES, SCATTER_FORWARD);
      MoFEMFunctionReturnHot(0);
    };
 
    auto setup_schur_block_solver = [&]() {
      CHKERR TSAppendOptionsPrefix(ts, "elastic_");
      CHKERR TSSetFromOptions(ts);
      CHKERR TSSetDM(ts, ep_ptr->dmElastic);
      // Adding field split solver
      boost::shared_ptr<EshelbianCore::SetUpSchur> schur_ptr;
      if constexpr (A == AssemblyType::BLOCK_MAT) {
        schur_ptr =
            EshelbianCore::SetUpSchur::createSetUpSchur(ep_ptr->mField, ep_ptr);
        CHK_THROW_MESSAGE(schur_ptr->setUp(ts), "setup schur");
      }
      return schur_ptr;
    };
 
    // Warning: sequence of construction is not guaranteed for tuple. You have
    // to enforce order by proper packaging.
 
#ifdef PYTHON_SDF
    return std::make_tuple(setup_sdf(), setup_ts_monitor(),
                           setup_snes_monitor(), setup_section(),
                           set_vector_on_mesh(), setup_schur_block_solver());
#else
    return std::make_tuple(setup_ts_monitor(), setup_snes_monitor(),
                           setup_section(), set_vector_on_mesh(),
                           setup_schur_block_solver());
#endif
  }

---

The documentation for this struct was generated from the following file:

• users_modules/eshelbian_plasticit/src/impl/EshelbianPlasticity.cpp