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

Inheritance diagram for BCs_RVELagrange_Periodic::OpRVEBCsPeriodicCalCU:

Collaboration diagram for BCs_RVELagrange_Periodic::OpRVEBCsPeriodicCalCU:

Public Member Functions
	OpRVEBCsPeriodicCalCU (const string field_name, RVEBC_Data_Periodic &data, Vec f, CommonDataPeriodic &common_data_periodic, CommonFunctionsPeriodic &common_functions_periodic, bool ho_geometry=false)

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

Public Attributes
Vec	F

RVEBC_Data_Periodic &	dAta

bool	hoGeometry

CommonDataPeriodic &	commonDataPeriodic

CommonFunctionsPeriodic &	commonFunctionsPeriodic

ublas::vector< int >	row_ind

VectorDouble	fe1

VectorDouble	fe2

MatrixDouble	N_mat

Detailed Description

Definition at line 825 of file BCs_RVELagrange_Periodic.hpp.

Constructor & Destructor Documentation

◆ OpRVEBCsPeriodicCalCU()

BCs_RVELagrange_Periodic::OpRVEBCsPeriodicCalCU::OpRVEBCsPeriodicCalCU	(	const string	field_name,
		RVEBC_Data_Periodic &	data,
		Vec	f,
		CommonDataPeriodic &	common_data_periodic,
		CommonFunctionsPeriodic &	common_functions_periodic,
		bool	ho_geometry = `false`
	)

inline

Definition at line 833 of file BCs_RVELagrange_Periodic.hpp.

     :
    FlatPrismElementForcesAndSurcesCore::UserDataOperator(
      field_name,UserDataOperator::OPROW
    ),
    F(f),
    dAta(data),
    hoGeometry(ho_geometry),
    commonDataPeriodic(common_data_periodic),
    commonFunctionsPeriodic(common_functions_periodic) {
    }

Member Function Documentation

◆ doWork()

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

inline

Definition at line 857 of file BCs_RVELagrange_Periodic.hpp.

                                                                                         {
      PetscFunctionBegin;
//      cout<<" OpRVEBCsPeriodicCalCU "<<endl;
      if(data.getIndices().size()==0) PetscFunctionReturn(0);
      if(type == MBQUAD) 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)
 
      int rank = data.getFieldDofs()[0]->getNbOfCoeffs();
 
      if(type == MBVERTEX) {
        int nb=data.getIndices().size()/2;  //for one face for nodes only  (total indieces are for both faces)
        row_ind.resize(nb); //1st face of prism -- nodes
        for(int ii=0; ii<nb; ii++){
          row_ind[ii]=data.getIndices()[ii];
        }
      }
 
      for(unsigned int gg = 0;gg<data.getN().size1();gg++) {
          double area;
          VectorDouble normal_f3 = getNormalsAtGaussPtF3(gg);
          area = cblas_dnrm2(3,&normal_f3[0],1)*0.5;
          double val = getGaussPts()(2,gg)*area;
 
          if(type == MBVERTEX) {
            ierr = commonFunctionsPeriodic.shapeMat(rank, gg, data, N_mat, 2); CHKERRQ(ierr);}
          else {
            ierr = commonFunctionsPeriodic.shapeMat(rank, gg, data, N_mat); CHKERRQ(ierr);}
 
          if(gg==0){
            fe1=-val*prod(trans(N_mat), commonDataPeriodic.DispAtGaussPts[0][gg]);
            fe2= val*prod(trans(N_mat), commonDataPeriodic.DispAtGaussPts[1][gg]);}
          else{
            fe1+=-val*prod(trans(N_mat), commonDataPeriodic.DispAtGaussPts[0][gg]);
            fe2+= val*prod(trans(N_mat), commonDataPeriodic.DispAtGaussPts[1][gg]);}
      }
 
      if(type == MBVERTEX) {
          ierr = VecSetValues(F,row_ind.size(),&row_ind[0],&fe1[0],ADD_VALUES); CHKERRQ(ierr);
          ierr = VecSetValues(F,row_ind.size(),&row_ind[0],&fe2[0],ADD_VALUES); CHKERRQ(ierr);}
      else{
          ierr = VecSetValues(F,data.getIndices().size(),&data.getIndices()[0],&fe1[0],ADD_VALUES); CHKERRQ(ierr);
          ierr = VecSetValues(F,data.getIndices().size(),&data.getIndices()[0],&fe2[0],ADD_VALUES); CHKERRQ(ierr);}
 
    PetscFunctionReturn(0);
    }