AdolCOps::AdolCHuHu< DIM > Struct Template Reference

Inheritance diagram for AdolCOps::AdolCHuHu< DIM >:

Collaboration diagram for AdolCOps::AdolCHuHu< DIM >:

Public Types
using	A = PhysicalEquations

Public Member Functions
ForcesAndSourcesCore::UserDataOperator *	createOp (boost::shared_ptr< PhysicalEquations > physical_ptr, bool eval_stress, bool eval_tangent) override
MoFEMErrorCode	getOptions (MoFEM::Interface *m_field_ptr=nullptr) override
MoFEMErrorCode	setParams (EntityHandle ent, int gg) override
MoFEMErrorCode	recordTape () override
ForcesAndSourcesCore::UserDataOperator *	createOp (boost::shared_ptr< PhysicalEquations > physical_ptr, bool eval_stress, bool eval_tangent)
ForcesAndSourcesCore::UserDataOperator *	createOp (boost::shared_ptr< PhysicalEquations > physical_ptr, bool eval_stress, bool eval_tangent)
	PhysicalEquations ()=delete
	PhysicalEquations (boost::shared_ptr< ADolCData > adolc_data_ptr, int tag, boost::shared_ptr< std::map< int, Range > > tag_vs_range_ptr=nullptr)
Public Member Functions inherited from AdolCOps::PhysicalEquations
	PhysicalEquations ()=delete
	PhysicalEquations (boost::shared_ptr< ADolCData > adolc_data_ptr, int tag, boost::shared_ptr< std::map< int, Range > > tag_vs_range_ptr=nullptr)
virtual	~PhysicalEquations ()=default
MoFEMErrorCode	setActiveContinuousVector ()
MoFEMErrorCode	setDependentContinuousVector ()
MoFEMErrorCode	setDependentDerivativesContinuousVector ()
MoFEMErrorCode	getActiveContinuousVector ()
MoFEMErrorCode	getDependentContinuousVector ()
MoFEMErrorCode	getDependentDerivativesContinuousVector ()

Protected Attributes
double	K = 1e-6
std::vector< double >	defaultMaterialParameters = {K}

Additional Inherited Members
Public Attributes inherited from AdolCOps::PhysicalEquations
int	tAg
boost::shared_ptr< std::map< int, Range > >	tagVsRangePtr
std::vector< std::pair< Range, std::vector< double > > >	paramVecByRange
boost::shared_ptr< ADolCData >	adolcDataPtr
std::vector< double >	activeVariables
std::vector< double >	dependentVariables
std::vector< double >	dependentVariablesDerivatives

Detailed Description

template<int DIM>
struct AdolCOps::AdolCHuHu< DIM >

Definition at line 21 of file AdolCHuHu.cpp.

Member Typedef Documentation

A

template<int DIM>

using AdolCOps::AdolCHuHu< DIM >::A = PhysicalEquations

Definition at line 24 of file AdolCHuHu.cpp.

Member Function Documentation

createOp() [1/3]

ForcesAndSourcesCore::UserDataOperator * AdolCOps::AdolCHuHu< 3 >::createOp ( boost::shared_ptr< PhysicalEquations >  physical_ptr, bool  eval_stress, bool  eval_tangent  )

virtual

Implements AdolCOps::PhysicalEquations.

Definition at line 307 of file AdolCHuHu.cpp.

                                                            {
  using EleOp = ForcesAndSourcesCore::UserDataOperator;
  return createOpImpl<3, EleOp>(physical_ptr, eval_stress, eval_tangent);
}

createOp() [2/3]

ForcesAndSourcesCore::UserDataOperator * AdolCOps::AdolCHuHu< 2 >::createOp ( boost::shared_ptr< PhysicalEquations >  physical_ptr, bool  eval_stress, bool  eval_tangent  )

virtual

Implements AdolCOps::PhysicalEquations.

Definition at line 315 of file AdolCHuHu.cpp.

                                                            {
  using EleOp = ForcesAndSourcesCore::UserDataOperator;
  return createOpImpl<2, EleOp>(physical_ptr, eval_stress, eval_tangent);
}

createOp() [3/3]

template<int DIM>

ForcesAndSourcesCore::UserDataOperator * AdolCOps::AdolCHuHu< DIM >::createOp ( boost::shared_ptr< PhysicalEquations >  physical_ptr, bool  eval_stress, bool  eval_tangent  )

overridevirtual

Implements AdolCOps::PhysicalEquations.

getOptions()

template<int DIM>

MoFEMErrorCode AdolCOps::AdolCHuHu< DIM >::getOptions ( MoFEM::Interface *  m_field_ptr = nullptr )

inlineoverridevirtual

Implements AdolCOps::PhysicalEquations.

Definition at line 30 of file AdolCHuHu.cpp.

                                                                            {
    MoFEMFunctionBegin;
    MOFEM_LOG_CHANNEL("WORLD");
 
    PetscOptionsBegin(PETSC_COMM_WORLD, "huhu_", "", "none");
    CHKERR PetscOptionsScalar("-K", "Bulk modulus K", "", K, &K, PETSC_NULLPTR);
    PetscOptionsEnd();
 
    MOFEM_TAG_AND_LOG("WORLD", Sev::inform, "Default HuHu parameters:")
        << " K = " << K;
    defaultMaterialParameters = {K};
 
    std::string block_name = "MAT_HUHU";
 
    for (auto &m :
 
         m_field_ptr->getInterface<MeshsetsManager>()->
 
         getCubitMeshsetPtr(
             std::regex((boost::format("%s(.*)") % block_name).str()))
 
    ) {
 
      std::vector<double> block_data;
      CHKERR m->getAttributes(block_data);
      if (block_data.size() < 1) {
        block_data.push_back(
            K); // Use default K if not provided in block attributes
      }
      auto get_block_ents = [&]() {
        Range ents;
        CHK_MOAB_THROW(m_field_ptr->get_moab().get_entities_by_handle(
                           m->meshset, ents, true),
                       "can not get block entities");
        return ents;
      };
 
      A::paramVecByRange.push_back({get_block_ents(), {block_data[0]}});
 
      MOFEM_TAG_AND_LOG("WORLD", Sev::inform, "MatBlock for HuHu")
          << *m << " K = " << block_data[0];
    }
 
    MoFEMFunctionReturn(0);
  };

PhysicalEquations() [1/2]

template<int DIM>

AdolCOps::PhysicalEquations::PhysicalEquations ( )

delete

PhysicalEquations() [2/2]

template<int DIM>

AdolCOps::PhysicalEquations::PhysicalEquations ( boost::shared_ptr< ADolCData >  adolc_data_ptr, int  tag, boost::shared_ptr< std::map< int, Range > >  tag_vs_range_ptr = nullptr  )

inline

Definition at line 126 of file AdolCOps.hpp.

      : adolcDataPtr(adolc_data_ptr), tAg(tag),
        tagVsRangePtr(tag_vs_range_ptr) {}

recordTape()

template<int DIM>

MoFEMErrorCode AdolCOps::AdolCHuHu< DIM >::recordTape ( )

inlineoverridevirtual

Implements AdolCOps::PhysicalEquations.

Definition at line 82 of file AdolCHuHu.cpp.

                                       {
    MoFEMFunctionBegin;
 
    A::adolcDataPtr->insertCommonData("grad", MatrixDouble(DIM * DIM, 1));
    A::adolcDataPtr->insertCommonData("k", MatrixDouble(1, 1));
    A::adolcDataPtr->insertCommonData("k_dF", MatrixDouble(DIM * DIM, 1));
 
    A::adolcDataPtr->insertActiveData("F", MatrixDouble(DIM, DIM));
    A::adolcDataPtr->insertDependentData("k", MatrixDouble(1, 1));
    A::adolcDataPtr->insertDependentDerivativesData("k_dF",
                                                    MatrixDouble(DIM, DIM));
 
    auto t_F = getFTensor2FromPtr<DIM, DIM>(
        A::adolcDataPtr->getActiveDataPtr("F")->data().data());
    auto k = (*A::adolcDataPtr->getDependentDataPtr("k"))(0, 0);
 
    auto t_kd = FTensor::Kronecker_Delta<int>();
 
    FTENSOR_INDEX(DIM, i);
    FTENSOR_INDEX(DIM, j);
    FTENSOR_INDEX(DIM, I);
    FTENSOR_INDEX(DIM, J);
 
    t_F(i, J) = t_kd(i, J);
 
    FTensor::Tensor2<adouble, DIM, DIM> ta_F;
    adouble ta_k;
 
    adouble det_aF;
    FTensor::Tensor2<adouble, DIM, DIM> ta_invF;
 
    trace_on(A::tAg);
    auto p_K = mkparam(K);
    ta_F(i, J) <<= t_F(i, J);
    ta_k = p_K;
    ta_k >>= k;
 
    trace_off();
 
    MoFEMFunctionReturn(0);
  }

setParams()

template<int DIM>

MoFEMErrorCode AdolCOps::AdolCHuHu< DIM >::setParams ( EntityHandle  ent, int  gg  )

inlineoverridevirtual

Implements AdolCOps::PhysicalEquations.

Definition at line 76 of file AdolCHuHu.cpp.

                                                              {
    set_param_vec(A::tAg, defaultMaterialParameters.size(),
                  defaultMaterialParameters.data());
    return 0;
  }

Member Data Documentation

defaultMaterialParameters

template<int DIM>

std::vector<double> AdolCOps::AdolCHuHu< DIM >::defaultMaterialParameters = {K}

protected

Definition at line 126 of file AdolCHuHu.cpp.

{K};

K

template<int DIM>

double AdolCOps::AdolCHuHu< DIM >::K = 1e-6

protected

Definition at line 125 of file AdolCHuHu.cpp.

---

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

• src/adolc/impl/AdolCHuHu.cpp