ConvectiveMassElement::OpMassLhs_dM_dX Struct Reference

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

Inheritance diagram for ConvectiveMassElement::OpMassLhs_dM_dX:

Collaboration diagram for ConvectiveMassElement::OpMassLhs_dM_dX:

Public Member Functions
	OpMassLhs_dM_dX (const std::string field_name, const std::string col_field, BlockData &data, CommonData &common_data)
MoFEMErrorCode	getJac (EntitiesFieldData::EntData &col_data, int gg)
Public Member Functions inherited from ConvectiveMassElement::OpMassLhs_dM_dv
	OpMassLhs_dM_dv (const std::string vel_field, const std::string field_name, BlockData &data, CommonData &common_data, Range *forcesonlyonentities_ptr=NULL)
MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Additional Inherited Members
Public Attributes inherited from ConvectiveMassElement::OpMassLhs_dM_dv
BlockData &	dAta
CommonData &	commonData
Range	forcesOnlyOnEntities
MatrixDouble	k
MatrixDouble	jac

Detailed Description

Definition at line 253 of file ConvectiveMassElement.hpp.

Constructor & Destructor Documentation

OpMassLhs_dM_dX()

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

Definition at line 664 of file ConvectiveMassElement.cpp.

    : OpMassLhs_dM_dv(field_name, col_field, data, common_data) {}

Member Function Documentation

getJac()

MoFEMErrorCode ConvectiveMassElement::OpMassLhs_dM_dX::getJac ( EntitiesFieldData::EntData &  col_data, int  gg  )

virtual

Reimplemented from ConvectiveMassElement::OpMassLhs_dM_dv.

Definition at line 669 of file ConvectiveMassElement.cpp.

                                                {
  MoFEMFunctionBeginHot;
  int nb_col = col_data.getIndices().size();
  jac.clear();
  double *base_ptr = const_cast<double *>(&col_data.getN(gg)[0]);
  FTensor::Tensor0<double *> base(base_ptr, 1);
  double *diff_ptr =
      const_cast<double *>(&(col_data.getDiffN(gg, nb_col / 3)(0, 0)));
  FTensor::Tensor1<double *, 3> diff(diff_ptr, &diff_ptr[1], &diff_ptr[2], 3);
  FTensor::Tensor2<double *, 3, 3> t_jac(&jac(0, 0), &jac(0, 1), &jac(0, 2),
                                         &jac(1, 0), &jac(1, 1), &jac(1, 2),
                                         &jac(2, 0), &jac(2, 1), &jac(2, 2), 3);
  const int u = 3 + 9 + 9;
  FTensor::Tensor2<double *, 3, 3> t_mass1(
      &commonData.jacMass[gg](0, u + 0), &commonData.jacMass[gg](0, u + 1),
      &commonData.jacMass[gg](0, u + 2), &commonData.jacMass[gg](1, u + 0),
      &commonData.jacMass[gg](1, u + 1), &commonData.jacMass[gg](1, u + 2),
      &commonData.jacMass[gg](2, u + 0), &commonData.jacMass[gg](2, u + 1),
      &commonData.jacMass[gg](2, u + 2));
  const int s = 3 + 9 + 9 + 3;
  FTensor::Tensor3<double *, 3, 3, 3> t_mass3(
      // T* d000, T* d001, T* d002,
      // T* d010, T* d011, T* d012,
      // T* d020, T* d021, T* d022,
      // T* d100, T* d101, T* d102,
      // T* d110, T* d111, T* d112,
      // T* d120, T* d121, T* d122,
      // T* d200, T* d201, T* d202,
      // T* d210, T* d211, T* d212,
      // T* d220, T* d221, T* d222,
      &commonData.jacMass[gg](0, s + 0), &commonData.jacMass[gg](0, s + 1),
      &commonData.jacMass[gg](0, s + 2), &commonData.jacMass[gg](0, s + 3),
      &commonData.jacMass[gg](0, s + 4), &commonData.jacMass[gg](0, s + 5),
      &commonData.jacMass[gg](0, s + 6), &commonData.jacMass[gg](0, s + 7),
      &commonData.jacMass[gg](0, s + 8), &commonData.jacMass[gg](1, s + 0),
      &commonData.jacMass[gg](1, s + 1), &commonData.jacMass[gg](1, s + 2),
      &commonData.jacMass[gg](1, s + 3), &commonData.jacMass[gg](1, s + 4),
      &commonData.jacMass[gg](1, s + 5), &commonData.jacMass[gg](1, s + 6),
      &commonData.jacMass[gg](1, s + 7), &commonData.jacMass[gg](1, s + 8),
      &commonData.jacMass[gg](2, s + 0), &commonData.jacMass[gg](2, s + 1),
      &commonData.jacMass[gg](2, s + 2), &commonData.jacMass[gg](2, s + 3),
      &commonData.jacMass[gg](2, s + 4), &commonData.jacMass[gg](2, s + 5),
      &commonData.jacMass[gg](2, s + 6), &commonData.jacMass[gg](2, s + 7),
      &commonData.jacMass[gg](2, s + 8));
  FTensor::Index<'i', 3> i;
  FTensor::Index<'j', 3> j;
  FTensor::Index<'k', 3> k;
  for (int dd = 0; dd < nb_col / 3; dd++) {
    t_jac(i, j) += t_mass1(i, j) * base * getFEMethod()->ts_a;
    t_jac(i, j) += t_mass3(i, j, k) * diff(k);
    ++base_ptr;
    ++diff_ptr;
    ++t_jac;
  }
  MoFEMFunctionReturnHot(0);
}

---

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