Fix quadrilateral boundary segments in GmshReader
- Add support for quadrilaterals as boundary segments
- Throw an exception if unsupported boundary segment type is used. Otherwise segment insertion indices will be messed up even for supported segment types.
- Add test case with unit cube and quadrilateral boundary segments.