#include <src/finite_elements/LinearFormsIntegratorsImpl.hpp>

Inheritance diagram for MoFEM::OpGradTimesSymTensorImpl< 1, SPACE_DIM, SPACE_DIM, S, GAUSS, OpBase >:

Collaboration diagram for MoFEM::OpGradTimesSymTensorImpl< 1, SPACE_DIM, SPACE_DIM, S, GAUSS, OpBase >:

Public Member Functions
	OpGradTimesSymTensorImpl (const std::string field_name, boost::shared_ptr< MatrixDouble > mat_vals, ScalarFun beta_coeff=[](double, double, double) constexpr { return 1;})

Protected Member Functions
MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &data)

Protected Attributes
boost::shared_ptr< MatrixDouble >	matVals

FTensor::Index< 'i', SPACE_DIM >	i

FTensor::Index< 'j', SPACE_DIM >	j

ScalarFun	betaCoeff

Detailed Description

template<int SPACE_DIM, int S, typename OpBase>
struct MoFEM::OpGradTimesSymTensorImpl< 1, SPACE_DIM, SPACE_DIM, S, GAUSS, OpBase >

Definition at line 224 of file LinearFormsIntegratorsImpl.hpp.

Constructor & Destructor Documentation

◆ OpGradTimesSymTensorImpl()

template<int SPACE_DIM, int S, typename OpBase >

MoFEM::OpGradTimesSymTensorImpl< 1, SPACE_DIM, SPACE_DIM, S, GAUSS, OpBase >::OpGradTimesSymTensorImpl	(	const std::string	field_name,
		boost::shared_ptr< MatrixDouble >	mat_vals,
		ScalarFun	beta_coeff = `[](double, double, double) constexpr { return 1; }`
	)

inline

Definition at line 227 of file LinearFormsIntegratorsImpl.hpp.

                                                                   { return 1; })
       : OpBase(field_name, field_name, OpBase::OPROW), matVals(mat_vals),
         betaCoeff(beta_coeff) {}

Member Function Documentation

◆ iNtegrate()

template<int SPACE_DIM, int S, typename OpBase >

MoFEMErrorCode MoFEM::OpGradTimesSymTensorImpl< 1, SPACE_DIM, SPACE_DIM, S, GAUSS, OpBase >::iNtegrate ( EntitiesFieldData::EntData & data )

protected

Definition at line 797 of file LinearFormsIntegratorsImpl.hpp.

                                         {
   MoFEMFunctionBegin;
   // get element volume
   const double vol = OpBase::getMeasure();
   // get coords
   auto t_coords = OpBase::getFTensor1CoordsAtGaussPts();
   // get integration weights
   auto t_w = OpBase::getFTensor0IntegrationWeight();
   // get base function gradient on rows
   auto t_row_grad = row_data.getFTensor1DiffN<SPACE_DIM>();
   // get filed gradient values
   auto t_val_mat = getFTensor2SymmetricFromMat<SPACE_DIM, S>(*(matVals));
   // loop over integration points
   for (int gg = 0; gg != OpBase::nbIntegrationPts; gg++) {
     // take into account Jacobian
     const double alpha =
         t_w * vol * betaCoeff(t_coords(0), t_coords(1), t_coords(2));
     // get rhs vector
     auto t_nf = OpBase::template getNf<SPACE_DIM>();
     // loop over rows base functions
     int rr = 0;
     for (; rr != OpBase::nbRows / SPACE_DIM; rr++) {
       // calculate element of local matrix
       t_nf(j) += alpha * (t_row_grad(i) * t_val_mat(i, j));
       ++t_row_grad; // move to another element of gradient of base
                     // function on row
       ++t_nf;
     }
     for (; rr < OpBase::nbRowBaseFunctions; ++rr)
       ++t_row_grad;
     ++t_val_mat;
     ++t_coords;
     ++t_w; // move to another integration weight
   }
   MoFEMFunctionReturn(0);
 }