\biref operator to calculate the RHS of the constrain for the RVE boundary conditions More...

#include "users_modules/mofem_um_homogenisation/src/BCs_RVELagrange_Trac.hpp"

Inheritance diagram for BCs_RVELagrange_Trac::OpRVEBCsRhs_Assemble:

Collaboration diagram for BCs_RVELagrange_Trac::OpRVEBCsRhs_Assemble:

Public Member Functions
	OpRVEBCsRhs_Assemble (const string lagrang_field_name, vector< Vec > &f, RVEBC_Data &data, CommonData &common_data)

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

Public Attributes
RVEBC_Data &	dAta

CommonData &	commonData

vector< Vec > &	F

VectorDouble	applied_strain

VectorDouble	f

Detailed Description

\biref operator to calculate the RHS of the constrain for the RVE boundary conditions

Definition at line 337 of file BCs_RVELagrange_Trac.hpp.

Constructor & Destructor Documentation

◆ OpRVEBCsRhs_Assemble()

BCs_RVELagrange_Trac::OpRVEBCsRhs_Assemble::OpRVEBCsRhs_Assemble	(	const string	lagrang_field_name,
		vector< Vec > &	f,
		RVEBC_Data &	data,
		CommonData &	common_data
	)

inline

Definition at line 343 of file BCs_RVELagrange_Trac.hpp.

                                                                                                                 :
    FaceElementForcesAndSourcesCore::UserDataOperator(lagrang_field_name,UserDataOperator::OPROW),
    dAta(data),commonData(common_data),F(f) {
    }

Member Function Documentation

◆ doWork()

PetscErrorCode BCs_RVELagrange_Trac::OpRVEBCsRhs_Assemble::doWork	(	int	side,
		EntityType	type,
		DataForcesAndSurcesCore::EntData &	data
	)

inline

Definition at line 351 of file BCs_RVELagrange_Trac.hpp.

                                                                                         {
      PetscFunctionBegin;
 
      if(data.getIndices().size()==0) PetscFunctionReturn(0);
 
      MatrixDouble &D_mat=commonData.D_mat;
      applied_strain.resize(D_mat.size2(),false);
 
      for(int ii = 0;ii<F.size();ii++) {
        applied_strain.clear();
        applied_strain(ii) = 1.0;
        f.resize(D_mat.size1(),false);
        noalias(f) = prod(D_mat,applied_strain);
        ierr = VecSetValues(F[ii],data.getIndices().size(),&data.getIndices()[0],&f[0],ADD_VALUES); CHKERRQ(ierr);
      }
 
      PetscFunctionReturn(0);
    }

Member Data Documentation

◆ applied_strain

VectorDouble BCs_RVELagrange_Trac::OpRVEBCsRhs_Assemble::applied_strain

Definition at line 349 of file BCs_RVELagrange_Trac.hpp.

◆ commonData

CommonData& BCs_RVELagrange_Trac::OpRVEBCsRhs_Assemble::commonData

Definition at line 340 of file BCs_RVELagrange_Trac.hpp.

◆ dAta

RVEBC_Data& BCs_RVELagrange_Trac::OpRVEBCsRhs_Assemble::dAta

Definition at line 339 of file BCs_RVELagrange_Trac.hpp.

◆ F

vector<Vec>& BCs_RVELagrange_Trac::OpRVEBCsRhs_Assemble::F

Definition at line 341 of file BCs_RVELagrange_Trac.hpp.

◆ f

VectorDouble BCs_RVELagrange_Trac::OpRVEBCsRhs_Assemble::f