#include <tutorials/adv-0/src/PlasticOpsGeneric.hpp>

Inheritance diagram for PlasticOps::OpCalculateConstraintsRhsImpl< GAUSS, AssemblyDomainEleOp >:

Collaboration diagram for PlasticOps::OpCalculateConstraintsRhsImpl< GAUSS, AssemblyDomainEleOp >:

Public Member Functions
	OpCalculateConstraintsRhsImpl (const std::string field_name, boost::shared_ptr< CommonData > common_data_ptr, boost::shared_ptr< MatrixDouble > m_D_ptr)

MoFEMErrorCode	iNtegrate (EntitiesFieldData::EntData &data)

Private Attributes
boost::shared_ptr< CommonData >	commonDataPtr

boost::shared_ptr< MatrixDouble >	mDPtr

Detailed Description

template<typename AssemblyDomainEleOp>
struct PlasticOps::OpCalculateConstraintsRhsImpl< GAUSS, AssemblyDomainEleOp >

Definition at line 877 of file PlasticOpsGeneric.hpp.

Constructor & Destructor Documentation

◆ OpCalculateConstraintsRhsImpl()

template<typename AssemblyDomainEleOp >

PlasticOps::OpCalculateConstraintsRhsImpl< GAUSS, AssemblyDomainEleOp >::OpCalculateConstraintsRhsImpl	(	const std::string	field_name,
		boost::shared_ptr< CommonData >	common_data_ptr,
		boost::shared_ptr< MatrixDouble >	m_D_ptr
	)

Definition at line 891 of file PlasticOpsGeneric.hpp.

894 : AssemblyDomainEleOp(field_name, field_name, AssemblyDomainEleOp::OPROW),

895 commonDataPtr(common_data_ptr), mDPtr(m_D_ptr) {}

Member Function Documentation

◆ iNtegrate()

template<typename AssemblyDomainEleOp >

MoFEMErrorCode PlasticOps::OpCalculateConstraintsRhsImpl< GAUSS, AssemblyDomainEleOp >::iNtegrate ( EntitiesFieldData::EntData & data )

Examples: PlasticOpsGeneric.hpp.

Definition at line 899 of file PlasticOpsGeneric.hpp.

                                     {
   MoFEMFunctionBegin;
  
   const size_t nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
   const size_t nb_base_functions = data.getN().size2();
  
   auto t_res_c = getFTensor0FromVec(commonDataPtr->resC);
  
   auto next = [&]() { ++t_res_c; };
  
   auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
   auto &nf = AssemblyDomainEleOp::locF;
   auto t_base = data.getFTensor0N();
   for (size_t gg = 0; gg != nb_integration_pts; ++gg) {
     const double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
     ++t_w;
     const auto res = alpha * t_res_c;
     next();
  
     size_t bb = 0;
     for (; bb != AssemblyDomainEleOp::nbRows; ++bb) {
       nf[bb] += t_base * res;
       ++t_base;
     }
     for (; bb < nb_base_functions; ++bb)
       ++t_base;
   }
  
   MoFEMFunctionReturn(0);
 }