BCs_RVELagrange_Trac::OpRVEBCsLhs Struct Reference

\biref operator to calculate the LHS for the RVE bounary conditions More...

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

Inheritance diagram for BCs_RVELagrange_Trac::OpRVEBCsLhs:

Collaboration diagram for BCs_RVELagrange_Trac::OpRVEBCsLhs:

Public Member Functions
	OpRVEBCsLhs (const string field_name, const string lagrang_field_name, Mat aij, RVEBC_Data &data, CommonFunctions &common_functions, bool ho_geometry=false)
PetscErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, DataForcesAndSurcesCore::EntData &row_data, DataForcesAndSurcesCore::EntData &col_data)

Public Attributes
RVEBC_Data &	dAta
CommonFunctions &	commonFunctions
bool	hoGeometry
Mat	Aij
MatrixDouble	K
MatrixDouble	transK
MatrixDouble	NTN
ublas::vector< MatrixDouble >	N_mat_row
MatrixDouble	N_mat_col
MatrixDouble	H_mat

Detailed Description

\biref operator to calculate the LHS for the RVE bounary conditions

Definition at line 164 of file BCs_RVELagrange_Trac.hpp.

Constructor & Destructor Documentation

OpRVEBCsLhs()

BCs_RVELagrange_Trac::OpRVEBCsLhs::OpRVEBCsLhs ( const string  field_name, const string  lagrang_field_name, Mat  aij, RVEBC_Data &  data, CommonFunctions &  common_functions, bool  ho_geometry = false  )

inline

Definition at line 171 of file BCs_RVELagrange_Trac.hpp.

     :
    FaceElementForcesAndSourcesCore::UserDataOperator(lagrang_field_name, field_name, UserDataOperator::OPROWCOL),
    dAta(data),commonFunctions(common_functions),hoGeometry(ho_geometry),Aij(aij) {
      //This will make sure to loop over all entities
      // (e.g. for order=2 it will make doWork to loop 16 time)
      sYmm = false;
    }

Member Function Documentation

doWork()

PetscErrorCode BCs_RVELagrange_Trac::OpRVEBCsLhs::doWork ( int  row_side, int  col_side, EntityType  row_type, EntityType  col_type, DataForcesAndSurcesCore::EntData &  row_data, DataForcesAndSurcesCore::EntData &  col_data  )

inline

Definition at line 185 of file BCs_RVELagrange_Trac.hpp.

      {
      PetscFunctionBegin;
 
      try {
        if(row_data.getIndices().size()==0) PetscFunctionReturn(0);
        if(col_data.getIndices().size()==0) PetscFunctionReturn(0);
 
        int nb_row = row_data.getIndices().size();
        int nb_col = col_data.getIndices().size();
 
        
        const FENumeredDofEntity *dof_ptr;
        ierr = getNumeredEntFiniteElementPtr()->getColDofsByPetscGlobalDofIdx(col_data.getIndices()[0],&dof_ptr); CHKERRQ(ierr);
        int rank = dof_ptr->getNbOfCoeffs();
 
        K.resize(nb_row,nb_col);
        K.clear();
 
        N_mat_row.resize(col_data.getN().size1());
        for(unsigned int gg = 0;gg<col_data.getN().size1();gg++) {
          double area;
          if(hoGeometry) {
            area = norm_2(getNormalsAtGaussPt(gg))*0.5;
          }   else {
            area = getArea();
          }
          double val = getGaussPts()(2,gg)*area;
 
          if(col_type==MBVERTEX) {
            ierr = commonFunctions.shapeMat(rank,gg,col_data,N_mat_row[gg]); CHKERRQ(ierr);
          }
          ierr = commonFunctions.shapeMat(rank,gg,col_data,N_mat_col); CHKERRQ(ierr);
          ierr = commonFunctions.hMat(getNormalsAtGaussPt(gg),rank,N_mat_row(gg),H_mat); CHKERRQ(ierr);
 
          if(gg==0){
            NTN = prod(trans(N_mat_row(gg)),N_mat_col);  //we don't need to define size NTN
            K = val*prod(H_mat,NTN); //we don't need to defiend size of K
          } else {
            noalias(NTN) = prod(trans(N_mat_row(gg)),N_mat_col);
            K += val*prod(H_mat, NTN);
          }
        }
 
        // Matrix C
        int nb_rows=row_data.getIndices().size();
        int nb_cols=col_data.getIndices().size();
        ierr = MatSetValues(
          Aij,nb_rows,&row_data.getIndices()[0],nb_cols,&col_data.getIndices()[0],&K(0,0),ADD_VALUES
        ); CHKERRQ(ierr);
 
        // Matrix CT
        transK = trans(K);
        ierr = MatSetValues(
          Aij,nb_cols,&col_data.getIndices()[0],nb_rows,&row_data.getIndices()[0],&transK(0,0),ADD_VALUES
        ); CHKERRQ(ierr);
 
      } catch (const std::exception& ex) {
        ostringstream ss;
        ss << "throw in method: " << ex.what() << endl;
        SETERRQ(PETSC_COMM_SELF,MOFEM_NOT_IMPLEMENTED,ss.str().c_str());
      }
      PetscFunctionReturn(0);
    }

Member Data Documentation

Aij

Mat BCs_RVELagrange_Trac::OpRVEBCsLhs::Aij

Definition at line 169 of file BCs_RVELagrange_Trac.hpp.

commonFunctions

CommonFunctions& BCs_RVELagrange_Trac::OpRVEBCsLhs::commonFunctions

Definition at line 167 of file BCs_RVELagrange_Trac.hpp.

dAta

RVEBC_Data& BCs_RVELagrange_Trac::OpRVEBCsLhs::dAta

Definition at line 166 of file BCs_RVELagrange_Trac.hpp.

H_mat

MatrixDouble BCs_RVELagrange_Trac::OpRVEBCsLhs::H_mat

Definition at line 183 of file BCs_RVELagrange_Trac.hpp.

hoGeometry

bool BCs_RVELagrange_Trac::OpRVEBCsLhs::hoGeometry

Definition at line 168 of file BCs_RVELagrange_Trac.hpp.

K

MatrixDouble BCs_RVELagrange_Trac::OpRVEBCsLhs::K

Definition at line 181 of file BCs_RVELagrange_Trac.hpp.

N_mat_col

MatrixDouble BCs_RVELagrange_Trac::OpRVEBCsLhs::N_mat_col

Definition at line 183 of file BCs_RVELagrange_Trac.hpp.

N_mat_row

ublas::vector<MatrixDouble > BCs_RVELagrange_Trac::OpRVEBCsLhs::N_mat_row

Definition at line 182 of file BCs_RVELagrange_Trac.hpp.

NTN

MatrixDouble BCs_RVELagrange_Trac::OpRVEBCsLhs::NTN

Definition at line 181 of file BCs_RVELagrange_Trac.hpp.

transK

MatrixDouble BCs_RVELagrange_Trac::OpRVEBCsLhs::transK

Definition at line 181 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