Operator calculate right hand side vector. More...

#include <users_modules/basic_finite_elements/cohesive_interface/src/CohesiveInterfaceElement.hpp>

Inheritance diagram for CohesiveElement::CohesiveInterfaceElement::OpRhs:

Collaboration diagram for CohesiveElement::CohesiveInterfaceElement::OpRhs:

Public Member Functions
	OpRhs (const std::string field_name, CommonData &common_data, PhysicalEquation &physical_eqations)

MoFEMErrorCode	doWork (int side, EntityType type, EntitiesFieldData::EntData &data)

Public Attributes
CommonData &	commonData

PhysicalEquation &	physicalEqations

VectorDouble	traction

VectorDouble	Nf

Detailed Description

Operator calculate right hand side vector.

Definition at line 431 of file CohesiveInterfaceElement.hpp.

Constructor & Destructor Documentation

◆ OpRhs()

CohesiveElement::CohesiveInterfaceElement::OpRhs::OpRhs	(	const std::string	field_name,
		CommonData &	common_data,
		PhysicalEquation &	physical_eqations
	)

inline

Definition at line 435 of file CohesiveInterfaceElement.hpp.

                                                                                                  :
       FlatPrismElementForcesAndSourcesCore::UserDataOperator(field_name,ForcesAndSourcesCore::UserDataOperator::OPROW),
       commonData(common_data),physicalEqations(physical_eqations) {}

Member Function Documentation

◆ doWork()

MoFEMErrorCode CohesiveElement::CohesiveInterfaceElement::OpRhs::doWork	(	int	side,
		EntityType	type,
		EntitiesFieldData::EntData &	data
	)

inline

Definition at line 440 of file CohesiveInterfaceElement.hpp.

                                                                                    {
       MoFEMFunctionBeginHot;
  
       try {
         int nb_dofs = data.getIndices().size();
         if(nb_dofs == 0) MoFEMFunctionReturnHot(0);
         if(physicalEqations.pRisms.find(getNumeredEntFiniteElementPtr()->getEnt()) == physicalEqations.pRisms.end()) {
           MoFEMFunctionReturnHot(0);
         }
         Nf.resize(nb_dofs);
         Nf.clear();
         int nb_gauss_pts = data.getN().size1();
         for(int gg = 0;gg<nb_gauss_pts;gg++) {
           ierr = physicalEqations.calculateTraction(traction,gg,commonData,getFEMethod()); CHKERRG(ierr);
           double w = getGaussPts()(2,gg)*cblas_dnrm2(3,&getNormalsAtGaussPtsF3()(gg,0),1)*0.5;
           for(int nn = 0;nn<nb_dofs/3;nn++) {
             for(int dd = 0;dd<3;dd++) {
               Nf[3*nn+dd] += w*data.getN(gg)[nn]*traction[dd];
             }
           }
         }
         ierr = VecSetValues(getFEMethod()->snes_f,
         data.getIndices().size(),&data.getIndices()[0],&Nf[0],ADD_VALUES); CHKERRG(ierr);
       } catch (const std::exception& ex) {
         std::ostringstream ss;
         ss << "throw in method: " << ex.what() << std::endl;
         SETERRQ(PETSC_COMM_SELF,1,ss.str().c_str());
       }
       MoFEMFunctionReturnHot(0);
     }