PeriodicNitscheConstrains::PeriodicVolume Struct Reference

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

Inheritance diagram for PeriodicNitscheConstrains::PeriodicVolume:

Collaboration diagram for PeriodicNitscheConstrains::PeriodicVolume:

Public Member Functions
	PeriodicVolume (MoFEM::Interface &m_field, CommonData &common_data)
int	getRule (int order)
PetscErrorCode	setGaussPts (int order)

Public Attributes
EntityHandle	forFace
CommonData &	commonData

Detailed Description

Definition at line 235 of file NitschePeriodicMethod.hpp.

Constructor & Destructor Documentation

PeriodicVolume()

PeriodicNitscheConstrains::PeriodicVolume::PeriodicVolume ( MoFEM::Interface &  m_field, CommonData &  common_data  )

inline

Definition at line 240 of file NitschePeriodicMethod.hpp.

                                                                   :
    VolumeElementForcesAndSourcesCore(m_field),
    forFace(0),
    commonData(common_data) {}

Member Function Documentation

getRule()

int PeriodicNitscheConstrains::PeriodicVolume::getRule ( int  order )

inline

Definition at line 244 of file NitschePeriodicMethod.hpp.

{ return -1; }

setGaussPts()

PetscErrorCode PeriodicNitscheConstrains::PeriodicVolume::setGaussPts ( int  order )

inline

Definition at line 245 of file NitschePeriodicMethod.hpp.

                                          {
 
      PetscFunctionBegin;
      try {
        gaussPts.resize(4,0,false);
        EntityHandle tet = numeredEntFiniteElementPtr->getEnt();
        int vffgg = 0;
        for(int ff = 0;ff<4;ff++) {
          EntityHandle face;
          rval = mField.get_moab().side_element(tet,2,ff,face); CHKERRQ_MOAB(rval);
          if(face!=forFace) continue;
          if(commonData.localCoordsMap.find(face)!=commonData.localCoordsMap.end()) {
            const double coords_tet[12] = { 0,0,0, 1,0,0, 0,1,0, 0,0,1 };
            int nb_face_gauss_pts = commonData.localCoordsMap[face].size();
            gaussPts.resize(4,nb_face_gauss_pts);
            commonData.inTetTetGaussPtsNumber[tet].resize(nb_face_gauss_pts);
            for(int ffgg = 0;ffgg!=nb_face_gauss_pts;ffgg++,vffgg++) {
              commonData.inTetTetGaussPtsNumber[tet][vffgg] = commonData.inTetFaceGaussPtsNumber[face][ffgg];
              double N[3];
              VectorDouble &local_coords = commonData.localCoordsMap[face][ffgg];
              ierr = ShapeMBTRI(N,&local_coords[0],&local_coords[1],1); CHKERRQ(ierr);
              for(int dd = 0;dd<3;dd++) {
                gaussPts(dd,vffgg) =
                N[0]*coords_tet[3*dataH1.facesNodes(ff,0)+dd]+
                N[1]*coords_tet[3*dataH1.facesNodes(ff,1)+dd]+
                N[2]*coords_tet[3*dataH1.facesNodes(ff,2)+dd];
              }
              gaussPts(3,vffgg) = 0;
            }
          } else {
            SETERRQ(PETSC_COMM_SELF,MOFEM_DATA_INCONSISTENCY,"data inconstancy");
          }
        }
        //cerr << gaussPts << endl;
      } catch (const std::exception& ex) {
        ostringstream ss;
        ss << "throw in method: " << ex.what() << endl;
        SETERRQ(PETSC_COMM_SELF,1,ss.str().c_str());
      }
      PetscFunctionReturn(0);
    }

Member Data Documentation

commonData

CommonData& PeriodicNitscheConstrains::PeriodicVolume::commonData

Definition at line 238 of file NitschePeriodicMethod.hpp.

forFace

EntityHandle PeriodicNitscheConstrains::PeriodicVolume::forFace

Definition at line 237 of file NitschePeriodicMethod.hpp.

---

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

• users_modules/homogenisation/src/NitschePeriodicMethod.hpp