auxiliary_logarithmic_stress_space_atom.cpp File Reference

Weak mesh-stress closure, observability and coefficient elimination. More...

#include <MoFEM.hpp>
#include <EshelbianPlasticity.hpp>
#include <NeoHookeanLogarithmicMaterial.hpp>
#include <IntegrationRules.hpp>

Go to the source code of this file.

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

Variables
static char	help []

Detailed Description

Weak mesh-stress closure, observability and coefficient elimination.

Definition in file auxiliary_logarithmic_stress_space_atom.cpp.

Function Documentation

main()

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

Definition at line 480 of file auxiliary_logarithmic_stress_space_atom.cpp.

                                 {
  MoFEM::Core::Initialize(&argc, &argv, nullptr, help);
  try {
    moab::Core moab;
    MoFEM::Core core(moab);
    MoFEM::Interface &m_field = core;
    auto *json = m_field.getInterface<JsonConfigManager>();
    auto *meshsets = m_field.getInterface<MeshsetsManager>();
    CHKERR meshsets->addMeshset(BLOCKSET, 1001, "MAT_NEOHOOKEAN");
    CHKERR meshsets->setMeshsetFromFile();
    const auto block = json->getParamsFromBlockset("MAT_NEOHOOKEAN", 1001);
    if (block.size() != 2 || !block.count("c10") || !block.count("k"))
      SETERRQ(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
              "Space atom requires MAT_NEOHOOKEAN JSON block 1001 with c10,k");
    const Material::Parameters parameters{block.at("c10"), block.at("k")};
    CHKERR Material::validateParameters(parameters);
    PetscInt quadrature_order = 0;
    PetscReal tolerance = 0., newton_tolerance = 0.;
    CHKERR PetscOptionsGetInt(nullptr, nullptr, "-auxiliary_space_quadrature_order",
                              &quadrature_order, nullptr);
    CHKERR PetscOptionsGetReal(nullptr, nullptr, "-auxiliary_space_tolerance",
                               &tolerance, nullptr);
    CHKERR PetscOptionsGetReal(nullptr, nullptr,
                               "-auxiliary_space_newton_tolerance",
                               &newton_tolerance, nullptr);
    if (quadrature_order < 4 || !std::isfinite(tolerance) || tolerance <= 0. ||
        !std::isfinite(newton_tolerance) || newton_tolerance <= 0. ||
        newton_tolerance >= tolerance)
      SETERRQ(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
              "Space atom requires positive finite JSON tolerances and "
              "quadrature order at least four");
    CHKERR runSpaceAtom(parameters, quadrature_order, tolerance, newton_tolerance);
  }
  CATCH_ERRORS;
  MoFEM::Core::Finalize();
  return 0;
}

Variable Documentation

help

char help[]

static

Initial value:

=
    "Check weak auxiliary stress spaces and coefficient Schur elimination.\n"

Definition at line 477 of file auxiliary_logarithmic_stress_space_atom.cpp.