#include <XFEMElementPairLocator.h>

Inheritance diagram for XFEMElementPairLocator:

Public Types
typedef std::list< std::pair< const Elem , const Elem > >	ElementPairList

Public Member Functions
	XFEMElementPairLocator (std::shared_ptr< XFEM > xfem, unsigned int interface_id, bool use_displaced_mesh=false)

virtual void	reinit ()

virtual void	update ()

const ElementPairList &	getElemPairs () const

const ElementPairInfo &	getElemPairInfo (std::pair< const Elem , const Elem > elem_pair) const

Protected Attributes
std::shared_ptr< XFEM >	_xfem

bool	_use_displaced_mesh

const ElementPairList *	_elem_pairs

std::map< std::pair< const Elem , const Elem >, ElementPairInfo >	_element_pair_info

unsigned int	_interface_id

Detailed Description

Definition at line 15 of file XFEMElementPairLocator.h.

Constructor & Destructor Documentation

◆ XFEMElementPairLocator()

XFEMElementPairLocator::XFEMElementPairLocator	(	std::shared_ptr< XFEM >	xfem,
		unsigned int	interface_id,
		bool	use_displaced_mesh = `false`
	)

Definition at line 12 of file XFEMElementPairLocator.C.

  : ElementPairLocator(interface_id), _xfem(xfem), _use_displaced_mesh(use_displaced_mesh)
{
  if (_use_displaced_mesh)
    _elem_pairs = _xfem->getXFEMDisplacedCutElemPairs(interface_id);
  else
    _elem_pairs = _xfem->getXFEMCutElemPairs(interface_id);
}

Member Function Documentation

◆ reinit()

void XFEMElementPairLocator::reinit ( )

virtual

Reimplemented from ElementPairLocator.

Definition at line 24 of file XFEMElementPairLocator.C.

{
  // Does not support secondary cut yet.
  if (_xfem->has_secondary_cut())
    return;
 
  _element_pair_info.clear();
 
  for (std::list<std::pair<const Elem *, const Elem *>>::const_iterator it = _elem_pairs->begin();
       it != _elem_pairs->end();
       ++it)
  {
    const Elem * elem1 = it->first;
    const Elem * elem2 = it->second;
 
    std::vector<Point> intersectionPoints1;
    Point normal1;
    std::vector<Point> q_points1;
    std::vector<Real> weights1;
 
    unsigned int plane_id = 0; // Only support one cut plane for the time being
 
    _xfem->getXFEMIntersectionInfo(
        elem1, plane_id, normal1, intersectionPoints1, _use_displaced_mesh);
 
    if (intersectionPoints1.size() == 2)
      _xfem->getXFEMqRuleOnLine(intersectionPoints1, q_points1, weights1);
    else if (intersectionPoints1.size() > 2)
      _xfem->getXFEMqRuleOnSurface(intersectionPoints1, q_points1, weights1);
 
    if (!_use_displaced_mesh)
    {
      ElementPairInfo new_elem_info(
          elem1, elem2, q_points1, q_points1, weights1, weights1, normal1, -normal1);
      _element_pair_info.insert(
          std::pair<std::pair<const Elem *, const Elem *>, ElementPairInfo>(*it, new_elem_info));
    }
    else
    {
      std::vector<Point> intersectionPoints2;
      Point normal2;
      std::vector<Point> q_points2;
      std::vector<Real> weights2;
 
      _xfem->getXFEMIntersectionInfo(
          elem2, plane_id, normal2, intersectionPoints2, _use_displaced_mesh);
 
      // reverse the order of intersectionPoints2
      std::reverse(std::begin(intersectionPoints2), std::end(intersectionPoints2));
 
      if (intersectionPoints2.size() == 2)
        _xfem->getXFEMqRuleOnLine(intersectionPoints2, q_points2, weights2);
      else if (intersectionPoints2.size() > 2)
        _xfem->getXFEMqRuleOnSurface(intersectionPoints2, q_points2, weights2);
 
      ElementPairInfo new_elem_info(
          elem1, elem2, q_points1, q_points2, weights1, weights2, normal1, normal2);
      _element_pair_info.insert(
          std::pair<std::pair<const Elem *, const Elem *>, ElementPairInfo>(*it, new_elem_info));
    }
  }
}

Referenced by update().

◆ update()

void XFEMElementPairLocator::update ( )

virtual

Reimplemented from ElementPairLocator.

Definition at line 88 of file XFEMElementPairLocator.C.

{
  reinit();
}

Member Data Documentation

◆ _use_displaced_mesh

bool XFEMElementPairLocator::_use_displaced_mesh

protected

Definition at line 26 of file XFEMElementPairLocator.h.

Referenced by reinit(), and XFEMElementPairLocator().

◆ _xfem

std::shared_ptr<XFEM> XFEMElementPairLocator::_xfem

protected

Definition at line 25 of file XFEMElementPairLocator.h.