SimpleContactProblem::SimpleContactElement Struct Reference

#include <users_modules/mortar_contact/src/SimpleContact.hpp>

Inheritance diagram for SimpleContactProblem::SimpleContactElement:

Collaboration diagram for SimpleContactProblem::SimpleContactElement:

Public Member Functions
	SimpleContactElement (MoFEM::Interface &m_field, bool newton_cotes=false)
MoFEMErrorCode	preProcess ()
MoFEMErrorCode	postProcess ()
int	getRule (int order)
virtual MoFEMErrorCode	setGaussPts (int order)

Public Attributes
MoFEM::Interface &	mField
bool	newtonCotes
SmartPetscObj< Vec >	contactStateVec
friend	ConvectSlaveIntegrationPts

Detailed Description

Definition at line 76 of file SimpleContact.hpp.

Constructor & Destructor Documentation

SimpleContactElement()

SimpleContactProblem::SimpleContactElement::SimpleContactElement ( MoFEM::Interface &  m_field, bool  newton_cotes = false  )

inline

Definition at line 82 of file SimpleContact.hpp.

        : ContactEle(m_field), mField(m_field), newtonCotes(newton_cotes),
          contactStateVec() {}

Member Function Documentation

getRule()

int SimpleContactProblem::SimpleContactElement::getRule ( int  order )

inline

Definition at line 119 of file SimpleContact.hpp.

                           {
      if (newtonCotes)
        return -1;
      else
        return 2 * order;
    }

postProcess()

MoFEMErrorCode SimpleContactProblem::SimpleContactElement::postProcess ( )

inline

Definition at line 97 of file SimpleContact.hpp.

                                 {
      MoFEMFunctionBegin;
 
      if (snes_ctx != CTX_SNESSETFUNCTION || !contactStateVec)
        MoFEMFunctionReturnHot(0);
 
      CHKERR VecAssemblyBegin(contactStateVec);
      CHKERR VecAssemblyEnd(contactStateVec);
 
      const double *array;
      CHKERR VecGetArrayRead(contactStateVec, &array);
      if (mField.get_comm_rank() == 0)
 
        MOFEM_LOG_C("WORLD", Sev::verbose,
                    "  Active Gauss pts: %d out of %d", (int)array[0],
                    (int)array[1]);
 
      CHKERR VecRestoreArrayRead(contactStateVec, &array);
 
      MoFEMFunctionReturn(0);
    }

preProcess()

MoFEMErrorCode SimpleContactProblem::SimpleContactElement::preProcess ( )

inline

Definition at line 86 of file SimpleContact.hpp.

                                {
      MoFEMFunctionBegin;
      if (snes_ctx == CTX_SNESSETFUNCTION && contactStateVec) {
        CHKERR VecAssemblyBegin(contactStateVec);
        CHKERR VecAssemblyEnd(contactStateVec);
 
        CHKERR VecZeroEntries(contactStateVec);
      }
      MoFEMFunctionReturn(0);
    }

setGaussPts()

MoFEMErrorCode SimpleContactProblem::SimpleContactElement::setGaussPts ( int  order )

virtual

Reimplemented in MortarContactProblem::MortarConvectSlaveContactElement, MortarContactProblem::MortarConvectMasterContactElement, SimpleContactProblem::ConvectSlaveContactElement, SimpleContactProblem::ConvectMasterContactElement, and MortarContactProblem::MortarContactElement.

Definition at line 29 of file SimpleContact.cpp.

                                                               {
  MoFEMFunctionBegin;
  if (newtonCotes) {
    int rule = order + 2;
    int nb_gauss_pts = IntRules::NCC::triangle_ncc_order_num(rule);
    gaussPtsMaster.resize(3, nb_gauss_pts, false);
    gaussPtsSlave.resize(3, nb_gauss_pts, false);
    double xy_coords[2 * nb_gauss_pts];
    double w_array[nb_gauss_pts];
    IntRules::NCC::triangle_ncc_rule(rule, nb_gauss_pts, xy_coords, w_array);
 
    for (int gg = 0; gg != nb_gauss_pts; ++gg) {
      gaussPtsMaster(0, gg) = xy_coords[gg * 2];
      gaussPtsMaster(1, gg) = xy_coords[gg * 2 + 1];
      gaussPtsMaster(2, gg) = w_array[gg];
      gaussPtsSlave(0, gg) = xy_coords[gg * 2];
      gaussPtsSlave(1, gg) = xy_coords[gg * 2 + 1];
      gaussPtsSlave(2, gg) = w_array[gg];
    }
  } else {
    CHKERR ContactEle::setDefaultGaussPts(2 * order);
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

contactStateVec

SmartPetscObj<Vec> SimpleContactProblem::SimpleContactElement::contactStateVec

Definition at line 80 of file SimpleContact.hpp.

ConvectSlaveIntegrationPts

friend SimpleContactProblem::SimpleContactElement::ConvectSlaveIntegrationPts

Definition at line 128 of file SimpleContact.hpp.

mField

MoFEM::Interface& SimpleContactProblem::SimpleContactElement::mField

Definition at line 78 of file SimpleContact.hpp.

newtonCotes

bool SimpleContactProblem::SimpleContactElement::newtonCotes

Definition at line 79 of file SimpleContact.hpp.

---

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

• users_modules/mortar_contact/src/SimpleContact.hpp
• users_modules/mortar_contact/src/impl/SimpleContact.cpp