BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU Struct Reference

\biref operator to calculate and assemble CU More...

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

Inheritance diagram for BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU:

Collaboration diagram for BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU:

Public Member Functions
	OpRVEBCsRhsForceExternal_CU (const string field_name, const string lagrang_field_name, Vec f, CommonFunctions &_common_functions, RVEBC_Data &data, 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
bool	hoGeometry
Vec	F
RVEBC_Data &	dAta
CommonFunctions	commonFunctions
MatrixDouble	K
MatrixDouble	transK
MatrixDouble	N_mat_row
MatrixDouble	N_mat_col
VectorDouble	CU
VectorDouble	CTLam

Detailed Description

\biref operator to calculate and assemble CU

Definition at line 420 of file BCs_RVELagrange_Disp.hpp.

Constructor & Destructor Documentation

OpRVEBCsRhsForceExternal_CU()

BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::OpRVEBCsRhsForceExternal_CU ( const string  field_name, const string  lagrang_field_name, Vec  f, CommonFunctions &  _common_functions, RVEBC_Data &  data, bool  ho_geometry = false  )

inline

Definition at line 428 of file BCs_RVELagrange_Disp.hpp.

     :
    FaceElementForcesAndSourcesCore::UserDataOperator(
      lagrang_field_name, field_name, UserDataOperator::OPROWCOL
    ),
    dAta(data),hoGeometry(ho_geometry),F(f), commonFunctions(_common_functions){
      // 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_Disp::OpRVEBCsRhsForceExternal_CU::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 444 of file BCs_RVELagrange_Disp.hpp.

      {
      PetscFunctionBegin;
 
      try {
        int nb_row = row_data.getIndices().size();
        int nb_col = col_data.getIndices().size();
 
        if(nb_row == 0) PetscFunctionReturn(0);
        if(nb_col == 0) PetscFunctionReturn(0);
 
 
        
        const FENumeredDofEntity *dof_ptr;
        ierr = getNumeredEntFiniteElementPtr()->getRowDofsByPetscGlobalDofIdx(row_data.getIndices()[0],&dof_ptr); CHKERRQ(ierr);
        int rank = dof_ptr->getNbOfCoeffs();
        for(unsigned int gg = 0;gg<row_data.getN().size1();gg++) {
          double area;
          if(hoGeometry) {
            area = norm_2(getNormalsAtGaussPt(gg))*0.5;
          }   else {
            area = getArea();
          }
          double val = getGaussPts()(2,gg)*area;
          ierr = commonFunctions.shapeMat(rank, gg,  row_data, N_mat_row); CHKERRQ(ierr);
          ierr = commonFunctions.shapeMat(rank, gg,  col_data, N_mat_col); CHKERRQ(ierr);
 
//        cout<<"row_data.getN() "<<row_data.getN()<<endl;
//        cout<<"N_mat_row() "<<N_mat_row<<endl;
//
//        cout<<"col_data.getN() "<<col_data.getN()<<endl;
//        cout<<"N_mat_col "<<N_mat_col<<endl;
 
          if(gg==0){
            K =val* prod(trans(N_mat_row),N_mat_col);  //we don't need to define size K
          } else {
            K += val*prod(trans(N_mat_row),N_mat_col);
          }
        }
        CU=prod(K,       col_data.getFieldData());
        CTLam=prod(trans(K),row_data.getFieldData());
 
 
        Vec myF = F;
        if(F == PETSC_NULL) {
          myF = getFEMethod()->snes_f;
        }
        
        ierr = VecSetValues(myF,row_data.getIndices().size(),&row_data.getIndices()[0],&CU[0],ADD_VALUES); CHKERRQ(ierr);
        ierr = VecSetValues(myF,col_data.getIndices().size(),&col_data.getIndices()[0],&CTLam[0],ADD_VALUES); CHKERRQ(ierr);
 
      } catch (const std::exception& ex) {
        ostringstream ss;
        ss << "throw in method: " << ex.what() << endl;
        SETERRQ(PETSC_COMM_SELF,1,ss.str().c_str());
      }
 
      PetscFunctionReturn(0);
    }

Member Data Documentation

commonFunctions

CommonFunctions BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::commonFunctions

Definition at line 425 of file BCs_RVELagrange_Disp.hpp.

CTLam

VectorDouble BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::CTLam

Definition at line 442 of file BCs_RVELagrange_Disp.hpp.

CU

VectorDouble BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::CU

Definition at line 442 of file BCs_RVELagrange_Disp.hpp.

dAta

RVEBC_Data& BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::dAta

Definition at line 424 of file BCs_RVELagrange_Disp.hpp.

F

Vec BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::F

Definition at line 423 of file BCs_RVELagrange_Disp.hpp.

hoGeometry

bool BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::hoGeometry

Definition at line 422 of file BCs_RVELagrange_Disp.hpp.

K

MatrixDouble BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::K

Definition at line 441 of file BCs_RVELagrange_Disp.hpp.

N_mat_col

MatrixDouble BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::N_mat_col

Definition at line 441 of file BCs_RVELagrange_Disp.hpp.

N_mat_row

MatrixDouble BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::N_mat_row

Definition at line 441 of file BCs_RVELagrange_Disp.hpp.

transK

MatrixDouble BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::transK

Definition at line 441 of file BCs_RVELagrange_Disp.hpp.

---

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

• users_modules/mofem_um_homogenisation/src/BCs_RVELagrange_Disp.hpp