anonymous_namespace{prenotch.cpp} Namespace Reference

Functions
template<typename Predicate >
std::size_t	breakBondsIf (geometry::Fracture &fracture, const std::vector< util::Point > &xRef, const std::vector< std::vector< std::size_t > > &neighbors, const std::vector< std::size_t > &nodeParticleId, std::size_t particleId, Predicate crosses)

Function Documentation

breakBondsIf()

template<typename Predicate >

std::size_t anonymous_namespace{prenotch.cpp}::breakBondsIf	(	geometry::Fracture &	fracture,
		const std::vector< util::Point > &	xRef,
		const std::vector< std::vector< std::size_t > > &	neighbors,
		const std::vector< std::size_t > &	nodeParticleId,
		std::size_t	particleId,
		Predicate	crosses
	)

Breaks the bonds of one particle for which the predicate is true.

Definition at line 19 of file prenotch.cpp.

                       {
  if (neighbors.size() > xRef.size() ||
      nodeParticleId.size() < neighbors.size())
    throw std::invalid_argument(
        "prenotch: node arrays are inconsistent with the neighbour list");
 
  std::size_t broken = 0;
  for (std::size_t i = 0; i < neighbors.size(); ++i) {
    if (nodeParticleId[i] != particleId)
      continue;
    const auto &xi = xRef[i];
    for (std::size_t k = 0; k < neighbors[i].size(); ++k) {
      const std::size_t j = neighbors[i][k];
      if (j >= nodeParticleId.size() || nodeParticleId[j] != particleId)
        continue;
      const auto &xj = xRef[j];
      if (crosses(xi.d_x, xi.d_y, xj.d_x, xj.d_y)) {
        fracture.setBondState(i, k, true);
        ++broken;
      }
    }
  }
  return broken;
}

References geometry::Fracture::setBondState().

Here is the call graph for this function: