NitscheMethod::OpGetFaceData Struct Reference

Get integration pts data on face. More...

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

Inheritance diagram for NitscheMethod::OpGetFaceData:

Collaboration diagram for NitscheMethod::OpGetFaceData:

Public Member Functions
	OpGetFaceData (CommonData &common_data)
MoFEMErrorCode	doWork (int side, EntityType type, DataForcesAndSourcesCore::EntData &data)

Public Attributes
CommonData &	commonData

Detailed Description

Get integration pts data on face.

Definition at line 153 of file NitscheMethod.hpp.

Constructor & Destructor Documentation

OpGetFaceData()

NitscheMethod::OpGetFaceData::OpGetFaceData

(

CommonData &

common_data

)

Definition at line 158 of file NitscheMethod.hpp.

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

Member Function Documentation

doWork()

MoFEMErrorCode NitscheMethod::OpGetFaceData::doWork	(	int	side,
		EntityType	type,
		DataForcesAndSourcesCore::EntData &	data
	)

Definition at line 163 of file NitscheMethod.hpp.

                                                                 {
      MoFEMFunctionBeginHot;
 
      int faceInRespectToTet = getFEMethod()->nInTheLoop;
      int nb_face_gauss_pts = getGaussPts().size2();
 
      try {
        if (type == MBVERTEX) {
          commonData.faceGaussPts.resize(4);
          commonData.faceGaussPts[faceInRespectToTet] = getGaussPts();
          for (int fgg = 0; fgg < nb_face_gauss_pts; fgg++) {
            int gg = commonData.nbTetGaussPts;
            commonData.inTetFaceGaussPtsNumber[faceInRespectToTet].push_back(
                gg);
            commonData.nbTetGaussPts++;
          }
          commonData.faceNormals.resize(4);
          commonData.faceNormals[faceInRespectToTet] =
              0.5 * getNormalsAtGaussPtss();
          commonData.hoCoordsAtGaussPts.resize(4);
          commonData.hoCoordsAtGaussPts[faceInRespectToTet] =
              getCoordsAtGaussPts();
          commonData.cOords.resize(4);
          commonData.cOords[faceInRespectToTet] = getCoords();
          commonData.coordsAtGaussPts.resize(4);
          commonData.coordsAtGaussPts[faceInRespectToTet] =
              getCoordsAtGaussPts();
          commonData.rAy.resize(4);
          commonData.rAy[faceInRespectToTet] = -getNormal();
          commonData.rAy[faceInRespectToTet] /= norm_2(getNormal());
        }
      } catch (const std::exception &ex) {
        std::ostringstream ss;
        ss << "throw in method: " << ex.what() << std::endl;
        SETERRQ(PETSC_COMM_SELF, 1, ss.str().c_str());
      }
 
      try {
        for (int fgg = 0; fgg < nb_face_gauss_pts; fgg++) {
          int gg = commonData.inTetFaceGaussPtsNumber[faceInRespectToTet][fgg];
          // std::cerr << fgg << " " << gg << " " << side << " " << type <<
          // std::endl;
          CommonData::MultiIndexData gauss_pt_data(gg, side, type);
          std::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);
          // std::cerr << "nb_shape_fun " << nb_shape_fun << std::endl;
          if (nb_shape_fun) {
            cblas_dcopy(nb_shape_fun, &data.getN()(fgg, 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();
          }
          // std::cerr << shape_fun << std::endl;
        }
      } catch (const std::exception &ex) {
        std::ostringstream ss;
        ss << "throw in method: " << ex.what() << std::endl;
        SETERRQ(PETSC_COMM_SELF, 1, ss.str().c_str());
      }
      MoFEMFunctionReturnHot(0);
    }

Member Data Documentation

commonData

CommonData& NitscheMethod::OpGetFaceData::commonData

Definition at line 156 of file NitscheMethod.hpp.

---

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

• users_modules/homogenisation/src/NitscheMethod.hpp