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

Inheritance diagram for BCs_RVEVolume::OpVolumeCal:

Collaboration diagram for BCs_RVEVolume::OpVolumeCal:

Public Member Functions
	OpVolumeCal (MoFEM::Interface &m_field, string field_name, Vec _RVE_volume_Vec, NonlinearElasticElement::BlockData &data)

PetscErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, DataForcesAndSurcesCore::EntData &row_data, DataForcesAndSurcesCore::EntData &col_data)

Public Attributes
Vec	RVE_volume_Vec

NonlinearElasticElement::BlockData &	dAta

MoFEM::Interface &	mField

Detailed Description

\biref operator to calculate left hand side of heat conductivity terms \infroup mofem_thermal_elem

Definition at line 41 of file BCs_RVEVolume.hpp.

Constructor & Destructor Documentation

◆ OpVolumeCal()

BCs_RVEVolume::OpVolumeCal::OpVolumeCal	(	MoFEM::Interface &	m_field,
		string	field_name,
		Vec	_RVE_volume_Vec,
		NonlinearElasticElement::BlockData &	data
	)

inline

Definition at line 46 of file BCs_RVEVolume.hpp.

                                                                                                                         :
     VolumeElementForcesAndSourcesCore::UserDataOperator(field_name, UserDataOperator::OPROWCOL),
     mField(m_field), RVE_volume_Vec(_RVE_volume_Vec), dAta(data) { }

Member Function Documentation

◆ doWork()

PetscErrorCode BCs_RVEVolume::OpVolumeCal::doWork	(	int	row_side,
		int	col_side,
		EntityType	row_type,
		EntityType	col_type,
		DataForcesAndSurcesCore::EntData &	row_data,
		DataForcesAndSurcesCore::EntData &	col_data
	)

inline

Definition at line 50 of file BCs_RVEVolume.hpp.

                                                 {
         PetscFunctionBegin;
         if(dAta.tEts.find(getNumeredEntFiniteElementPtr()->getEnt()) == dAta.tEts.end()) {
           PetscFunctionReturn(0);
         }
  
         //      cout<<"Hi from OpVolumeCal "<<endl;
         if(row_type == MBVERTEX && col_type==MBVERTEX) {
           //        cout<<"Hi from MBVERTEX "<<endl;
           ParallelComm* pcomm = ParallelComm::get_pcomm(&mField.get_moab(),MYPCOMM_INDEX);
           int Indices[1];  Indices[0]=pcomm->rank();
           double Vol_elm[1];  Vol_elm[0]=0;
           for(unsigned int gg = 0;gg<row_data.getN().size1();gg++) {
             if(getHoGaussPtsDetJac().size()>0) {
               //              cout<<"getHoGaussPtsDetJac()[gg] "<<getHoGaussPtsDetJac()[gg]<<endl;
               //              cout<<"High order geometry "<<endl;
               Vol_elm[0]+=getVolume()*getGaussPts()(3,gg)*getHoGaussPtsDetJac()[gg];
             }else{
               //              cout<<"Low order geometry "<<endl;
               Vol_elm[0]+=getVolume()*getGaussPts()(3,gg);
             }
           }
           
           ierr = VecSetValues(RVE_volume_Vec,1,Indices,Vol_elm,ADD_VALUES); CHKERRQ(ierr);
           //          cout<<"Indices[0] "<<Indices[0] << endl;
           //          cout<<"Vol_elm[0] "<<Vol_elm[0] << endl;
         }
         PetscFunctionReturn(0);
       }

Member Data Documentation

◆ dAta

NonlinearElasticElement::BlockData& BCs_RVEVolume::OpVolumeCal::dAta

Definition at line 43 of file BCs_RVEVolume.hpp.

◆ mField

MoFEM::Interface& BCs_RVEVolume::OpVolumeCal::mField

Definition at line 44 of file BCs_RVEVolume.hpp.

◆ RVE_volume_Vec

Vec BCs_RVEVolume::OpVolumeCal::RVE_volume_Vec

Definition at line 42 of file BCs_RVEVolume.hpp.

The documentation for this struct was generated from the following file:

users_modules/homogenisation/src/BCs_RVEVolume.hpp

Public Member Functions

Public Attributes