CSGUtils Namespace Reference

Functions
CSG::CSGRegion	getInnerRegion (const std::vector< std::reference_wrapper< const CSG::CSGSurface >> &surfaces, const libMesh::Point &origin)
	Get inner region of given surfaces, defined as the intersection of halfspaces of each surface. More...
void	checkValidCSGName (const std::string &name)
	Check name of CSG component for disallowed characters and symbols. More...

Function Documentation

checkValidCSGName()

void CSGUtils::checkValidCSGName

(

const std::string &

name

)

Check name of CSG component for disallowed characters and symbols.

Parameters

name	name of CSG component to check

Definition at line 36 of file CSGUtils.C.

Referenced by CSG::CSGCell::CSGCell(), CSG::CSGLattice::CSGLattice(), CSG::CSGSurface::CSGSurface(), and CSG::CSGUniverse::CSGUniverse().

{
  // Check if invalid symbol is present in the name. These include whitespaces and symbols
  // for halfspaces and region operators
  const std::regex invalid_symbols(R"([\+\-~|& ])");
  std::smatch matches;

  if (std::regex_search(name, matches, invalid_symbols))
  {
    char matched_char = name[matches.position(0)];
    if (matched_char == ' ')
      mooseError("Detected whitespace in CSG component with name ", name, ". This is not allowed.");
    else
      mooseError("Invalid symbol in CSG component with name ", name, ": ", matched_char);
  }
}

getInnerRegion()

CSGRegion CSGUtils::getInnerRegion	(	const std::vector< std::reference_wrapper< const CSG::CSGSurface >> &	surfaces,
		const libMesh::Point &	origin
	)

Get inner region of given surfaces, defined as the intersection of halfspaces of each surface.

Here, the halfspace direction is determined based on the origin.

Parameters

surfaces	List of references to surfaces used to define inner region. This should ideally be defined with the minimum number of surfaces needed to enclose the region
origin	Point used to determine halfspace direction when defining intersected region

Returns

inner region defined by surfaces

Definition at line 19 of file CSGUtils.C.

{
  CSGRegion inner_region;
  for (const auto & surf_ref : surfaces)
  {
    const auto & surf = surf_ref.get();
    const auto direction = surf.getHalfspaceFromPoint(origin);
    auto halfspace = (direction == CSGSurface::Halfspace::POSITIVE) ? +surf : -surf;
    inner_region = (inner_region.getRegionType() == CSGRegion::RegionType::EMPTY)
                       ? halfspace
                       : inner_region & halfspace;
  }
  return inner_region;
}