BCs_RVELagrange_Trac::OpRVEHomoStress_Assemble Struct Reference

#include "users_modules/mofem_um_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, DataForcesAndSurcesCore::EntData &data)

Public Attributes
RVEBC_Data &	dAta
CommonData &	commonData
Vec	stressHomo
VectorDouble	stressVector

Detailed Description

Definition at line 396 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 402 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, DataForcesAndSurcesCore::EntData &  data  )

inline

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

dAta

RVEBC_Data& BCs_RVELagrange_Trac::OpRVEHomoStress_Assemble::dAta

Definition at line 398 of file BCs_RVELagrange_Trac.hpp.

stressHomo

Vec BCs_RVELagrange_Trac::OpRVEHomoStress_Assemble::stressHomo

Definition at line 400 of file BCs_RVELagrange_Trac.hpp.

stressVector

VectorDouble BCs_RVELagrange_Trac::OpRVEHomoStress_Assemble::stressVector

Definition at line 410 of file BCs_RVELagrange_Trac.hpp.

---

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

• users_modules/mofem_um_homogenisation/src/BCs_RVELagrange_Trac.hpp