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

Inheritance diagram for BCs_RVELagrange_Periodic::OpRVEBCsPeriodicCalDispAtGaussPts:

Collaboration diagram for BCs_RVELagrange_Periodic::OpRVEBCsPeriodicCalDispAtGaussPts:

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

Constructor & Destructor Documentation

◆ OpRVEBCsPeriodicCalDispAtGaussPts()

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

inline

Definition at line 747 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::OpRVEBCsPeriodicCalDispAtGaussPts::doWork	(	int	side,
		EntityType	type,
		DataForcesAndSourcesCore::EntData &	data
	)

inline

Definition at line 768 of file BCs_RVELagrange_Periodic.hpp.

                                                                                           {
       PetscFunctionBegin;
 //      cout<<" OpRVEBCsPeriodicCalDispAtGaussPts "<<endl;
       if(data.getIndices().size()==0) PetscFunctionReturn(0);
       int nb_gauss_pts = data.getN().size1();
       int rank = data.getFieldDofs()[0]->getNbOfCoeffs();
  
        //initialize
       if(type == MBVERTEX) {
         commonDataPeriodic.DispAtGaussPts.resize(2);
         field_data_nodes.resize(2);
  
         commonDataPeriodic.DispAtGaussPts[0].resize(nb_gauss_pts);
         commonDataPeriodic.DispAtGaussPts[1].resize(nb_gauss_pts);
         for(int gg = 0;gg<nb_gauss_pts;gg++) {
           commonDataPeriodic.DispAtGaussPts[0][gg].resize(rank);  commonDataPeriodic.DispAtGaussPts[0][gg].clear();
           commonDataPeriodic.DispAtGaussPts[1][gg].resize(rank);  commonDataPeriodic.DispAtGaussPts[1][gg].clear();
         }
  
         int nb=data.getFieldData().size()/2;  //for one face for nodes only  (total indieces are for both faces)
         field_data_nodes[0].resize(nb); //1st face of prism -- nodes
         field_data_nodes[1].resize(nb); //2nd face of prism -- nodes
  
         for(int ii=0; ii<nb; ii++){
           field_data_nodes[0][ii]=data.getFieldData()[ii];
           field_data_nodes[1][ii]=data.getFieldData()[ii+nb];
         }
 //        cout<<"field_data_nodes[0] "<<field_data_nodes[0]<<endl;
 //        cout<<"field_data_nodes[1] "<<field_data_nodes[1]<<endl;
       }
  
       for(unsigned int gg = 0; gg<nb_gauss_pts;gg++) {
           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(type == MBVERTEX) {
             commonDataPeriodic.DispAtGaussPts[0][gg]=prod(N_mat, field_data_nodes[0]);
             commonDataPeriodic.DispAtGaussPts[1][gg]=prod(N_mat, field_data_nodes[1]);}
           else if(type == MBEDGE && side < 3){
             commonDataPeriodic.DispAtGaussPts[0][gg]+=prod(N_mat, data.getFieldData());}
           else if(type == MBEDGE && side  >= 6){
             commonDataPeriodic.DispAtGaussPts[1][gg]+=prod(N_mat, data.getFieldData());}
           else if(type == MBTRI && side  == 3){
             commonDataPeriodic.DispAtGaussPts[0][gg]+=prod(N_mat, data.getFieldData());}
           else if(type == MBTRI && side  == 4){
             commonDataPeriodic.DispAtGaussPts[1][gg]+=prod(N_mat, data.getFieldData());}
       }
  
     PetscFunctionReturn(0);
     }