Specialization for TemperatureCubitBcData. More...

#include <src/boundary_conditions/EssentialTemperatureCubitBcData.hpp>

Collaboration diagram for MoFEM::EssentialPreProc< TemperatureCubitBcData >:

Public Member Functions
	EssentialPreProc (MoFEM::Interface &m_field, boost::shared_ptr< FEMethod > fe_ptr, std::vector< boost::shared_ptr< ScalingMethod >> smv)

MoFEMErrorCode	operator() ()

Protected Attributes
MoFEM::Interface &	mField

boost::weak_ptr< FEMethod >	fePtr

VecOfTimeScalingMethods	vecOfTimeScalingMethods

Detailed Description

Specialization for TemperatureCubitBcData.

Specialization to enforce blocksets which TemperatureCubitBcData ptr. That is to enforce constrains on temperature. set

Template Parameters

Examples: poisson_2d_homogeneous.cpp.

Definition at line 91 of file EssentialTemperatureCubitBcData.hpp.

Constructor & Destructor Documentation

◆ EssentialPreProc()

MoFEM::EssentialPreProc< TemperatureCubitBcData >::EssentialPreProc	(	MoFEM::Interface &	m_field,
		boost::shared_ptr< FEMethod >	fe_ptr,
		std::vector< boost::shared_ptr< ScalingMethod >>	smv
	)

Definition at line 15 of file EssentialTemperatureCubitBcData.cpp.

18 : mField(m_field), fePtr(fe_ptr), vecOfTimeScalingMethods(smv) {}

Member Function Documentation

◆ operator()()

MoFEMErrorCode MoFEM::EssentialPreProc< TemperatureCubitBcData >::operator() ( )

Definition at line 20 of file EssentialTemperatureCubitBcData.cpp.

                                                                     {
   MoFEMFunctionBegin;
  
   if (auto fe_method_ptr = fePtr.lock()) {
  
     auto bc_mng = mField.getInterface<BcManager>();
     auto fb = mField.getInterface<FieldBlas>();
     const auto problem_name = fe_method_ptr->problemPtr->getName();
  
     for (auto bc : bc_mng->getBcMapByBlockName()) {
       if (auto temp_bc = bc.second->tempBcPtr) {
  
         auto &bc_id = bc.first;
  
         auto regex_str = (boost::format("%s_(.*)") % problem_name).str();
         if (std::regex_match(bc_id, std::regex(regex_str))) {
  
           auto [field_name, block_name] =
               BcManager::extractStringFromBlockId(bc_id, problem_name);
  
           MOFEM_LOG("WORLD", Sev::noisy)
               << "Apply EssentialPreProc<TemperatureCubitBcData>: "
               << problem_name << "_" << field_name << "_" << block_name;
  
           auto verts = bc.second->bcEnts.subset_by_type(MBVERTEX);
           auto v = temp_bc->data.value1;
           for (auto s : vecOfTimeScalingMethods) {
             v *= s->getScale(fe_method_ptr->ts_t);
           }
  
           auto lambda = [&](boost::shared_ptr<FieldEntity> field_entity_ptr) {
             MoFEMFunctionBegin;
             for (auto &d : field_entity_ptr->getEntFieldData())
               d = v;
             MoFEMFunctionReturn(0);
           };
           CHKERR fb->fieldLambdaOnEntities(lambda, field_name, &verts);
         }
       }
     }
  
   } else {
     SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "Can not lock shared pointer");
   }
  
   MoFEMFunctionReturn(0);
 }