#include "src/finite_elements/FormsBrokenSpaceConstraintImpl.hpp"

Inheritance diagram for MoFEM::OpBrokenTopoBase< OpBrokenBase >:

Collaboration diagram for MoFEM::OpBrokenTopoBase< OpBrokenBase >:

Public Types
using	OP = OpBrokenBase

Protected Member Functions
	OpBrokenTopoBase (const std::string row_field, boost::shared_ptr< MatrixDouble > tangent1_diff_ptr, boost::shared_ptr< MatrixDouble > tangent2_diff_ptr, SmartPetscObj< Vec > assemble_vec, Tag th, boost::shared_ptr< Range > ents_ptr=nullptr)

MoFEMErrorCode	aSsemble (EntitiesFieldData::EntData &row_data)

Protected Attributes
boost::shared_ptr< MatrixDouble >	tangent1DiffPtr

boost::shared_ptr< MatrixDouble >	tangent2DiffPtr

SmartPetscObj< Vec >	assembleVec

Tag	thGradTag

Detailed Description

template<typename OpBrokenBase>
struct MoFEM::OpBrokenTopoBase< OpBrokenBase >

Definition at line 668 of file FormsBrokenSpaceConstraintImpl.hpp.

Member Typedef Documentation

◆ OP

template<typename OpBrokenBase >

using MoFEM::OpBrokenTopoBase< OpBrokenBase >::OP = OpBrokenBase

Definition at line 669 of file FormsBrokenSpaceConstraintImpl.hpp.

Constructor & Destructor Documentation

◆ OpBrokenTopoBase()

template<typename OpBrokenBase >

MoFEM::OpBrokenTopoBase< OpBrokenBase >::OpBrokenTopoBase	(	const std::string	row_field,
		boost::shared_ptr< MatrixDouble >	tangent1_diff_ptr,
		boost::shared_ptr< MatrixDouble >	tangent2_diff_ptr,
		SmartPetscObj< Vec >	assemble_vec,
		Tag	th,
		boost::shared_ptr< Range >	ents_ptr = `nullptr`
	)

inlineprotected

Definition at line 673 of file FormsBrokenSpaceConstraintImpl.hpp.

      : OP(row_field, row_field, OP::OPROW, ents_ptr),
        tangent1DiffPtr(tangent1_diff_ptr),
        tangent2DiffPtr(tangent2_diff_ptr), assembleVec(assemble_vec),
        thGradTag(th) {}

Member Function Documentation

◆ aSsemble()

template<typename OpBrokenBase >

MoFEMErrorCode MoFEM::OpBrokenTopoBase< OpBrokenBase >::aSsemble ( EntitiesFieldData::EntData & row_data )

inlineprotected

Definition at line 684 of file FormsBrokenSpaceConstraintImpl.hpp.

                                                              {
    MoFEMFunctionBegin;
    if (!this->timeScalingFun.empty())
      this->locF *= this->timeScalingFun(this->getFEMethod()->ts_t);
    if (!this->feScalingFun.empty())
      this->locF *= this->feScalingFun(this->getFEMethod());
    if (assembleVec) {
      auto *vec_ptr = this->locF.data().data();
      const auto nb_dofs = row_data.getIndices().size();
      auto *ind_ptr = row_data.getIndices().data().data();
      return VecSetValues(assembleVec, nb_dofs, ind_ptr, vec_ptr, ADD_VALUES);
    }
    if (thGradTag) {
      const auto field_ents = row_data.getFieldEntities();
      std::vector<EntityHandle> ents(field_ents.size());
      std::transform(field_ents.begin(), field_ents.end(), ents.begin(),
                     [](const auto *fe) { return fe->getEnt(); });
      if (field_ents.empty())
        MoFEMFunctionReturnHot(0);
      if (type_from_handle(ents[0]) != MBVERTEX)
        MoFEMFunctionReturnHot(0);
      auto &moab = this->getMoab();
      VectorDouble topo_values(this->locF.size());
      CHKERR moab.tag_get_data(thGradTag, ents.data(), ents.size(),
                               topo_values.data().data());
      topo_values += this->locF;
      CHKERR moab.tag_set_data(thGradTag, ents.data(), ents.size(),
                               topo_values.data().data());
    }
    MoFEMFunctionReturn(0);
  }