SimpleContactProblem::OpGetGaussPtsState Struct Reference

#include <users_modules/mortar_contact/src/SimpleContact.hpp>

Inheritance diagram for SimpleContactProblem::OpGetGaussPtsState:

Collaboration diagram for SimpleContactProblem::OpGetGaussPtsState:

Public Member Functions
	OpGetGaussPtsState (const string lagrange_field_name, boost::shared_ptr< CommonDataSimpleContact > common_data_contact, const double cn, const bool alm_flag=false)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Private Attributes
boost::shared_ptr< CommonDataSimpleContact >	commonDataSimpleContact
const double	cN
const bool	almFlag
VectorDouble	vecR

Detailed Description

Definition at line 3523 of file SimpleContact.hpp.

Constructor & Destructor Documentation

OpGetGaussPtsState()

SimpleContactProblem::OpGetGaussPtsState::OpGetGaussPtsState ( const string  lagrange_field_name, boost::shared_ptr< CommonDataSimpleContact >  common_data_contact, const double  cn, const bool  alm_flag = false  )

inline

Definition at line 3525 of file SimpleContact.hpp.

        : ContactOp(lagrange_field_name, UserDataOperator::OPCOL,
                    ContactOp::FACESLAVE),
          commonDataSimpleContact(common_data_contact), cN(cn),
          almFlag(alm_flag) {}

Member Function Documentation

doWork()

MoFEMErrorCode SimpleContactProblem::OpGetGaussPtsState::doWork	(	int	side,
		EntityType	type,
		EntData &	data
	)

Definition at line 3869 of file SimpleContact.cpp.

                                                                               {
  MoFEMFunctionBegin;
 
  if (data.getIndices().size() == 0)
    MoFEMFunctionReturnHot(0);
 
  const int nb_gauss_pts = data.getN().size1();
 
  vecR.resize(CommonDataSimpleContact::LAST_ELEMENT, false);
  vecR.clear();
 
  commonDataSimpleContact->gaussPtsStatePtr->resize(nb_gauss_pts, false);
  commonDataSimpleContact->gaussPtsStatePtr->clear();
 
  auto t_state_gp =
      getFTensor0FromVec(*commonDataSimpleContact->gaussPtsStatePtr);
 
  auto t_lagrange_slave =
      getFTensor0FromVec(*commonDataSimpleContact->lagMultAtGaussPtsPtr);
  auto t_gap_gp = getFTensor0FromVec(*commonDataSimpleContact->gapPtr);
 
  for (int gg = 0; gg != nb_gauss_pts; gg++) {
    vecR[CommonDataSimpleContact::TOTAL] += 1;
 
    if (!almFlag &&
        SimpleContactProblem::State(cN, t_gap_gp, t_lagrange_slave)) {
      vecR[CommonDataSimpleContact::ACTIVE] += 1;
      t_state_gp = 1;
    }
 
    if (almFlag &&
        SimpleContactProblem::StateALM(cN, t_gap_gp, t_lagrange_slave)) {
      vecR[CommonDataSimpleContact::ACTIVE] += 1;
      t_state_gp = 1;
    }
 
    ++t_lagrange_slave;
    ++t_gap_gp;
    ++t_state_gp;
  } // for gauss points
 
  constexpr std::array<int, 2> indices = {CommonDataSimpleContact::ACTIVE,
                                          CommonDataSimpleContact::TOTAL};
  CHKERR VecSetValues(commonDataSimpleContact->gaussPtsStateVec, 2,
                      indices.data(), &vecR[0], ADD_VALUES);
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

almFlag

const bool SimpleContactProblem::OpGetGaussPtsState::almFlag

private

Definition at line 3539 of file SimpleContact.hpp.

cN

const double SimpleContactProblem::OpGetGaussPtsState::cN

private

Definition at line 3538 of file SimpleContact.hpp.

commonDataSimpleContact

boost::shared_ptr<CommonDataSimpleContact> SimpleContactProblem::OpGetGaussPtsState::commonDataSimpleContact

private

Definition at line 3537 of file SimpleContact.hpp.

vecR

VectorDouble SimpleContactProblem::OpGetGaussPtsState::vecR

private

Definition at line 3540 of file SimpleContact.hpp.

---

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

• users_modules/mortar_contact/src/SimpleContact.hpp
• users_modules/mortar_contact/src/impl/SimpleContact.cpp