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

#include "src/materials/MatHuHu.hpp"

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

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

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

Public Member Functions
	OpLhsHuGrad (const std::string field_name, boost::shared_ptr< MatrixDouble > mat_vals, boost::shared_ptr< MatrixDouble > mat_diff_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 >	matVals
boost::shared_ptr< MatrixDouble >	matDiffK

Detailed Description

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

Definition at line 205 of file MatHuHu.hpp.

Member Typedef Documentation

OP

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

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

Definition at line 207 of file MatHuHu.hpp.

Constructor & Destructor Documentation

OpLhsHuGrad()

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

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

inline

Definition at line 208 of file MatHuHu.hpp.

      : OP(field_name, field_name, OP::OPROWCOL), matVals(mat_vals),
        matDiffK(mat_diff_K) {
    OP::sYmm = false;
  }

Member Function Documentation

iNtegrate()

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

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

protected

Definition at line 339 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);
  FTENSOR_INDEX(SPACE_DIM, M);
 
  auto get_val_grad_at_pts = MatrixSizeHelper<
      GetFTensor3FromMatType<FIELD_DIM, SPACE_DIM, SPACE_DIM, -1, DL>,
      DL>::get(*matVals, OP::nbIntegrationPts);
  auto get_diff_K_at_pts =
      MatrixSizeHelper<GetFTensor2FromMatType<FIELD_DIM, SPACE_DIM, -1, DL>,
                       DL>::get(*matDiffK, OP::nbIntegrationPts);
 
  // 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 field gradient values
  auto t_val_grad_at_pts =
      get_val_grad_at_pts(); // tensor of second derivatives
  // get K values
  auto t_diff_K_at_pts = get_diff_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;
    // 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.getFTensor1DiffN<SPACE_DIM>(gg, 0);
      for (int bb = 0; bb != OP::nbCols / FIELD_DIM; ++bb) {
        t_mat(i, j) += alpha * (t_row_grad(J, K) * t_val_grad_at_pts(i, J, K)) *
                       (t_diff_K_at_pts(j, M) * t_col_grad(M));
 
        ++t_mat;      // move to another element of local matrix
        ++t_col_grad; // move to another element of gradient of base
                      // function on column
      }
 
      ++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_val_grad_at_pts;
    ++t_diff_K_at_pts;
    ++t_w; // move to another integration weight
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

matDiffK

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

boost::shared_ptr<MatrixDouble> MatOps::OpLhsHuGrad< FIELD_DIM, SPACE_DIM, A, GAUSS, OpBase >::matDiffK

protected

Definition at line 219 of file MatHuHu.hpp.

matVals

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

boost::shared_ptr<MatrixDouble> MatOps::OpLhsHuGrad< FIELD_DIM, SPACE_DIM, A, GAUSS, OpBase >::matVals

protected

Definition at line 218 of file MatHuHu.hpp.

---

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

• src/materials/MatHuHu.hpp