adjoint.cpp File Reference

#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>

Go to the source code of this file.

Classes
struct	DomainBCs
	[Define entities] More...
struct	BoundaryBCs
struct	PostProcEleByDim< 2 >
struct	PostProcEleByDim< 3 >
struct	ObjectiveFunctionData
struct	Example
	[Example] More...
struct	OpAdJointTestOp< SPACE_DIM, IntegrationType::GAUSS, DomainBaseOp >
struct	OpAdJointGradTimesSymTensor< SPACE_DIM, IntegrationType::GAUSS, DomainBaseOp >
struct	OpAdJointObjective
struct	ObjectiveFunctionDataImpl

Typedefs
using	DomainEle = PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::DomainEle
using	BoundaryEle
using	DomainRhsBCs = NaturalBC<DomainEleOp>::Assembly<A>::LinearForm<I>
using	OpDomainRhsBCs = DomainRhsBCs::OpFlux<DomainBCs, 1, SPACE_DIM>
using	BoundaryRhsBCs = NaturalBC<BoundaryEleOp>::Assembly<A>::LinearForm<I>
using	OpBoundaryRhsBCs = BoundaryRhsBCs::OpFlux<BoundaryBCs, 1, SPACE_DIM>
using	BoundaryLhsBCs = NaturalBC<BoundaryEleOp>::Assembly<A>::BiLinearForm<I>
using	OpBoundaryLhsBCs = BoundaryLhsBCs::OpFlux<BoundaryBCs, 1, SPACE_DIM>
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]

Functions
boost::shared_ptr< ObjectiveFunctionData >	create_python_objective_function (std::string py_file)
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)
template<int DIM>
auto	diff_symmetrize (FTensor::Number< DIM >)
int	main (int argc, char *argv[])

Variables
constexpr int	BASE_DIM = 1
constexpr int	SPACE_DIM
	[Define dimension]
constexpr AssemblyType	A = AssemblyType::PETSC
	[Define dimension]
constexpr IntegrationType	I
constexpr double	young_modulus = 1
constexpr double	poisson_ratio = 0.3
constexpr double	bulk_modulus_K = young_modulus / (3 * (1 - 2 * poisson_ratio))
constexpr double	shear_modulus_G = young_modulus / (2 * (1 + poisson_ratio))
PetscBool	is_plane_strain = PETSC_FALSE
static char	help [] = "...\n\n"
	[calculateGradient]

Typedef Documentation

BoundaryEle

using BoundaryEle

Initial value:

 
    PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::BoundaryEle

Definition at line 33 of file adjoint.cpp.

BoundaryLhsBCs

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

Definition at line 46 of file adjoint.cpp.

BoundaryRhsBCs

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

Definition at line 44 of file adjoint.cpp.

DomainBaseOp

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

[Postprocess results]

[calculateGradient]

Examples

mofem/tutorials/vec-7/adjoint.cpp.

Definition at line 1355 of file adjoint.cpp.

DomainEle

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

Definition at line 32 of file adjoint.cpp.

DomainRhsBCs

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

Definition at line 42 of file adjoint.cpp.

OpBoundaryLhsBCs

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

Definition at line 47 of file adjoint.cpp.

OpBoundaryRhsBCs

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

Definition at line 45 of file adjoint.cpp.

OpDomainRhsBCs

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

Definition at line 43 of file adjoint.cpp.

PostProcEleBdy

using PostProcEleBdy = PostProcEleByDim<SPACE_DIM>::PostProcEleBdy

Definition at line 65 of file adjoint.cpp.

PostProcEleDomain

using PostProcEleDomain = PostProcEleByDim<SPACE_DIM>::PostProcEleDomain

Definition at line 63 of file adjoint.cpp.

SideEle

using SideEle = PostProcEleByDim<SPACE_DIM>::SideEle

Definition at line 64 of file adjoint.cpp.

Function Documentation

create_python_objective_function()

boost::shared_ptr< ObjectiveFunctionData > create_python_objective_function

(

std::string

py_file

)

Examples

mofem/tutorials/vec-7/adjoint.cpp.

Definition at line 2185 of file adjoint.cpp.

                                                    {
  auto ptr = boost::make_shared<ObjectiveFunctionDataImpl>();
  CHK_THROW_MESSAGE(ptr->initPython(py_file), "init python");
  return ptr;
}

diff_symmetrize()

template<int DIM>

auto diff_symmetrize ( FTensor::Number< DIM > )

inline

Examples

mofem/tutorials/vec-7/adjoint.cpp.

Definition at line 1441 of file adjoint.cpp.

                                                                   {
 
  FTensor::Index<'i', DIM> i;
  FTensor::Index<'j', DIM> j;
  FTensor::Index<'k', DIM> k;
  FTensor::Index<'l', DIM> l;
 
  FTensor::Tensor4<double, DIM, DIM, DIM, DIM> t_diff;
 
  t_diff(i, j, k, l) = 0;
  t_diff(0, 0, 0, 0) = 1;
  t_diff(1, 1, 1, 1) = 1;
 
  t_diff(1, 0, 1, 0) = 0.5;
  t_diff(1, 0, 0, 1) = 0.5;
 
  t_diff(0, 1, 0, 1) = 0.5;
  t_diff(0, 1, 1, 0) = 0.5;
 
  if constexpr (DIM == 3) {
    t_diff(2, 2, 2, 2) = 1;
 
    t_diff(2, 0, 2, 0) = 0.5;
    t_diff(2, 0, 0, 2) = 0.5;
    t_diff(0, 2, 0, 2) = 0.5;
    t_diff(0, 2, 2, 0) = 0.5;
 
    t_diff(2, 1, 2, 1) = 0.5;
    t_diff(2, 1, 1, 2) = 0.5;
    t_diff(1, 2, 1, 2) = 0.5;
    t_diff(1, 2, 2, 1) = 0.5;
  }
 
  return t_diff;
};

get_range_from_block()

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

Examples

EshelbianPlasticity.cpp, and mofem/tutorials/vec-7/adjoint.cpp.

Definition at line 2579 of file adjoint.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]

[Example]

Definition at line 2036 of file adjoint.cpp.

                                 {
 
  Py_Initialize();
  np::initialize();
 
  // Initialisation of 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 2618 of file adjoint.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

AssemblyType A = AssemblyType::PETSC

constexpr

[Define dimension]

Definition at line 26 of file adjoint.cpp.

BASE_DIM

int BASE_DIM = 1

constexpr

Definition at line 19 of file adjoint.cpp.

bulk_modulus_K

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

constexpr

Examples

dynamic_first_order_con_law.cpp, and eigen_elastic.cpp.

Definition at line 75 of file adjoint.cpp.

help

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

static

[calculateGradient]

Definition at line 2034 of file adjoint.cpp.

I

IntegrationType I

constexpr

Initial value:

=
    IntegrationType::GAUSS

Examples

mofem/tutorials/vec-7/adjoint.cpp.

Definition at line 27 of file adjoint.cpp.

is_plane_strain

PetscBool is_plane_strain = PETSC_FALSE

Examples

mofem/tutorials/vec-7/adjoint.cpp.

Definition at line 78 of file adjoint.cpp.

poisson_ratio

double poisson_ratio = 0.3

constexpr

Examples

eigen_elastic.cpp, and plastic.cpp.

Definition at line 74 of file adjoint.cpp.

shear_modulus_G

double shear_modulus_G = young_modulus / (2 * (1 + poisson_ratio))

constexpr

Examples

dynamic_first_order_con_law.cpp, and eigen_elastic.cpp.

Definition at line 76 of file adjoint.cpp.

SPACE_DIM

int SPACE_DIM

constexpr

Initial value:

=
    EXECUTABLE_DIMENSION

[Define dimension]

Examples

dynamic_first_order_con_law.cpp, eigen_elastic.cpp, mofem/tutorials/vec-7/adjoint.cpp, plastic.cpp, and plot_base.cpp.

Definition at line 22 of file adjoint.cpp.

young_modulus

double young_modulus = 1

constexpr

Examples

eigen_elastic.cpp, and plastic.cpp.

Definition at line 73 of file adjoint.cpp.