ConvectiveMassElement::OpMassRhs Struct Reference

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

Inheritance diagram for ConvectiveMassElement::OpMassRhs:

Collaboration diagram for ConvectiveMassElement::OpMassRhs:

Public Member Functions
	OpMassRhs (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 207 of file ConvectiveMassElement.hpp.

Constructor & Destructor Documentation

OpMassRhs()

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

Definition at line 383 of file ConvectiveMassElement.cpp.

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

Member Function Documentation

doWork()

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

Definition at line 391 of file ConvectiveMassElement.cpp.

                                                                             {
  MoFEMFunctionBeginHot;
 
  if (dAta.tEts.find(getNumeredEntFiniteElementPtr()->getEnt()) ==
      dAta.tEts.end()) {
    MoFEMFunctionReturnHot(0);
  }
  if (row_data.getIndices().size() == 0)
    MoFEMFunctionReturnHot(0);
  int nb_dofs = row_data.getIndices().size();
 
  auto base = row_data.getFTensor0N();
  int nb_base_functions = row_data.getN().size2();
 
  {
 
    nf.resize(nb_dofs);
    nf.clear();
 
    FTensor::Index<'i', 3> i;
 
    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.valMass[gg][0],
                                        &commonData.valMass[gg][1],
                                        &commonData.valMass[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 ((unsigned int)nb_dofs > 3 * row_data.getN().size2()) {
      SETERRQ(PETSC_COMM_SELF, 1, "data inconsistency");
    }
    CHKERR VecSetValues(getFEMethod()->ts_F, nb_dofs, &row_data.getIndices()[0],
                        &nf[0], ADD_VALUES);
  }
 
  MoFEMFunctionReturnHot(0);
}

Member Data Documentation

commonData

CommonData& ConvectiveMassElement::OpMassRhs::commonData

Definition at line 211 of file ConvectiveMassElement.hpp.

dAta

BlockData& ConvectiveMassElement::OpMassRhs::dAta

Definition at line 210 of file ConvectiveMassElement.hpp.

nf

VectorDouble ConvectiveMassElement::OpMassRhs::nf

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