MatOps::OpLhsHuHu< FIELD_DIM, SPACE_DIM, A, GAUSS, OpBase > Struct Template Reference

#include "src/materials/MatHuHu.hpp"

Inheritance diagram for MatOps::OpLhsHuHu< FIELD_DIM, SPACE_DIM, A, GAUSS, OpBase >:

Collaboration diagram for MatOps::OpLhsHuHu< FIELD_DIM, SPACE_DIM, A, GAUSS, OpBase >:

Public Types
using	OP = typename FormsIntegrators< OpBase >::template Assembly< A >::OpBase

Public Member Functions
	OpLhsHuHu (const std::string field_name, boost::shared_ptr< MatrixDouble > mat_K, boost::shared_ptr< Range > ents_ptr=nullptr)

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Protected Attributes
boost::shared_ptr< MatrixDouble >	matK

Detailed Description

template<int FIELD_DIM, int SPACE_DIM, AssemblyType A, typename OpBase>
struct MatOps::OpLhsHuHu< FIELD_DIM, SPACE_DIM, A, GAUSS, OpBase >

Definition at line 190 of file MatHuHu.hpp.

Member Typedef Documentation

OP

template<int FIELD_DIM, int SPACE_DIM, AssemblyType A, typename OpBase >

using MatOps::OpLhsHuHu< FIELD_DIM, SPACE_DIM, A, GAUSS, OpBase >::OP = typename FormsIntegrators<OpBase>::template Assembly<A>::OpBase

Definition at line 193 of file MatHuHu.hpp.

Constructor & Destructor Documentation

OpLhsHuHu()

template<int FIELD_DIM, int SPACE_DIM, AssemblyType A, typename OpBase >

MatOps::OpLhsHuHu< FIELD_DIM, SPACE_DIM, A, GAUSS, OpBase >::OpLhsHuHu ( const std::string  field_name, boost::shared_ptr< MatrixDouble >  mat_K, boost::shared_ptr< Range >  ents_ptr = nullptr  )

inline

Definition at line 194 of file MatHuHu.hpp.

      : OP(field_name, field_name, OP::OPROWCOL), matK(mat_K) {}

Member Function Documentation

iNtegrate()

template<int FIELD_DIM, int SPACE_DIM, AssemblyType A, typename OpBase >

MoFEMErrorCode MatOps::OpLhsHuHu< FIELD_DIM, SPACE_DIM, A, GAUSS, OpBase >::iNtegrate ( EntitiesFieldData::EntData &  row_data, EntitiesFieldData::EntData &  col_data  )

protected

Definition at line 279 of file MatHuHu.hpp.

                                        {
  MoFEMFunctionBegin;
 
  FTENSOR_INDEX(FIELD_DIM, i);
  FTENSOR_INDEX(FIELD_DIM, j);
  FTENSOR_INDEX(SPACE_DIM, J);
  FTENSOR_INDEX(SPACE_DIM, K);
 
  auto get_K_at_pts =
      MatrixSizeHelper<GetFTensor1FromMatType<1, -1, DL>, DL>::get(
          *matK, OP::nbIntegrationPts);
 
  constexpr auto t_kd = FTensor::Kronecker_Delta<int>();
 
  // get element volume
  const double vol = OP::getMeasure();
  // get integration weights
  auto t_w = OP::getFTensor0IntegrationWeight();
  // get base function gradient on rows
  auto t_row_grad = row_data.getFTensor2DiffN2<SPACE_DIM>();
  // get k values
  auto t_K_at_pts = get_K_at_pts(); // tensor of second derivatives
  // get coordinate at integration points
  auto t_coords = OP::getFTensor1CoordsAtGaussPts();
  // loop over integration points
  for (int gg = 0; gg != OP::nbIntegrationPts; gg++) {
    // take into account Jacobian
    const double alpha = t_w * vol * t_K_at_pts(0);
    // loop over rows base functions
    int rr = 0;
    for (; rr != OP::nbRows / FIELD_DIM; rr++) {
 
      auto t_mat = getFTensor2FromArray<FIELD_DIM, FIELD_DIM, FIELD_DIM>(
          OP::locMat, rr * FIELD_DIM, 0);
      auto t_col_grad = col_data.getFTensor2DiffN2<SPACE_DIM>(gg, 0);
 
      // calculate element of local matrix
      for (int cc = 0; cc != OP::nbCols / FIELD_DIM; cc++) {
        t_mat(i, j) +=
            alpha * (t_row_grad(J, K) * t_col_grad(J, K)) * t_kd(i, j);
        ++t_col_grad; // move to another element of gradient of base
                      // function on column
        ++t_mat;      // move to another element of local matrix
      }
      ++t_row_grad; // move to another element of gradient of base
                    // function on row
    }
    for (; rr < OP::nbRowBaseFunctions; ++rr)
      ++t_row_grad;
 
    ++t_coords;
    ++t_K_at_pts;
    ++t_w; // move to another integration weight
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

matK

template<int FIELD_DIM, int SPACE_DIM, AssemblyType A, typename OpBase >

boost::shared_ptr<MatrixDouble> MatOps::OpLhsHuHu< FIELD_DIM, SPACE_DIM, A, GAUSS, OpBase >::matK

protected

Definition at line 199 of file MatHuHu.hpp.

---

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

• src/materials/MatHuHu.hpp