#include "src/finite_elements/BiLinearFormsIntegratorsImpl.hpp"

Inheritance diagram for MoFEM::OpMassCacheImpl< 1, FIELD_DIM, I, OpBase >:

Collaboration diagram for MoFEM::OpMassCacheImpl< 1, FIELD_DIM, I, OpBase >:

Public Member Functions
	OpMassCacheImpl (CacheMatsTypeType cache, const std::string row_field_name, const std::string col_field_name, const double beta, boost::shared_ptr< Range > ents_ptr=nullptr)

MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Protected Attributes
double	scalarBeta

CacheMatsTypeType	cacheLocMats

Detailed Description

template<int FIELD_DIM, IntegrationType I, typename OpBase>
struct MoFEM::OpMassCacheImpl< 1, FIELD_DIM, I, OpBase >

Definition at line 185 of file BiLinearFormsIntegratorsImpl.hpp.

Constructor & Destructor Documentation

◆ OpMassCacheImpl()

template<int FIELD_DIM, IntegrationType I, typename OpBase >

MoFEM::OpMassCacheImpl< 1, FIELD_DIM, I, OpBase >::OpMassCacheImpl	(	CacheMatsTypeType	cache,
		const std::string	row_field_name,
		const std::string	col_field_name,
		const double	beta,
		boost::shared_ptr< Range >	ents_ptr = `nullptr`
	)

inline

Definition at line 187 of file BiLinearFormsIntegratorsImpl.hpp.

      : OpMassImpl<1, FIELD_DIM, I, OpBase>(
            row_field_name, col_field_name,
            [](double, double, double) constexpr { return 1; }, ents_ptr),
        cacheLocMats(cache), scalarBeta(beta) {}

Member Function Documentation

◆ iNtegrate()

template<int FIELD_DIM, IntegrationType I, typename OpBase >

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

Definition at line 887 of file BiLinearFormsIntegratorsImpl.hpp.

                                        {
  MoFEMFunctionBegin;
 
  const auto vol = this->getMeasure();
  const auto row_type = this->rowType;
  const auto col_type = this->colType;
  auto &loc_mat = this->locMat;
 
  auto p = std::make_pair(row_type, col_type);
 
  if (cacheLocMats[p]) {
    if (cacheLocMats[p]->size1() != loc_mat.size1() &&
        cacheLocMats[p]->size2() != loc_mat.size2()) {
      cacheLocMats[p]->resize(loc_mat.size1(), loc_mat.size2());
      CHKERR this->integrateImpl(row_data, col_data, 1);
      *(cacheLocMats[p]) = loc_mat;
    } else {
      loc_mat = *(cacheLocMats[p]);
    }
    loc_mat *= scalarBeta * this->getMeasure();
  } else {
    CHKERR this->integrateImpl(row_data, col_data,
                               scalarBeta * this->getMeasure());
  }
  MoFEMFunctionReturn(0);
}