FractureMechanics::ConstantArea::OpTangentC Struct Reference

#include <users_modules/fracture_mechanics/src/ConstantArea.hpp>

Inheritance diagram for FractureMechanics::ConstantArea::OpTangentC:

Collaboration diagram for FractureMechanics::ConstantArea::OpTangentC:

Public Member Functions
	OpTangentC (int tag, CommonData &common_data)
PetscErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::EntData &col_data)
Public Member Functions inherited from FractureMechanics::ConstantArea::AuxOp
	AuxOp (int tag, CommonData &common_data)
PetscErrorCode	setVariables (FaceElementForcesAndSourcesCore::UserDataOperator *fe_ptr, int side, EntityType type, DataForcesAndSourcesCore::EntData &data)

Public Attributes
MatrixDouble	cMat
Public Attributes inherited from FractureMechanics::ConstantArea::AuxOp
int	tAg
CommonData &	commonData
ublas::vector< int >	rowIndices
VectorDouble	activeVariables

Detailed Description

Definition at line 700 of file ConstantArea.hpp.

Constructor & Destructor Documentation

OpTangentC()

FractureMechanics::ConstantArea::OpTangentC::OpTangentC ( int  tag, CommonData &  common_data  )

inline

Definition at line 703 of file ConstantArea.hpp.

        : FaceElementForcesAndSourcesCore::UserDataOperator(
              "LAMBDA_CRACKFRONT_AREA_TANGENT", "MESH_NODE_POSITIONS",
              UserDataOperator::OPROWCOL),
          AuxOp(tag, common_data) {
      sYmm = false;
    }

Member Function Documentation

doWork()

PetscErrorCode FractureMechanics::ConstantArea::OpTangentC::doWork ( int  row_side, int  col_side, EntityType  row_type, EntityType  col_type, DataForcesAndSourcesCore::EntData &  row_data, DataForcesAndSourcesCore::EntData &  col_data  )

inline

Definition at line 713 of file ConstantArea.hpp.

                                                                     {
      MoFEMFunctionBegin;
 
      if (row_type != MBVERTEX) {
        MoFEMFunctionReturnHot(0);
      }
 
      VectorInt &row_indices = row_data.getIndices();
      const int nb_rows = row_indices.size();
      if (!nb_rows)
        MoFEMFunctionReturnHot(0);
      VectorInt &col_indices = col_data.getIndices();
      const int nb_cols = col_indices.size();
      if (!nb_cols)
        MoFEMFunctionReturnHot(0);
      cMat.resize(3, 9, false);
      cMat.clear();
 
      CHKERR setVariables(this, col_side, col_type, col_data);
 
      int nb_dofs = col_data.getFieldData().size();
      commonData.crackFrontTangent.resize(nb_dofs, false);
      int r;
      // play recorder for values
      r = ::function(tAg, nb_dofs, 18, &activeVariables[0],
                     &commonData.crackFrontTangent[0]);
      if (r < 3) {
        SETERRQ1(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
                 "ADOL-C function evaluation with error r = %d", r);
      }
 
      // Set face orientation
      {
        // Get orientation in respect to adjacent tet
        const EntityHandle face = getFEEntityHandle();
        const BitRefLevel bit = getFEMethod()->problemPtr->getBitRefLevel();
        Range adj_side_elems;
        CHKERR commonData.mField.getInterface<BitRefManager>()->getAdjacencies(
            bit, &face, 1, 3, adj_side_elems);
        adj_side_elems = adj_side_elems.subset_by_type(MBTET);
        if (adj_side_elems.size() != 1) {
          SETERRQ1(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                   "expect 1 tet but is %u", adj_side_elems.size());
        }
        EntityHandle side_elem = *adj_side_elems.begin();
        int side_number, sense, offset;
        CHKERR commonData.mField.get_moab().side_number(
            side_elem, face, side_number, sense, offset);
        if (sense == -1) {
          commonData.crackFrontTangent *= -1;
        }
        // Get orientaton of face in respect to face in the prism
        int side;
        CHKERR commonData.mField.get_moab().tag_get_data(
            commonData.thInterfaceSide, &face, 1, &side);
        if (side == 1) {
          commonData.crackFrontTangent *= -1;
        }
      }
 
      auto check_isnormal = [](double v) {
        MoFEMFunctionBeginHot;
        if (v != 0 && !std::isnormal(v)) {
          SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "Not a number");
        }
        MoFEMFunctionReturnHot(0);
      };
 
      for (int rr = 0; rr != nb_rows; ++rr) {
        if (row_indices[rr] != -1) {
          // cerr << commonData.crackFrontTangent(3*rr+0) << " "
          // << commonData.crackFrontTangent(3*rr+1) << " "
          // << commonData.crackFrontTangent(3*rr+2) << endl;
          CHKERR check_isnormal(commonData.crackFrontTangent[3 * rr + 0]);
          CHKERR check_isnormal(commonData.crackFrontTangent[3 * rr + 1]);
          CHKERR check_isnormal(commonData.crackFrontTangent[3 * rr + 2]);
          cMat(rr, 3 * rr + 0) = commonData.crackFrontTangent[3 * rr + 0];
          cMat(rr, 3 * rr + 1) = commonData.crackFrontTangent[3 * rr + 1];
          cMat(rr, 3 * rr + 2) = commonData.crackFrontTangent[3 * rr + 2];
        }
      }
 
      if (!commonData.setMapV) {
        CHKERR MatSetValues(getFEMethod()->snes_B, nb_rows, &row_indices[0],
                            nb_cols, &col_indices[0], &cMat(0, 0), ADD_VALUES);
      } else {
        for (int rr = 0; rr != nb_rows; ++rr) {
          if (row_indices[rr] != -1) {
            if (commonData.mapV.find(row_indices[rr]) ==
                commonData.mapV.end()) {
              SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                      "Vector for given DOF in map not found");
            }
            for (int cc = 0; cc != nb_cols; ++cc) {
              if (col_indices[cc] != -1) {
                CHKERR check_isnormal(cMat(rr, cc));
              }
            }
            CHKERR VecSetValues(commonData.mapV[row_indices[rr]], nb_cols,
                                &col_indices[0], &cMat(rr, 0), ADD_VALUES);
          }
        }
      }
 
      MoFEMFunctionReturn(0);
    }

Member Data Documentation

cMat

MatrixDouble FractureMechanics::ConstantArea::OpTangentC::cMat

Definition at line 711 of file ConstantArea.hpp.

---

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

• users_modules/fracture_mechanics/src/ConstantArea.hpp