AdolCMetaElastic.cpp

Go to the documentation of this file.

/**
 * @file AdolCMetaElastic.cpp
 * @brief 
 * @version 0.1
 * @date 2026-03-28
 * 
 * @copyright Copyright (c) 2026
 * 
 */
 
 namespace AdolCOps {
 
template <int DIM> struct AdolCMetaElastic : public AdolCElasticImpl<DIM> {
  using AdolCElasticImpl<DIM>::AdolCElasticImpl;
 
  using A = AdolCElasticImpl<DIM>;
 
  MoFEMErrorCode getOptions(MoFEM::Interface *m_field_ptr = nullptr) override {
    MoFEMFunctionBegin;
    for (auto &m : metaElasticAtPts) {
      CHKERR m.second->getOptions(m_field_ptr);
    }
    A::tagVsRangePtr = boost::make_shared<std::map<int, Range>>();
    for (const auto &m : metaElasticAtPts) {
      A::tagVsRangePtr->insert(
          {m.first, Range()}); // Initialize with empty range
      for (auto &rp : A::paramVecByRange) {
        A::tagVsRangePtr->at(m.first).merge(
            rp.first); // Merge ranges for the same tag
      }
    }
    MoFEMFunctionReturn(0);
  }
 
  MoFEMErrorCode setParams(EntityHandle ent, int gg) override {
    MoFEMFunctionBegin;
    for(auto &m : metaElasticAtPts) {
      CHKERR m.second->setParams(ent, gg);
    }
    MoFEMFunctionReturn(0);
  }
 
  MoFEMErrorCode recordTape() override {
    MoFEMFunctionBegin;
    for (auto &m : metaElasticAtPts) {
      CHKERR m.second->recordTape();
    }
    MoFEMFunctionReturn(0);
  }
 
protected:
  static inline constexpr int modelType =
      (DIM == 3) ? AdolCOps::MODEL_3D : AdolCOps::MODEL_2D_PLANE_STRAIN;
 
  std::map<int, boost::shared_ptr<PhysicalEquations>> metaElasticAtPts = {
 
      {
 
          // First tag is always MetaElastic tag
          AdolCOps::AdolCTagsRegistry::setTagName("MetaElastic"),
          createAdolCPhysicalEquationsPtr<ELASTICITY::NEOHOOKEAN, modelType>(
              this->adolcDataPtr,
              AdolCOps::AdolCTagsRegistry::setTagName("MetaElastic"))
 
      }
 
  };
};
 
template <>
boost::shared_ptr<PhysicalEquations>
createAdolCPhysicalEquationsPtr<ELASTICITY::META, MODEL_3D>(
    boost::shared_ptr<ADolCData> adolc_data_ptr, int tag) {
  return boost::make_shared<AdolCMetaElastic<3>>(adolc_data_ptr, tag);
}
 
template <>
boost::shared_ptr<PhysicalEquations>
createAdolCPhysicalEquationsPtr<ELASTICITY::META, MODEL_2D_PLANE_STRAIN>(
    boost::shared_ptr<ADolCData> adolc_data_ptr, int tag) {
  return boost::make_shared<AdolCMetaElastic<2>>(adolc_data_ptr, tag);
}
}