MoFEM::OpMassImpl< 1, FIELD_DIM, GAUSS, OpBase > Struct Template Reference

#include <src/finite_elements/BiLinearFormsIntegratorsImpl.hpp>

Inheritance diagram for MoFEM::OpMassImpl< 1, FIELD_DIM, GAUSS, OpBase >:

Collaboration diagram for MoFEM::OpMassImpl< 1, FIELD_DIM, GAUSS, OpBase >:

Protected Member Functions
MoFEMErrorCode	integrateImpl (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data, double vol)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)
Protected Member Functions inherited from MoFEM::OpMassImpl< 1, 1, GAUSS, OpBase >
MoFEMErrorCode	integrateImpl (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data, double vol)
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Additional Inherited Members
Public Member Functions inherited from MoFEM::OpMassImpl< 1, 1, GAUSS, OpBase >
	OpMassImpl (const std::string row_field_name, const std::string col_field_name, ScalarFun beta=[](double, double, double) constexpr { return 1;}, boost::shared_ptr< Range > ents_ptr=nullptr, boost::shared_ptr< MatrixDouble > cache_mat=nullptr)
Protected Attributes inherited from MoFEM::OpMassImpl< 1, 1, GAUSS, OpBase >
ScalarFun	betaCoeff = [](double, double, double) constexpr { return 1; }

Detailed Description

template<int FIELD_DIM, typename OpBase>
struct MoFEM::OpMassImpl< 1, FIELD_DIM, GAUSS, OpBase >

Definition at line 92 of file BiLinearFormsIntegratorsImpl.hpp.

Member Function Documentation

iNtegrate()

template<int FIELD_DIM, typename OpBase >

MoFEMErrorCode MoFEM::OpMassImpl< 1, FIELD_DIM, GAUSS, OpBase >::iNtegrate ( EntitiesFieldData::EntData &  row_data, EntitiesFieldData::EntData &  col_data  )

inlineprotected

Definition at line 101 of file BiLinearFormsIntegratorsImpl.hpp.

                                                                      {
    return integrateImpl(row_data, col_data, OpBase::getMeasure());
  }

integrateImpl()

template<int FIELD_DIM, typename OpBase >

MoFEMErrorCode MoFEM::OpMassImpl< 1, FIELD_DIM, GAUSS, OpBase >::integrateImpl ( EntitiesFieldData::EntData &  row_data, EntitiesFieldData::EntData &  col_data, double  vol  )

protected

Definition at line 602 of file BiLinearFormsIntegratorsImpl.hpp.

                {
  MoFEMFunctionBegin;
 
  auto integrate = [&](auto &mat) {
    MoFEMFunctionBeginHot;
    // get integration weights
    auto t_w = OpBase::getFTensor0IntegrationWeight();
    // get base function gradient on rows
    auto t_row_base = row_data.getFTensor0N();
    // get coordinate at integration points
    auto t_coords = OpBase::getFTensor1CoordsAtGaussPts();
 
    FTensor::Index<'i', FIELD_DIM> i;
    auto get_t_vec = [&](const int rr) {
      std::array<double *, FIELD_DIM> ptrs;
      for (auto i = 0; i != FIELD_DIM; ++i)
        ptrs[i] = &mat(rr + i, i);
      return FTensor::Tensor1<FTensor::PackPtr<double *, FIELD_DIM>, FIELD_DIM>(
          ptrs);
    };
 
    // loop over integration points
    for (int gg = 0; gg != OpBase::nbIntegrationPts; gg++) {
      const double beta = betaCoeff(t_coords(0), t_coords(1), t_coords(2));
      // take into account Jacobian
      const double alpha = t_w * beta;
      // loop over rows base functions
      int rr = 0;
      for (; rr != OpBase::nbRows / FIELD_DIM; rr++) {
        // get column base functions gradient at gauss point gg
        auto t_col_base = col_data.getFTensor0N(gg, 0);
        // get mat vec
        auto t_vec = get_t_vec(FIELD_DIM * rr);
        // loop over columns
        for (int cc = 0; cc != OpBase::nbCols / FIELD_DIM; cc++) {
          // calculate element of local matrix
          t_vec(i) += alpha * (t_row_base * t_col_base);
          ++t_col_base;
          ++t_vec;
        }
        ++t_row_base;
      }
      for (; rr < OpBase::nbRowBaseFunctions; ++rr)
        ++t_row_base;
      ++t_coords;
      ++t_w; // move to another integration weight
    }
    MoFEMFunctionReturnHot(0);
  };
 
  CHKERR integrate(OpBase::locMat);
  OpBase::locMat *= vol;
 
  MoFEMFunctionReturn(0);
};

---

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

• src/finite_elements/BiLinearFormsIntegratorsImpl.hpp