FractureMechanics::OpGetCrackFrontDataGradientAtGaussPts Struct Reference

#include <users_modules/fracture_mechanics/src/CrackFrontElement.hpp>

Inheritance diagram for FractureMechanics::OpGetCrackFrontDataGradientAtGaussPts:

Collaboration diagram for FractureMechanics::OpGetCrackFrontDataGradientAtGaussPts:

Public Member Functions
	OpGetCrackFrontDataGradientAtGaussPts (const std::string field_name, boost::shared_ptr< MatrixDouble > data_at_pts, bool &singular_element, MatrixDouble &inv_s_jac)
MoFEMErrorCode	doWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)

Public Attributes
bool &	singularElement
MatrixDouble &	invSJac

Detailed Description

Definition at line 478 of file CrackFrontElement.hpp.

Constructor & Destructor Documentation

OpGetCrackFrontDataGradientAtGaussPts()

FractureMechanics::OpGetCrackFrontDataGradientAtGaussPts::OpGetCrackFrontDataGradientAtGaussPts ( const std::string  field_name, boost::shared_ptr< MatrixDouble >  data_at_pts, bool &  singular_element, MatrixDouble &  inv_s_jac  )

inline

Definition at line 484 of file CrackFrontElement.hpp.

      : OpCalculateVectorFieldGradient<3, 3>(field_name, data_at_pts),
        singularElement(singular_element), invSJac(inv_s_jac) {}

Member Function Documentation

doWork()

MoFEMErrorCode FractureMechanics::OpGetCrackFrontDataGradientAtGaussPts::doWork	(	int	side,
		EntityType	type,
		DataForcesAndSourcesCore::EntData &	data
	)

Definition at line 263 of file CrackFrontElement.cpp.

                                                                      {
  MoFEMFunctionBegin;
 
  if (!singularElement) {
    return OpCalculateVectorFieldGradient<3, 3>::doWork(side, type, data);
  }
 
  const int nb_dofs = data.getFieldData().size();
  if (!nb_dofs && type == this->zeroType) {
    this->dataPtr->resize(9, 0, false);
    MoFEMFunctionReturnHot(0);
  }
  if (!nb_dofs) {
    MoFEMFunctionReturnHot(0);
  }
  const int nb_gauss_pts = data.getN().size1();
  const int nb_base_functions = data.getN().size2();
  ublas::matrix<double, ublas::row_major, DoubleAllocator> &mat =
      *this->dataPtr;
  if (type == this->zeroType) {
    mat.resize(9, nb_gauss_pts, false);
    mat.clear();
  }
  auto diff_base_function = data.getFTensor1DiffN<3>();
  auto gradients_at_gauss_pts = getFTensor2FromMat<3, 3>(mat);
  FTensor::Index<'I', 3> I;
  FTensor::Index<'J', 3> J;
  FTensor::Index<'K', 3> K;
  int size = nb_dofs / 3;
  if (nb_dofs % 3) {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "Data inconsistency");
  }
  double *inv_jac_ptr = &*invSJac.data().begin();
  FTensor::Tensor2<double *, 3, 3> t_inv_singular_jacobian(
      inv_jac_ptr, &inv_jac_ptr[1], &inv_jac_ptr[2], &inv_jac_ptr[3],
      &inv_jac_ptr[4], &inv_jac_ptr[5], &inv_jac_ptr[6], &inv_jac_ptr[7],
      &inv_jac_ptr[8], 9);
 
  for (int gg = 0; gg < nb_gauss_pts; ++gg) {
    auto field_data = data.getFTensor1FieldData<3>();
    int bb = 0;
    for (; bb < size; ++bb) {
      gradients_at_gauss_pts(I, J) +=
          field_data(I) *
          (t_inv_singular_jacobian(K, J) * diff_base_function(K));
      ++field_data;
      ++diff_base_function;
    }
    // Number of dofs can be smaller than number of Tensor_Dim0 x base functions
    for (; bb != nb_base_functions; ++bb)
      ++diff_base_function;
    ++gradients_at_gauss_pts;
    ++t_inv_singular_jacobian;
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

invSJac

MatrixDouble& FractureMechanics::OpGetCrackFrontDataGradientAtGaussPts::invSJac

Definition at line 482 of file CrackFrontElement.hpp.

singularElement

bool& FractureMechanics::OpGetCrackFrontDataGradientAtGaussPts::singularElement

Definition at line 481 of file CrackFrontElement.hpp.

---

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

• users_modules/fracture_mechanics/src/CrackFrontElement.hpp
• users_modules/fracture_mechanics/src/impl/CrackFrontElement.cpp