PeriodicNitscheConstrains::OpGetFaceData Struct Reference

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

Inheritance diagram for PeriodicNitscheConstrains::OpGetFaceData:

Collaboration diagram for PeriodicNitscheConstrains::OpGetFaceData:

Public Member Functions
	OpGetFaceData (CommonData &common_data, bool field_disp=true)
PetscErrorCode	doWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)

Public Attributes
CommonData &	commonData
bool	fieldDisp

Detailed Description

Definition at line 83 of file NitschePeriodicMethod.hpp.

Constructor & Destructor Documentation

OpGetFaceData()

PeriodicNitscheConstrains::OpGetFaceData::OpGetFaceData ( CommonData &  common_data, bool  field_disp = true  )

inline

Definition at line 88 of file NitschePeriodicMethod.hpp.

                                                                 :
    FaceElementForcesAndSourcesCore::UserDataOperator("DISPLACEMENT",OPROW),
    commonData(common_data),
    fieldDisp(field_disp) {
    }

Member Function Documentation

doWork()

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

inline

Definition at line 94 of file NitschePeriodicMethod.hpp.

                                                                                          {
      PetscFunctionBegin;
 
      if(data.getIndices().size()==0) PetscFunctionReturn(0);
      try {
        EntityHandle this_face = getNumeredEntFiniteElementPtr()->getEnt();
        int nb_gauss_pts_on_this_face = data.getN().size1();
        for(int ffgg = 0;ffgg<nb_gauss_pts_on_this_face;ffgg++) {
          int gg = commonData.inTetFaceGaussPtsNumber[this_face][ffgg];
          CommonData::MultiIndexData gauss_pt_data(gg,side,type);
          pair<CommonData::Container::iterator,bool> p;
          p = commonData.facesContainer.insert(gauss_pt_data);
          if(!p.second) {
            SETERRQ(PETSC_COMM_SELF,MOFEM_DATA_INCONSISTENCY,"data not inserted");
          }
          CommonData::MultiIndexData &p_data = const_cast<CommonData::MultiIndexData&>(*p.first);
          VectorDouble &shape_fun = p_data.shapeFunctions;
          int nb_shape_fun = data.getN().size2();
          shape_fun.resize(nb_shape_fun);
          cblas_dcopy(nb_shape_fun,&data.getN()(ffgg,0),1,&shape_fun[0],1);
          p_data.iNdices = data.getIndices();
          p_data.dofOrders.resize(data.getFieldDofs().size(),false);
          for(unsigned int dd = 0;dd<data.getFieldDofs().size();dd++) {
            p_data.dofOrders[dd] = data.getFieldDofs()[dd]->getDofOrder();
          }
          int nb_dofs = data.getFieldData().size();
          if(type == MBVERTEX) {
            commonData.dIsplacements[gg].resize(3,false);
            commonData.coordsAtGaussPts[gg].resize(3,false);
            commonData.hoCoordsAtGaussPts[gg].resize(3,false);
            for(int rr = 0;rr<3;rr++) {
              (commonData.dIsplacements[gg])[rr] = cblas_ddot(
                nb_dofs/3,&data.getFieldData()[rr],3,&data.getN()(ffgg,0),1
              );
              if(!fieldDisp) {
                (commonData.dIsplacements[gg])[rr] -= getCoordsAtGaussPts()(ffgg,rr);
              }
              commonData.coordsAtGaussPts[gg][rr] = getCoordsAtGaussPts()(ffgg,rr);
              commonData.hoCoordsAtGaussPts[gg][rr] = getCoordsAtGaussPts()(ffgg,rr);
            }
          } else {
            for(int rr = 0;rr<3;rr++) {
              commonData.dIsplacements[gg][rr] += cblas_ddot(
                nb_dofs/3,&data.getFieldData()[rr],3,&data.getN()(ffgg,0),1
              );
            }
          }
        }
      } catch (const std::exception& ex) {
        ostringstream ss;
        ss << "throw in method: " << ex.what() << endl;
        SETERRQ(PETSC_COMM_SELF,1,ss.str().c_str());
      }
 
      //cerr << "done\n";
      PetscFunctionReturn(0);
    }

Member Data Documentation

commonData

CommonData& PeriodicNitscheConstrains::OpGetFaceData::commonData

Definition at line 85 of file NitschePeriodicMethod.hpp.

fieldDisp

bool PeriodicNitscheConstrains::OpGetFaceData::fieldDisp

Definition at line 86 of file NitschePeriodicMethod.hpp.

---

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

• users_modules/homogenisation/src/NitschePeriodicMethod.hpp