Test all coeffs cases

Author: zfergus

tests/src/tests/candidates/test_coefficients.cpp

@@ -7,47 +7,52 @@
 
 using namespace ipc;
 
-TEST_CASE("Face-Vertex collision stencil coeffs.", "[fv][stencil][coeffs]")
+TEST_CASE("Vertex-vertex collision stencil coeffs.", "[vv][stencil][coeffs]")
 {
-    Eigen::MatrixXd V(4, 3);
-    V << 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1;
+    Eigen::MatrixXd V(2, 3);
+    V << -1, 0, 0, /**/ 1, 0, 0;
 
-    SECTION("default") { }
-    SECTION("random") { V.row(3).setRandom(); }
-    SECTION("off triangle 0") { V.row(3) << 2, 0.5, 1; }
-    SECTION("off triangle 1") { V.row(3) << 2, 2, 1; }
-    SECTION("off triangle 2") { V.row(3) << 0.5, 2, 1; }
-    SECTION("off triangle 3") { V.row(3) << -1, 2, 1; }
-    SECTION("off triangle 4") { V.row(3) << -1, 0.5, 1; }
-    SECTION("off triangle 5") { V.row(3) << -1, -1, 1; }
-    SECTION("off triangle 6") { V.row(3) << 0.5, -1, 1; }
-    SECTION("off triangle 7") { V.row(3) << 2, -0.5, 1; }
-    SECTION("failure case 1") { V.row(3) << 0.680375, -0.211234, 0.566198; }
-    SECTION("failure case 1") { V.row(3) << -0.997497, 0.127171, -0.613392; }
+    Eigen::MatrixXi E, F;
 
-    Eigen::MatrixXi F(1, 3);
-    F << 0, 1, 2;
+    VertexVertexCandidate vv(0, 1);
 
-    Eigen::MatrixXi E;
-    igl::edges(F, E);
+    VectorMax4d coeffs = vv.compute_coefficients(V, E, F);
 
-    FaceVertexCandidate fv(0, 3);
+    std::array<VectorMax3d, 4> vertices = vv.vertices(V, E, F);
+    Eigen::Vector3d n = Eigen::Vector3d::Zero();
+    for (int i = 0; i < vv.num_vertices(); i++) {
+        n += coeffs[i] * vertices[i];
+    }
 
-    VectorMax4d coeffs = fv.compute_coefficients(V, E, F);
-    CAPTURE(
-        point_triangle_distance_type(V.row(3), V.row(0), V.row(1), V.row(2)));
-    CAPTURE(V.row(3), coeffs.transpose());
+    CHECK(n.squaredNorm() == Catch::Approx(vv.compute_distance(V, E, F)));
+    CHECK(coeffs[0] == 1);
+    CHECK(coeffs[1] == -1);
+}
 
-    CHECK(-coeffs[1] == Catch::Approx(1 + coeffs[2] + coeffs[3]));
+TEST_CASE("Edge-vertex collision stencil coeffs.", "[ev][stencil][coeffs]")
+{
+    Eigen::MatrixXd V(3, 3);
+    V << -1, 0, 0, /**/ 1, 0, 0, /**/ 0, 1, 0;
 
-    std::array<VectorMax3d, 4> vertices = fv.vertices(V, E, F);
+    SECTION("default") { }
+    SECTION("random") { V.row(2).setRandom(); }
+    SECTION("e0") { V(0, 2) = -2; }
+    SECTION("e1") { V(0, 2) = 2; }
+
+    Eigen::MatrixXi E(1, 2), F;
+    E << 0, 1;
+
+    EdgeVertexCandidate ev(0, 2);
+
+    VectorMax4d coeffs = ev.compute_coefficients(V, E, F);
+
+    std::array<VectorMax3d, 4> vertices = ev.vertices(V, E, F);
     Eigen::Vector3d n = Eigen::Vector3d::Zero();
-    for (int i = 0; i < fv.num_vertices(); i++) {
+    for (int i = 0; i < ev.num_vertices(); i++) {
         n += coeffs[i] * vertices[i];
     }
-    CAPTURE(n.transpose());
 
-    CHECK(n.squaredNorm() == Catch::Approx(fv.compute_distance(V, E, F)));
+    CHECK(n.squaredNorm() == Catch::Approx(ev.compute_distance(V, E, F)));
 }
 
 TEST_CASE("Edge-edge collision stencil coeffs.", "[ee][stencil][coeffs]")
@@ -85,4 +90,47 @@ TEST_CASE("Edge-edge collision stencil coeffs.", "[ee][stencil][coeffs]")
     }
 
     CHECK(n.squaredNorm() == Catch::Approx(ee.compute_distance(V, E, F)));
-}
+}
+
+TEST_CASE("Face-Vertex collision stencil coeffs.", "[fv][stencil][coeffs]")
+{
+    Eigen::MatrixXd V(4, 3);
+    V << 0, 0, 0, /**/ 1, 0, 0, /**/ 0, 1, 0, /**/ 0.333, 0.333, 1;
+
+    SECTION("default") { }
+    SECTION("random") { V.row(3).setRandom(); }
+    SECTION("off triangle 0") { V.row(3) << 2, 0.5, 1; }
+    SECTION("off triangle 1") { V.row(3) << 2, 2, 1; }
+    SECTION("off triangle 2") { V.row(3) << 0.5, 2, 1; }
+    SECTION("off triangle 3") { V.row(3) << -1, 2, 1; }
+    SECTION("off triangle 4") { V.row(3) << -1, 0.5, 1; }
+    SECTION("off triangle 5") { V.row(3) << -1, -1, 1; }
+    SECTION("off triangle 6") { V.row(3) << 0.5, -1, 1; }
+    SECTION("off triangle 7") { V.row(3) << 2, -0.5, 1; }
+    SECTION("failure case 1") { V.row(3) << 0.680375, -0.211234, 0.566198; }
+    SECTION("failure case 1") { V.row(3) << -0.997497, 0.127171, -0.613392; }
+
+    Eigen::MatrixXi F(1, 3);
+    F << 0, 1, 2;
+
+    Eigen::MatrixXi E;
+    igl::edges(F, E);
+
+    FaceVertexCandidate fv(0, 3);
+
+    VectorMax4d coeffs = fv.compute_coefficients(V, E, F);
+    CAPTURE(
+        point_triangle_distance_type(V.row(3), V.row(0), V.row(1), V.row(2)));
+    CAPTURE(V.row(3), coeffs.transpose());
+
+    CHECK(-coeffs[1] == Catch::Approx(1 + coeffs[2] + coeffs[3]));
+
+    std::array<VectorMax3d, 4> vertices = fv.vertices(V, E, F);
+    Eigen::Vector3d n = Eigen::Vector3d::Zero();
+    for (int i = 0; i < fv.num_vertices(); i++) {
+        n += coeffs[i] * vertices[i];
+    }
+    CAPTURE(n.transpose());
+
+    CHECK(n.squaredNorm() == Catch::Approx(fv.compute_distance(V, E, F)));
+}