BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU Struct Reference

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

#include <users_modules/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, DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::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 438 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 446 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, DataForcesAndSourcesCore::EntData &  row_data, DataForcesAndSourcesCore::EntData &  col_data  )

inline

Definition at line 462 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);
 
        auto weak_ptr_dof =
            getNumeredEntFiniteElementPtr()->getRowDofsByPetscGlobalDofIdx(
                row_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");
 
        int rank = dof_ptr->getNbOfCoeffs();
 
        for(unsigned int gg = 0;gg<row_data.getN().size1();gg++) {
          double area;
          if(hoGeometry) {
            area = norm_2(getNormalsAtGaussPts(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 443 of file BCs_RVELagrange_Disp.hpp.

CTLam

VectorDouble BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::CTLam

Definition at line 460 of file BCs_RVELagrange_Disp.hpp.

CU

VectorDouble BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::CU

Definition at line 460 of file BCs_RVELagrange_Disp.hpp.

dAta

RVEBC_Data& BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::dAta

Definition at line 442 of file BCs_RVELagrange_Disp.hpp.

F

Vec BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::F

Definition at line 441 of file BCs_RVELagrange_Disp.hpp.

hoGeometry

bool BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::hoGeometry

Definition at line 440 of file BCs_RVELagrange_Disp.hpp.

K

MatrixDouble BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::K

Definition at line 459 of file BCs_RVELagrange_Disp.hpp.

N_mat_col

MatrixDouble BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::N_mat_col

Definition at line 459 of file BCs_RVELagrange_Disp.hpp.

N_mat_row

MatrixDouble BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::N_mat_row

Definition at line 459 of file BCs_RVELagrange_Disp.hpp.

transK

MatrixDouble BCs_RVELagrange_Disp::OpRVEBCsRhsForceExternal_CU::transK

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