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

Inheritance diagram for BCs_RVELagrange_Periodic::OpRVEBCsPeriodicColInd:

Collaboration diagram for BCs_RVELagrange_Periodic::OpRVEBCsPeriodicColInd:

Public Member Functions
	OpRVEBCsPeriodicColInd (const string field_name, 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 Attributes
RVEBC_Data_Periodic &	dAta

bool	hoGeometry

CommonDataPeriodic &	commonDataPeriodic

CommonFunctionsPeriodic &	commonFunctionsPeriodic

ublas::vector< VectorDouble >	field_data_nodes

MatrixDouble	N_mat

Detailed Description

Definition at line 161 of file BCs_RVELagrange_Periodic.hpp.

Constructor & Destructor Documentation

◆ OpRVEBCsPeriodicColInd()

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

inline

Definition at line 168 of file BCs_RVELagrange_Periodic.hpp.

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

Member Function Documentation

◆ doWork()

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

inline

Definition at line 189 of file BCs_RVELagrange_Periodic.hpp.

                                                                                         {
      PetscFunctionBegin;
 
        if(type==MBVERTEX){
          commonDataPeriodic.ColInd.resize(2);
          commonDataPeriodic.ColN.resize(2);
        }
        int nb;
        //Global indices Columns
        if(data.getIndices().size()==0) PetscFunctionReturn(0);
 
        switch(type){
          case MBVERTEX:
            nb=data.getIndices().size()/2;  //for one face for nodes only  (total indieces are for both faces)
            commonDataPeriodic.ColInd[0][type][side].resize(nb); //1st face of prism -- nodes
            commonDataPeriodic.ColInd[1][type][side].resize(nb); //2nd face of prism -- nodes
            //different indices for cols
            for(int ii=0; ii<nb; ii++){
              commonDataPeriodic.ColInd[0][type][side][ii]=data.getIndices()[ii];
              commonDataPeriodic.ColInd[1][type][side][ii]=data.getIndices()[ii+nb];
            }
//              cout<<"commonDataPeriodic.ColInd[0] "<<"type "<<type << " side "<< side << "     = "<<  commonDataPeriodic.ColInd[0][type][side]<<endl;
//              cout<<"commonDataPeriodic.ColInd[0] "<<"type "<<type << " side "<< side << "     = "<<  commonDataPeriodic.ColInd[1][type][side]<<endl;
//              cout<<"type "<<type << " side "<< side <<endl;
            commonDataPeriodic.ColN[0][type][side]=data.getN();
            commonDataPeriodic.ColN[1][type][side]=data.getN();
            break;
          case MBEDGE:
            if(side<3){ //1st face eges 0, 1, 2 (cananical numbering)-- same indices for both faces
              commonDataPeriodic.ColInd[0][type][side]=data.getIndices();
              commonDataPeriodic.ColN[0][type][side]=data.getN();
//              cout<<"col_type "<<col_type << " col_side "<< col_side <<endl;
//              cout<<"col_data.getIndices "<<col_data.getIndices()<<endl;
              }
            else if (side >= 6){ //2nd face edges 6,7,8
              commonDataPeriodic.ColInd[1][type][side-6]=data.getIndices();
              commonDataPeriodic.ColN[1][type][side-6]=data.getN();
//              cout<<"col_type "<<col_type << " col_side "<< col_side <<endl;
              }
            break;
          case MBTRI:
            if(side==3){//1st face face 3
              commonDataPeriodic.ColInd[0][type][side]=data.getIndices();
              commonDataPeriodic.ColN[0][type][side]=data.getN();
              }
            else{ //2nd face face 4
              commonDataPeriodic.ColInd[1][type][side-1]=data.getIndices();
              commonDataPeriodic.ColN[1][type][side-1]=data.getN();
              }
            break;
          default:
            SETERRQ(PETSC_COMM_SELF,1,"data inconsitency");
        }
      PetscFunctionReturn(0);
    }