FractureMechanics::ConstantArea::OpAreaC Struct Reference

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

Inheritance diagram for FractureMechanics::ConstantArea::OpAreaC:

Collaboration diagram for FractureMechanics::ConstantArea::OpAreaC:

Public Member Functions
	OpAreaC (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 506 of file ConstantArea.hpp.

Constructor & Destructor Documentation

OpAreaC()

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

inline

Definition at line 509 of file ConstantArea.hpp.

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

Member Function Documentation

doWork()

PetscErrorCode FractureMechanics::ConstantArea::OpAreaC::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 519 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);
 
      const EntityHandle *conn = getConn();
      for (int nn = 0; nn != 3; ++nn) {
        if (row_indices[nn] != -1) {
          if (row_data.getFieldDofs()[nn]->getEnt() != conn[nn]) {
            SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                    "Data inconsistency");
          }
        }
        for (int dd = 0; dd != 3; ++dd) {
          if (col_indices[3 * nn + dd] != -1) {
            if (col_data.getFieldDofs()[3 * nn + dd]->getEnt() != conn[nn]) {
              SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                      "Data inconsistency");
            }
          }
        }
      }
 
      auto row_dofs = getFEMethod()->getRowDofsPtr();
      for (auto it = row_dofs->begin(); it != row_dofs->end(); ++it) {
        int side = it->get()->getSideNumber();
        if (conn[side] != it->get()->getEnt()) {
          SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                  "Data inconsistency");
        }
      }
 
      auto col_dofs = getFEMethod()->getColDofsPtr();
      for (auto it = col_dofs->begin(); it != col_dofs->end(); ++it) {
        int side = it->get()->getSideNumber();
        if (conn[side] != it->get()->getEnt()) {
          SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                  "Data inconsistency");
        }
      }
 
      cMat.resize(3, 9, false);
      cMat.clear();
 
      CHKERR setVariables(this, col_side, col_type, col_data);
 
      int nb_dofs = col_data.getFieldData().size();
      commonData.crackAreaIncrement.resize(nb_dofs, false);
      int r;
      // play recorder for values
      r = ::function(tAg, nb_dofs, 18, &activeVariables[0],
                     &commonData.crackAreaIncrement[0]);
      if (r < 3) {
        SETERRQ1(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
                 "ADOL-C function evaluation with error r = %d", r);
      }
 
      auto check_isnormal = [](double v) {
        MoFEMFunctionBeginHot;
        if (v != 0 && !std::isnormal(v)) {
          SETERRQ1(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                   "Not a number %3.4e", v);
        }
        MoFEMFunctionReturnHot(0);
      };
 
      for (int rr = 0; rr != nb_rows; rr++) {
        if (row_indices[rr] != -1) {
          CHKERR check_isnormal(commonData.crackAreaIncrement[3 * rr + 0]);
          CHKERR check_isnormal(commonData.crackAreaIncrement[3 * rr + 1]);
          CHKERR check_isnormal(commonData.crackAreaIncrement[3 * rr + 2]);
          cMat(rr, 3 * rr + 0) = commonData.crackAreaIncrement[3 * rr + 0];
          cMat(rr, 3 * rr + 1) = commonData.crackAreaIncrement[3 * rr + 1];
          cMat(rr, 3 * rr + 2) = commonData.crackAreaIncrement[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::OpAreaC::cMat

Definition at line 517 of file ConstantArea.hpp.

---

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

• users_modules/fracture_mechanics/src/ConstantArea.hpp