ConvectiveMassElement::OpEshelbyDynamicMaterialMomentumRhs Struct Reference

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

Inheritance diagram for ConvectiveMassElement::OpEshelbyDynamicMaterialMomentumRhs:

Collaboration diagram for ConvectiveMassElement::OpEshelbyDynamicMaterialMomentumRhs:

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

Public Attributes
BlockData &	dAta
CommonData &	commonData
Range	forcesOnlyOnEntities
VectorDouble	nf

Detailed Description

Definition at line 370 of file ConvectiveMassElement.hpp.

Constructor & Destructor Documentation

OpEshelbyDynamicMaterialMomentumRhs()

ConvectiveMassElement::OpEshelbyDynamicMaterialMomentumRhs::OpEshelbyDynamicMaterialMomentumRhs	(	const std::string	field_name,
		BlockData &	data,
		CommonData &	common_data,
		Range *	forcesonlyonentities_ptr
	)

Definition at line 1456 of file ConvectiveMassElement.cpp.

    : VolumeElementForcesAndSourcesCore::UserDataOperator(
          field_name, ForcesAndSourcesCore::UserDataOperator::OPROW),
      dAta(data), commonData(common_data) {
  if (forcesonlyonentities_ptr != NULL) {
    forcesOnlyOnEntities = *forcesonlyonentities_ptr;
  }
}

Member Function Documentation

doWork()

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

Definition at line 1470 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);
 
  try {
 
    nf.resize(nb_dofs);
    nf.clear();
 
    auto base = row_data.getFTensor0N();
    int nb_base_functions = row_data.getN().size2();
    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.valT[gg][0],
                                        &commonData.valT[gg][1],
                                        &commonData.valT[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");
    }
    if (!forcesOnlyOnEntities.empty()) {
      VectorInt indices = row_data.getIndices();
      VectorDofs &dofs = row_data.getFieldDofs();
      VectorDofs::iterator dit = dofs.begin();
      for (int ii = 0; dit != dofs.end(); dit++, ii++) {
        if (forcesOnlyOnEntities.find((*dit)->getEnt()) ==
            forcesOnlyOnEntities.end()) {
          // std::cerr << **dit << std::endl;
          indices[ii] = -1;
        }
      }
      // std::cerr << indices << std::endl;
      CHKERR VecSetValues(getFEMethod()->ts_F, indices.size(), &indices[0],
                          &nf[0], ADD_VALUES);
    } else {
      CHKERR VecSetValues(getFEMethod()->ts_F, row_data.getIndices().size(),
                          &row_data.getIndices()[0], &nf[0], ADD_VALUES);
    }
 
  } catch (const std::exception &ex) {
    std::ostringstream ss;
    ss << "throw in method: " << ex.what() << std::endl;
    SETERRQ(PETSC_COMM_SELF, 1, ss.str().c_str());
  }
 
  MoFEMFunctionReturnHot(0);
}

Member Data Documentation

commonData

CommonData& ConvectiveMassElement::OpEshelbyDynamicMaterialMomentumRhs::commonData

Definition at line 375 of file ConvectiveMassElement.hpp.

dAta

BlockData& ConvectiveMassElement::OpEshelbyDynamicMaterialMomentumRhs::dAta

Definition at line 374 of file ConvectiveMassElement.hpp.

forcesOnlyOnEntities

Range ConvectiveMassElement::OpEshelbyDynamicMaterialMomentumRhs::forcesOnlyOnEntities

Definition at line 376 of file ConvectiveMassElement.hpp.

nf

VectorDouble ConvectiveMassElement::OpEshelbyDynamicMaterialMomentumRhs::nf

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