snappyHexMesh with empty regions in surfaces

Summary

shm can potentially snap to a particular region. This is done by constructing a per-region search tree (triSurfaceRegionSearch). triSurfaceRegionSearch fails if there are some zero-sized regions.

Steps to reproduce

Have surfaces with zero-sized regions. Run with

snapControls:multiRegionFeatureSnap true

Example case

What is the current bug behaviour?

out-of-bounds access in Map

What is the expected correct behavior?

Relevant logs and/or images

Environment information

OpenFOAM version : Operating system : Compiler :

Possible fixes