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

Inheritance diagram for BCs_RVELagrange_Periodic::OpRVEBCsPeriodicRhs:

Collaboration diagram for BCs_RVELagrange_Periodic::OpRVEBCsPeriodicRhs:

Public Member Functions
	OpRVEBCsPeriodicRhs (const string field_name, const string lagrang_field_name, vector< Vec > &f, RVEBC_Data_Periodic &data, CommonDataPeriodic &common_data_periodic, CommonFunctionsPeriodic &common_functions_periodic, bool ho_geometry=false)

PetscErrorCode	doWork (int side, EntityType type, DataForcesAndSurcesCore::EntData &data)

Public Member Functions inherited from BCs_RVELagrange_Periodic::OpDmatRhs
	OpDmatRhs (const string field_name, const string lagrang_field_name, RVEBC_Data_Periodic &data, CommonDataPeriodic &common_data_periodic, CommonFunctionsPeriodic &common_functions_periodic, bool ho_geometry=false)

PetscErrorCode	calculateDmat (int side, EntityType type, DataForcesAndSurcesCore::EntData &data)

Public Attributes
vector< Vec > &	F

VectorDouble	appliedStrain

VectorDouble	nF

Public Attributes inherited from BCs_RVELagrange_Periodic::OpDmatRhs
RVEBC_Data_Periodic &	dAta

bool	hoGeometry

CommonDataPeriodic &	commonDataPeriodic

CommonFunctionsPeriodic &	commonFunctionsPeriodic

ublas::vector< MatrixDouble >	X_mat

MatrixDouble	N_mat

Detailed Description

Definition at line 597 of file BCs_RVELagrange_Periodic.hpp.

Constructor & Destructor Documentation

◆ OpRVEBCsPeriodicRhs()

BCs_RVELagrange_Periodic::OpRVEBCsPeriodicRhs::OpRVEBCsPeriodicRhs	(	const string	field_name,
		const string	lagrang_field_name,
		vector< Vec > &	f,
		RVEBC_Data_Periodic &	data,
		CommonDataPeriodic &	common_data_periodic,
		CommonFunctionsPeriodic &	common_functions_periodic,
		bool	ho_geometry = `false`
	)

inline

Definition at line 601 of file BCs_RVELagrange_Periodic.hpp.

     :
    OpDmatRhs(
      field_name,lagrang_field_name,data,common_data_periodic,common_functions_periodic,ho_geometry
    ),
    F(f) {
    }

Member Function Documentation

◆ doWork()

PetscErrorCode BCs_RVELagrange_Periodic::OpRVEBCsPeriodicRhs::doWork	(	int	side,
		EntityType	type,
		DataForcesAndSurcesCore::EntData &	data
	)

inline

Definition at line 619 of file BCs_RVELagrange_Periodic.hpp.

                                                                                         {
      PetscFunctionBegin;
      if(data.getIndices().size()==0) PetscFunctionReturn(0);
      if(type == MBEDGE && side >= 3) PetscFunctionReturn(0); //ignore second triangel edges (cananical number 6 to 8)
      if(type == MBTRI && side == 4) PetscFunctionReturn(0);  //ignore second triangel face (cananical number 4)
 
      
      ierr = calculateDmat(side,type,data); CHKERRQ(ierr);
 
//      if(type == MBVERTEX) {
//      cout<<"\n\nNew element "<<endl;
//      }
//      cout<<"type "<<type << " side "<< side <<endl;
//      cout<< "commonDataPeriodic.D_mat[0] "<<commonDataPeriodic.D_mat[0]<<endl;
//      cout<< "commonDataPeriodic.D_mat[1] "<<commonDataPeriodic.D_mat[1]<<endl<<endl;
//      string aaa;
//      cin >> aaa;
 
 
      ublas::vector<int> rowvec;
      if(type == MBVERTEX) {
        int nb=data.getIndices().size()/2;  //for one face
        rowvec.resize(nb);
        for(int ii=0; ii<nb; ii++){
          rowvec[ii]=data.getIndices()[ii];
        }
      } else{
        rowvec=data.getIndices();
      }
 
      appliedStrain.resize(commonDataPeriodic.D_mat[0].size2(),false);
 
      int rank = data.getFieldDofs()[0]->getNbOfCoeffs();
      int size = (rank == 3) ? 6 : 3;
      for(int ii = 0;ii<size;ii++) {
        appliedStrain.clear();
        appliedStrain[ii] = 1.0;
        nF = prod(commonDataPeriodic.D_mat[0], appliedStrain);
        ierr = VecSetValues(F[ii],rowvec.size(),&rowvec[0],&nF[0],ADD_VALUES); CHKERRQ(ierr);
        nF = prod(commonDataPeriodic.D_mat[1], appliedStrain);
        ierr = VecSetValues(F[ii],rowvec.size(),&rowvec[0],&nF[0],ADD_VALUES); CHKERRQ(ierr);
      }
 
      PetscFunctionReturn(0);
    }