\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