ElementFragmentAlgorithm Class Reference

#include <ElementFragmentAlgorithm.h>

Public Member Functions
	ElementFragmentAlgorithm (std::ostream &os)
	Constructor.
	~ElementFragmentAlgorithm ()
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 ()
const std::vector< EFAElement * > &	getParentElements ()
const std::vector< EFANode * > &	getNewNodes ()
const std::set< EFAElement * > &	getCrackTipElements ()
const std::map< unsigned int, EFANode * > &	getPermanentNodes ()
const std::map< unsigned int, EFANode * > &	getTempNodes ()
const std::map< unsigned int, EFANode * > &	getEmbeddedNodes ()
EFAElement *	getElemByID (unsigned int id)
unsigned int	getElemIdByNodes (unsigned int *node_id)
void	clearPotentialIsolatedNodes ()

Private Attributes
std::ostream &	_ostream
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

Detailed Description

Definition at line 19 of file ElementFragmentAlgorithm.h.

Constructor & Destructor Documentation

ElementFragmentAlgorithm()

ElementFragmentAlgorithm::ElementFragmentAlgorithm

(

std::ostream &

os

)

Constructor.

Definition at line 30 of file ElementFragmentAlgorithm.C.

: _ostream(os) {}

~ElementFragmentAlgorithm()

ElementFragmentAlgorithm::~ElementFragmentAlgorithm

(

)

Definition at line 32 of file ElementFragmentAlgorithm.C.

{
  std::map<unsigned int, EFANode *>::iterator mit;
  for (mit = _permanent_nodes.begin(); mit != _permanent_nodes.end(); ++mit)
  {
    delete mit->second;
    mit->second = nullptr;
  }
  for (mit = _embedded_nodes.begin(); mit != _embedded_nodes.end(); ++mit)
  {
    delete mit->second;
    mit->second = nullptr;
  }
  for (mit = _embedded_permanent_nodes.begin(); mit != _embedded_permanent_nodes.end(); ++mit)
  {
    delete mit->second;
    mit->second = nullptr;
  }
  for (mit = _temp_nodes.begin(); mit != _temp_nodes.end(); ++mit)
  {
    delete mit->second;
    mit->second = nullptr;
  }
  std::map<unsigned int, EFAElement *>::iterator eit;
  for (eit = _elements.begin(); eit != _elements.end(); ++eit)
  {
    delete eit->second;
    eit->second = nullptr;
  }
}

Member Function Documentation

add2DElement()

EFAElement * ElementFragmentAlgorithm::add2DElement	(	const std::vector< unsigned int > &	quad,
		unsigned int	id
	)

Definition at line 102 of file ElementFragmentAlgorithm.C.

{
  unsigned int num_nodes = quad.size();
 
  std::map<unsigned int, EFAElement *>::iterator mit = _elements.find(id);
  if (mit != _elements.end())
    EFAError("In add2DElement element with id: ", id, " already exists");
 
  EFAElement2D * newElem = new EFAElement2D(id, num_nodes);
  _elements.insert(std::make_pair(id, newElem));
 
  for (unsigned int j = 0; j < num_nodes; ++j)
  {
    EFANode * currNode = nullptr;
    std::map<unsigned int, EFANode *>::iterator mit = _permanent_nodes.find(quad[j]);
    if (mit == _permanent_nodes.end())
    {
      currNode = new EFANode(quad[j], EFANode::N_CATEGORY_PERMANENT);
      _permanent_nodes.insert(std::make_pair(quad[j], currNode));
    }
    else
      currNode = mit->second;
 
    newElem->setNode(j, currNode);
    _inverse_connectivity[currNode].insert(newElem);
  }
  newElem->createEdges();
  return newElem;
}

Referenced by XFEM::buildEFAMesh(), and TEST().

add2DElements()

unsigned int ElementFragmentAlgorithm::add2DElements

(

std::vector< std::vector< unsigned int > > &

quads

)

Definition at line 64 of file ElementFragmentAlgorithm.C.

{
  unsigned int first_id = 0;
  unsigned int num_nodes = quads[0].size();
 
  if (quads.size() == 0)
    EFAError("add2DElements called with empty vector of quads");
 
  for (unsigned int i = 0; i < quads.size(); ++i)
  {
    unsigned int new_elem_id = Efa::getNewID(_elements);
    EFAElement2D * newElem = new EFAElement2D(new_elem_id, num_nodes);
    _elements.insert(std::make_pair(new_elem_id, newElem));
 
    if (i == 0)
      first_id = new_elem_id;
 
    for (unsigned int j = 0; j < num_nodes; ++j)
    {
      EFANode * currNode = nullptr;
      std::map<unsigned int, EFANode *>::iterator mit = _permanent_nodes.find(quads[i][j]);
      if (mit == _permanent_nodes.end())
      {
        currNode = new EFANode(quads[i][j], EFANode::N_CATEGORY_PERMANENT);
        _permanent_nodes.insert(std::make_pair(quads[i][j], currNode));
      }
      else
        currNode = mit->second;
 
      newElem->setNode(j, currNode);
      _inverse_connectivity[currNode].insert(newElem);
    }
    newElem->createEdges();
  }
  return first_id;
}

Referenced by case1Common(), case2Mesh(), case5Mesh(), and TEST().

add3DElement()

EFAElement * ElementFragmentAlgorithm::add3DElement	(	const std::vector< unsigned int > &	quad,
		unsigned int	id
	)

Definition at line 133 of file ElementFragmentAlgorithm.C.

{
  unsigned int num_nodes = quad.size();
  unsigned int num_faces = 0;
  if (num_nodes == 8 || num_nodes == 20 || num_nodes == 27)
    num_faces = 6;
  else if (num_nodes == 4 || num_nodes == 10 || num_nodes == 14)
    num_faces = 4;
  else
    EFAError("In add3DElement element with id: ", id, " has invalid num_nodes");
 
  std::map<unsigned int, EFAElement *>::iterator mit = _elements.find(id);
  if (mit != _elements.end())
    EFAError("In add3DElement element with id: ", id, " already exists");
 
  EFAElement3D * newElem = new EFAElement3D(id, num_nodes, num_faces);
  _elements.insert(std::make_pair(id, newElem));
 
  for (unsigned int j = 0; j < num_nodes; ++j)
  {
    EFANode * currNode = nullptr;
    std::map<unsigned int, EFANode *>::iterator mit = _permanent_nodes.find(quad[j]);
    if (mit == _permanent_nodes.end())
    {
      currNode = new EFANode(quad[j], EFANode::N_CATEGORY_PERMANENT);
      _permanent_nodes.insert(std::make_pair(quad[j], currNode));
    }
    else
      currNode = mit->second;
 
    newElem->setNode(j, currNode);
    _inverse_connectivity[currNode].insert(newElem);
  }
  newElem->createFaces();
  return newElem;
}

Referenced by XFEM::buildEFAMesh(), case6Mesh(), and case7Mesh().

addElemEdgeIntersection()

void ElementFragmentAlgorithm::addElemEdgeIntersection	(	unsigned int	elemid,
		unsigned int	edgeid,
		double	position
	)

Definition at line 202 of file ElementFragmentAlgorithm.C.

{
  // this method is called when we are marking cut edges
  auto eit = _elements.find(elemid);
  if (eit == _elements.end())
    EFAError("Could not find element with id: ", elemid, " in addEdgeIntersection");
 
  EFAElement2D * curr_elem = dynamic_cast<EFAElement2D *>(eit->second);
  if (!curr_elem)
    EFAError("addElemEdgeIntersection: elem ", elemid, " is not of type EFAelement2D");
  curr_elem->addEdgeCut(edgeid, position, nullptr, _embedded_nodes, true);
}

Referenced by case1Common(), case2Intersections(), XFEM::markCutEdgesByGeometry(), XFEM::markCutEdgesByState(), TEST(), TEST(), TEST(), TEST(), TEST(), and TEST().

addElemFaceIntersection()

void ElementFragmentAlgorithm::addElemFaceIntersection	(	unsigned int	elemid,
		unsigned int	faceid,
		const std::vector< unsigned int > &	edgeid,
		const std::vector< double > &	position
	)

Definition at line 251 of file ElementFragmentAlgorithm.C.

{
  // this method is called when we are marking cut edges
  auto eit = _elements.find(elemid);
  if (eit == _elements.end())
    EFAError("Could not find element with id: ", elemid, " in addEdgeIntersection");
 
  EFAElement3D * curr_elem = dynamic_cast<EFAElement3D *>(eit->second);
  if (!curr_elem)
    EFAError("addElemEdgeIntersection: elem ", elemid, " is not of type EFAelement2D");
 
  // add cuts to two face edges at the same time
  curr_elem->addFaceEdgeCut(faceid, edgeid[0], position[0], nullptr, _embedded_nodes, true, true);
  curr_elem->addFaceEdgeCut(faceid, edgeid[1], position[1], nullptr, _embedded_nodes, true, true);
}

Referenced by XFEM::markCutFacesByGeometry(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), and TEST().

addElemNodeIntersection()

void ElementFragmentAlgorithm::addElemNodeIntersection	(	unsigned int	elemid,
		unsigned int	nodeid
	)

Definition at line 218 of file ElementFragmentAlgorithm.C.

{
  // this method is called when we are marking cut nodes
  auto eit = _elements.find(elemid);
  if (eit == _elements.end())
    EFAError("Could not find element with id: ", elemid, " in addElemNodeIntersection");
 
  EFAElement2D * curr_elem = dynamic_cast<EFAElement2D *>(eit->second);
  if (!curr_elem)
    EFAError("addElemNodeIntersection: elem ", elemid, " is not of type EFAelement2D");
 
  // Only add cut node when the curr_elem does not have any fragment
  if (curr_elem->numFragments() == 0)
    curr_elem->addNodeCut(nodeid, nullptr, _permanent_nodes, _embedded_permanent_nodes);
}

Referenced by XFEM::markCutEdgesByGeometry().

addFragEdgeIntersection()

bool ElementFragmentAlgorithm::addFragEdgeIntersection	(	unsigned int	elemid,
		unsigned int	frag_edge_id,
		double	position
	)

Definition at line 235 of file ElementFragmentAlgorithm.C.

{
  // N.B. this method must be called after addEdgeIntersection
  auto eit = _elements.find(elemid);
  if (eit == _elements.end())
    EFAError("Could not find element with id: ", elemid, " in addFragEdgeIntersection");
 
  EFAElement2D * elem = dynamic_cast<EFAElement2D *>(eit->second);
  if (!elem)
    EFAError("addFragEdgeIntersection: elem ", elemid, " is not of type EFAelement2D");
  return elem->addFragmentEdgeCut(frag_edge_id, position, _embedded_nodes);
}

Referenced by XFEM::markCutEdgesByGeometry(), XFEM::markCutEdgesByState(), TEST(), and TEST().

addFragFaceIntersection()

void ElementFragmentAlgorithm::addFragFaceIntersection	(	unsigned int	ElemID,
		unsigned int	FragFaceID,
		const std::vector< unsigned int > &	FragFaceEdgeID,
		const std::vector< double > &	position
	)

Definition at line 271 of file ElementFragmentAlgorithm.C.

{
  // TODO: need to finish this for 3D problems
}

Referenced by XFEM::markCutFacesByGeometry().

clearAncestry()

void ElementFragmentAlgorithm::clearAncestry

(

)

Definition at line 350 of file ElementFragmentAlgorithm.C.

{
  _inverse_connectivity.clear();
  for (unsigned int i = 0; i < _parent_elements.size(); ++i)
  {
    if (!Efa::deleteFromMap(_elements, _parent_elements[i]))
      EFAError("Attempted to delete parent element: ",
               _parent_elements[i]->id(),
               " from _elements, but couldn't find it");
  }
  _parent_elements.clear();
 
  std::map<unsigned int, EFAElement *>::iterator eit;
  for (eit = _elements.begin(); eit != _elements.end(); ++eit)
  {
    EFAElement * curr_elem = eit->second;
    curr_elem->clearParentAndChildren();
    for (unsigned int j = 0; j < curr_elem->numNodes(); j++)
    {
      EFANode * curr_node = curr_elem->getNode(j);
      _inverse_connectivity[curr_node].insert(curr_elem);
    }
  }
 
  std::map<unsigned int, EFANode *>::iterator mit;
  for (mit = _permanent_nodes.begin(); mit != _permanent_nodes.end(); ++mit)
    mit->second->removeParent();
 
  for (mit = _temp_nodes.begin(); mit != _temp_nodes.end(); ++mit)
  {
    delete mit->second;
    mit->second = nullptr;
  }
  _temp_nodes.clear();
 
  _new_nodes.clear();
  _child_elements.clear();
 
  // TODO: Sanity check to make sure that there are no nodes that are not connected
  //      to an element -- there shouldn't be any
}

Referenced by TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), and TEST().

clearPotentialIsolatedNodes()

void ElementFragmentAlgorithm::clearPotentialIsolatedNodes

(

)

Definition at line 600 of file ElementFragmentAlgorithm.C.

{
  // compile a set of all child element nodes
  std::set<EFANode *> child_nodes;
  for (const auto & child_element : _child_elements)
  {
    const auto & nodes = child_element->getNodes();
    child_nodes.insert(nodes.begin(), nodes.end());
  }
 
  // Collect all parent nodes that will be isolated
  std::map<EFANode *, EFANode *> isolate_parent_to_child;
  for (unsigned int i = 0; i < _new_nodes.size(); ++i)
  {
    EFANode * parent_node = _new_nodes[i]->parent();
    if (!parent_node)
      EFAError("a new permanent node must have a parent node!");
 
    auto it = isolate_parent_to_child.lower_bound(parent_node);
    if (it != isolate_parent_to_child.end() && it->first == parent_node)
      continue;
 
    if (child_nodes.count(parent_node) == 0)
      isolate_parent_to_child.emplace_hint(it, parent_node, _new_nodes[i]);
  }
 
  // For each isolated parent node, pick one of its child new node
  // Then, switch that child with its parent for all new elems
  for (const auto [parent_node, child_node] : isolate_parent_to_child)
  {
    for (unsigned int i = 0; i < _child_elements.size(); ++i)
      if (_child_elements[i]->containsNode(child_node))
        _child_elements[i]->switchNode(parent_node, child_node, true);
 
    _new_nodes.erase(std::remove(_new_nodes.begin(), _new_nodes.end(), child_node),
                     _new_nodes.end());
    Efa::deleteFromMap(_permanent_nodes, child_node);
  }
}

Referenced by updateTopology().

connectFragments()

void ElementFragmentAlgorithm::connectFragments

(

bool

mergeUncutVirtualEdges

)

Definition at line 423 of file ElementFragmentAlgorithm.C.

{
  // now perform the comparison on the children
  for (unsigned int elem_iter = 0; elem_iter < _child_elements.size(); elem_iter++)
  {
    EFAElement * childElem = _child_elements[elem_iter];
    childElem->connectNeighbors(
        _permanent_nodes, _temp_nodes, _inverse_connectivity, mergeUncutVirtualEdges);
    childElem->updateFragmentNode();
  }
 
  // remove all deleted children
  _child_elements.erase(std::remove(_child_elements.begin(), _child_elements.end(), nullptr),
                        _child_elements.end());
}

Referenced by updateTopology().

createChildElements()

void ElementFragmentAlgorithm::createChildElements

(

)

Definition at line 400 of file ElementFragmentAlgorithm.C.

{
  // temporary container for new elements -- will be merged with Elements
  std::map<unsigned int, EFAElement *> newChildElements;
 
  // loop over the original elements in the mesh
  std::map<unsigned int, EFAElement *>::iterator eit;
  std::map<unsigned int, EFAElement *>::iterator ElementsEnd = _elements.end();
  for (eit = _elements.begin(); eit != ElementsEnd; ++eit)
  {
    EFAElement * curr_elem = eit->second;
    curr_elem->createChild(_crack_tip_elements,
                           _elements,
                           newChildElements,
                           _child_elements,
                           _parent_elements,
                           _temp_nodes);
  } // loop over elements
  // Merge newChildElements back in with Elements
  _elements.insert(newChildElements.begin(), newChildElements.end());
}

Referenced by updateTopology().

error()

void ElementFragmentAlgorithm::error

(

const std::string &

error_string

)

getChildElements()

const std::vector< EFAElement * > & ElementFragmentAlgorithm::getChildElements ( )

inline

Definition at line 78 of file ElementFragmentAlgorithm.h.

{ return _child_elements; };

Referenced by XFEM::cutMeshWithEFA(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), and TEST().

getCrackTipElements()

const std::set< EFAElement * > & ElementFragmentAlgorithm::getCrackTipElements ( )

inline

Definition at line 81 of file ElementFragmentAlgorithm.h.

{ return _crack_tip_elements; };

Referenced by XFEM::cutMeshWithEFA(), XFEM::storeCrackTipOriginAndDirection(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), and TEST().

getElemByID()

EFAElement * ElementFragmentAlgorithm::getElemByID

(

unsigned int

id

)

Definition at line 568 of file ElementFragmentAlgorithm.C.

{
  std::map<unsigned int, EFAElement *>::iterator mit = _elements.find(id);
  if (mit == _elements.end())
    EFAError("in getElemByID() could not find element: ", id);
  return mit->second;
}

Referenced by XFEM::buildEFAMesh(), XFEM::getEFAElem2D(), and XFEM::getEFAElem3D().

getElemIdByNodes()

unsigned int ElementFragmentAlgorithm::getElemIdByNodes

(

unsigned int *

node_id

)

Definition at line 577 of file ElementFragmentAlgorithm.C.

{
  unsigned int elem_id = std::numeric_limits<unsigned int>::max();
  std::map<unsigned int, EFAElement *>::iterator eit;
  for (eit = _elements.begin(); eit != _elements.end(); ++eit)
  {
    EFAElement * curr_elem = eit->second;
    unsigned int counter = 0;
    for (unsigned int i = 0; i < curr_elem->numNodes(); ++i)
    {
      if (curr_elem->getNode(i)->id() == node_id[i])
        counter += 1;
    }
    if (counter == curr_elem->numNodes())
    {
      elem_id = curr_elem->id();
      break;
    }
  }
  return elem_id;
}

getEmbeddedNodes()

const std::map< unsigned int, EFANode * > & ElementFragmentAlgorithm::getEmbeddedNodes ( )

inline

Definition at line 84 of file ElementFragmentAlgorithm.h.

{ return _embedded_nodes; }

Referenced by TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), and TEST().

getNewNodes()

const std::vector< EFANode * > & ElementFragmentAlgorithm::getNewNodes ( )

inline

Definition at line 80 of file ElementFragmentAlgorithm.h.

{ return _new_nodes; };

Referenced by XFEM::cutMeshWithEFA().

getParentElements()

const std::vector< EFAElement * > & ElementFragmentAlgorithm::getParentElements ( )

inline

Definition at line 79 of file ElementFragmentAlgorithm.h.

{ return _parent_elements; };

Referenced by XFEM::cutMeshWithEFA(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), and TEST().

getPermanentNodes()

const std::map< unsigned int, EFANode * > & ElementFragmentAlgorithm::getPermanentNodes ( )

inline

Definition at line 82 of file ElementFragmentAlgorithm.h.

{ return _permanent_nodes; }

Referenced by TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), and TEST().

getTempNodes()

const std::map< unsigned int, EFANode * > & ElementFragmentAlgorithm::getTempNodes ( )

inline

Definition at line 83 of file ElementFragmentAlgorithm.h.

{ return _temp_nodes; }

Referenced by TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), and TEST().

initCrackTipTopology()

void ElementFragmentAlgorithm::initCrackTipTopology

(

)

Definition at line 194 of file ElementFragmentAlgorithm.C.

{
  _crack_tip_elements.clear(); // re-build CrackTipElements!
  for (auto pair : _elements)
    pair.second->initCrackTip(_crack_tip_elements); // CrackTipElements changed here
}

Referenced by XFEM::buildEFAMesh(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), and TEST().

printMesh()

void ElementFragmentAlgorithm::printMesh

(

)

Definition at line 482 of file ElementFragmentAlgorithm.C.

{
  _ostream << "============================================================"
           << "==================================================" << std::endl;
  _ostream << "                                            CutElemMesh Data" << std::endl;
  _ostream << "============================================================"
           << "==================================================" << std::endl;
  _ostream << "Permanent Nodes:" << std::endl;
  std::map<unsigned int, EFANode *>::iterator mit;
  unsigned int counter = 0;
  for (mit = _permanent_nodes.begin(); mit != _permanent_nodes.end(); ++mit)
  {
    _ostream << "  " << mit->second->id();
    counter += 1;
    if (counter % 10 == 0)
      _ostream << std::endl;
  }
  _ostream << std::endl;
  _ostream << "Temp Nodes:" << std::endl;
  counter = 0;
  for (mit = _temp_nodes.begin(); mit != _temp_nodes.end(); ++mit)
  {
    _ostream << "  " << mit->second->id();
    counter += 1;
    if (counter % 10 == 0)
      _ostream << std::endl;
  }
  _ostream << std::endl;
  _ostream << "Embedded Nodes:" << std::endl;
  counter = 0;
  for (mit = _embedded_nodes.begin(); mit != _embedded_nodes.end(); ++mit)
  {
    _ostream << "  " << mit->second->id();
    counter += 1;
    if (counter % 10 == 0)
      _ostream << std::endl;
  }
  _ostream << std::endl;
  _ostream << "Embedded Permanent Nodes:" << std::endl;
  counter = 0;
  for (mit = _embedded_permanent_nodes.begin(); mit != _embedded_permanent_nodes.end(); ++mit)
  {
    _ostream << "  " << mit->second->id();
    counter += 1;
    if (counter % 10 == 0)
      _ostream << std::endl;
  }
  _ostream << std::endl;
  _ostream << "Parent Elements:" << std::endl;
  counter = 0;
  for (unsigned int i = 0; i < _parent_elements.size(); ++i)
  {
    _ostream << " " << _parent_elements[i]->id();
    counter += 1;
    if (counter % 10 == 0)
      _ostream << std::endl;
  }
  _ostream << std::endl;
  _ostream << "Child Elements:" << std::endl;
  counter = 0;
  for (unsigned int i = 0; i < _child_elements.size(); ++i)
  {
    _ostream << " " << _child_elements[i]->id();
    counter += 1;
    if (counter % 10 == 0)
      _ostream << std::endl;
  }
  _ostream << std::endl;
  _ostream << "Elements:" << std::endl;
  _ostream << "  id "
           << "|  nodes                "
           << "|  embedded nodes       "
           << "|  edge neighbors       "
           << "|  frag "
           << "|  frag link      ...   " << std::endl;
  _ostream << "------------------------------------------------------------"
           << "--------------------------------------------------" << std::endl;
  std::map<unsigned int, EFAElement *>::iterator eit;
  for (eit = _elements.begin(); eit != _elements.end(); ++eit)
  {
    EFAElement * currElem = eit->second;
    currElem->printElement(_ostream);
  }
}

Referenced by XFEM::cutMeshWithEFA(), sanityCheck(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), and TEST().

reset()

void ElementFragmentAlgorithm::reset

(

)

Definition at line 318 of file ElementFragmentAlgorithm.C.

{
  _new_nodes.clear();
  _child_elements.clear();
  _parent_elements.clear();
  _crack_tip_elements.clear();
  _inverse_connectivity.clear();
 
  std::map<unsigned int, EFANode *>::iterator mit;
  for (mit = _permanent_nodes.begin(); mit != _permanent_nodes.end(); ++mit)
  {
    delete mit->second;
    mit->second = nullptr;
  }
  _permanent_nodes.clear();
 
  for (mit = _temp_nodes.begin(); mit != _temp_nodes.end(); ++mit)
  {
    delete mit->second;
    mit->second = nullptr;
  }
  _temp_nodes.clear();
  std::map<unsigned int, EFAElement *>::iterator eit;
  for (eit = _elements.begin(); eit != _elements.end(); ++eit)
  {
    delete eit->second;
    eit->second = nullptr;
  }
  _elements.clear();
}

Referenced by XFEM::buildEFAMesh().

restoreFragmentInfo()

void ElementFragmentAlgorithm::restoreFragmentInfo	(	EFAElement *const	elem,
		const EFAElement *const	from_elem
	)

Definition at line 393 of file ElementFragmentAlgorithm.C.

{
  elem->restoreFragment(from_elem);
}

Referenced by XFEM::buildEFAMesh().

sanityCheck()

void ElementFragmentAlgorithm::sanityCheck

(

)

Definition at line 440 of file ElementFragmentAlgorithm.C.

{
  // Make sure there are no remaining TempNodes
  if (_temp_nodes.size() > 0)
  {
    _ostream << "_temp_nodes size > 0.  size=" << _temp_nodes.size() << std::endl;
    printMesh();
    throw std::runtime_error("_temp_nodes size > 0");
  }
}

Referenced by updateTopology().

updateCrackTipElements()

void ElementFragmentAlgorithm::updateCrackTipElements

(

)

Definition at line 452 of file ElementFragmentAlgorithm.C.

{
  std::set<EFAElement *>::iterator sit;
  // Delete all elements that were previously flagged as crack tip elements if they have
  // been split (and hence appear in ParentElements).
  for (unsigned int i = 0; i < _parent_elements.size(); ++i)
  {
    sit = _crack_tip_elements.find(_parent_elements[i]);
    if (sit != _crack_tip_elements.end())
      _crack_tip_elements.erase(sit);
  }
 
  // Go through new child elements to find elements that are newly at the crack tip due to
  // crack growth.
  for (unsigned int elem_iter = 0; elem_iter < _child_elements.size(); elem_iter++)
  {
    EFAElement * childElem = _child_elements[elem_iter];
    if (childElem->isCrackTipElement())
      _crack_tip_elements.insert(childElem);
  } // loop over (new) child elements
 
  //_ostream << "Crack tip elements: ";
  // for (sit=CrackTipElements.begin(); sit!=CrackTipElements.end(); ++sit)
  //{
  //  _ostream << (*sit)->id<<" ";
  //}
  //_ostream << std::endl;
}

Referenced by updateTopology().

updateEdgeNeighbors()

void ElementFragmentAlgorithm::updateEdgeNeighbors

(

)

Definition at line 171 of file ElementFragmentAlgorithm.C.

{
  std::map<unsigned int, EFAElement *>::iterator eit;
  for (eit = _elements.begin(); eit != _elements.end(); ++eit)
  {
    EFAElement * elem = eit->second;
    elem->clearNeighbors();
  }
 
  for (eit = _elements.begin(); eit != _elements.end(); ++eit)
  {
    EFAElement * curr_elem = eit->second;
    curr_elem->setupNeighbors(_inverse_connectivity);
  } // loop over all elements
 
  for (eit = _elements.begin(); eit != _elements.end(); ++eit)
  {
    EFAElement * curr_elem = eit->second;
    curr_elem->neighborSanityCheck();
  }
}

Referenced by XFEM::buildEFAMesh(), case1Common(), case2Mesh(), case5Mesh(), case6Mesh(), case7Mesh(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), and TEST().

updatePhysicalLinksAndFragments()

void ElementFragmentAlgorithm::updatePhysicalLinksAndFragments

(

)

Definition at line 281 of file ElementFragmentAlgorithm.C.

{
  // loop over the elements in the mesh
  std::map<unsigned int, EFAElement *>::iterator eit;
  for (eit = _elements.begin(); eit != _elements.end(); ++eit)
  {
    EFAElement * curr_elem = eit->second;
    curr_elem->updateFragments(_crack_tip_elements, _embedded_nodes);
  } // loop over all elements
}

Referenced by case1Common(), case2Intersections(), XFEM::cutMeshWithEFA(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), and TEST().

updateTopology()

void ElementFragmentAlgorithm::updateTopology

(

bool

mergeUncutVirtualEdges = true

)

Definition at line 293 of file ElementFragmentAlgorithm.C.

{
  // If mergeUncutVirtualEdges=true, this algorithm replicates the
  // behavior of classical XFEM.  If false, it gives the behavior of
  // the Richardson et. al. (2011) paper
 
  _new_nodes.clear();
  _child_elements.clear();
  _parent_elements.clear();
 
  unsigned int first_new_node_id = Efa::getNewID(_permanent_nodes);
 
  createChildElements();
  connectFragments(mergeUncutVirtualEdges);
  sanityCheck();
  updateCrackTipElements();
 
  for (const auto & [id, node] : _permanent_nodes)
    if (id >= first_new_node_id)
      _new_nodes.push_back(node);
 
  clearPotentialIsolatedNodes(); // _new_nodes and _permanent_nodes may change here
}

Referenced by XFEM::cutMeshWithEFA(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), and TEST().

Member Data Documentation

_child_elements

std::vector<EFAElement *> ElementFragmentAlgorithm::_child_elements

private

Definition at line 40 of file ElementFragmentAlgorithm.h.

Referenced by clearAncestry(), clearPotentialIsolatedNodes(), connectFragments(), createChildElements(), getChildElements(), printMesh(), reset(), updateCrackTipElements(), and updateTopology().

_crack_tip_elements

std::set<EFAElement *> ElementFragmentAlgorithm::_crack_tip_elements

private

Definition at line 37 of file ElementFragmentAlgorithm.h.

Referenced by createChildElements(), getCrackTipElements(), initCrackTipTopology(), reset(), updateCrackTipElements(), and updatePhysicalLinksAndFragments().

_deleted_nodes

std::vector<EFANode *> ElementFragmentAlgorithm::_deleted_nodes

private

Definition at line 39 of file ElementFragmentAlgorithm.h.

_elements

std::map<unsigned int, EFAElement *> ElementFragmentAlgorithm::_elements

private

Definition at line 36 of file ElementFragmentAlgorithm.h.

Referenced by add2DElement(), add2DElements(), add3DElement(), addElemEdgeIntersection(), addElemFaceIntersection(), addElemNodeIntersection(), addFragEdgeIntersection(), clearAncestry(), createChildElements(), getElemByID(), getElemIdByNodes(), initCrackTipTopology(), printMesh(), reset(), updateEdgeNeighbors(), updatePhysicalLinksAndFragments(), and ~ElementFragmentAlgorithm().

_embedded_nodes

std::map<unsigned int, EFANode *> ElementFragmentAlgorithm::_embedded_nodes

private

Definition at line 33 of file ElementFragmentAlgorithm.h.

Referenced by addElemEdgeIntersection(), addElemFaceIntersection(), addFragEdgeIntersection(), getEmbeddedNodes(), printMesh(), updatePhysicalLinksAndFragments(), and ~ElementFragmentAlgorithm().

_embedded_permanent_nodes

std::map<unsigned int, EFANode *> ElementFragmentAlgorithm::_embedded_permanent_nodes

private

Definition at line 35 of file ElementFragmentAlgorithm.h.

Referenced by addElemNodeIntersection(), printMesh(), and ~ElementFragmentAlgorithm().

_inverse_connectivity

std::map<EFANode *, std::set<EFAElement *> > ElementFragmentAlgorithm::_inverse_connectivity

private

Definition at line 42 of file ElementFragmentAlgorithm.h.

Referenced by add2DElement(), add2DElements(), add3DElement(), clearAncestry(), connectFragments(), reset(), and updateEdgeNeighbors().

_new_nodes

std::vector<EFANode *> ElementFragmentAlgorithm::_new_nodes

private

Definition at line 38 of file ElementFragmentAlgorithm.h.

Referenced by clearAncestry(), clearPotentialIsolatedNodes(), getNewNodes(), reset(), and updateTopology().

_ostream

std::ostream& ElementFragmentAlgorithm::_ostream

private

Definition at line 30 of file ElementFragmentAlgorithm.h.

Referenced by printMesh(), and sanityCheck().

_parent_elements

std::vector<EFAElement *> ElementFragmentAlgorithm::_parent_elements

private

Definition at line 41 of file ElementFragmentAlgorithm.h.

Referenced by clearAncestry(), createChildElements(), getParentElements(), printMesh(), reset(), updateCrackTipElements(), and updateTopology().

_permanent_nodes

std::map<unsigned int, EFANode *> ElementFragmentAlgorithm::_permanent_nodes

private

Definition at line 32 of file ElementFragmentAlgorithm.h.

Referenced by add2DElement(), add2DElements(), add3DElement(), addElemNodeIntersection(), clearAncestry(), clearPotentialIsolatedNodes(), connectFragments(), getPermanentNodes(), printMesh(), reset(), updateTopology(), and ~ElementFragmentAlgorithm().

_temp_nodes

std::map<unsigned int, EFANode *> ElementFragmentAlgorithm::_temp_nodes

private

Definition at line 34 of file ElementFragmentAlgorithm.h.

Referenced by clearAncestry(), connectFragments(), createChildElements(), getTempNodes(), printMesh(), reset(), sanityCheck(), and ~ElementFragmentAlgorithm().

---

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

• xfem/include/efa/ElementFragmentAlgorithm.h
• xfem/src/efa/ElementFragmentAlgorithm.C