#include <users_modules/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, DataForcesAndSourcesCore::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,
		DataForcesAndSourcesCore::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);
     }