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

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

Inheritance diagram for BCs_RVELagrange_Trac_Rigid_Trans::OpRVEBCsLhs:

Collaboration diagram for BCs_RVELagrange_Trac_Rigid_Trans::OpRVEBCsLhs:

Public Member Functions
	OpRVEBCsLhs (const string field_name, const string lagrang_field_name, Mat aij, RVEBC_Data &data)

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

Mat	Aij

MatrixDouble	Mat_face

MatrixDouble	Mat_face_Tran

Detailed Description

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

Definition at line 24 of file BCs_RVELagrange_Trac_Rigid_Trans.hpp.

Constructor & Destructor Documentation

◆ OpRVEBCsLhs()

BCs_RVELagrange_Trac_Rigid_Trans::OpRVEBCsLhs::OpRVEBCsLhs	(	const string	field_name,
		const string	lagrang_field_name,
		Mat	aij,
		RVEBC_Data &	data
	)

inline

Definition at line 28 of file BCs_RVELagrange_Trac_Rigid_Trans.hpp.

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

Member Function Documentation

◆ doWork()

PetscErrorCode BCs_RVELagrange_Trac_Rigid_Trans::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 40 of file BCs_RVELagrange_Trac_Rigid_Trans.hpp.

      {
      PetscFunctionBegin;
 
      try {
        if(row_data.getIndices().size()==0) PetscFunctionReturn(0);
        if(col_data.getIndices().size()==0) PetscFunctionReturn(0);
        if(col_type==MBVERTEX) {
 
          
          const FENumeredDofEntity *dof_ptr;
          ierr = getNumeredEntFiniteElementPtr()->getColDofsByPetscGlobalDofIdx(col_data.getIndices()[0],&dof_ptr); CHKERRQ(ierr);
          int rank = dof_ptr->getNbOfCoeffs();
 
          Mat_face.resize(rank,3*rank);
          Mat_face.clear();
          Mat_face_Tran.resize(3*rank,rank);
          Mat_face_Tran.clear();
 
          switch(rank) {
            case 3:
            for(int nn=0; nn<3; nn++){
              Mat_face(0,3*nn+0) = 1.0;
              Mat_face(1,3*nn+1) = 1.0;
              Mat_face(2,3*nn+2) = 1.0;
            }
            break;
            case 1:
            Mat_face(0,0) = 1.0;
            Mat_face(0,1) = 1.0;
            Mat_face(0,2) = 1.0;
            break;
            default:
            SETERRQ(PETSC_COMM_SELF,1,"not implemented");
          }
 
          // Matrix C1
          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],&Mat_face(0,0),ADD_VALUES
          ); CHKERRQ(ierr);
 
          // Matrix C1T
          noalias(Mat_face_Tran) = trans(Mat_face);
          ierr = MatSetValues(
            Aij,nb_cols,&col_data.getIndices()[0],nb_rows,&row_data.getIndices()[0],&Mat_face_Tran(0,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);
    }