#include <boost/python.hpp>
#include <boost/python/def.hpp>
#include <boost/python/numpy.hpp>
#include <MoFEM.hpp>
#include <ElasticSpring.hpp>
#include <FluidLevel.hpp>
#include <CalculateTraction.hpp>
#include <NaturalDomainBC.hpp>
#include <NaturalBoundaryBC.hpp>
#include <HookeOps.hpp>
#include <ElasticPostProc.hpp>
#include <ObjectiveFunctionData.hpp>

Classes
struct	DomainBCs
	[Define entities] More...

struct	BoundaryBCs
	Boundary conditions marker. More...

struct	PostProcEleByDim< 2 >

struct	PostProcEleByDim< 3 >

struct	Example
	[Example] More...

struct	OpStateSensitivity

struct	OpObjective

struct	OpAdJointObjective

Typedefs
using	EntData = EntitiesFieldData::EntData
	Entity data for field operations.

using	DomainEle = PipelineManager::ElementsAndOpsByDim< SPACE_DIM >::DomainEle
	Domain finite elements.

using	BoundaryEle = PipelineManager::ElementsAndOpsByDim< SPACE_DIM >::BoundaryEle
	Boundary finite elements.

using	DomainEleOp = DomainEle::UserDataOperator
	Domain element operators.

using	BoundaryEleOp = BoundaryEle::UserDataOperator

using	DomainRhsBCs = NaturalBC< DomainEleOp >::Assembly< A >::LinearForm< I >
	Domain RHS natural BCs.

using	OpDomainRhsBCs = DomainRhsBCs::OpFlux< DomainBCs, 1, SPACE_DIM >
	Domain flux operator.

using	BoundaryRhsBCs = NaturalBC< BoundaryEleOp >::Assembly< A >::LinearForm< I >
	Boundary RHS natural BCs.

using	OpBoundaryRhsBCs = BoundaryRhsBCs::OpFlux< BoundaryBCs, 1, SPACE_DIM >
	Boundary flux operator.

using	BoundaryLhsBCs = NaturalBC< BoundaryEleOp >::Assembly< A >::BiLinearForm< I >
	Boundary LHS natural BCs.

using	OpBoundaryLhsBCs = BoundaryLhsBCs::OpFlux< BoundaryBCs, 1, SPACE_DIM >
	Boundary LHS flux operator.

using	PostProcEleDomain = PostProcEleByDim< SPACE_DIM >::PostProcEleDomain

using	SideEle = PostProcEleByDim< SPACE_DIM >::SideEle

using	PostProcEleBdy = PostProcEleByDim< SPACE_DIM >::PostProcEleBdy

using	DomainBaseOp = FormsIntegrators< DomainEleOp >::Assembly< A >::OpBase
	[Postprocess results]

Enumerations
enum	SensitivityMethod { DIRECT , ADJOINT }

Functions
Range	get_range_from_block (MoFEM::Interface &m_field, const std::string block_name, int dim)

MoFEMErrorCode	save_range (moab::Interface &moab, const std::string name, const Range r)

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

Variables
constexpr int	BASE_DIM = 1
	[Constants and material properties]

constexpr int	SPACE_DIM
	[Define dimension]

constexpr AssemblyType	A
	[Define dimension]

constexpr IntegrationType	I
	Use Gauss quadrature for integration.

constexpr double	young_modulus = 1
	[Material properties for linear elasticity]

constexpr double	poisson_ratio = 0.3
	Poisson's ratio ν

constexpr double	bulk_modulus_K
	Bulk modulus K = E/(3(1-2ν))

constexpr double	shear_modulus_G
	Shear modulus G = E/(2(1+ν))

PetscBool	is_plane_strain

SensitivityMethod	derivative_type = ADJOINT

static char	help [] = "...\n\n"
	[Finite difference check]

Typedef Documentation

◆ BoundaryEle

using BoundaryEle = PipelineManager::ElementsAndOpsByDim< SPACE_DIM>::BoundaryEle

Boundary finite elements.

Definition at line 46 of file gradient.cpp.

◆ BoundaryEleOp

using BoundaryEleOp = BoundaryEle::UserDataOperator

Boundary element operators

Definition at line 49 of file gradient.cpp.

◆ BoundaryLhsBCs

using BoundaryLhsBCs = NaturalBC<BoundaryEleOp>::Assembly<A>::BiLinearForm< I>

Boundary LHS natural BCs.

Definition at line 67 of file gradient.cpp.

◆ BoundaryRhsBCs

using BoundaryRhsBCs = NaturalBC<BoundaryEleOp>::Assembly<A>::LinearForm< I>

Boundary RHS natural BCs.

Definition at line 62 of file gradient.cpp.

◆ DomainBaseOp

using DomainBaseOp = FormsIntegrators<DomainEleOp>::Assembly<A>::OpBase

[Postprocess results]

Examples: mofem/tutorials/vec-7_shape_optimisation/adjoint.cpp.

Definition at line 703 of file gradient.cpp.

◆ DomainEle

using DomainEle = PipelineManager::ElementsAndOpsByDim< SPACE_DIM>::DomainEle

Domain finite elements.

Definition at line 44 of file gradient.cpp.

◆ DomainEleOp

using DomainEleOp = DomainEle::UserDataOperator

Domain element operators.

Definition at line 48 of file gradient.cpp.

◆ DomainRhsBCs

using DomainRhsBCs = NaturalBC<DomainEleOp>::Assembly<A>::LinearForm< I>

Domain RHS natural BCs.

Definition at line 58 of file gradient.cpp.

◆ EntData

using EntData = EntitiesFieldData::EntData

Entity data for field operations.

Definition at line 42 of file gradient.cpp.

◆ OpBoundaryLhsBCs

using OpBoundaryLhsBCs = BoundaryLhsBCs::OpFlux<BoundaryBCs, 1, SPACE_DIM>

Boundary LHS flux operator.

Definition at line 69 of file gradient.cpp.

◆ OpBoundaryRhsBCs

using OpBoundaryRhsBCs = BoundaryRhsBCs::OpFlux<BoundaryBCs, 1, SPACE_DIM>

Boundary flux operator.

Definition at line 64 of file gradient.cpp.

◆ OpDomainRhsBCs

using OpDomainRhsBCs = DomainRhsBCs::OpFlux<DomainBCs, 1, SPACE_DIM>

Domain flux operator.

Definition at line 60 of file gradient.cpp.

◆ PostProcEleBdy

using PostProcEleBdy = PostProcEleByDim<SPACE_DIM>::PostProcEleBdy

Definition at line 89 of file gradient.cpp.

◆ PostProcEleDomain

using PostProcEleDomain = PostProcEleByDim<SPACE_DIM>::PostProcEleDomain

Definition at line 87 of file gradient.cpp.

◆ SideEle

using SideEle = PostProcEleByDim<SPACE_DIM>::SideEle

Definition at line 88 of file gradient.cpp.

Enumeration Type Documentation

◆ SensitivityMethod

enum SensitivityMethod

Enumerator
DIRECT
ADJOINT

Definition at line 94 of file gradient.cpp.

94{ DIRECT, ADJOINT };

DIRECT

@ DIRECT

Definition gradient.cpp:94

ADJOINT

@ ADJOINT

Definition gradient.cpp:94

Function Documentation

◆ get_range_from_block()

Range get_range_from_block	(	MoFEM::Interface &	m_field,
		const std::string	block_name,
		int	dim
	)

Definition at line 1347 of file gradient.cpp.

                                                                {
  Range r;
 
  auto mesh_mng = m_field.getInterface<MeshsetsManager>();
  auto bcs = mesh_mng->getCubitMeshsetPtr(
 
      std::regex((boost::format("%s(.*)") % block_name).str())
 
  );
 
  for (auto bc : bcs) {
    Range faces;
    CHK_MOAB_THROW(bc->getMeshsetIdEntitiesByDimension(m_field.get_moab(), dim,
                                                       faces, true),
                   "get meshset ents");
    r.merge(faces);
  }
 
  for (auto dd = dim - 1; dd >= 0; --dd) {
    if (dd >= 0) {
      Range ents;
      CHK_MOAB_THROW(m_field.get_moab().get_adjacencies(r, dd, false, ents,
                                                        moab::Interface::UNION),
                     "get adjs");
      r.merge(ents);
    } else {
      Range verts;
      CHK_MOAB_THROW(m_field.get_moab().get_connectivity(r, verts),
                     "get verts");
      r.merge(verts);
    }
    CHK_MOAB_THROW(
        m_field.getInterface<CommInterface>()->synchroniseEntities(r), "comm");
  }
 
  return r;
};

◆ main()

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

[Register MoFEM discrete manager in PETSc]

[Register MoFEM discrete manager in PETSc

[Create MoAB]

< mesh database

< mesh database interface

[Create MoAB]

[Create MoFEM]

< finite element database

< finite element database interface

[Create MoFEM]

[Example]

Definition at line 1294 of file gradient.cpp.

                                 {
 
  // Initialize Python environment for objective function interface
  Py_Initialize();
  np::initialize();
 
  // Initialize MoFEM/PETSc and MOAB data structures
  const char param_file[] = "param_file.petsc";
  MoFEM::Core::Initialize(&argc, &argv, param_file, help);
 
  auto core_log = logging::core::get();
  core_log->add_sink(
      LogManager::createSink(LogManager::getStrmWorld(), "TIMER"));
 
  core_log->add_sink(
      LogManager::createSink(LogManager::getStrmSync(), "FieldEvaluator"));
  LogManager::setLog("FieldEvaluator");
  MOFEM_LOG_TAG("FieldEvaluator", "field_eval");
 
  try {
 
    //! [Register MoFEM discrete manager in PETSc]
    DMType dm_name = "DMMOFEM";
    CHKERR DMRegister_MoFEM(dm_name);
    DMType dm_name_mg = "DMMOFEM_MG";
    CHKERR DMRegister_MGViaApproxOrders(dm_name_mg);
    //! [Register MoFEM discrete manager in PETSc
 
    //! [Create MoAB]
    moab::Core mb_instance;              ///< mesh database
    moab::Interface &moab = mb_instance; ///< mesh database interface
    //! [Create MoAB]
 
    //! [Create MoFEM]
    MoFEM::Core core(moab);           ///< finite element database
    MoFEM::Interface &m_field = core; ///< finite element database interface
    //! [Create MoFEM]
 
    //! [Example]
 
    Example ex(m_field);
    CHKERR ex.runProblem();
    //! [Example]
  }
  CATCH_ERRORS;
 
  CHKERR MoFEM::Core::Finalize();
 
  if (Py_FinalizeEx() < 0) {
    exit(120);
  }
}

◆ save_range()

MoFEMErrorCode save_range	(	moab::Interface &	moab,
		const std::string	name,
		const Range	r
	)

Definition at line 1386 of file gradient.cpp.

                                         {
  MoFEMFunctionBegin;
  auto out_meshset = get_temp_meshset_ptr(moab);
  CHKERR moab.add_entities(*out_meshset, r);
  CHKERR moab.write_file(name.c_str(), "VTK", "", out_meshset->get_ptr(), 1);
  MoFEMFunctionReturn(0);
};

Variable Documentation

◆ A

constexpr AssemblyType A

constexpr

Initial value:

=

AssemblyType::PETSC

[Define dimension]

Use PETSc for matrix/vector assembly

Definition at line 23 of file gradient.cpp.

◆ BASE_DIM

constexpr int BASE_DIM = 1

constexpr

[Constants and material properties]

Dimension of the base functions

Definition at line 16 of file gradient.cpp.

◆ bulk_modulus_K

constexpr double bulk_modulus_K

constexpr

Initial value:

=
    young_modulus /
    (3 * (1 - 2 * poisson_ratio))

Bulk modulus K = E/(3(1-2ν))

Definition at line 31 of file gradient.cpp.

◆ derivative_type

SensitivityMethod derivative_type = ADJOINT

Definition at line 96 of file gradient.cpp.

◆ help

char help[] = "...\n\n"

static

[Finite difference check]

Definition at line 1292 of file gradient.cpp.

◆ I

constexpr IntegrationType I

constexpr

Initial value:

=

IntegrationType::GAUSS

Use Gauss quadrature for integration.

Examples: mofem/tutorials/vec-7_shape_optimisation/adjoint.cpp.

Definition at line 25 of file gradient.cpp.

◆ is_plane_strain

PetscBool is_plane_strain

Initial value:

=

PETSC_FALSE

Flag for plane strain vs plane stress in 2D

Definition at line 37 of file gradient.cpp.

◆ poisson_ratio

constexpr double poisson_ratio = 0.3

constexpr

Poisson's ratio ν

Definition at line 30 of file gradient.cpp.

◆ shear_modulus_G

constexpr double shear_modulus_G

constexpr

Initial value:

=

young_modulus / (2 * (1 + poisson_ratio))

Shear modulus G = E/(2(1+ν))

Definition at line 34 of file gradient.cpp.

◆ SPACE_DIM

constexpr int SPACE_DIM

constexpr

Initial value:

=

EXECUTABLE_DIMENSION

#define EXECUTABLE_DIMENSION

Definition plastic.cpp:13

[Define dimension]

Space dimension of problem (2D or 3D), set at compile time

Examples: mofem/tutorials/vec-7_shape_optimisation/adjoint.cpp.

Definition at line 19 of file gradient.cpp.

◆ young_modulus

constexpr double young_modulus = 1

constexpr

[Material properties for linear elasticity]

Young's modulus E

Definition at line 29 of file gradient.cpp.

Classes

Typedefs

Enumerations

Functions

Variables

Typedef Documentation

◆ BoundaryEle

◆ BoundaryEleOp

◆ BoundaryLhsBCs

◆ BoundaryRhsBCs

◆ DomainBaseOp

◆ DomainEle

◆ DomainEleOp

◆ DomainRhsBCs

◆ EntData

◆ OpBoundaryLhsBCs

◆ OpBoundaryRhsBCs

◆ OpDomainRhsBCs

◆ PostProcEleBdy

◆ PostProcEleDomain

◆ SideEle

Enumeration Type Documentation

◆ SensitivityMethod

Function Documentation

◆ get_range_from_block()

◆ main()

◆ save_range()

Variable Documentation

◆ A

◆ BASE_DIM

◆ bulk_modulus_K

◆ derivative_type

◆ help

◆ I

◆ is_plane_strain

◆ poisson_ratio

◆ shear_modulus_G

◆ SPACE_DIM

◆ young_modulus