#include <src/finite_elements/BiLinearFormsIntegratorsImpl.hpp>

Inheritance diagram for MoFEM::OpMassImpl< 3, 9, GAUSS, OpBase >:

Collaboration diagram for MoFEM::OpMassImpl< 3, 9, GAUSS, OpBase >:

Public Member Functions
	OpMassImpl (const std::string row_field_name, const std::string col_field_name, ScalarFun beta=scalar_fun_one, boost::shared_ptr< Range > ents_ptr=nullptr)

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Protected Attributes
ScalarFun	betaCoeff

Detailed Description

template<typename OpBase>
struct MoFEM::OpMassImpl< 3, 9, GAUSS, OpBase >

Definition at line 143 of file BiLinearFormsIntegratorsImpl.hpp.

Constructor & Destructor Documentation

◆ OpMassImpl()

template<typename OpBase >

MoFEM::OpMassImpl< 3, 9, GAUSS, OpBase >::OpMassImpl	(	const std::string	row_field_name,
		const std::string	col_field_name,
		ScalarFun	beta = `scalar_fun_one`,
		boost::shared_ptr< Range >	ents_ptr = `nullptr`
	)

inline

Definition at line 144 of file BiLinearFormsIntegratorsImpl.hpp.

       : 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< 3, 9, GAUSS, OpBase >::iNtegrate	(	EntitiesFieldData::EntData &	row_data,
		EntitiesFieldData::EntData &	col_data
	)

protected

Definition at line 752 of file BiLinearFormsIntegratorsImpl.hpp.

                                         {
   MoFEMFunctionBegin;
   FTensor::Index<'i', 3> i;
   FTensor::Index<'k', 3> k;
   auto get_t_vec = [&](const int rr) {
     return FTensor::Tensor1<FTensor::PackPtr<double *, 3>, 3>{
         &OpBase::locMat(rr + 0, 0), &OpBase::locMat(rr + 1, 1),
         &OpBase::locMat(rr + 2, 2)};
   };
   size_t nb_base_functions = row_data.getN().size2() / 3;
   // // get element volume
   const double vol = OpBase::getMeasure();
   // get integration weights
   auto t_w = OpBase::getFTensor0IntegrationWeight();
   // get base function gradient on rows
   auto t_row_base = row_data.getFTensor1N<3>();
   // 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 = vol * 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 / 3; rr++) {
       // get column base functions gradient at gauss point gg
       auto t_col_base = col_data.getFTensor1N<3>(gg, 0);
       auto t_vec = get_t_vec(3 * rr);
       // loop over columns
       for (int cc = 0; cc != OpBase::nbCols / 3; cc++) {
         // calculate element of local matrix
         t_vec(i) += alpha * (t_row_base(k) * t_col_base(k));
         ++t_col_base;
         ++t_vec;
       }
       ++t_row_base;
     }
     for (; rr < nb_base_functions; ++rr)
       ++t_row_base;
     ++t_coords;
     ++t_w; // move to another integration weight
   }
   MoFEMFunctionReturn(0);
 };