#include <users_modules/homogenisation/src/BCs_RVELagrange_Trac.hpp>

Inheritance diagram for BCs_RVELagrange_Trac::OpRVEHomoStress_Assemble:

Collaboration diagram for BCs_RVELagrange_Trac::OpRVEHomoStress_Assemble:

Public Member Functions
	OpRVEHomoStress_Assemble (const string lagrang_field_name, Vec stress_homo, RVEBC_Data &data, CommonData &common_data)

PetscErrorCode	doWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)

Public Attributes
RVEBC_Data &	dAta

CommonData &	commonData

Vec	stressHomo

VectorDouble	stressVector

Detailed Description

Definition at line 408 of file BCs_RVELagrange_Trac.hpp.

Constructor & Destructor Documentation

◆ OpRVEHomoStress_Assemble()

BCs_RVELagrange_Trac::OpRVEHomoStress_Assemble::OpRVEHomoStress_Assemble	(	const string	lagrang_field_name,
		Vec	stress_homo,
		RVEBC_Data &	data,
		CommonData &	common_data
	)

inline

Definition at line 414 of file BCs_RVELagrange_Trac.hpp.

      :
     FaceElementForcesAndSourcesCore::UserDataOperator(lagrang_field_name,UserDataOperator::OPROW),
     dAta(data),commonData(common_data),stressHomo(stress_homo) {
     }

Member Function Documentation

◆ doWork()

PetscErrorCode BCs_RVELagrange_Trac::OpRVEHomoStress_Assemble::doWork	(	int	side,
		EntityType	type,
		DataForcesAndSourcesCore::EntData &	data
	)

inline

Definition at line 424 of file BCs_RVELagrange_Trac.hpp.

                                                                                           {
       PetscFunctionBegin;
       if(data.getIndices().size()==0) PetscFunctionReturn(0);
       if(dAta.tRis.find(getNumeredEntFiniteElementPtr()->getEnt())==dAta.tRis.end()) PetscFunctionReturn(0);
       stressVector.resize(commonData.D_mat.size2());
       //Lamda=data.getFieldData() is reaction force (so multiply for -1 to get the force)
       noalias(stressVector) = prod(trans(commonData.D_mat),-data.getFieldData());
       const int indices_6[6] = {0, 1, 2, 3, 4, 5};
       const int indices_3[3] = {0, 1, 2};
       switch(stressVector.size()) {
         case 6:
           ierr = VecSetValues(stressHomo,6,indices_6,&stressVector[0],ADD_VALUES); CHKERRQ(ierr);
           break;
         case 3:
           ierr = VecSetValues(stressHomo,3,indices_3,&stressVector[0],ADD_VALUES); CHKERRQ(ierr);
           break;
         default:
           SETERRQ(PETSC_COMM_SELF,MOFEM_NOT_IMPLEMENTED,"not implemented");
       }
       PetscFunctionReturn(0);
     }

Member Data Documentation

◆ commonData

CommonData& BCs_RVELagrange_Trac::OpRVEHomoStress_Assemble::commonData

Definition at line 411 of file BCs_RVELagrange_Trac.hpp.

◆ dAta

RVEBC_Data& BCs_RVELagrange_Trac::OpRVEHomoStress_Assemble::dAta

Definition at line 410 of file BCs_RVELagrange_Trac.hpp.

◆ stressHomo

Vec BCs_RVELagrange_Trac::OpRVEHomoStress_Assemble::stressHomo

Definition at line 412 of file BCs_RVELagrange_Trac.hpp.

◆ stressVector

VectorDouble BCs_RVELagrange_Trac::OpRVEHomoStress_Assemble::stressVector

Definition at line 422 of file BCs_RVELagrange_Trac.hpp.

The documentation for this struct was generated from the following file:

users_modules/homogenisation/src/BCs_RVELagrange_Trac.hpp

Public Member Functions

Public Attributes