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

#include <users_modules/homogenisation/src/BCs_RVELagrange_Trac_Rigid_Rot.hpp>

Inheritance diagram for BCs_RVELagrange_Trac_Rigid_Rot::OpRVEBCsLhs:

Collaboration diagram for BCs_RVELagrange_Trac_Rigid_Rot::OpRVEBCsLhs:

Public Member Functions
	OpRVEBCsLhs (MoFEM::Interface &_mField, 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, DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::EntData &col_data)

Public Attributes
MoFEM::Interface &	mField

Mat	Aij

RVEBC_Data &	dAta

MatrixDouble	Mat_face

MatrixDouble	Mat_face_Tran

Detailed Description

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

Definition at line 27 of file BCs_RVELagrange_Trac_Rigid_Rot.hpp.

Constructor & Destructor Documentation

◆ OpRVEBCsLhs()

BCs_RVELagrange_Trac_Rigid_Rot::OpRVEBCsLhs::OpRVEBCsLhs	(	MoFEM::Interface &	_mField,
		const string	field_name,
		const string	lagrang_field_name,
		Mat	aij,
		RVEBC_Data &	data
	)

inline

Definition at line 33 of file BCs_RVELagrange_Trac_Rigid_Rot.hpp.

      :
     FaceElementForcesAndSourcesCore::UserDataOperator(
       lagrang_field_name, field_name, UserDataOperator::OPROWCOL),
       mField(_mField),Aij(aij), dAta(data){
       //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_Rigid_Rot::OpRVEBCsLhs::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 47 of file BCs_RVELagrange_Trac_Rigid_Rot.hpp.

       {
       PetscFunctionBegin;
  
       try {
         if(row_data.getIndices().size()==0) PetscFunctionReturn(0);
         if(col_data.getIndices().size()==0) PetscFunctionReturn(0);
         if(col_type==MBVERTEX){
           
  
           Mat_face.resize(3,9);
           Mat_face.clear();
           Mat_face_Tran.resize(9,3);
           Mat_face_Tran.clear();
  
           int num_nodes;
           const EntityHandle* conn_face;
           double coords_face[9];
           EntityHandle face_tri = getNumeredEntFiniteElementPtr()->getEnt(); //handle of finite element
  
           
           rval = mField.get_moab().get_connectivity(face_tri,conn_face,num_nodes,true); CHKERRQ_MOAB(rval);
           rval = mField.get_moab().get_coords(conn_face,num_nodes,coords_face); CHKERRQ_MOAB(rval);
  
           for(int nn=0; nn<3; nn++){
             Mat_face(0,3*nn+1)=-coords_face[3*nn+2];    Mat_face(0,3*nn+2)= coords_face[3*nn+1];
             Mat_face(1,3*nn+0)= coords_face[3*nn+2];    Mat_face(1,3*nn+2)=-coords_face[3*nn+0];
             Mat_face(2,3*nn+0)=-coords_face[3*nn+1];    Mat_face(2,3*nn+1)= coords_face[3*nn+0];
           }
  
           // 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);
     }