OpZeroNoSideBrokenDofs Struct Reference

Zeroes the DOFs which aro not at side for broken base field. More...

#include "users_modules/eshelbian_plasticity/src/EshelbianOperators.hpp"

Inheritance diagram for OpZeroNoSideBrokenDofs:

Collaboration diagram for OpZeroNoSideBrokenDofs:

Public Types
using	OP = SideEleOp

Public Member Functions
	OpZeroNoSideBrokenDofs (std::string field_name)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Detailed Description

Zeroes the DOFs which aro not at side for broken base field.

Definition at line 985 of file EshelbianOperators.hpp.

Member Typedef Documentation

OP

using OpZeroNoSideBrokenDofs::OP = SideEleOp

Definition at line 986 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

OpZeroNoSideBrokenDofs()

OpZeroNoSideBrokenDofs::OpZeroNoSideBrokenDofs ( std::string  field_name )

inline

Definition at line 987 of file EshelbianOperators.hpp.

: OP(field_name, OP::OPROW) {}

Member Function Documentation

doWork()

MoFEMErrorCode OpZeroNoSideBrokenDofs::doWork	(	int	side,
		EntityType	type,
		EntData &	data
	)

Examples

mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 3712 of file EshelbianOperators.cpp.

                                                             {
 
  MoFEMFunctionBegin;
 
  auto nb_dofs = data.getIndices().size();
  if (nb_dofs == 0)
    MoFEMFunctionReturnHot(0);
 
  auto continuity = data.getFieldEntities()[0]->getContinuity();
  if (continuity != DISCONTINUOUS)
    SETERRQ(PETSC_COMM_WORLD, MOFEM_DATA_INCONSISTENCY,
            "OpZeroNoSideBrokenDofs works only for discontinuous field");
 
  auto fe_type = OP::getFEType();
  BaseFunction::DofsSideMap &side_dof_map =
      data.getFieldEntities()[0]->getDofSideMap().at(fe_type);
  auto face_side_number = getFaceSideNumber();
 
  for (
 
      auto it = side_dof_map.get<1>().begin();
      it != side_dof_map.get<1>().end(); ++it
 
  ) {
    if ((3 * it->dof) >= data.getIndices().size())
      break;
    if (it->type != MBTRI || it->side != face_side_number) {
      for (auto dd = 0; dd < 3; ++dd) {
        data.getFieldData()[3 * it->dof + dd] = 0.;
      }
    }
  }
 
  MoFEMFunctionReturn(0);
}

---

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

• users_modules/eshelbian_plasticity/src/EshelbianOperators.hpp
• users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp