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

#include "src/materials/MatHuHu.hpp"

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

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

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

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

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &data)

Protected Attributes
boost::shared_ptr< MatrixDouble >	matVals
boost::shared_ptr< MatrixDouble >	matK

Detailed Description

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

Definition at line 172 of file MatHuHu.hpp.

Member Typedef Documentation

OP

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

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

Definition at line 175 of file MatHuHu.hpp.

Constructor & Destructor Documentation

OpRhsHu()

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

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

inline

Definition at line 177 of file MatHuHu.hpp.

      : OP(field_name, field_name, OP::OPROW), matVals(mat_vals), matK(mat_K) {}

Member Function Documentation

iNtegrate()

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

MoFEMErrorCode MatOps::OpRhsHu< FIELD_DIM, SPACE_DIM, A, GAUSS, OpBase >::iNtegrate ( EntitiesFieldData::EntData &  data )

protected

Definition at line 225 of file MatHuHu.hpp.

                                        {
  MoFEMFunctionBegin;
 
  FTENSOR_INDEX(FIELD_DIM, i);
  FTENSOR_INDEX(SPACE_DIM, J);
  FTENSOR_INDEX(SPACE_DIM, K);
 
  auto get_val_grad_at_pts = MatrixSizeHelper<
      GetFTensor3FromMatType<FIELD_DIM, SPACE_DIM, SPACE_DIM, -1, DL>,
      DL>::get(*matVals, OP::nbIntegrationPts);
  auto get_K_at_pts =
      MatrixSizeHelper<GetFTensor1FromMatType<1, -1, DL>, DL>::get(
          *matK, 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_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);
    // calculate rhs
    auto t_nf = getFTensor1FromArray<FIELD_DIM, FIELD_DIM>(OP::locF);
    // loop over rows base functions
    int rr = 0;
    for (; rr != OP::nbRows / FIELD_DIM; rr++) {
      // calculate element of local rhs vector
      t_nf(i) += alpha * (t_row_grad(J, K) * t_val_grad_at_pts(i, J, K));
      ++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_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::OpRhsHu< FIELD_DIM, SPACE_DIM, A, GAUSS, OpBase >::matK

protected

Definition at line 185 of file MatHuHu.hpp.

matVals

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

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

protected

Definition at line 184 of file MatHuHu.hpp.

---

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

• src/materials/MatHuHu.hpp