\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);
    }