OpRotatingFrameTools::OpVolumeSideCalculateEP Struct Reference

#include <users_modules/multifield_plasticity/src/RotatingFrameOperators.hpp>

Inheritance diagram for OpRotatingFrameTools::OpVolumeSideCalculateEP:

Collaboration diagram for OpRotatingFrameTools::OpVolumeSideCalculateEP:

Public Member Functions
	OpVolumeSideCalculateEP (const std::string field_name, boost::shared_ptr< CommonData > common_data_ptr)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Private Attributes
boost::shared_ptr< CommonData >	commonDataPtr

Detailed Description

Definition at line 101 of file RotatingFrameOperators.hpp.

Constructor & Destructor Documentation

OpVolumeSideCalculateEP()

OpRotatingFrameTools::OpVolumeSideCalculateEP::OpVolumeSideCalculateEP	(	const std::string	field_name,
		boost::shared_ptr< CommonData >	common_data_ptr
	)

Definition at line 470 of file RotatingFrameOperators.cpp.

    : DomainSideEleOp(field_name, field_name, DomainEleOp::OPROW),
      commonDataPtr(common_data_ptr) {
  std::fill(&doEntities[MBVERTEX], &doEntities[MBMAXTYPE], false);
  doEntities[MBTET] = true;
  doEntities[MBHEX] = true;
}

Member Function Documentation

doWork()

MoFEMErrorCode OpRotatingFrameTools::OpVolumeSideCalculateEP::doWork	(	int	side,
		EntityType	type,
		EntData &	data
	)

Definition at line 479 of file RotatingFrameOperators.cpp.

                                                              {
  MoFEMFunctionBegin;
  const int nb_gauss_pts = getGaussPts().size2();
  const int nb_base_functions = data.getN().size2();
  // commonDataPtr->mStressPtr->resize(6, nb_gauss_pts);
  // auto &t_D = commonDataPtr->tD;
  // auto t_strain =
  // getFTensor2SymmetricFromMat<3>(*(commonDataPtr->mStrainPtr)); auto t_stress
  // = getFTensor2SymmetricFromMat<3>(*(commonDataPtr->mStressPtr));
  // FIXME: this is wrong, resize only on zero type
  const size_t nb_dofs = data.getIndices().size();
  // data.getFieldData().size() == 0
  // if (side == getFaceSideNumber()) {
  if (nb_dofs) {
 
    // for testing //FIXME:
    MatrixDouble diff = getCoordsAtGaussPts() - getFaceCoordsAtGaussPts();
    const double eps = 1e-12;
    if (norm_inf(diff) > eps)
      SETERRQ(PETSC_COMM_WORLD, MOFEM_ATOM_TEST_INVALID,
              "coordinates at integration pts are different");
 
    const size_t nb_gauss_pts = data.getN().size1();
 
    auto base_function = data.getFTensor0N();
    // MatrixDouble *ep_elem_data = nullptr;
    int test_nb_in_loop = getFEMethod()->nInTheLoop;
    int test_face_sense = getFaceSense();
    MatrixDouble &mat_ep = commonDataPtr->plasticStrainSideMap[getFaceSense()];
    mat_ep.resize(6, nb_gauss_pts, false);
    mat_ep.clear();
 
    // if (getFEMethod()->nInTheLoop == 0)
    //   ep_elem_data = &(commonDataPtr->plasticStrainLeft);
    // else
    //   ep_elem_data = &(commonDataPtr->plasticStrainRight);
 
    // mat_ep.resize(6, nb_gauss_pts, false);
    // mat_ep.clear();
    auto values_at_gauss_pts = getFTensor2SymmetricFromMat<3>(mat_ep);
    for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
      auto field_data = data.getFTensor2SymmetricFieldData<3>();
      size_t bb = 0;
      for (; bb != nb_dofs / 6; ++bb) {
        values_at_gauss_pts(i, j) += field_data(i, j) * base_function;
        ++field_data;
        ++base_function;
      }
      for (; bb != nb_base_functions; ++bb)
        ++base_function;
      ++values_at_gauss_pts;
    }
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataPtr

boost::shared_ptr<CommonData> OpRotatingFrameTools::OpVolumeSideCalculateEP::commonDataPtr

private

Definition at line 107 of file RotatingFrameOperators.hpp.

---

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

• users_modules/multifield_plasticity/src/RotatingFrameOperators.hpp
• users_modules/multifield_plasticity/src/impl/RotatingFrameOperators.cpp