#include "users_modules/mofem_um_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/mofem_um_homogenisation/src/BCs_RVEVolume.hpp

Public Member Functions

Public Attributes