EdgeFlippingOps::OpLhsSymTensorLeastSquaresProjImpl< DIM, GAUSS, AssemblyDomainEleOp > Struct Template Reference

#include "users_modules/multifield-thermoplasticity-private/src/EdgeFlippingOps.hpp"

Inheritance diagram for EdgeFlippingOps::OpLhsSymTensorLeastSquaresProjImpl< DIM, GAUSS, AssemblyDomainEleOp >:

Collaboration diagram for EdgeFlippingOps::OpLhsSymTensorLeastSquaresProjImpl< DIM, GAUSS, AssemblyDomainEleOp >:

Public Member Functions
	OpLhsSymTensorLeastSquaresProjImpl (const std::string row_field_name, const std::string col_field_name)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Detailed Description

template<int DIM, typename AssemblyDomainEleOp>
struct EdgeFlippingOps::OpLhsSymTensorLeastSquaresProjImpl< DIM, GAUSS, AssemblyDomainEleOp >

Definition at line 359 of file EdgeFlippingOps.hpp.

Constructor & Destructor Documentation

OpLhsSymTensorLeastSquaresProjImpl()

template<int DIM, typename AssemblyDomainEleOp >

EdgeFlippingOps::OpLhsSymTensorLeastSquaresProjImpl< DIM, GAUSS, AssemblyDomainEleOp >::OpLhsSymTensorLeastSquaresProjImpl	(	const std::string	row_field_name,
		const std::string	col_field_name
	)

Definition at line 370 of file EdgeFlippingOps.hpp.

    : AssemblyDomainEleOp(row_field_name, col_field_name,
                          AssemblyDomainEleOp::OPROWCOL) {
  AssemblyDomainEleOp::sYmm = false;
}

Member Function Documentation

iNtegrate()

template<int DIM, typename AssemblyDomainEleOp >

MoFEMErrorCode EdgeFlippingOps::OpLhsSymTensorLeastSquaresProjImpl< DIM, GAUSS, AssemblyDomainEleOp >::iNtegrate	(	EntitiesFieldData::EntData &	row_data,
		EntitiesFieldData::EntData &	col_data
	)

Definition at line 380 of file EdgeFlippingOps.hpp.

                                        {
  MoFEMFunctionBegin;
 
  FTENSOR_INDEX(DIM, i);
  FTENSOR_INDEX(DIM, j);
  FTENSOR_INDEX(DIM, k);
  FTENSOR_INDEX(DIM, l);
  constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
  FTENSOR_INDEX(size_symm, L);
  FTENSOR_INDEX(size_symm, O);
 
  auto &locMat = AssemblyDomainEleOp::locMat;
 
  const auto nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
  const auto nb_row_base_functions = row_data.getN().size2();
 
  auto t_L = PlasticOps::symm_L_tensor(FTensor::Number<DIM>());
 
  auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
  auto t_row_base = row_data.getFTensor0N();
  for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
    double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
    ++t_w;
 
    FTensor::Tensor2<double, size_symm, size_symm> t_res_mat;
    t_res_mat(O, L) = alpha * (t_L(i, j, O) * t_L(i, j, L));
 
    size_t rr = 0;
    for (; rr != AssemblyDomainEleOp::nbRows / size_symm; ++rr) {
      auto t_mat = PlasticOps::get_mat_tensor_sym_dtensor_sym(
          rr, locMat, FTensor::Number<SPACE_DIM>());
      auto t_col_base = col_data.getFTensor0N(gg, 0);
      for (size_t cc = 0; cc != AssemblyDomainEleOp::nbCols / size_symm; ++cc) {
        t_mat(O, L) += ((t_row_base * t_col_base) * t_res_mat(O, L));
        ++t_mat;
        ++t_col_base;
      }
 
      ++t_row_base;
    }
 
    for (; rr < nb_row_base_functions; ++rr)
      ++t_row_base;
  }
 
  MoFEMFunctionReturn(0);
}

---

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

• users_modules/multifield-thermoplasticity-private/src/EdgeFlippingOps.hpp