ObjectiveFunctionDataImpl Struct Reference

Inheritance diagram for ObjectiveFunctionDataImpl:

Collaboration diagram for ObjectiveFunctionDataImpl:

Public Member Functions
	ObjectiveFunctionDataImpl ()=default
virtual	~ObjectiveFunctionDataImpl ()=default
MoFEMErrorCode	initPython (const std::string py_file)
MoFEMErrorCode	evalObjectiveFunction (MatrixDouble &coords, boost::shared_ptr< MatrixDouble > u_ptr, boost::shared_ptr< MatrixDouble > stress_ptr, boost::shared_ptr< MatrixDouble > strain_ptr, boost::shared_ptr< VectorDouble > o_ptr)
MoFEMErrorCode	evalObjectiveGradientStress (MatrixDouble &coords, boost::shared_ptr< MatrixDouble > u_ptr, boost::shared_ptr< MatrixDouble > stress_ptr, boost::shared_ptr< MatrixDouble > strain_ptr, boost::shared_ptr< MatrixDouble > o_ptr)
MoFEMErrorCode	evalObjectiveGradientStrain (MatrixDouble &coords, boost::shared_ptr< MatrixDouble > u_ptr, boost::shared_ptr< MatrixDouble > stress_ptr, boost::shared_ptr< MatrixDouble > strain_ptr, boost::shared_ptr< MatrixDouble > o_ptr)
MoFEMErrorCode	evalObjectiveGradientU (MatrixDouble &coords, boost::shared_ptr< MatrixDouble > u_ptr, boost::shared_ptr< MatrixDouble > stress_ptr, boost::shared_ptr< MatrixDouble > strain_ptr, boost::shared_ptr< MatrixDouble > o_ptr)
MoFEMErrorCode	numberOfModes (int block_id, int &modes)
MoFEMErrorCode	blockModes (int block_id, MatrixDouble &coords, std::array< double, 3 > &centroid, std::array< double, 6 > &bbodx, MatrixDouble &o_ptr)
Public Member Functions inherited from ObjectiveFunctionData
virtual	~ObjectiveFunctionData ()=default

Private Member Functions
MoFEMErrorCode	objectiveFunctionImpl (np::ndarray coords, np::ndarray u, np::ndarray stress, np::ndarray strain, np::ndarray &o)
MoFEMErrorCode	objectiveGradientStressImpl (np::ndarray coords, np::ndarray u, np::ndarray stress, np::ndarray strain, np::ndarray &o)
MoFEMErrorCode	objectiveGradientStrainImpl (np::ndarray coords, np::ndarray u, np::ndarray stress, np::ndarray strain, np::ndarray &o)
MoFEMErrorCode	objectiveGradientUImpl (np::ndarray coords, np::ndarray u, np::ndarray stress, np::ndarray strain, np::ndarray &o)
MoFEMErrorCode	blockModesImpl (int block_id, np::ndarray coords, np::ndarray centroid, np::ndarray bbodx, np::ndarray &o_ptr)
np::ndarray	convertToNumPy (std::vector< double > &data, int rows, int nb_gauss_pts)
	Converts a std::vector<double> to a NumPy ndarray.
np::ndarray	convertToNumPy (double *ptr, int s)
MatrixDouble	copyToFull (MatrixDouble &s)
void	copyToSymmetric (double *ptr, MatrixDouble &s)

Private Attributes
bp::object	mainNamespace
bp::object	objectiveFunction
bp::object	objectiveGradientStress
bp::object	objectiveGradientStrain
bp::object	objectiveGradientU
bp::object	objectNumberOfModes
bp::object	objectBlockModes

Detailed Description

Examples

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

Definition at line 2087 of file adjoint.cpp.

Constructor & Destructor Documentation

ObjectiveFunctionDataImpl()

ObjectiveFunctionDataImpl::ObjectiveFunctionDataImpl ( )

default

Examples

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

~ObjectiveFunctionDataImpl()

virtual ObjectiveFunctionDataImpl::~ObjectiveFunctionDataImpl ( )

virtualdefault

Examples

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

Member Function Documentation

blockModes()

MoFEMErrorCode ObjectiveFunctionDataImpl::blockModes ( int block_id, MatrixDouble & coords, std::array< double, 3 > & centroid, std::array< double, 6 > & bbodx, MatrixDouble & o_ptr )

virtual

Implements ObjectiveFunctionData.

Examples

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

Definition at line 2391 of file adjoint.cpp.

                                                     {
  MoFEMFunctionBegin;
  try {
 
    int nb_modes = bp::extract<int>(objectNumberOfModes(block_id));
 
    auto np_coords =
        convertToNumPy(coords.data(), coords.size1(), coords.size2());
 
    auto np_centroid = convertToNumPy(centroid.data(), 3);
    auto np_bbodx = convertToNumPy(bbodx.data(), 6);
 
    np::ndarray np_output =
        np::empty(bp::make_tuple(coords.size1(), nb_modes, 3),
                  np::dtype::get_builtin<double>());
    CHKERR blockModesImpl(block_id, np_coords, np_centroid, np_bbodx,
                          np_output);
 
    o_ptr.resize(nb_modes, coords.size1() * coords.size2(), false);
    double *val_ptr = reinterpret_cast<double *>(np_output.get_data());
    std::copy(val_ptr, val_ptr + coords.size1() * coords.size2() * nb_modes,
              o_ptr.data().begin());
 
  } catch (bp::error_already_set const &) {
    PyErr_Print();
    CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL, "Python error");
  }
  MoFEMFunctionReturn(0);
}

blockModesImpl()

MoFEMErrorCode ObjectiveFunctionDataImpl::blockModesImpl ( int block_id, np::ndarray coords, np::ndarray centroid, np::ndarray bbodx, np::ndarray & o_ptr )

private

Examples

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

Definition at line 2524 of file adjoint.cpp.

                                                                       {
  MoFEMFunctionBegin;
  try {
    // call python function
    o = bp::extract<np::ndarray>(
        objectBlockModes(block_id, coords, centroid, bbodx));
  } catch (bp::error_already_set const &) {
    // print all other errors to stderr
    PyErr_Print();
    CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL, "Python error");
  }
  MoFEMFunctionReturn(0);
}

convertToNumPy() [1/2]

np::ndarray ObjectiveFunctionDataImpl::convertToNumPy ( double * ptr, int s )

inlineprivate

Definition at line 2571 of file adjoint.cpp.

                                                                    {
  auto dtype = np::dtype::get_builtin<double>();
  auto size = bp::make_tuple(s);
  auto stride = bp::make_tuple(sizeof(double));
  return (np::from_data(ptr, dtype, size, stride, bp::object()));
}

convertToNumPy() [2/2]

np::ndarray ObjectiveFunctionDataImpl::convertToNumPy ( std::vector< double > & data, int rows, int nb_gauss_pts )

inlineprivate

Converts a std::vector<double> to a NumPy ndarray.

This function wraps the given vector data into a NumPy array with the specified number of rows and Gauss points. The resulting ndarray shares memory with the input vector, so changes to one will affect the other.

Parameters

data	Reference to the vector containing double values to be converted.
rows	Number of rows in the resulting NumPy array.
nb_gauss_pts	Number of Gauss points (columns) in the resulting NumPy array.

Returns

np::ndarray NumPy array view of the input data.

Note

• size specifies the shape of the resulting ndarray as a tuple (rows, nb_gauss_pts).
• stride specifies the step size in bytes to move to the next element in memory. Here, it is set to sizeof(double), indicating contiguous storage for each element.

Examples

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

Definition at line 2563 of file adjoint.cpp.

                                                            {
  auto dtype = np::dtype::get_builtin<double>();
  auto size = bp::make_tuple(rows, nb_gauss_pts);
  auto stride = bp::make_tuple(nb_gauss_pts * sizeof(double), sizeof(double));
  return (np::from_data(data.data(), dtype, size, stride, bp::object()));
}

copyToFull()

MatrixDouble ObjectiveFunctionDataImpl::copyToFull ( MatrixDouble & s )

private

Examples

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

Definition at line 2216 of file adjoint.cpp.

                                                                  {
  MatrixDouble f(9, s.size2());
  f.clear();
  FTENSOR_INDEX(SPACE_DIM, i);
  FTENSOR_INDEX(SPACE_DIM, j);
  auto t_f = getFTensor2FromMat<3, 3>(f);
  auto t_s = getFTensor2SymmetricFromMat<SPACE_DIM>(s);
  for (int ii = 0; ii != s.size2(); ++ii) {
    t_f(i, j) = t_s(i, j);
    ++t_f;
    ++t_s;
  }
  return f;
};

copyToSymmetric()

void ObjectiveFunctionDataImpl::copyToSymmetric ( double * ptr, MatrixDouble & s )

private

Examples

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

Definition at line 2231 of file adjoint.cpp.

                                                                            {
  FTENSOR_INDEX(SPACE_DIM, i);
  FTENSOR_INDEX(SPACE_DIM, j);
  auto t_f = FTensor::Tensor2<FTensor::PackPtr<double *, 1>, 3, 3>{
 
      ptr + 0 * s.size2(), ptr + 1 * s.size2(),
      ptr + 2 * s.size2(), ptr + 3 * s.size2(),
      ptr + 4 * s.size2(), ptr + 5 * s.size2(),
      ptr + 6 * s.size2(), ptr + 7 * s.size2(),
      ptr + 8 * s.size2()
 
  };
  auto t_s = getFTensor2SymmetricFromMat<SPACE_DIM>(s);
  for (int ii = 0; ii != s.size2(); ++ii) {
    t_s(i, j) = (t_f(i, j) || t_f(j, i)) / 2.0;
    ++t_f;
    ++t_s;
  }
}

evalObjectiveFunction()

MoFEMErrorCode ObjectiveFunctionDataImpl::evalObjectiveFunction ( MatrixDouble & coords, boost::shared_ptr< MatrixDouble > u_ptr, boost::shared_ptr< MatrixDouble > stress_ptr, boost::shared_ptr< MatrixDouble > strain_ptr, boost::shared_ptr< VectorDouble > o_ptr )

virtual

Implements ObjectiveFunctionData.

Examples

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

Definition at line 2251 of file adjoint.cpp.

                                         {
  MoFEMFunctionBegin;
  try {
 
    auto np_coords =
        convertToNumPy(coords.data(), coords.size1(), coords.size2());
    auto np_u = convertToNumPy(u_ptr->data(), u_ptr->size1(), u_ptr->size2());
 
    auto full_stress = copyToFull(*(stress_ptr));
    auto full_strain = copyToFull(*(strain_ptr));
 
    auto np_stress = convertToNumPy(full_stress.data(), full_stress.size1(),
                                    full_stress.size2());
    auto np_strain = convertToNumPy(full_strain.data(), full_strain.size1(),
                                    full_strain.size2());
    np::ndarray np_output = np::empty(bp::make_tuple(strain_ptr->size2()),
                                      np::dtype::get_builtin<double>());
    CHKERR objectiveFunctionImpl(np_coords, np_u, np_stress, np_strain,
                                 np_output);
    o_ptr->resize(stress_ptr->size2(), false);
    double *val_ptr = reinterpret_cast<double *>(np_output.get_data());
    std::copy(val_ptr, val_ptr + strain_ptr->size2(), o_ptr->data().begin());
 
  } catch (bp::error_already_set const &) {
    PyErr_Print();
    CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL, "Python error");
  }
  MoFEMFunctionReturn(0);
}

evalObjectiveGradientStrain()

MoFEMErrorCode ObjectiveFunctionDataImpl::evalObjectiveGradientStrain ( MatrixDouble & coords, boost::shared_ptr< MatrixDouble > u_ptr, boost::shared_ptr< MatrixDouble > stress_ptr, boost::shared_ptr< MatrixDouble > strain_ptr, boost::shared_ptr< MatrixDouble > o_ptr )

virtual

Implements ObjectiveFunctionData.

Examples

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

Definition at line 2320 of file adjoint.cpp.

                                         {
  MoFEMFunctionBegin;
  try {
 
    auto np_coords =
        convertToNumPy(coords.data(), coords.size1(), coords.size2());
    auto np_u = convertToNumPy(u_ptr->data(), u_ptr->size1(), u_ptr->size2());
 
    auto full_stress = copyToFull(*(stress_ptr));
    auto full_strain = copyToFull(*(strain_ptr));
    auto np_stress = convertToNumPy(full_stress.data(), full_stress.size1(),
                                    full_stress.size2());
    auto np_strain = convertToNumPy(full_strain.data(), full_strain.size1(),
                                    full_strain.size2());
 
    np::ndarray np_output =
        np::empty(bp::make_tuple(full_strain.size1(), full_strain.size2()),
                  np::dtype::get_builtin<double>());
    CHKERR objectiveGradientStrainImpl(np_coords, np_u, np_stress, np_strain,
                                       np_output);
    o_ptr->resize(strain_ptr->size1(), strain_ptr->size2(), false);
    double *val_ptr = reinterpret_cast<double *>(np_output.get_data());
    copyToSymmetric(val_ptr, *(o_ptr));
 
  } catch (bp::error_already_set const &) {
    PyErr_Print();
    CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL, "Python error");
  }
  MoFEMFunctionReturn(0);
}

evalObjectiveGradientStress()

MoFEMErrorCode ObjectiveFunctionDataImpl::evalObjectiveGradientStress ( MatrixDouble & coords, boost::shared_ptr< MatrixDouble > u_ptr, boost::shared_ptr< MatrixDouble > stress_ptr, boost::shared_ptr< MatrixDouble > strain_ptr, boost::shared_ptr< MatrixDouble > o_ptr )

virtual

Implements ObjectiveFunctionData.

Examples

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

Definition at line 2285 of file adjoint.cpp.

                                         {
  MoFEMFunctionBegin;
  try {
 
    auto np_coords =
        convertToNumPy(coords.data(), coords.size1(), coords.size2());
    auto np_u = convertToNumPy(u_ptr->data(), u_ptr->size1(), u_ptr->size2());
 
    auto full_stress = copyToFull(*(stress_ptr));
    auto full_strain = copyToFull(*(strain_ptr));
 
    auto np_stress = convertToNumPy(full_stress.data(), full_stress.size1(),
                                    full_stress.size2());
    auto np_strain = convertToNumPy(full_strain.data(), full_strain.size1(),
                                    full_strain.size2());
    np::ndarray np_output =
        np::empty(bp::make_tuple(full_strain.size1(), full_strain.size2()),
                  np::dtype::get_builtin<double>());
    CHKERR objectiveGradientStressImpl(np_coords, np_u, np_stress, np_strain,
                                       np_output);
    o_ptr->resize(stress_ptr->size1(), stress_ptr->size2(), false);
    double *val_ptr = reinterpret_cast<double *>(np_output.get_data());
    copyToSymmetric(val_ptr, *(o_ptr));
 
  } catch (bp::error_already_set const &) {
    PyErr_Print();
    CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL, "Python error");
  }
  MoFEMFunctionReturn(0);
}

evalObjectiveGradientU()

MoFEMErrorCode ObjectiveFunctionDataImpl::evalObjectiveGradientU	(	MatrixDouble &	coords,
		boost::shared_ptr< MatrixDouble >	u_ptr,
		boost::shared_ptr< MatrixDouble >	stress_ptr,
		boost::shared_ptr< MatrixDouble >	strain_ptr,
		boost::shared_ptr< MatrixDouble >	o_ptr )

Examples

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

Definition at line 2355 of file adjoint.cpp.

                                         {
  MoFEMFunctionBegin;
  try {
 
    auto np_coords =
        convertToNumPy(coords.data(), coords.size1(), coords.size2());
    auto np_u = convertToNumPy(u_ptr->data(), u_ptr->size1(), u_ptr->size2());
 
    auto full_stress = copyToFull(*(stress_ptr));
    auto full_strain = copyToFull(*(strain_ptr));
    auto np_stress = convertToNumPy(full_stress.data(), full_stress.size1(),
                                    full_stress.size2());
    auto np_strain = convertToNumPy(full_strain.data(), full_strain.size1(),
                                    full_strain.size2());
 
    np::ndarray np_output =
        np::empty(bp::make_tuple(u_ptr->size1(), u_ptr->size2()),
                  np::dtype::get_builtin<double>());
    CHKERR objectiveGradientStrainImpl(np_coords, np_u, np_stress, np_strain,
                                       np_output);
    o_ptr->resize(u_ptr->size1(), u_ptr->size2(), false);
    double *val_ptr = reinterpret_cast<double *>(np_output.get_data());
    std::copy(val_ptr, val_ptr + u_ptr->size1() * u_ptr->size2(),
              o_ptr->data().begin());
 
  } catch (bp::error_already_set const &) {
    PyErr_Print();
    CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL, "Python error");
  }
  MoFEMFunctionReturn(0);
}

initPython()

MoFEMErrorCode ObjectiveFunctionDataImpl::initPython

(

const std::string

py_file

)

Examples

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

Definition at line 2192 of file adjoint.cpp.

                                                             {
  MoFEMFunctionBegin;
  try {
 
    // create main module
    auto main_module = bp::import("__main__");
    mainNamespace = main_module.attr("__dict__");
    bp::exec_file(py_file.c_str(), mainNamespace, mainNamespace);
    // create a reference to python function
    objectiveFunction = mainNamespace["f"];
    objectiveGradientStress = mainNamespace["f_stress"];
    objectiveGradientStrain = mainNamespace["f_strain"];
    objectiveGradientU = mainNamespace["f_u"];
    objectNumberOfModes = mainNamespace["number_of_modes"];
    objectBlockModes = mainNamespace["block_modes"];
 
  } catch (bp::error_already_set const &) {
    // print all other errors to stderr
    PyErr_Print();
    CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL, "Python error");
  }
  MoFEMFunctionReturn(0);
}

numberOfModes()

MoFEMErrorCode ObjectiveFunctionDataImpl::numberOfModes ( int block_id, int & modes )

virtual

Implements ObjectiveFunctionData.

Examples

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

Definition at line 2509 of file adjoint.cpp.

                                                                    {
  MoFEMFunctionBegin;
  try {
 
    modes = bp::extract<int>(objectNumberOfModes(block_id));
 
  } catch (bp::error_already_set const &) {
    // print all other errors to stderr
    PyErr_Print();
    CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL, "Python error");
  }
  MoFEMFunctionReturn(0);
}

objectiveFunctionImpl()

MoFEMErrorCode ObjectiveFunctionDataImpl::objectiveFunctionImpl ( np::ndarray coords, np::ndarray u, np::ndarray stress, np::ndarray strain, np::ndarray & o )

private

Examples

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

Definition at line 2423 of file adjoint.cpp.

  {
  MoFEMFunctionBegin;
  try {
 
    // call python function
    o = bp::extract<np::ndarray>(objectiveFunction(coords, u, stress, strain));
 
  } catch (bp::error_already_set const &) {
    // print all other errors to stderr
    PyErr_Print();
    CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL, "Python error");
  }
  MoFEMFunctionReturn(0);
}

objectiveGradientStrainImpl()

MoFEMErrorCode ObjectiveFunctionDataImpl::objectiveGradientStrainImpl ( np::ndarray coords, np::ndarray u, np::ndarray stress, np::ndarray strain, np::ndarray & o )

private

Examples

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

Definition at line 2466 of file adjoint.cpp.

  {
  MoFEMFunctionBegin;
  try {
 
    // call python function
    o = bp::extract<np::ndarray>(
        objectiveGradientStrain(coords, u, stress, strain));
 
  } catch (bp::error_already_set const &) {
    // print all other errors to stderr
    PyErr_Print();
    CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL, "Python error");
  }
  MoFEMFunctionReturn(0);
}

objectiveGradientStressImpl()

MoFEMErrorCode ObjectiveFunctionDataImpl::objectiveGradientStressImpl ( np::ndarray coords, np::ndarray u, np::ndarray stress, np::ndarray strain, np::ndarray & o )

private

Examples

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

Definition at line 2444 of file adjoint.cpp.

  {
  MoFEMFunctionBegin;
  try {
 
    // call python function
    o = bp::extract<np::ndarray>(
        objectiveGradientStress(coords, u, stress, strain));
 
  } catch (bp::error_already_set const &) {
    // print all other errors to stderr
    PyErr_Print();
    CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL, "Python error");
  }
  MoFEMFunctionReturn(0);
}

objectiveGradientUImpl()

MoFEMErrorCode ObjectiveFunctionDataImpl::objectiveGradientUImpl ( np::ndarray coords, np::ndarray u, np::ndarray stress, np::ndarray strain, np::ndarray & o )

private

Examples

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

Definition at line 2488 of file adjoint.cpp.

  {
  MoFEMFunctionBegin;
  try {
 
    // call python function
    o = bp::extract<np::ndarray>(objectiveGradientU(coords, u, stress, strain));
 
  } catch (bp::error_already_set const &) {
    // print all other errors to stderr
    PyErr_Print();
    CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL, "Python error");
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

mainNamespace

bp::object ObjectiveFunctionDataImpl::mainNamespace

private

Examples

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

Definition at line 2128 of file adjoint.cpp.

objectBlockModes

bp::object ObjectiveFunctionDataImpl::objectBlockModes

private

Examples

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

Definition at line 2134 of file adjoint.cpp.

objectiveFunction

bp::object ObjectiveFunctionDataImpl::objectiveFunction

private

Examples

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

Definition at line 2129 of file adjoint.cpp.

objectiveGradientStrain

bp::object ObjectiveFunctionDataImpl::objectiveGradientStrain

private

Examples

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

Definition at line 2131 of file adjoint.cpp.

objectiveGradientStress

bp::object ObjectiveFunctionDataImpl::objectiveGradientStress

private

Examples

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

Definition at line 2130 of file adjoint.cpp.

objectiveGradientU

bp::object ObjectiveFunctionDataImpl::objectiveGradientU

private

Examples

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

Definition at line 2132 of file adjoint.cpp.

objectNumberOfModes

bp::object ObjectiveFunctionDataImpl::objectNumberOfModes

private

Examples

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

Definition at line 2133 of file adjoint.cpp.

---

The documentation for this struct was generated from the following file:

• tutorials/vec-7/adjoint.cpp