FractureMechanics::ConstantArea::MyTriangleFE Struct Reference

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

Inheritance diagram for FractureMechanics::ConstantArea::MyTriangleFE:

Collaboration diagram for FractureMechanics::ConstantArea::MyTriangleFE:

Public Member Functions
	MyTriangleFE (MoFEM::Interface &m_field, CommonData &common_data, std::string field_name="LAMBDA_CRACKFRONT_AREA")
int	getRule (int order)
PetscErrorCode	preProcess ()
PetscErrorCode	postProcess ()

Public Attributes
CommonData &	commonData
std::string	fieldName
Mat	petscB
Vec	petscF
Mat	petscQ

Detailed Description

Definition at line 52 of file ConstantArea.hpp.

Constructor & Destructor Documentation

MyTriangleFE()

FractureMechanics::ConstantArea::MyTriangleFE::MyTriangleFE ( MoFEM::Interface &  m_field, CommonData &  common_data, std::string  field_name = "LAMBDA_CRACKFRONT_AREA"  )

inline

Definition at line 61 of file ConstantArea.hpp.

        : FaceElementForcesAndSourcesCore(m_field), commonData(common_data),
          fieldName(field_name), petscB(PETSC_NULL), petscF(PETSC_NULL),
          petscQ(PETSC_NULL) {}

Member Function Documentation

getRule()

int FractureMechanics::ConstantArea::MyTriangleFE::getRule ( int  order )

inline

Definition at line 66 of file ConstantArea.hpp.

{ return order; };

postProcess()

PetscErrorCode FractureMechanics::ConstantArea::MyTriangleFE::postProcess ( )

inline

Definition at line 132 of file ConstantArea.hpp.

                                 {
      MoFEMFunctionBegin;
      CHKERR FaceElementForcesAndSourcesCore::postProcess();
      if (petscQ != PETSC_NULL) {
        if (petscB == PETSC_NULL) {
          SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                  "Matrix B should be set");
        }
 
        // Q vectors
        void *void_ctx;
        CHKERR MatShellGetContext(petscQ, &void_ctx);
        ConstrainMatrixCtx *q_mat_ctx = (ConstrainMatrixCtx *)void_ctx;
 
        Vec Qv;
        CHKERR mField.getInterface<VecManager>()->vecCreateGhost(
            problemPtr->getName(), COL, &Qv);
        CHKERR MatAssemblyBegin(petscB, MAT_FLUSH_ASSEMBLY);
        CHKERR MatAssemblyEnd(petscB, MAT_FLUSH_ASSEMBLY);
        auto project_vectors_and_assemble = [this, Qv](auto dit) {
          MoFEMFunctionBegin;
          if (commonData.mapV.find(dit->get()->getPetscGlobalDofIdx()) ==
              commonData.mapV.end()) {
            SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                    "Vector on row not found");
          }
          CHKERR VecAssemblyBegin(
              commonData.mapV[dit->get()->getPetscGlobalDofIdx()]);
          CHKERR VecAssemblyEnd(
              commonData.mapV[dit->get()->getPetscGlobalDofIdx()]);
          CHKERR VecGhostUpdateBegin(
              commonData.mapV[dit->get()->getPetscGlobalDofIdx()], ADD_VALUES,
              SCATTER_REVERSE);
          CHKERR VecGhostUpdateEnd(
              commonData.mapV[dit->get()->getPetscGlobalDofIdx()], ADD_VALUES,
              SCATTER_REVERSE);
          CHKERR VecGhostUpdateBegin(
              commonData.mapV[dit->get()->getPetscGlobalDofIdx()],
              INSERT_VALUES, SCATTER_FORWARD);
          CHKERR VecGhostUpdateEnd(
              commonData.mapV[dit->get()->getPetscGlobalDofIdx()],
              INSERT_VALUES, SCATTER_FORWARD);
 
          CHKERR MatMult(
              petscQ, commonData.mapV[dit->get()->getPetscGlobalDofIdx()], Qv);
          CHKERR VecGhostUpdateBegin(Qv, INSERT_VALUES, SCATTER_FORWARD);
          CHKERR VecGhostUpdateEnd(Qv, INSERT_VALUES, SCATTER_FORWARD);
 
          if (0) {
            auto save_vec = [](std::string name, Vec v) {
              PetscViewer viewer;
              PetscViewerASCIIOpen(PETSC_COMM_WORLD, name.c_str(), &viewer);
              PetscViewerPushFormat(viewer, PETSC_VIEWER_ASCII_MATLAB);
              VecView(v, viewer);
              PetscViewerPopFormat(viewer);
              PetscViewerDestroy(&viewer);
            };
            save_vec(
                "V_" + boost::lexical_cast<std::string>(dit->get()->getEnt()) +
                    ".m",
                commonData.mapV[dit->get()->getPetscGlobalDofIdx()]);
            save_vec(
                "Qv_" + boost::lexical_cast<std::string>(dit->get()->getEnt()) +
                    ".m",
                Qv);
          }
 
          if (dit->get()->getPart() == mField.get_comm_rank()) {
            double *array;
            CHKERR VecGetArray(Qv, &array);
            vector<int> glob_idx;
            vector<double> vals;
            int row = dit->get()->getPetscGlobalDofIdx();
            for (_IT_NUMEREDDOF_COL_BY_BITNUMBER_FOR_LOOP_(
                     problemPtr,
                     mField.get_field_bit_number("MESH_NODE_POSITIONS"),
                     diit)) {
              int idx = diit->get()->getPetscGlobalDofIdx();
              double val = array[diit->get()->getPetscGlobalDofIdx()];
              glob_idx.push_back(idx);
              vals.push_back(val);
            }
            CHKERR MatSetValues(petscB, 1, &row, glob_idx.size(),
                                &*glob_idx.begin(), &*vals.begin(), ADD_VALUES);
            CHKERR VecRestoreArray(Qv, &array);
          }
 
          MoFEMFunctionReturn(0);
        };
 
        for (_IT_NUMEREDDOF_ROW_BY_BITNUMBER_FOR_LOOP_(
                 problemPtr, mField.get_field_bit_number(fieldName), dit)) {
          CHKERR project_vectors_and_assemble(dit);
          CHKERR VecDestroy(
              &commonData.mapV[dit->get()->getPetscGlobalDofIdx()]);
        }
 
        CHKERR VecDestroy(&Qv);
        commonData.mapV.clear();
      }
 
      MoFEMFunctionReturn(0);
    }

preProcess()

PetscErrorCode FractureMechanics::ConstantArea::MyTriangleFE::preProcess ( )

inline

Definition at line 68 of file ConstantArea.hpp.

                                {
      MoFEMFunctionBegin;
      CHKERR FaceElementForcesAndSourcesCore::preProcess();
 
      if (petscB != PETSC_NULL) {
        snes_B = petscB;
      }
 
      if (petscF != PETSC_NULL) {
        snes_f = petscF;
      }
 
      switch (ts_ctx) {
      case CTX_TSSETIFUNCTION: {
        if (!petscF) {
          snes_ctx = CTX_SNESSETFUNCTION;
          snes_f = ts_F;
        }
        break;
      }
      case CTX_TSSETIJACOBIAN: {
        if (!petscB) {
          snes_ctx = CTX_SNESSETJACOBIAN;
          snes_B = ts_B;
        }
        break;
      }
      default:
        break;
      }
 
      // Tag to detect which side is it
      CHKERR mField.get_moab().tag_get_handle("INTERFACE_SIDE",
                                              commonData.thInterfaceSide);
 
      // If petscQ matrix is defined appply projection to each row
      if (petscQ != PETSC_NULL) {
        if (petscB == PETSC_NULL) {
          SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                  "Matrix B should be set");
        }
        const auto field_bit_number = mField.get_field_bit_number(fieldName);
        for (_IT_NUMEREDDOF_ROW_BY_BITNUMBER_FOR_LOOP_(problemPtr, field_bit_number,
                                                  dit)) {
          CHKERR mField.getInterface<VecManager>()->vecCreateGhost(
              problemPtr->getName(), COL,
              &commonData.mapV[dit->get()->getPetscGlobalDofIdx()]);
          CHKERR VecZeroEntries(
              commonData.mapV[dit->get()->getPetscGlobalDofIdx()]);
          CHKERR VecGhostUpdateBegin(
              commonData.mapV[dit->get()->getPetscGlobalDofIdx()],
              INSERT_VALUES, SCATTER_FORWARD);
          CHKERR VecGhostUpdateEnd(
              commonData.mapV[dit->get()->getPetscGlobalDofIdx()],
              INSERT_VALUES, SCATTER_FORWARD);
        }
        commonData.setMapV = true;
      } else {
        commonData.setMapV = false;
      }
 
      MoFEMFunctionReturn(0);
    }

Member Data Documentation

commonData

CommonData& FractureMechanics::ConstantArea::MyTriangleFE::commonData

Definition at line 54 of file ConstantArea.hpp.

fieldName

std::string FractureMechanics::ConstantArea::MyTriangleFE::fieldName

Definition at line 55 of file ConstantArea.hpp.

petscB

Mat FractureMechanics::ConstantArea::MyTriangleFE::petscB

Definition at line 57 of file ConstantArea.hpp.

petscF

Vec FractureMechanics::ConstantArea::MyTriangleFE::petscF

Definition at line 58 of file ConstantArea.hpp.

petscQ

Mat FractureMechanics::ConstantArea::MyTriangleFE::petscQ

Definition at line 59 of file ConstantArea.hpp.

---

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

• users_modules/fracture_mechanics/src/ConstantArea.hpp