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);
    }