#include "src/finite_elements/FormsBrokenSpaceConstraintImpl.hpp"

Inheritance diagram for MoFEM::OpBrokenTopoBase< OpBrokenBase >:

Collaboration diagram for MoFEM::OpBrokenTopoBase< OpBrokenBase >:

Public Types
using	OP = OpBrokenBase

Public Member Functions
	OpBrokenTopoBase (boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_base_side_data, boost::shared_ptr< MatrixDouble > adjoint_lambda_ptr, boost::shared_ptr< MatrixDouble > tangent1_diff_ptr, boost::shared_ptr< MatrixDouble > tangent2_diff_ptr, boost::shared_ptr< double > beta_ptr, SmartPetscObj< Vec > assemble_vec, Tag th, boost::shared_ptr< Range > ents_ptr=nullptr)

	OpBrokenTopoBase (const std::string row_field, boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_side_data_ptr, boost::shared_ptr< MatrixDouble > adjoint_lambda_ptr, boost::shared_ptr< MatrixDouble > tangent1_ptr, boost::shared_ptr< MatrixDouble > tangent2_ptr, SmartPetscObj< Vec > assemble_vec, Tag th, boost::shared_ptr< Range > ents_ptr=nullptr)

Protected Member Functions
MoFEMErrorCode	aSsemble (EntitiesFieldData::EntData &row_data)

Protected Attributes
boost::shared_ptr< MatrixDouble >	adjointLambdaPtr

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 593 of file FormsBrokenSpaceConstraintImpl.hpp.

Member Typedef Documentation

◆ OP

template<typename OpBrokenBase >

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

Definition at line 594 of file FormsBrokenSpaceConstraintImpl.hpp.

Constructor & Destructor Documentation

◆ OpBrokenTopoBase() [1/2]

template<typename OpBrokenBase >

MoFEM::OpBrokenTopoBase< OpBrokenBase >::OpBrokenTopoBase	(	boost::shared_ptr< std::vector< BrokenBaseSideData > >	broken_base_side_data,
		boost::shared_ptr< MatrixDouble >	adjoint_lambda_ptr,
		boost::shared_ptr< MatrixDouble >	tangent1_diff_ptr,
		boost::shared_ptr< MatrixDouble >	tangent2_diff_ptr,
		boost::shared_ptr< double >	beta_ptr,
		SmartPetscObj< Vec >	assemble_vec,
		Tag	th,
		boost::shared_ptr< Range >	ents_ptr = `nullptr`
	)

inline

Definition at line 596 of file FormsBrokenSpaceConstraintImpl.hpp.

      : OP(broken_base_side_data, ents_ptr),
        adjointLambdaPtr(adjoint_lambda_ptr),
        tangent1DiffPtr(tangent1_diff_ptr), tangent2DiffPtr(tangent2_diff_ptr),
        assembleVec(assemble_vec), thGradTag(th) {}

◆ OpBrokenTopoBase() [2/2]

template<typename OpBrokenBase >

MoFEM::OpBrokenTopoBase< OpBrokenBase >::OpBrokenTopoBase	(	const std::string	row_field,
		boost::shared_ptr< std::vector< BrokenBaseSideData > >	broken_side_data_ptr,
		boost::shared_ptr< MatrixDouble >	adjoint_lambda_ptr,
		boost::shared_ptr< MatrixDouble >	tangent1_ptr,
		boost::shared_ptr< MatrixDouble >	tangent2_ptr,
		SmartPetscObj< Vec >	assemble_vec,
		Tag	th,
		boost::shared_ptr< Range >	ents_ptr = `nullptr`
	)

inline

Definition at line 608 of file FormsBrokenSpaceConstraintImpl.hpp.

      : OpBase(row_field, row_field, OpBase::OPROW, ents_ptr),
        brokenSideDataPtr(broken_side_data_ptr),
        adjointLambdaPtr(adjoint_lambda_ptr), tangent1Ptr(tangent1_ptr),
        tangent2Ptr(tangent2_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 622 of file FormsBrokenSpaceConstraintImpl.hpp.

                                                              {
    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 = 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 = getMoab();
      VectorDouble topo_values(this->locF.size());
      CHKERR moab.tag_set_data(thGradTag, ents.data(), ents.size(),
                               topo_values.data().data());
      topo_values += this->locF;
      CHKERR moab.tag_set_data(thGradTag, ents.data(), ents.size(),
                               this->locF.data().data());
    }
  }