#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 163 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 165 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 801 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);
 }