storakers_stress_to_stretch_atom.cpp

Go to the documentation of this file.

#include <MoFEM.hpp>
 
using namespace MoFEM;
 
namespace EshelbianPlasticity {
MoFEMErrorCode testHMHStorakersStressToStretch();
} // namespace EshelbianPlasticity
 
static char help[] =
    "Validate Storakers stress-to-stretch inversion and compliance.\n"
    "Options:\n"
    "  -storakers_atom_fd_step <real>\n"
    "  -storakers_atom_stretch_tol <real>\n"
    "  -storakers_atom_compliance_tol <real>\n"
    "  -storakers_atom_solver_tol <real>\n"
    "Reproduction options:\n"
    "  -storakers_atom_eta <real>\n"
    "  -storakers_atom_mu <real>\n"
    "  -storakers_atom_beta <real>\n"
    "  -storakers_atom_q <real>\n"
    "  -storakers_atom_initial_y <real>\n"
    "  -storakers_atom_biot_stress <s00,s01,s02,s11,s12,s22>\n";
 
int main(int argc, char *argv[]) {
  MoFEM::Core::Initialize(&argc, &argv, (char *)0, help);
 
  auto core_log = logging::core::get();
  core_log->add_sink(LogManager::createSink(LogManager::getStrmSelf(), "EP"));
  LogManager::setLog("EP");
  MOFEM_LOG_TAG("EP", "storakers_stress_to_stretch_atom");
 
  try {
    CHKERR EshelbianPlasticity::testHMHStorakersStressToStretch();
  }
  CATCH_ERRORS;
 
  MoFEM::Core::Finalize();
  return 0;
}