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

#include <src/finite_elements/BiLinearFormsIntegratorsImpl.hpp>

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

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

Public Member Functions
	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 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 Attributes
ScalarFun	betaCoeff = [](double, double, double) constexpr { return 1; }

Detailed Description

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

Definition at line 67 of file BiLinearFormsIntegratorsImpl.hpp.

Constructor & Destructor Documentation

OpMassImpl()

template<typename OpBase >

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  )

inline

Definition at line 69 of file BiLinearFormsIntegratorsImpl.hpp.

                                                            { return 1; },
      boost::shared_ptr<Range> ents_ptr = nullptr,
      boost::shared_ptr<MatrixDouble> cache_mat = nullptr)
      : OpBase(row_field_name, col_field_name, OpBase::OPROWCOL, ents_ptr),
        betaCoeff(beta) {
    if (row_field_name == col_field_name)
      this->sYmm = true;
  }

Member Function Documentation

iNtegrate()

template<typename OpBase >

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

inlineprotected

Definition at line 85 of file BiLinearFormsIntegratorsImpl.hpp.

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

integrateImpl()

template<typename OpBase >

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

protected

Definition at line 525 of file BiLinearFormsIntegratorsImpl.hpp.

                                                    {
  MoFEMFunctionBegin;
 
#ifndef NDEBUG
  auto log_error = [&]() {
    MOFEM_LOG("SELF", Sev::error) << "Row side " << OpBase::rowSide << " "
                                  << CN::EntityTypeName(OpBase::rowType);
    MOFEM_LOG("SELF", Sev::error) << "Col side " << OpBase::colSide << " "
                                  << CN::EntityTypeName(OpBase::colType);
  };
 
  if (row_data.getN().size2() < OpBase::nbRows) {
    log_error();
    SETERRQ2(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "Wrong number of base functions on rows %d < %d",
             row_data.getN().size2(), OpBase::nbRows);
  }
  if (col_data.getN().size2() < OpBase::nbCols) {
    log_error();
    SETERRQ2(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "Wrong number of base functions on cols %d < %d",
             col_data.getN().size2(), OpBase::nbCols);
  }
  if (row_data.getN().size1() != OpBase::nbIntegrationPts) {
    log_error();
    SETERRQ2(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "Wrong number of integration points on rows %d != %d",
             row_data.getN().size1(), OpBase::nbIntegrationPts);
  }
  if (col_data.getN().size1() != OpBase::nbIntegrationPts) {
    log_error();
    SETERRQ2(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "Wrong number of integration points on cols %d != %d",
             col_data.getN().size1(), OpBase::nbIntegrationPts);
  }
#endif
 
  // 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();
  // 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
    auto a_mat_ptr = &*OpBase::locMat.data().begin();
    int rr = 0;
    for (; rr != OpBase::nbRows; rr++) {
      // get column base functions gradient at gauss point gg
      auto t_col_base = col_data.getFTensor0N(gg, 0);
      // loop over columns
      for (int cc = 0; cc != OpBase::nbCols; cc++) {
        // calculate element of local matrix
        *a_mat_ptr += alpha * (t_row_base * t_col_base);
        ++t_col_base;
        ++a_mat_ptr;
      }
      ++t_row_base;
    }
    for (; rr < OpBase::nbRowBaseFunctions; ++rr)
      ++t_row_base;
    ++t_coords;
    ++t_w; // move to another integration weight
  }
 
  OpBase::locMat *= vol;
 
  MoFEMFunctionReturn(0);
};

Member Data Documentation

betaCoeff

template<typename OpBase >

ScalarFun MoFEM::OpMassImpl< 1, 1, GAUSS, OpBase >::betaCoeff = [](double, double, double) constexpr { return 1; }

protected

Definition at line 81 of file BiLinearFormsIntegratorsImpl.hpp.

---

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

• src/finite_elements/BiLinearFormsIntegratorsImpl.hpp