ConvectiveMassElement::OpVelocityRhs Struct Reference

#include <users_modules/basic_finite_elements/src/ConvectiveMassElement.hpp>

Inheritance diagram for ConvectiveMassElement::OpVelocityRhs:

Collaboration diagram for ConvectiveMassElement::OpVelocityRhs:

Public Member Functions
	OpVelocityRhs (const std::string field_name, BlockData &data, CommonData &common_data)
MoFEMErrorCode	doWork (int row_side, EntityType row_type, EntitiesFieldData::EntData &row_data)

Public Attributes
BlockData &	dAta
CommonData &	commonData
VectorDouble	nf

Detailed Description

Definition at line 301 of file ConvectiveMassElement.hpp.

Constructor & Destructor Documentation

OpVelocityRhs()

ConvectiveMassElement::OpVelocityRhs::OpVelocityRhs	(	const std::string	field_name,
		BlockData &	data,
		CommonData &	common_data
	)

Definition at line 1009 of file ConvectiveMassElement.cpp.

    : VolumeElementForcesAndSourcesCore::UserDataOperator(
          field_name, ForcesAndSourcesCore::UserDataOperator::OPROW),
      dAta(data), commonData(common_data) {}

Member Function Documentation

doWork()

MoFEMErrorCode ConvectiveMassElement::OpVelocityRhs::doWork	(	int	row_side,
		EntityType	row_type,
		EntitiesFieldData::EntData &	row_data
	)

Definition at line 1015 of file ConvectiveMassElement.cpp.

                                        {
  MoFEMFunctionBeginHot;
 
  if (dAta.tEts.find(getNumeredEntFiniteElementPtr()->getEnt()) ==
      dAta.tEts.end()) {
    MoFEMFunctionReturnHot(0);
  }
  int nb_dofs = row_data.getIndices().size();
  if (nb_dofs == 0)
    MoFEMFunctionReturnHot(0);
 
  auto base = row_data.getFTensor0N();
  int nb_base_functions = row_data.getN().size2();
  FTensor::Index<'i', 3> i;
 
  {
 
    nf.resize(nb_dofs);
    nf.clear();
 
    for (unsigned int gg = 0; gg < row_data.getN().size1(); gg++) {
      FTensor::Tensor1<double *, 3> t_nf(&nf[0], &nf[1], &nf[2], 3);
      FTensor::Tensor1<double *, 3> res(&commonData.valVel[gg][0],
                                        &commonData.valVel[gg][1],
                                        &commonData.valVel[gg][2]);
      int dd = 0;
      for (; dd < nb_dofs / 3; dd++) {
        t_nf(i) += base * res(i);
        ++base;
        ++t_nf;
      }
      for (; dd != nb_base_functions; dd++) {
        ++base;
      }
    }
 
    if (row_data.getIndices().size() > 3 * row_data.getN().size2()) {
      SETERRQ(PETSC_COMM_SELF, 1, "data inconsistency");
    }
    CHKERR VecSetValues(getFEMethod()->ts_F, row_data.getIndices().size(),
                        &row_data.getIndices()[0], &nf[0], ADD_VALUES);
  }
  MoFEMFunctionReturnHot(0);
}

Member Data Documentation

commonData

CommonData& ConvectiveMassElement::OpVelocityRhs::commonData

Definition at line 306 of file ConvectiveMassElement.hpp.

dAta

BlockData& ConvectiveMassElement::OpVelocityRhs::dAta

Definition at line 305 of file ConvectiveMassElement.hpp.

nf

VectorDouble ConvectiveMassElement::OpVelocityRhs::nf

Definition at line 311 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