I have a set of polygons of which I want to find the union. In order to allow for some floating point inaccuracies for polygons that touch each other, I am using boost::geometry::buffer with a inflation strategy via boost::geometry::strategy::buffer::distance_symmetric<double>. For certain inputs, the buffer operation produces results that I do not understand: