Operator calculate right hand side vector. More...

#include "tutorials/cor-12_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 449 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 453 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 460 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++) {
          CHKERR physicalEqations.calculateTraction(traction, gg, commonData,
                                                    getFEMethod());
          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];
            }
          }
        }
        CHKERR VecSetValues(getFEMethod()->snes_f, data.getIndices().size(),
                            &data.getIndices()[0], &Nf[0], ADD_VALUES);
      } catch (const std::exception &ex) {
        std::ostringstream ss;
        ss << "throw in method: " << ex.what() << std::endl;
        SETERRQ(PETSC_COMM_SELF, MOFEM_STD_EXCEPTION_THROW, "%s",
                ss.str().c_str());
      }
      MoFEMFunctionReturnHot(0);
    }