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

Inheritance diagram for BCs_RVELagrange_Periodic::OpRVEBCsPeriodicCalCTLam:

Collaboration diagram for BCs_RVELagrange_Periodic::OpRVEBCsPeriodicCalCTLam:

Public Member Functions
	OpRVEBCsPeriodicCalCTLam (const string field_name, RVEBC_Data_Periodic &data, Vec f, CommonDataPeriodic &common_data_periodic, CommonFunctionsPeriodic &common_functions_periodic, bool ho_geometry=false)

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

Public Attributes
Vec	F

RVEBC_Data_Periodic &	dAta

bool	hoGeometry

CommonDataPeriodic &	commonDataPeriodic

CommonFunctionsPeriodic &	commonFunctionsPeriodic

ublas::vector< ublas::vector< int > >	col_ind

VectorDouble	fe1

VectorDouble	fe2

MatrixDouble	N_mat

Detailed Description

Definition at line 981 of file BCs_RVELagrange_Periodic.hpp.

Constructor & Destructor Documentation

◆ OpRVEBCsPeriodicCalCTLam()

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

inline

Definition at line 989 of file BCs_RVELagrange_Periodic.hpp.

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

Member Function Documentation

◆ doWork()

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

inline

Definition at line 1013 of file BCs_RVELagrange_Periodic.hpp.

                                                                                           {
       PetscFunctionBegin;
 //      cout<<" OpRVEBCsPeriodicCalCTLam "<<endl;
       if(data.getIndices().size()==0) PetscFunctionReturn(0);
       int rank = data.getFieldDofs()[0]->getNbOfCoeffs();
       if(type == MBVERTEX) {
         col_ind.resize(2);
         int nb=data.getIndices().size()/2;  //for one face for nodes only  (total indieces are for both faces)
         col_ind[0].resize(nb); //1st face of prism -- nodes
         col_ind[1].resize(nb); //2nd face of prism -- nodes
         for(int ii=0; ii<nb; ii++){
           col_ind[0][ii]=data.getIndices()[ii];
           col_ind[1][ii]=data.getIndices()[ii+nb];
         }
       }
  
  
 //      cout<<"data.getN().size1() "<<data.getN().size1()<<endl;
       for(unsigned int gg = 0;gg<data.getN().size1();gg++) {
 //         cout<<"gg "<<gg<<endl;
           double area;
           VectorDouble normal_f3 = getNormalsAtGaussPtsF3(gg);
           area = cblas_dnrm2(3,&normal_f3[0],1)*0.5;
           double val = getGaussPts()(2,gg)*area;
  
           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(gg==0){
               if(type == MBVERTEX) {
                 fe1=-val*prod(trans(N_mat), commonDataPeriodic.LagMulAtGaussPts[gg]);
                 fe2= val*prod(trans(N_mat), commonDataPeriodic.LagMulAtGaussPts[gg]);}
               else if(type == MBEDGE && side < 3){
                 fe1=-val*prod(trans(N_mat), commonDataPeriodic.LagMulAtGaussPts[gg]);}
               else if(type == MBEDGE && side  >= 6){
                 fe2= val*prod(trans(N_mat), commonDataPeriodic.LagMulAtGaussPts[gg]);}
               else if(type == MBTRI && side  == 3){
                 fe1=-val*prod(trans(N_mat), commonDataPeriodic.LagMulAtGaussPts[gg]);}
               else if(type == MBTRI && side  == 4){
                 fe2= val*prod(trans(N_mat), commonDataPeriodic.LagMulAtGaussPts[gg]);}
           }else{
           if(type == MBVERTEX) {
             fe1+=-val*prod(trans(N_mat), commonDataPeriodic.LagMulAtGaussPts[gg]);
             fe2+= val*prod(trans(N_mat), commonDataPeriodic.LagMulAtGaussPts[gg]);}
           else if(type == MBEDGE && side < 3){
             fe1+=-val*prod(trans(N_mat), commonDataPeriodic.LagMulAtGaussPts[gg]);}
           else if(type == MBEDGE && side  >= 6){
             fe2+= val*prod(trans(N_mat), commonDataPeriodic.LagMulAtGaussPts[gg]);}
           else if(type == MBTRI && side  == 3){
             fe1+=-val*prod(trans(N_mat), commonDataPeriodic.LagMulAtGaussPts[gg]);}
           else if(type == MBTRI && side  == 4){
             fe2+= val*prod(trans(N_mat), commonDataPeriodic.LagMulAtGaussPts[gg]);}
         }
  
       }
 //      cout<<" Before assembley of F"<<endl;
  
       if(type == MBVERTEX) {
           ierr = VecSetValues(F,col_ind[0].size(),&col_ind[0][0],&fe1[0],ADD_VALUES); CHKERRQ(ierr);
           ierr = VecSetValues(F,col_ind[1].size(),&col_ind[1][0],&fe2[0],ADD_VALUES); CHKERRQ(ierr);}
       else if(type == MBEDGE && side < 3){
           ierr = VecSetValues(F,data.getIndices().size(),&data.getIndices()[0],&fe1[0],ADD_VALUES); CHKERRQ(ierr);}
       else if(type == MBEDGE && side  >= 6){
           ierr = VecSetValues(F,data.getIndices().size(),&data.getIndices()[0],&fe2[0],ADD_VALUES); CHKERRQ(ierr);}
       else if(type == MBTRI && side  == 3){
           ierr = VecSetValues(F,data.getIndices().size(),&data.getIndices()[0],&fe1[0],ADD_VALUES); CHKERRQ(ierr);}
       else if(type == MBTRI && side  == 4){
           ierr = VecSetValues(F,data.getIndices().size(),&data.getIndices()[0],&fe2[0],ADD_VALUES); CHKERRQ(ierr);}
     PetscFunctionReturn(0);
     }