Nativepolytope: Fixed a bug in quickhull when invoked on just a single point.

src/storm/storage/geometry/nativepolytopeconversion/QuickHull.cpp

@@ -128,11 +128,17 @@ namespace storm {
                     // get one hyperplane that holds all points
                     const uint_fast64_t dimension = points.front().rows();
                     EigenVector refPoint = points.front();
-                    EigenMatrix constraints(points.size() - 1, dimension);
-                    for (unsigned row = 1; row < points.size(); ++row) {
-                        constraints.row(row - 1) = points[row] - refPoint;
+                    EigenVector normal;
+                    if (points.size() == 1) {
+                        normal.resize(dimension);
+                        normal(0) = storm::utility::one<ValueType>();
+                    } else {
+                        EigenMatrix constraints(points.size() - 1, dimension);
+                        for (unsigned row = 1; row < points.size(); ++row) {
+                            constraints.row(row - 1) = points[row] - refPoint;
+                        }
+                        normal = constraints.fullPivLu().kernel().col(0);
                     }
-                    EigenVector normal = constraints.fullPivLu().kernel().col(0);
 
                     // Eigen returns the column vector 0...0 if the kernel is empty (i.e., there is no such hyperplane)
                     if (normal.isZero()) {