\biref operator to calculate the RHS for the calculation of the homgoensied stiffness matrix More...

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

Inheritance diagram for BCs_RVELagrange_Disp::OpRVEBCsRhsHomoC:

Collaboration diagram for BCs_RVELagrange_Disp::OpRVEBCsRhsHomoC:

Public Member Functions
	OpRVEBCsRhsHomoC (const string field_name, const string lagrang_field_name, vector< Vec > &f, RVEBC_Data &data, bool ho_geometry=false)

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

Public Member Functions inherited from BCs_RVELagrange_Disp::OpDmatRhs
	OpDmatRhs (const string field_name, const string lagrang_field_name, RVEBC_Data &data, bool ho_geometry=false)

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

Public Attributes
vector< Vec > &	F

Public Attributes inherited from BCs_RVELagrange_Disp::OpDmatRhs
RVEBC_Data &	dAta

bool	hoGeometry

VectorDouble	f

VectorDouble	applied_strain

MatrixDouble	X_mat

MatrixDouble	N_mat

MatrixDouble	D_mat

int	rank

int	nb_row_dofs

Detailed Description

\biref operator to calculate the RHS for the calculation of the homgoensied stiffness matrix

Definition at line 365 of file BCs_RVELagrange_Disp.hpp.

Constructor & Destructor Documentation

◆ OpRVEBCsRhsHomoC()

BCs_RVELagrange_Disp::OpRVEBCsRhsHomoC::OpRVEBCsRhsHomoC	(	const string	field_name,
		const string	lagrang_field_name,
		vector< Vec > &	f,
		RVEBC_Data &	data,
		bool	ho_geometry = `false`
	)

inline

Definition at line 369 of file BCs_RVELagrange_Disp.hpp.

                  :
     OpDmatRhs(field_name,lagrang_field_name,data,ho_geometry),F(f) {
     }

Member Function Documentation

◆ doWork()

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

inline

Definition at line 376 of file BCs_RVELagrange_Disp.hpp.

                                                                                           {
       PetscFunctionBegin;
       
       nb_row_dofs = data.getIndices().size();
       
       if(nb_row_dofs ==0) PetscFunctionReturn(0);
       if(dAta.tRis.find(getNumeredEntFiniteElementPtr()->getEnt())==dAta.tRis.end()) PetscFunctionReturn(0);
       
       auto weak_ptr_dof =
           getNumeredEntFiniteElementPtr()->getRowDofsByPetscGlobalDofIdx(
               data.getIndices()[0]);
       const FENumeredDofEntity *dof_ptr;
       if (auto ptr = weak_ptr_dof.lock())
         dof_ptr = ptr.get();
       else
         SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "Dof not found");
  
       rank = dof_ptr->getNbOfCoeffs();
       nb_row_dofs /= rank;
       
       ierr = calculateDmat(side,type,data); CHKERRQ(ierr);
       
       applied_strain.resize(D_mat.size2(),false);
       f.resize(D_mat.size1(),false);
       
       int size = (rank == 3) ? 6 : 3;
       for(int ii = 0;ii<size;ii++) {
         applied_strain.clear();
         applied_strain[ii] = 1.0;
         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

◆ F

vector<Vec>& BCs_RVELagrange_Disp::OpRVEBCsRhsHomoC::F

Definition at line 367 of file BCs_RVELagrange_Disp.hpp.

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

users_modules/homogenisation/src/BCs_RVELagrange_Disp.hpp

Public Member Functions