Merge pull request #1409 from vissarion/feat/polysurf_is_valid

feat: Implement is_valid algorithm for polyhedral surfaces

Author: vissarion

doc/concept/polyhedral_surface.qbk

@@ -25,6 +25,7 @@ The PolyhedralSurface Concept is defined as following:
 * It must behave like a Boost.Range Random Access Range
 * The type defined by the metafunction `range_value<...>::type` must fulfill
     the [link geometry.reference.concepts.concept_polygon Polygon Concept]
+* It must be 3-dimensional and cartesian
 
 [heading Rules]
 

include/boost/geometry/algorithms/detail/is_valid/implementation.hpp

@@ -14,6 +14,7 @@
 #include <boost/geometry/algorithms/detail/is_valid/pointlike.hpp>
 #include <boost/geometry/algorithms/detail/is_valid/linear.hpp>
 #include <boost/geometry/algorithms/detail/is_valid/polygon.hpp>
+#include <boost/geometry/algorithms/detail/is_valid/polyhedral_surface.hpp>
 #include <boost/geometry/algorithms/detail/is_valid/multipolygon.hpp>
 #include <boost/geometry/algorithms/detail/is_valid/ring.hpp>
 #include <boost/geometry/algorithms/detail/is_valid/segment.hpp>

include/boost/geometry/algorithms/detail/is_valid/interface.hpp

@@ -26,7 +26,8 @@
 #include <boost/geometry/strategies/default_strategy.hpp>
 #include <boost/geometry/strategies/detail.hpp>
 #include <boost/geometry/strategies/relate/services.hpp>
-
+#include <boost/geometry/strategies/relate3/services.hpp>
+#include <boost/geometry/strategies/relate3/cartesian.hpp>
 
 namespace boost { namespace geometry
 {
@@ -77,13 +78,15 @@ struct is_valid<default_strategy, false>
                              default_strategy)
     {
         // NOTE: Currently the strategy is only used for Areal geometries
-        typedef typename strategies::relate::services::default_strategy
+        // Select strategy_type based on Geometry's dimension
+        using strategy_type = typename std::conditional_t
             <
-                Geometry, Geometry
-            >::type strategy_type;
+                geometry::dimension<Geometry>::value == 3,
+                typename strategies::relate3::services::default_strategy<Geometry, Geometry>,
+                typename strategies::relate::services::default_strategy<Geometry, Geometry>
+            >::type;
 
-        return dispatch::is_valid<Geometry>::apply(geometry, visitor,
-                                                   strategy_type());
+        return dispatch::is_valid<Geometry>::apply(geometry, visitor, strategy_type());
     }
 };