neohookean_stretch_gradient_atom.cpp

Go to the documentation of this file.

#include <MoFEM.hpp>
 
using namespace MoFEM;
 
namespace EshelbianPlasticity {
MoFEMErrorCode testHMHNeohookeanStretchGradient();
} // namespace EshelbianPlasticity
 
static char help[] =
    "Finite-difference check of HMHNeohookean dHenckyStretch/dBiotStress.\n"
    "Reproduction options:\n"
    "  -nh_stretch_atom_c10 <real>\n"
    "  -nh_stretch_atom_k <real>\n"
    "  -nh_stretch_atom_q <real>\n"
    "  -nh_stretch_atom_max_it <int>\n"
    "  -nh_stretch_atom_tol <real>\n"
    "  -nh_stretch_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", "neohookean_stretch_gradient_atom");
 
  try {
    CHKERR EshelbianPlasticity::testHMHNeohookeanStretchGradient();
  }
  CATCH_ERRORS;
 
  MoFEM::Core::Finalize();
  return 0;
}