ConvectiveMassElement::OpGetDataAtGaussPts Struct Reference

#include "users_modules/basic_finite_elements/src/ConvectiveMassElement.hpp"

Inheritance diagram for ConvectiveMassElement::OpGetDataAtGaussPts:

Collaboration diagram for ConvectiveMassElement::OpGetDataAtGaussPts:

Public Member Functions
	OpGetDataAtGaussPts (const std::string field_name, std::vector< VectorDouble > &values_at_gauss_pts, std::vector< MatrixDouble > &gardient_at_gauss_pts)
MoFEMErrorCode	doWork (int side, EntityType type, EntitiesFieldData::EntData &data)
	operator calculating deformation gradient

Public Attributes
std::vector< VectorDouble > &	valuesAtGaussPts
std::vector< MatrixDouble > &	gradientAtGaussPts
const EntityType	zeroAtType

Detailed Description

Definition at line 152 of file ConvectiveMassElement.hpp.

Constructor & Destructor Documentation

OpGetDataAtGaussPts()

ConvectiveMassElement::OpGetDataAtGaussPts::OpGetDataAtGaussPts	(	const std::string	field_name,
		std::vector< VectorDouble > &	values_at_gauss_pts,
		std::vector< MatrixDouble > &	gardient_at_gauss_pts
	)

Definition at line 80 of file ConvectiveMassElement.cpp.

    : VolumeElementForcesAndSourcesCore::UserDataOperator(
          field_name, ForcesAndSourcesCore::UserDataOperator::OPROW),
      valuesAtGaussPts(values_at_gauss_pts),
      gradientAtGaussPts(gardient_at_gauss_pts), zeroAtType(MBVERTEX) {}

Member Function Documentation

doWork()

MoFEMErrorCode ConvectiveMassElement::OpGetDataAtGaussPts::doWork	(	int	side,
		EntityType	type,
		EntitiesFieldData::EntData &	data
	)

operator calculating deformation gradient

Definition at line 89 of file ConvectiveMassElement.cpp.

                                                               {
  MoFEMFunctionBeginHot;
 
  int nb_dofs = data.getFieldData().size();
  if (nb_dofs == 0) {
    MoFEMFunctionReturnHot(0);
  }
  int nb_gauss_pts = data.getN().size1();
  int nb_base_functions = data.getN().size2();
 
  // initialize
  // VectorDouble& values = data.getFieldData();
  valuesAtGaussPts.resize(nb_gauss_pts);
  gradientAtGaussPts.resize(nb_gauss_pts);
  for (int gg = 0; gg < nb_gauss_pts; gg++) {
    valuesAtGaussPts[gg].resize(3);
    gradientAtGaussPts[gg].resize(3, 3);
  }
 
  if (type == zeroAtType) {
    for (int gg = 0; gg < nb_gauss_pts; gg++) {
      valuesAtGaussPts[gg].clear();
      gradientAtGaussPts[gg].clear();
    }
  }
 
  auto base_function = data.getFTensor0N();
  auto diff_base_functions = data.getFTensor1DiffN<3>();
  FTensor::Index<'i', 3> i;
  FTensor::Index<'j', 3> j;
 
  for (int gg = 0; gg != nb_gauss_pts; gg++) {
    auto field_data = data.getFTensor1FieldData<3>();
    FTensor::Tensor1<double *, 3> values(&valuesAtGaussPts[gg][0],
                                         &valuesAtGaussPts[gg][1],
                                         &valuesAtGaussPts[gg][2]);
    FTensor::Tensor2<double *, 3, 3> gradient(
        &gradientAtGaussPts[gg](0, 0), &gradientAtGaussPts[gg](0, 1),
        &gradientAtGaussPts[gg](0, 2), &gradientAtGaussPts[gg](1, 0),
        &gradientAtGaussPts[gg](1, 1), &gradientAtGaussPts[gg](1, 2),
        &gradientAtGaussPts[gg](2, 0), &gradientAtGaussPts[gg](2, 1),
        &gradientAtGaussPts[gg](2, 2));
    int bb = 0;
    for (; bb != nb_dofs / 3; bb++) {
      values(i) += base_function * field_data(i);
      gradient(i, j) += field_data(i) * diff_base_functions(j);
      ++diff_base_functions;
      ++base_function;
      ++field_data;
    }
    for (; bb != nb_base_functions; bb++) {
      ++diff_base_functions;
      ++base_function;
    }
  }
  MoFEMFunctionReturnHot(0);
}

Member Data Documentation

gradientAtGaussPts

std::vector<MatrixDouble>& ConvectiveMassElement::OpGetDataAtGaussPts::gradientAtGaussPts

Definition at line 156 of file ConvectiveMassElement.hpp.

valuesAtGaussPts

std::vector<VectorDouble>& ConvectiveMassElement::OpGetDataAtGaussPts::valuesAtGaussPts

Definition at line 155 of file ConvectiveMassElement.hpp.

zeroAtType

const EntityType ConvectiveMassElement::OpGetDataAtGaussPts::zeroAtType

Definition at line 157 of file ConvectiveMassElement.hpp.

---

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

• users_modules/basic_finite_elements/src/ConvectiveMassElement.hpp
• users_modules/basic_finite_elements/src/impl/ConvectiveMassElement.cpp