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

#include <users_modules/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, DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::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,
		DataForcesAndSourcesCore::EntData &	row_data,
		DataForcesAndSourcesCore::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) {
  
           auto weak_ptr_dof =
               getNumeredEntFiniteElementPtr()->getRowDofsByPetscGlobalDofIdx(
                   col_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();
  
           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);
     }