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

#include <users_modules/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, DataForcesAndSourcesCore::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 349 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 355 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,
		DataForcesAndSourcesCore::EntData &	data
	)

inline

Definition at line 363 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 361 of file BCs_RVELagrange_Trac.hpp.

◆ commonData

CommonData& BCs_RVELagrange_Trac::OpRVEBCsRhs_Assemble::commonData

Definition at line 352 of file BCs_RVELagrange_Trac.hpp.

◆ dAta

RVEBC_Data& BCs_RVELagrange_Trac::OpRVEBCsRhs_Assemble::dAta

Definition at line 351 of file BCs_RVELagrange_Trac.hpp.

◆ F

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

Definition at line 353 of file BCs_RVELagrange_Trac.hpp.

◆ f

VectorDouble BCs_RVELagrange_Trac::OpRVEBCsRhs_Assemble::f