libigl/tests/include/igl/copyleft/cgal/outer_facet.cpp

#include <test_common.h>

#include <igl/copyleft/cgal/outer_facet.h>

TEST_CASE("OuterFacet: Simple", "[igl/copyleft/cgal]")
{
    Eigen::MatrixXd V;
    Eigen::MatrixXi F;
    igl::read_triangle_mesh(test_common::data_path("cube.obj"), V, F);

    const size_t num_faces = F.rows();

    Eigen::VectorXi I(num_faces);
    I.setLinSpaced(num_faces, 0, num_faces-1);

    Eigen::Index fid = num_faces + 1;
    bool flipped;
    igl::copyleft::cgal::outer_facet(V, F, I, fid, flipped);

    REQUIRE (num_faces > fid);
    REQUIRE (!flipped);
}

TEST_CASE("OuterFacet: DuplicatedOppositeFaces", "[igl/copyleft/cgal]")
{
    Eigen::MatrixXd V;
    Eigen::MatrixXi F1;
    igl::read_triangle_mesh(test_common::data_path("cube.obj"), V, F1);

    Eigen::MatrixXi F2 = F1;
    F2.col(0).swap(F2.col(1));

    Eigen::MatrixXi F(F1.rows()*2, F1.cols());
    F << F1, F2;

    Eigen::VectorXi I(F.rows());
    I.setLinSpaced(F.rows(), 0, F.rows()-1);

    Eigen::Index fid = F.rows() + 1;
    bool flipped;
    igl::copyleft::cgal::outer_facet(V, F, I, fid, flipped);

    REQUIRE (F.rows() > fid);
    REQUIRE (!flipped);
}

TEST_CASE("OuterFacet: FullyDegnerated", "[igl/copyleft/cgal]")
{
    Eigen::MatrixXd V;
    Eigen::MatrixXi F;
    igl::read_triangle_mesh(test_common::data_path("degenerated.obj"), V, F);

    Eigen::VectorXi I(F.rows());
    I.setLinSpaced(F.rows(), 0, F.rows()-1);

    Eigen::Index fid = F.rows() + 1;
    bool flipped;
    igl::copyleft::cgal::outer_facet(V, F, I, fid, flipped);

    REQUIRE (F.rows() > fid);
    REQUIRE (!flipped);
}

TEST_CASE("OuterFacet: InvertedNormal", "[igl/copyleft/cgal]")
{
    Eigen::MatrixXd V;
    Eigen::MatrixXi F;
    igl::read_triangle_mesh(test_common::data_path("cube.obj"), V, F);
    F.col(0).swap(F.col(1));

    Eigen::VectorXi I(F.rows());
    I.setLinSpaced(F.rows(), 0, F.rows()-1);

    Eigen::Index fid = F.rows() + 1;
    bool flipped;
    igl::copyleft::cgal::outer_facet(V, F, I, fid, flipped);

    REQUIRE (F.rows() > fid);
    REQUIRE (flipped);
}

TEST_CASE("OuterFacet: SliverTet", "[igl/copyleft/cgal]")
{
    Eigen::MatrixXd V;
    Eigen::MatrixXi F;
    igl::read_triangle_mesh(test_common::data_path("sliver_tet.ply"), V, F);

    Eigen::VectorXi I(F.rows());
    I.setLinSpaced(F.rows(), 0, F.rows()-1);

    Eigen::Index fid = F.rows() + 1;
    bool flipped;
    igl::copyleft::cgal::outer_facet(V, F, I, fid, flipped);

    REQUIRE (F.rows() > fid);
    REQUIRE (!flipped);
}