ElementFragmentAlgorithm.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#pragma once

#include <vector>
#include <map>
#include <set>

#include "EFANode.h"
#include "EFAElement.h"

class ElementFragmentAlgorithm
{
public:
  ElementFragmentAlgorithm(std::ostream & os);

  ~ElementFragmentAlgorithm();

private:
  std::ostream & _ostream;
  // unsigned int MaxElemId;
  std::map<unsigned int, EFANode *> _permanent_nodes;
  std::map<unsigned int, EFANode *> _embedded_nodes;
  std::map<unsigned int, EFANode *> _temp_nodes;
  std::map<unsigned int, EFANode *> _embedded_permanent_nodes;
  std::map<unsigned int, EFAElement *> _elements;
  std::set<EFAElement *> _crack_tip_elements;
  std::vector<EFANode *> _new_nodes;
  std::vector<EFANode *> _deleted_nodes;
  std::vector<EFAElement *> _child_elements;
  std::vector<EFAElement *> _parent_elements;
  std::map<EFANode *, std::set<EFAElement *>> _inverse_connectivity;

public:
  unsigned int add2DElements(std::vector<std::vector<unsigned int>> & quads);
  EFAElement * add2DElement(const std::vector<unsigned int> & quad, unsigned int id);
  EFAElement * add3DElement(const std::vector<unsigned int> & quad, unsigned int id);

  void updateEdgeNeighbors();
  void initCrackTipTopology();
  void addElemEdgeIntersection(unsigned int elemid, unsigned int edgeid, double position);
  void addElemNodeIntersection(unsigned int elemid, unsigned int nodeid);
  bool addFragEdgeIntersection(unsigned int elemid, unsigned int frag_edge_id, double position);
  void addElemFaceIntersection(unsigned int elemid,
                               unsigned int faceid,
                               const std::vector<unsigned int> & edgeid,
                               const std::vector<double> & position);
  void addFragFaceIntersection(unsigned int ElemID,
                               unsigned int FragFaceID,
                               const std::vector<unsigned int> & FragFaceEdgeID,
                               const std::vector<double> & position);

  void updatePhysicalLinksAndFragments();

  void updateTopology(bool mergeUncutVirtualEdges = true);
  void reset();
  void clearAncestry();
  void restoreFragmentInfo(EFAElement * const elem, const EFAElement * const from_elem);

  void createChildElements();
  void connectFragments(bool mergeUncutVirtualEdges);

  void sanityCheck();
  void updateCrackTipElements();
  void printMesh();
  void error(const std::string & error_string);

  const std::vector<EFAElement *> & getChildElements() { return _child_elements; };
  const std::vector<EFAElement *> & getParentElements() { return _parent_elements; };
  const std::vector<EFANode *> & getNewNodes() { return _new_nodes; };
  const std::set<EFAElement *> & getCrackTipElements() { return _crack_tip_elements; };
  const std::map<unsigned int, EFANode *> & getPermanentNodes() { return _permanent_nodes; }
  const std::map<unsigned int, EFANode *> & getTempNodes() { return _temp_nodes; }
  const std::map<unsigned int, EFANode *> & getEmbeddedNodes() { return _embedded_nodes; }
  EFAElement * getElemByID(unsigned int id);
  unsigned int getElemIdByNodes(unsigned int * node_id);
  void clearPotentialIsolatedNodes();
};