123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177 |
- #include <test_common.h>
- #include <igl/cotmatrix.h>
- TEST_CASE("cotmatrix: constant_in_null_space", "[igl]")
- {
- const auto test_case = [](const std::string ¶m)
- {
- Eigen::MatrixXd V;
- Eigen::MatrixXi F;
- Eigen::SparseMatrix<double> L;
-
- test_common::load_mesh(param, V, F);
- igl::cotmatrix(V,F,L);
- REQUIRE (L.rows() == V.rows());
- REQUIRE (L.cols() == L.rows());
- Eigen::VectorXd C = Eigen::VectorXd::Ones(L.rows());
- Eigen::VectorXd Z = Eigen::VectorXd::Zero(L.rows());
-
-
-
- REQUIRE (1e-12 > ((L*C)-(Z)).norm());
- };
- test_common::run_test_cases(test_common::all_meshes(), test_case);
- }
- TEST_CASE("cotmatrix: cube", "[igl]")
- {
-
- const double epsilon = 1e-15;
- Eigen::MatrixXd V;
- Eigen::MatrixXi F;
-
- test_common::load_mesh("cube.obj", V, F);
-
- Eigen::MatrixXd V_huge = V * 1.0e8;
-
- Eigen::MatrixXd V_tiny = V * 1.0e-8;
-
-
-
-
-
- Eigen::SparseMatrix<double> L1;
- igl::cotmatrix(V,F,L1);
- REQUIRE (L1.rows() == V.rows());
- REQUIRE (L1.cols() == V.rows());
- for(int f = 0;f<L1.rows();f++)
- {
- #ifdef IGL_EDGE_LENGTHS_SQUARED_H
-
- REQUIRE (L1.coeff(f,f) == -3.0);
- REQUIRE (L1.row(f).sum() == 0.0);
- REQUIRE (L1.col(f).sum() == 0.0);
- #else
-
- REQUIRE (L1.coeff(f,f) == Approx (-3.0).margin( epsilon));
- REQUIRE (L1.row(f).sum() == Approx (0.0).margin( epsilon));
- REQUIRE (L1.col(f).sum() == Approx (0.0).margin( epsilon));
- #endif
- }
-
- igl::cotmatrix(V_huge,F,L1);
- REQUIRE (L1.rows() == V.rows());
- REQUIRE (L1.cols() == V.rows());
- for(int f = 0;f<L1.rows();f++)
- {
- REQUIRE (L1.coeff(f,f) == Approx (-3.0).margin( epsilon));
- REQUIRE (L1.row(f).sum() == Approx (0.0).margin( epsilon));
- REQUIRE (L1.col(f).sum() == Approx (0.0).margin( epsilon));
- }
-
- igl::cotmatrix(V_tiny,F,L1);
- REQUIRE (L1.rows() == V.rows());
- REQUIRE (L1.cols() == V.rows());
- for(int f = 0;f<L1.rows();f++)
- {
- REQUIRE (L1.coeff(f,f) == Approx (-3.0).margin( epsilon));
- REQUIRE (L1.row(f).sum() == Approx (0.0).margin( epsilon));
- REQUIRE (L1.col(f).sum() == Approx (0.0).margin( epsilon));
- }
- }
- TEST_CASE("cotmatrix: tetrahedron", "[igl]")
- {
-
- const double epsilon = 1e-15;
- Eigen::MatrixXd V;
- Eigen::MatrixXi F;
-
- test_common::load_mesh("cube.obj", V, F);
-
-
- Eigen::MatrixXi F_equi(4,3);
- F_equi << 4,6,1,
- 6,4,3,
- 4,1,3,
- 1,6,3;
-
- Eigen::MatrixXd V_huge = V * 1.0e8;
-
- Eigen::MatrixXd V_tiny = V * 1.0e-8;
-
-
-
-
-
- Eigen::SparseMatrix<double> L1;
-
- igl::cotmatrix(V,F_equi,L1);
- REQUIRE (L1.rows() == V.rows());
- REQUIRE (L1.cols() == V.rows());
- for(int f = 0;f<L1.rows();f++)
- {
-
- if (L1.coeff(f,f) < -0.1)
- REQUIRE (L1.coeff(f,f) == Approx (-3 / tan(M_PI / 3.0)).margin( epsilon));
- else
- REQUIRE (L1.coeff(f,f) == Approx (0.0).margin( epsilon));
- #ifdef IGL_EDGE_LENGTHS_SQUARED_H
-
- REQUIRE (L1.row(f).sum() == 0.0);
- REQUIRE (L1.col(f).sum() == 0.0);
- #else
-
- REQUIRE (L1.row(f).sum() == Approx (0.0).margin( epsilon));
- REQUIRE (L1.col(f).sum() == Approx (0.0).margin( epsilon));
- #endif
- }
-
- igl::cotmatrix(V_huge,F_equi,L1);
- REQUIRE (L1.rows() == V.rows());
- REQUIRE (L1.cols() == V.rows());
- for(int f = 0;f<L1.rows();f++)
- {
-
- if (L1.coeff(f,f) < -0.1)
- REQUIRE (L1.coeff(f,f) == Approx (-3 / tan(M_PI / 3.0)).margin( epsilon));
- else
- REQUIRE (L1.coeff(f,f) == Approx (0.0).margin( epsilon));
- REQUIRE (L1.row(f).sum() == Approx (0.0).margin( epsilon));
- REQUIRE (L1.col(f).sum() == Approx (0.0).margin( epsilon));
- }
-
- igl::cotmatrix(V_tiny,F_equi,L1);
- REQUIRE (L1.rows() == V.rows());
- REQUIRE (L1.cols() == V.rows());
- for(int f = 0;f<L1.rows();f++)
- {
-
- if (L1.coeff(f,f) < -0.1)
- REQUIRE (L1.coeff(f,f) == Approx (-3 / tan(M_PI / 3.0)).margin( epsilon));
- else
- REQUIRE (L1.coeff(f,f) == Approx (0.0).margin( epsilon));
- REQUIRE (L1.row(f).sum() == Approx (0.0).margin( epsilon));
- REQUIRE (L1.col(f).sum() == Approx (0.0).margin( epsilon));
- }
- }
|