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

Inheritance diagram for BCs_RVELagrange_Periodic::OpRVEBCsPeriodicRhs_givenStrain:

Collaboration diagram for BCs_RVELagrange_Periodic::OpRVEBCsPeriodicRhs_givenStrain:

Public Member Functions
	OpRVEBCsPeriodicRhs_givenStrain (const string field_name, const string lagrang_field_name, Vec f, VectorDouble given_strain, boost::ptr_vector< MethodForForceScaling > &methods_op, 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 Member Functions inherited from BCs_RVELagrange_Periodic::OpDmatRhs
	OpDmatRhs (const string field_name, const string lagrang_field_name, RVEBC_Data_Periodic &data, CommonDataPeriodic &common_data_periodic, CommonFunctionsPeriodic &common_functions_periodic, bool ho_geometry=false)

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

Public Attributes
Vec	F

VectorDouble	givenStrain

boost::ptr_vector< MethodForForceScaling > &	methodsOp

VectorDouble	nF

Public Attributes inherited from BCs_RVELagrange_Periodic::OpDmatRhs
RVEBC_Data_Periodic &	dAta

bool	hoGeometry

CommonDataPeriodic &	commonDataPeriodic

CommonFunctionsPeriodic &	commonFunctionsPeriodic

ublas::vector< MatrixDouble >	X_mat

MatrixDouble	N_mat

Detailed Description

Definition at line 669 of file BCs_RVELagrange_Periodic.hpp.

Constructor & Destructor Documentation

◆ OpRVEBCsPeriodicRhs_givenStrain()

BCs_RVELagrange_Periodic::OpRVEBCsPeriodicRhs_givenStrain::OpRVEBCsPeriodicRhs_givenStrain	(	const string	field_name,
		const string	lagrang_field_name,
		Vec	f,
		VectorDouble	given_strain,
		boost::ptr_vector< MethodForForceScaling > &	methods_op,
		RVEBC_Data_Periodic &	data,
		CommonDataPeriodic &	common_data_periodic,
		CommonFunctionsPeriodic &	common_functions_periodic,
		bool	ho_geometry = `false`
	)

inline

Definition at line 675 of file BCs_RVELagrange_Periodic.hpp.

      :
     OpDmatRhs(
       field_name,lagrang_field_name,data,common_data_periodic,common_functions_periodic,ho_geometry
     ),
     F(f), givenStrain(given_strain),methodsOp(methods_op){
     }

Member Function Documentation

◆ doWork()

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

inline

Definition at line 694 of file BCs_RVELagrange_Periodic.hpp.

                                                                                           {
       PetscFunctionBegin;
       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)
  
       
       ierr = calculateDmat(side,type,data); CHKERRQ(ierr);
  
       ublas::vector<int> rowvec;
       if(type == MBVERTEX) {
         int nb=data.getIndices().size()/2;  //for one face
         rowvec.resize(nb);
         for(int ii=0; ii<nb; ii++){
           rowvec[ii]=data.getIndices()[ii];
         }
       } else{
         rowvec=data.getIndices();
       }
  
       //This is to scale givenStrain vector to apply the strian in steps
       VectorDouble scaled_given_strain;
       scaled_given_strain.resize(6);
       scaled_given_strain=givenStrain;
 //      cout <<"givenStrain Before = "<<scaled_given_strain<<endl;
       ierr = MethodForForceScaling::applyScale(getFEMethod(),methodsOp,scaled_given_strain); CHKERRQ(ierr);
 //      cout <<"givenStrain  After = "<<scaled_given_strain<<endl;
  
 //      SNES snes = getFEMethod()->snes;
 //      int iter;
 //      SNESGetIterationNumber(snes,&iter);
 //      if(iter>0)  PetscFunctionReturn(0);
  
       nF = prod(commonDataPeriodic.D_mat[0], scaled_given_strain);
       nF*=-1;
       ierr = VecSetValues(F,rowvec.size(),&rowvec[0],&nF[0],ADD_VALUES); CHKERRQ(ierr);
  
       nF = prod(commonDataPeriodic.D_mat[1], scaled_given_strain);
       nF*=-1;
       ierr = VecSetValues(F,rowvec.size(),&rowvec[0],&nF[0],ADD_VALUES); CHKERRQ(ierr);
  
       PetscFunctionReturn(0);
     }