#include <users_modules/homogenisation/src/BCs_RVELagrange_Periodic.hpp>

Inheritance diagram for BCs_RVELagrange_Periodic::OpRVEBCsPeriodicRowInd:

Collaboration diagram for BCs_RVELagrange_Periodic::OpRVEBCsPeriodicRowInd:

Public Member Functions
	OpRVEBCsPeriodicRowInd (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 249 of file BCs_RVELagrange_Periodic.hpp.

Constructor & Destructor Documentation

◆ OpRVEBCsPeriodicRowInd()

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

inline

Definition at line 256 of file BCs_RVELagrange_Periodic.hpp.

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

Member Function Documentation

◆ doWork()

PetscErrorCode BCs_RVELagrange_Periodic::OpRVEBCsPeriodicRowInd::doWork	(	int	side,
		EntityType	type,
		DataForcesAndSourcesCore::EntData &	data
	)

inline

Definition at line 277 of file BCs_RVELagrange_Periodic.hpp.

                                                                                           {
       PetscFunctionBegin;
 //        cout<<"type "<<type << " side "<< side <<endl;
         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)
  
  
         if(type==MBVERTEX){
           commonDataPeriodic.RowInd.resize(2);
           commonDataPeriodic.RowN.resize(2);
         }
  
         int nb;
         //Global indices Rows
         switch(type){
           case MBVERTEX:
             nb=data.getIndices().size()/2;
             commonDataPeriodic.RowInd[0][type][side].resize(nb); //1st face of prism -- nodes
             commonDataPeriodic.RowInd[1][type][side].resize(nb); //2nd face of prism -- nodes
  
             //same indices for rows
             for(int ii=0; ii<nb; ii++){
               commonDataPeriodic.RowInd[0][type][side][ii]=data.getIndices()[ii];
               commonDataPeriodic.RowInd[1][type][side][ii]=data.getIndices()[ii];
             }
 //            cout<<"type "<<type << " side "<< side <<endl;
 //            cout<<"commonDataPeriodic.RowInd[0] "<<"type "<<type << " side "<< side << "     = "<<  commonDataPeriodic.RowInd[0][type][side]<<endl;
             commonDataPeriodic.RowN[0][type][side]=data.getN();
             commonDataPeriodic.RowN[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.RowInd[0][type][side]=data.getIndices();
               commonDataPeriodic.RowInd[1][type][side]=data.getIndices();
               commonDataPeriodic.RowN[0][type][side]=data.getN();
               commonDataPeriodic.RowN[1][type][side]=data.getN();
 //              cout<<"commonDataPeriodic.RowInd[0] "<<"type "<<type << " side "<< side << "     = "<<  commonDataPeriodic.RowInd[0][type][side]<<endl;
 //            cout<<"type "<<type << " side "<< side <<endl;
             }
             break;
           case MBTRI:
             if(side==3){//1st face face 3
               commonDataPeriodic.RowInd[0][type][side]=data.getIndices();
               commonDataPeriodic.RowInd[1][type][side]=data.getIndices();
               commonDataPeriodic.RowN[0][type][side]=data.getN();
               commonDataPeriodic.RowN[1][type][side]=data.getN();
               }
             break;
           default:
             SETERRQ(PETSC_COMM_SELF,1,"data inconsitency");
         }
  
 //        string aaa;
 //        cin>>aaa;
       PetscFunctionReturn(0);
     }