MatOps::MatOpsDataImpl Struct Reference

Inheritance diagram for MatOps::MatOpsDataImpl:

Collaboration diagram for MatOps::MatOpsDataImpl:

Classes
struct	Data

Public Types
using	MatrixPtr = MatOpsData::MatrixPtr
using	StateMatrixPtr = MatOpsData::StateMatrixPtr
using	DataContainer = boost::multi_index::multi_index_container< Data, boost::multi_index::indexed_by< boost::multi_index::sequenced< boost::multi_index::tag< DataSequence_mi_tag > >, boost::multi_index::hashed_unique< boost::multi_index::tag< DataName_mi_tag >, boost::multi_index::member< Data, std::string, &Data::name > > > >
Public Types inherited from MatOps::MatOpsData
using	MatrixPtr = boost::shared_ptr< MatrixDouble >
using	StateMatrixPtr = boost::shared_ptr< MatrixAdaptor >

Public Member Functions
auto &	getDataCommonByName ()
auto &	getDataActiveByName ()
auto &	getDataDependentByName ()
auto &	getDataDependentDerivativesByName ()
auto &	getStateDataByName ()
auto &	getStateDataBySequence ()
std::pair< bool, MatrixPtr >	insertCommonData (const std::string &name, const MatrixDouble &data=MatrixDouble(), int shift=0) override
std::pair< bool, MatrixPtr >	insertActiveData (const std::string &name, const MatrixDouble &data=MatrixDouble(), int shift=0) override
std::pair< bool, MatrixPtr >	insertDependentData (const std::string &name, const MatrixDouble &data=MatrixDouble(), int shift=0) override
std::pair< bool, MatrixPtr >	insertDependentDerivativesData (const std::string &name, const MatrixDouble &data=MatrixDouble(), int shift=0) override
MatrixPtr	getCommonDataPtr (const std::string &name) override
MatrixPtr	getActiveDataPtr (const std::string &name) override
MatrixPtr	getDependentDataPtr (const std::string &name) override
MatrixPtr	getDependentDerivativesDataPtr (const std::string &name) override
StateMatrixPtr	getStateDataPtr (const std::string &name, EntityHandle ent, int gg) override
StateTags	getStateTags () const override
MoFEMErrorCode	bindStateTag (MoFEM::Interface &m_field, const std::string &name, int num_components) override
MoFEMErrorCode	setupStateData () override
MoFEMErrorCode	setActiveContinuousVector (std::vector< double > &active_variables) override
MoFEMErrorCode	setDependentContinuousVector (std::vector< double > &dependent_variables) override
MoFEMErrorCode	setDependentDerivativesContinuousVector (std::vector< double > &dependent_variables_derivatives) override
MoFEMErrorCode	getActiveContinuousVector (const std::vector< double > &active_variables) override
MoFEMErrorCode	getDependentContinuousVector (const std::vector< double > &dependent_variables) override
MoFEMErrorCode	getDependentDerivativesContinuousVector (const std::vector< double > &dependent_variables_derivatives) override
Public Member Functions inherited from MatOps::MatOpsData
virtual	~MatOpsData ()=default

Protected Member Functions
MoFEMErrorCode	setContinuousVector (DataContainer &data_container, std::vector< double > &variables)
MoFEMErrorCode	getContinuousVector (DataContainer &data_container, const std::vector< double > &variables)

Protected Attributes
DataContainer	dataCommon
DataContainer	dataActive
DataContainer	dataDependent
DataContainer	dataDependentDerivatives
DataContainer	dataState
MoFEM::Interface *	mFieldPtr = nullptr

Detailed Description

Definition at line 92 of file MatOps.cpp.

Member Typedef Documentation

DataContainer

using MatOps::MatOpsDataImpl::DataContainer = boost::multi_index::multi_index_container< Data, boost::multi_index::indexed_by< boost::multi_index::sequenced< boost::multi_index::tag<DataSequence_mi_tag> >, boost::multi_index::hashed_unique< boost::multi_index::tag<DataName_mi_tag>, boost::multi_index::member<Data, std::string, &Data::name> >> >

Definition at line 115 of file MatOps.cpp.

MatrixPtr

using MatOps::MatOpsDataImpl::MatrixPtr = MatOpsData::MatrixPtr

Definition at line 93 of file MatOps.cpp.

StateMatrixPtr

using MatOps::MatOpsDataImpl::StateMatrixPtr = MatOpsData::StateMatrixPtr

Definition at line 94 of file MatOps.cpp.

Member Function Documentation

bindStateTag()

MoFEMErrorCode MatOps::MatOpsDataImpl::bindStateTag ( MoFEM::Interface &  m_field, const std::string &  name, int  num_components  )

inlineoverridevirtual

Implements MatOps::MatOpsData.

Definition at line 324 of file MatOps.cpp.

                                                           {
    MoFEMFunctionBegin;
    mFieldPtr = &m_field;
 
    if (name.empty()) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
              "State data name can not be empty");
    }
 
    if (num_components <= 0) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
              "State data %s has non-positive number of components %d",
              name.c_str(), num_components);
    }
 
    auto &state_by_name = getStateDataByName();
    auto it = state_by_name.find(name);
    if (it != state_by_name.end()) {
      if (it->numComponents != num_components) {
        SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
                "State data %s already exists with incompatible size",
                name.c_str());
      }
    } else {
      state_by_name.insert(Data{name, num_components});
    }
 
    MoFEMFunctionReturn(0);
  }

getActiveContinuousVector()

MoFEMErrorCode MatOps::MatOpsDataImpl::getActiveContinuousVector ( const std::vector< double > &  active_variables )

inlineoverridevirtual

Implements MatOps::MatOpsData.

Definition at line 395 of file MatOps.cpp.

                                                          {
    return getContinuousVector(dataActive, active_variables);
  }

getActiveDataPtr()

MatrixPtr MatOps::MatOpsDataImpl::getActiveDataPtr ( const std::string &  name )

inlineoverridevirtual

Implements MatOps::MatOpsData.

Definition at line 204 of file MatOps.cpp.

                                                             {
    auto &data_by_name = getDataActiveByName();
    auto it = data_by_name.find(name);
    if (it == data_by_name.end()) {
      CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL,
                        "Active data with name " + name + " not found");
      return MatrixPtr();
    }
    return boost::shared_ptr<MatrixDouble>(this->shared_from_this(), &it->data);
  }

getCommonDataPtr()

MatrixPtr MatOps::MatOpsDataImpl::getCommonDataPtr ( const std::string &  name )

inlineoverridevirtual

Implements MatOps::MatOpsData.

Definition at line 193 of file MatOps.cpp.

                                                             {
    auto &data_by_name = getDataCommonByName();
    auto it = data_by_name.find(name);
    if (it == data_by_name.end()) {
      CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL,
                        "Common data with name " + name + " not found");
      return MatrixPtr();
    }
    return boost::shared_ptr<MatrixDouble>(this->shared_from_this(), &it->data);
  }

getContinuousVector()

MoFEMErrorCode MatOps::MatOpsDataImpl::getContinuousVector ( DataContainer &  data_container, const std::vector< double > &  variables  )

inlineprotected

Definition at line 436 of file MatOps.cpp.

                                                                         {
    MoFEMFunctionBeginHot;
 
    auto &data_by_sequence =
        boost::multi_index::get<DataSequence_mi_tag>(data_container);
 
    int nb_variables = 0;
    for (const auto &d : data_by_sequence)
      nb_variables += d.data.data().size();
 
    if (variables.size() != static_cast<size_t>(nb_variables)) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
              "Inconsistent size of continuous vector");
    }
 
    int shift = 0;
    for (const auto &d : data_by_sequence) {
      d.shift = shift;
      for (auto it = d.data.data().begin(); it != d.data.data().end();
           ++it, ++shift)
        *it = variables[shift];
    }
 
    MoFEMFunctionReturnHot(0);
  }

getDataActiveByName()

auto & MatOps::MatOpsDataImpl::getDataActiveByName ( )

inline

Definition at line 132 of file MatOps.cpp.

                              {
    return boost::multi_index::get<DataName_mi_tag>(dataActive);
  }

getDataCommonByName()

auto & MatOps::MatOpsDataImpl::getDataCommonByName ( )

inline

Definition at line 128 of file MatOps.cpp.

                              {
    return boost::multi_index::get<DataName_mi_tag>(dataCommon);
  }

getDataDependentByName()

auto & MatOps::MatOpsDataImpl::getDataDependentByName ( )

inline

Definition at line 136 of file MatOps.cpp.

                                 {
    return boost::multi_index::get<DataName_mi_tag>(dataDependent);
  }

getDataDependentDerivativesByName()

auto & MatOps::MatOpsDataImpl::getDataDependentDerivativesByName ( )

inline

Definition at line 140 of file MatOps.cpp.

                                            {
    return boost::multi_index::get<DataName_mi_tag>(dataDependentDerivatives);
  }

getDependentContinuousVector()

MoFEMErrorCode MatOps::MatOpsDataImpl::getDependentContinuousVector ( const std::vector< double > &  dependent_variables )

inlineoverridevirtual

Implements MatOps::MatOpsData.

Definition at line 400 of file MatOps.cpp.

                                                             {
    return getContinuousVector(dataDependent, dependent_variables);
  }

getDependentDataPtr()

MatrixPtr MatOps::MatOpsDataImpl::getDependentDataPtr ( const std::string &  name )

inlineoverridevirtual

Implements MatOps::MatOpsData.

Definition at line 215 of file MatOps.cpp.

                                                                {
    auto &data_by_name = getDataDependentByName();
    auto it = data_by_name.find(name);
    if (it == data_by_name.end()) {
      CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL,
                        "Dependent data with name " + name + " not found");
      return MatrixPtr();
    }
    return boost::shared_ptr<MatrixDouble>(this->shared_from_this(), &it->data);
  }

getDependentDerivativesContinuousVector()

MoFEMErrorCode MatOps::MatOpsDataImpl::getDependentDerivativesContinuousVector ( const std::vector< double > &  dependent_variables_derivatives )

inlineoverridevirtual

Implements MatOps::MatOpsData.

Definition at line 405 of file MatOps.cpp.

                                                                         {
    return getContinuousVector(dataDependentDerivatives,
                               dependent_variables_derivatives);
  }

getDependentDerivativesDataPtr()

MatrixPtr MatOps::MatOpsDataImpl::getDependentDerivativesDataPtr ( const std::string &  name )

inlineoverridevirtual

Implements MatOps::MatOpsData.

Definition at line 226 of file MatOps.cpp.

                                                                           {
    auto &data_by_name = getDataDependentDerivativesByName();
    auto it = data_by_name.find(name);
    if (it == data_by_name.end()) {
      CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL,
                        "Dependent derivatives data with name " + name +
                            " not found");
      return MatrixPtr();
    }
    return boost::shared_ptr<MatrixDouble>(this->shared_from_this(), &it->data);
  }

getStateDataByName()

auto & MatOps::MatOpsDataImpl::getStateDataByName ( )

inline

Definition at line 144 of file MatOps.cpp.

                             {
    return boost::multi_index::get<DataName_mi_tag>(dataState);
  }

getStateDataBySequence()

auto & MatOps::MatOpsDataImpl::getStateDataBySequence ( )

inline

Definition at line 148 of file MatOps.cpp.

                                 {
    return boost::multi_index::get<DataSequence_mi_tag>(dataState);
  }

getStateDataPtr()

StateMatrixPtr MatOps::MatOpsDataImpl::getStateDataPtr ( const std::string &  name, EntityHandle  ent, int  gg  )

inlineoverridevirtual

Implements MatOps::MatOpsData.

Definition at line 238 of file MatOps.cpp.

                                                  {
    auto &state_by_name = getStateDataByName();
    auto it = state_by_name.find(name);
    if (it == state_by_name.end()) {
      CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL,
                        "State data with name " + name + " not found");
      return StateMatrixPtr();
    }
    if (!mFieldPtr) {
      CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL,
                        "MoFEM interface for state data is not set");
      return StateMatrixPtr();
    }
    if (!ent) {
      CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL,
                        "State entity is not set");
      return StateMatrixPtr();
    }
    if (gg < 0) {
      CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL,
                        "State gauss point out of range");
      return StateMatrixPtr();
    }
    if (!it->tag) {
      CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL,
                        "State data " + name +
                            " is not initialised; call setupStateData() after binding");
      return StateMatrixPtr();
    }
 
    double *tag_data = nullptr;
    int tag_size = 0;
    auto &moab = mFieldPtr->get_moab();
    const auto required_size = it->numComponents * (gg + 1);
    auto rval = moab.tag_get_by_ptr(it->tag, &ent, 1,
                                    (const void **)&tag_data, &tag_size);
    if (tag_data && tag_size > 0 && tag_size % it->numComponents != 0) {
      CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL,
                        "Inconsistent state tag size for " + name);
      return StateMatrixPtr();
    }
    if (rval != MB_SUCCESS || !tag_data || tag_size < required_size) {
      std::vector<double> storage(required_size, 0);
      if (rval == MB_SUCCESS && tag_data && tag_size > 0) {
        std::copy_n(tag_data, std::min(tag_size, required_size), storage.data());
      }
      void const *storage_ptr[] = {storage.data()};
      int storage_size = required_size;
      CHK_THROW_MESSAGE(moab.tag_set_by_ptr(it->tag, &ent, 1, storage_ptr,
                                            &storage_size),
                        "Failed to set state data " + name);
      CHK_THROW_MESSAGE(moab.tag_get_by_ptr(it->tag, &ent, 1,
                                            (const void **)&tag_data, &tag_size),
                        "Failed to get state data " + name);
    }
    if (!tag_data || tag_size < required_size ||
        (tag_size > 0 && tag_size % it->numComponents != 0)) {
      CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL,
                        "Inconsistent state tag size for " + name);
      return StateMatrixPtr();
    }
 
    auto *ptr = tag_data + gg * it->numComponents;
    return boost::make_shared<MatrixAdaptor>(
        it->numComponents, 1,
        ublas::shallow_array_adaptor<double>(it->numComponents, ptr));
  }

getStateTags()

StateTags MatOps::MatOpsDataImpl::getStateTags ( ) const

inlineoverridevirtual

Implements MatOps::MatOpsData.

Definition at line 307 of file MatOps.cpp.

                                          {
    StateTags state_tags;
    const auto &state_by_sequence =
        boost::multi_index::get<DataSequence_mi_tag>(dataState);
 
    for (const auto &state_data : state_by_sequence) {
      if (!state_data.tag || state_data.name.empty() ||
          state_data.numComponents <= 0)
        continue;
      state_tags.push_back(
          {state_data.name, state_data.tag, state_data.numComponents});
    }
 
    sortStateTags(state_tags);
    return state_tags;
  }

insertActiveData()

std::pair< bool, MatrixPtr > MatOps::MatOpsDataImpl::insertActiveData ( const std::string &  name, const MatrixDouble &  data = MatrixDouble(), int  shift = 0  )

inlineoverridevirtual

Implements MatOps::MatOpsData.

Definition at line 163 of file MatOps.cpp.

                                           {
    auto insert_ret = getDataActiveByName().insert(Data{name, data, shift});
    return {insert_ret.second,
            boost::shared_ptr<MatrixDouble>(this->shared_from_this(),
                                            &insert_ret.first->data)};
  }

insertCommonData()

std::pair< bool, MatrixPtr > MatOps::MatOpsDataImpl::insertCommonData ( const std::string &  name, const MatrixDouble &  data = MatrixDouble(), int  shift = 0  )

inlineoverridevirtual

Implements MatOps::MatOpsData.

Definition at line 153 of file MatOps.cpp.

                                           {
    auto insert_ret = getDataCommonByName().insert(Data{name, data, shift});
    return {insert_ret.second,
            boost::shared_ptr<MatrixDouble>(this->shared_from_this(),
                                            &insert_ret.first->data)};
  }

insertDependentData()

std::pair< bool, MatrixPtr > MatOps::MatOpsDataImpl::insertDependentData ( const std::string &  name, const MatrixDouble &  data = MatrixDouble(), int  shift = 0  )

inlineoverridevirtual

Implements MatOps::MatOpsData.

Definition at line 173 of file MatOps.cpp.

                                              {
    auto insert_ret = getDataDependentByName().insert(Data{name, data, shift});
    return {insert_ret.second,
            boost::shared_ptr<MatrixDouble>(this->shared_from_this(),
                                            &insert_ret.first->data)};
  }

insertDependentDerivativesData()

std::pair< bool, MatrixPtr > MatOps::MatOpsDataImpl::insertDependentDerivativesData ( const std::string &  name, const MatrixDouble &  data = MatrixDouble(), int  shift = 0  )

inlineoverridevirtual

Implements MatOps::MatOpsData.

Definition at line 183 of file MatOps.cpp.

                                                         {
    auto insert_ret =
        getDataDependentDerivativesByName().insert(Data{name, data, shift});
    return {insert_ret.second,
            boost::shared_ptr<MatrixDouble>(this->shared_from_this(),
                                            &insert_ret.first->data)};
  }

setActiveContinuousVector()

MoFEMErrorCode MatOps::MatOpsDataImpl::setActiveContinuousVector ( std::vector< double > &  active_variables )

inlineoverridevirtual

Implements MatOps::MatOpsData.

Definition at line 380 of file MatOps.cpp.

                                                                          {
    return setContinuousVector(dataActive, active_variables);
  }

setContinuousVector()

MoFEMErrorCode MatOps::MatOpsDataImpl::setContinuousVector ( DataContainer &  data_container, std::vector< double > &  variables  )

inlineprotected

Definition at line 412 of file MatOps.cpp.

                                                                   {
    MoFEMFunctionBeginHot;
 
    auto &data_by_sequence =
        boost::multi_index::get<DataSequence_mi_tag>(data_container);
 
    int nb_variables = 0;
    for (const auto &d : data_by_sequence)
      nb_variables += d.data.data().size();
 
    variables.resize(nb_variables);
 
    int shift = 0;
    for (const auto &d : data_by_sequence) {
      d.shift = shift;
      for (auto it = d.data.data().begin(); it != d.data.data().end();
           ++it, ++shift)
        variables[shift] = *it;
    }
 
    MoFEMFunctionReturnHot(0);
  }

setDependentContinuousVector()

MoFEMErrorCode MatOps::MatOpsDataImpl::setDependentContinuousVector ( std::vector< double > &  dependent_variables )

inlineoverridevirtual

Implements MatOps::MatOpsData.

Definition at line 384 of file MatOps.cpp.

                                                       {
    return setContinuousVector(dataDependent, dependent_variables);
  }

setDependentDerivativesContinuousVector()

MoFEMErrorCode MatOps::MatOpsDataImpl::setDependentDerivativesContinuousVector ( std::vector< double > &  dependent_variables_derivatives )

inlineoverridevirtual

Implements MatOps::MatOpsData.

Definition at line 389 of file MatOps.cpp.

                                                                   {
    return setContinuousVector(dataDependentDerivatives,
                               dependent_variables_derivatives);
  }

setupStateData()

MoFEMErrorCode MatOps::MatOpsDataImpl::setupStateData ( )

inlineoverridevirtual

Implements MatOps::MatOpsData.

Definition at line 356 of file MatOps.cpp.

                                           {
    MoFEMFunctionBegin;
    if (!mFieldPtr) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
              "MoFEM interface for state data is not set");
    }
 
    for (auto &state_data : getStateDataBySequence()) {
      Tag tag = 0;
      const auto state_tag_name = "_ADOLC_STATE_" +
                                  std::to_string(state_data.numComponents) +
                                  "_" + state_data.name;
      CHKERR mFieldPtr->get_moab().tag_get_handle(
          state_tag_name.c_str(), 0, MB_TYPE_DOUBLE, tag,
          MB_TAG_CREAT | MB_TAG_VARLEN | MB_TAG_SPARSE, PETSC_NULLPTR);
      getStateDataBySequence().modify(getStateDataBySequence().iterator_to(
                                          state_data),
                                      [&](auto &data) { data.tag = tag; });
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

dataActive

DataContainer MatOps::MatOpsDataImpl::dataActive

protected

Definition at line 464 of file MatOps.cpp.

dataCommon

DataContainer MatOps::MatOpsDataImpl::dataCommon

protected

Definition at line 463 of file MatOps.cpp.

dataDependent

DataContainer MatOps::MatOpsDataImpl::dataDependent

protected

Definition at line 465 of file MatOps.cpp.

dataDependentDerivatives

DataContainer MatOps::MatOpsDataImpl::dataDependentDerivatives

protected

Definition at line 466 of file MatOps.cpp.

dataState

DataContainer MatOps::MatOpsDataImpl::dataState

protected

Definition at line 467 of file MatOps.cpp.

mFieldPtr

MoFEM::Interface* MatOps::MatOpsDataImpl::mFieldPtr = nullptr

protected

Definition at line 468 of file MatOps.cpp.

---

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

• src/materials/impl/MatOps.cpp