[#27] QHull DOD refactoring

Epona/include/Epona/QuickhullConvexHull.hpp

@@ -25,6 +25,23 @@
 
 namespace epona
 {
+struct ConvexHull
+{
+    HalfEdgeDataStructure heds;
+    std::vector<glm::uvec3> faces;
+    std::vector<uint32_t> vertices;
+};
+
+inline ConvexHull CalculateConvexHull(std::vector<glm::vec3> const& vertices)
+{
+    return {};
+}
+
+inline bool InsertVertexConvexHull(ConvexHull& qhull, std::vector<glm::vec3> const& vertices, uint32_t index)
+{
+    return false;
+}
+
 /**
  * @brief Quickhull convex hull calculation algorithm
  * @tparam VerticesType STL compatible random access glm::vec3 container

Epona/sources/Analysis.cpp

@@ -3,6 +3,7 @@
 * This code is licensed under the MIT license (MIT)
 * (http://opensource.org/licenses/MIT)
 */
+#include <Epona/FloatingPoint.hpp>
 #include <Epona/Analysis.hpp>
 #include <glm/glm.hpp>
 
@@ -29,10 +30,17 @@ glm::vec3 epona::CalculateBarycentricCoordinates(glm::vec3 p, glm::vec3 a, glm::
     float const d21 = glm::dot(v2, v1);
     float const denominator = d00 * d11 - d01 * d01;
 
-    glm::vec3 coordinates;
-    coordinates.y = (d11 * d20 - d01 * d21) / denominator;
-    coordinates.z = (d00 * d21 - d01 * d20) / denominator;
-    coordinates.x = 1.0f - coordinates.y - coordinates.z;
+    glm::vec3 coordinates{ 1, 0, 0 };
+
+    if (!fp::IsZero(denominator))
+    {
+        coordinates.y = (d11 * d20 - d01 * d21) / denominator;
+        coordinates.z = (d00 * d21 - d01 * d20) / denominator;
+        coordinates.x = 1.0f - coordinates.y - coordinates.z;
+    }
+
+    assert(!isnan(coordinates.x) && !isnan(coordinates.y) && !isnan(coordinates.z));
+    assert(!isinf(coordinates.x) && !isinf(coordinates.y) && !isinf(coordinates.z));
 
     return coordinates;
 }