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