Reactions Struct Reference

calculate reactions from vector of internal forces on meshsets More...

#include "users_modules/basic_finite_elements/src/DirichletBC.hpp"

Collaboration diagram for Reactions:

Public Types
typedef std::map< int, VectorDouble >	ReactionsMap

Public Member Functions
	Reactions (MoFEM::Interface &m_field, string problem_name, string field_name)
const ReactionsMap &	getReactionsMap () const
	Get the Reactions Map.
const VectorDouble &	getReactionsFromSet (const int &id) const
	Get the Reactions at specified meshset id.
MoFEMErrorCode	calculateReactions (Vec &internal)
	calculate reactions from a given vector

Public Attributes
MoFEM::Interface &	mField

Private Attributes
std::string	problemName
std::string	fieldName
ReactionsMap	reactionsMap

Detailed Description

calculate reactions from vector of internal forces on meshsets

example usage

Vec F_int;
DMCreateGlobalVector_MoFEM(dm, &F_int);
 
feRhs->snes_ctx = FEMethod::CTX_SNESSETFUNCTION;
feRhs->snes_f = F_int;
DMoFEMLoopFiniteElements(dm, "ELASTIC", feRhs);
 
VecAssemblyBegin(F_int);
VecAssemblyEnd(F_int);
VecGhostUpdateBegin(F_int, INSERT_VALUES, SCATTER_FORWARD);
VecGhostUpdateEnd(F_int, INSERT_VALUES, SCATTER_FORWARD);
 
Reactions my_react(m_field, "DM_ELASTIC", "U");
my_react.calculateReactions(F_int);
int fix_nodes_meshset_id = 1;
cout << my_react.getReactionsFromSet(fix_nodes_meshset_id) << endl;

Definition at line 407 of file DirichletBC.hpp.

Member Typedef Documentation

ReactionsMap

std::map<int, VectorDouble> Reactions::ReactionsMap

Definition at line 412 of file DirichletBC.hpp.

Constructor & Destructor Documentation

Reactions()

Reactions::Reactions ( MoFEM::Interface & m_field, string problem_name, string field_name )

inline

Definition at line 409 of file DirichletBC.hpp.

      : mField(m_field), problemName(problem_name), fieldName(field_name) {}

Member Function Documentation

calculateReactions()

MoFEMErrorCode Reactions::calculateReactions

(

Vec &

internal

)

calculate reactions from a given vector

Parameters

internal

forces vector

Returns

MoFEMErrorCode

Definition at line 782 of file DirichletBC.cpp.

                                                          {
 
  MoFEMFunctionBegin;
 
  const Problem *problem_ptr;
  CHKERR mField.get_problem(problemName.c_str(), &problem_ptr);
  const double *array;
  CHKERR VecGetArrayRead(internal, &array);
 
  auto field_ptr = mField.get_field_structure(fieldName);
  const int nb_coefficients = field_ptr->getNbOfCoeffs();
 
  std::vector<int> ghosts(nb_coefficients);
  for (int g = 0; g != nb_coefficients; ++g)
    ghosts[g] = g;
 
  Vec v;
  CHKERR VecCreateGhost(
      mField.get_comm(), (mField.get_comm_rank() ? 0 : nb_coefficients),
      nb_coefficients, (mField.get_comm_rank() ? nb_coefficients : 0),
      &*ghosts.begin(), &v);
 
  for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(
           mField, NODESET | DISPLACEMENTSET, it)) {
 
    const int id = it->getMeshsetId();
    VectorDouble &reaction_vec = reactionsMap[id];
    reaction_vec.resize(nb_coefficients);
    reaction_vec.clear();
 
    Range verts;
    for (int dim = 0; dim != 3; ++dim) {
      Range ents;
      CHKERR it->getMeshsetIdEntitiesByDimension(mField.get_moab(), dim, ents,
                                                 true);
      Range nodes;
      CHKERR mField.get_moab().get_connectivity(ents, nodes, true);
      verts.insert(nodes.begin(), nodes.end());
    }
 
    auto for_each_dof = [&](auto &dof) {
      MoFEMFunctionBeginHot;
      reaction_vec[dof->getDofCoeffIdx()] += array[dof->getPetscLocalDofIdx()];
      MoFEMFunctionReturnHot(0);
    };
 
    CHKERR set_numered_dofs_on_ents(problem_ptr,
                                    mField.get_field_bit_number(fieldName),
                                    verts, for_each_dof);
 
    double *res_array;
 
    CHKERR VecGetArray(v, &res_array);
    for (int dd = 0; dd != reaction_vec.size(); ++dd)
      res_array[dd] = reaction_vec[dd];
    CHKERR VecRestoreArray(v, &res_array);
 
    CHKERR VecGetArray(v, &res_array);
    for (int dd = 0; dd != reaction_vec.size(); ++dd)
      reaction_vec[dd] = res_array[dd];
    CHKERR VecRestoreArray(v, &res_array);
  }
 
  CHKERR VecGhostUpdateBegin(v, ADD_VALUES, SCATTER_REVERSE);
  CHKERR VecGhostUpdateEnd(v, ADD_VALUES, SCATTER_REVERSE);
  CHKERR VecGhostUpdateBegin(v, INSERT_VALUES, SCATTER_FORWARD);
  CHKERR VecGhostUpdateEnd(v, INSERT_VALUES, SCATTER_FORWARD);
 
  CHKERR VecDestroy(&v);
  CHKERR VecRestoreArrayRead(internal, &array);
  MoFEMFunctionReturn(0);
}

getReactionsFromSet()

const VectorDouble & Reactions::getReactionsFromSet ( const int & id ) const

inline

Get the Reactions at specified meshset id.

Parameters

id	meshset id (from Cubit)

Returns

const VectorDouble&

Definition at line 426 of file DirichletBC.hpp.

                                                                      {
    return reactionsMap.at(id);
  }

getReactionsMap()

const ReactionsMap & Reactions::getReactionsMap ( ) const

inline

Get the Reactions Map.

Returns

const ReactionsMap&

Definition at line 419 of file DirichletBC.hpp.

{ return reactionsMap; }

Member Data Documentation

fieldName

std::string Reactions::fieldName

private

Definition at line 439 of file DirichletBC.hpp.

mField

MoFEM::Interface& Reactions::mField

Definition at line 413 of file DirichletBC.hpp.

problemName

std::string Reactions::problemName

private

Definition at line 438 of file DirichletBC.hpp.

reactionsMap

ReactionsMap Reactions::reactionsMap

private

Definition at line 440 of file DirichletBC.hpp.

---

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

• users_modules/basic_finite_elements/src/DirichletBC.hpp
• users_modules/basic_finite_elements/src/impl/DirichletBC.cpp